R Under development (unstable) (2025-09-04 r88794 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(capushe)
> data(datacapushe)
> data(datavalidcapushe)
> data(datapartialcapushe)
> 
> dataDDSE=data(datacapushe)
> print(dataDDSE)
[1] "datacapushe"
> summary(dataDDSE)
   Length     Class      Mode 
        1 character character 
> # plot(dataDDSE,newwindow=FALSE)
> 
> # DDSE(datacapushe,100)
> # DDSE(datacapushe,pct=10)
> # DDSE(datacapushe,pct=-0.5)
> # DDSE(datacapushe,pct="a")
> 
> 
> DDSE(datacapushe,point=10)
[1] "K=21"
> # DDSE(datacapushe,point=45)
> # DDSE(datacapushe,point=51)
> # DDSE(datacapushe,point=1.5)
> # DDSE(datacapushe,point="a")
> # DDSE(datacapushe,point=-3)
> 
> library(MASS)
> DDSE(datacapushe,psi.rlm=psi.bisquare)
[1] "K=21"
> DDSE(datacapushe,psi.rlm="lm")
[1] "K=21"
Warning message:
In DDSE(datacapushe, psi.rlm = "lm") :
  The function lm is used instead of rlm
> DDSE(datacapushe,psi.rlm=1)
[1] "K=21"
Warning message:
In DDSE(datacapushe, psi.rlm = 1) : The function lm is used instead of rlm
> # DDSE(datacapushe,psi.rlm=psi.huber)
> # DDSE(datacapushe,psi.rlm=psi.hampel)
> 
> DDSE(datacapushe,scoef=3)
[1] "K=21"
> # DDSE(datacapushe,scoef="a")
> DDSE(datacapushe,scoef=-1)
[1] "K=50"
> DDSE(datacapushe,scoef=Inf)
[1] "K=1"
> 
> # plot(dataDDSE,newwindow=FALSE)
> # plot(dataDDSE,newwindow=3)
> 
> # DDSE(datacapushe[1:5,])
> 
> 
> dataNan=datacapushe
> dataNan[2,3]=NaN
> DDSE(dataNan)
[1] "K=21"
Warning message:
In DDSE(dataNan) : 1 line has been removed by DDSE
> dataNan[15,1]=NaN
> DDSE(dataNan)
[1] "K=21"
Warning message:
In DDSE(dataNan) : 1 line has been removed by DDSE
> dataNan[15,3]=NaN
> DDSE(dataNan)
[1] "K=21"
Warning message:
In DDSE(dataNan) : 2 lines have been removed by DDSE
> dataNan[8:50,4]=rep(NaN,43)
> # DDSE(dataNan)
> 
> datacomplex=datacapushe
> datacomplex[4,3]=-3
> # DDSE(datacomplex)
> datacomplex[4,3]=0.001
> # DDSE(datacomplex)
> 
> datapen=datacapushe
> datapen[6,]=datapen[5,]
> datapen[6,4]=1
> plot(DDSE(datapen),newwindow=FALSE)
Warning messages:
1: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
2: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
3: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
4: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
> datapen[6,4]=0.005
> plot(DDSE(datapen),newwindow=FALSE)
Warning messages:
1: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
2: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
3: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
4: In title(...) :
  font metrics unknown for character 0x0a in encoding cp1252
> 
> DDSEpartial=DDSE(datapartialcapushe)
> validation(DDSEpartial,datavalidcapushe,newwindow=FALSE)
> 
> dataNan=datavalidcapushe
> dataNan[2,3]=NaN
> validation(DDSEpartial,dataNan,newwindow=FALSE)
Warning message:
In .local(x, data2, ...) : 1 line in data2 has been removed by validation
> dataNan[,3]=rep(NaN,3)
> validation(DDSEpartial,dataNan,newwindow=FALSE)
Warning message:
In .local(x, data2, ...) : 3 lines in data2 have been removed by validation
> 
> datacontrast=datacapushe
> datacontrast[6,4]=-3
> DDSE(datacontrast)
[1] "K=6"
> 
> proc.time()
   user  system elapsed 
   2.53    0.15    2.67