Last updated on 2026-02-01 01:49:20 CET.
| Package | ERROR | NOTE | OK |
|---|---|---|---|
| Bessel | 13 | ||
| bitops | 13 | ||
| CLA | 13 | ||
| classGraph | 13 | ||
| cluster | 5 | 8 | |
| cobs | 13 | ||
| copula | 5 | 8 | |
| diptest | 13 | ||
| DPQ | 13 | ||
| DPQmpfr | 13 | ||
| expm | 13 | ||
| fracdiff | 13 | ||
| lokern | 13 | ||
| longmemo | 13 | ||
| lpridge | 13 | ||
| nor1mix | 13 | ||
| plugdensity | 2 | 11 | |
| Rmpfr | 13 | ||
| robustbase | 1 | 12 | |
| robustX | 13 | ||
| round | 13 | ||
| sca | 2 | 11 | |
| sfsmisc | 13 | ||
| simest | 3 | 10 | |
| stabledist | 13 | ||
| supclust | 2 | 11 | |
| VLMC | 13 |
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: NOTE: 5, OK: 8
Version: 2.1.8.1
Check: tests
Result: NOTE
Running ‘agnes-ex.R’ [2s/3s]
Comparing ‘agnes-ex.Rout’ to ‘agnes-ex.Rout.save’ ... OK
Running ‘clara-NAs.R’ [0s/1s]
Comparing ‘clara-NAs.Rout’ to ‘clara-NAs.Rout.save’ ...
47c47
< Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
---
> Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
Running ‘clara-ex.R’ [2s/3s]
Comparing ‘clara-ex.Rout’ to ‘clara-ex.Rout.save’ ... OK
Running ‘clara-gower.R’ [0s/0s]
Running ‘clara.R’ [3s/4s]
Comparing ‘clara.Rout’ to ‘clara.Rout.save’ ... OK
Running ‘clusplot-out.R’ [1s/1s]
Comparing ‘clusplot-out.Rout’ to ‘clusplot-out.Rout.save’ ... OK
Running ‘daisy-ex.R’ [1s/1s]
Comparing ‘daisy-ex.Rout’ to ‘daisy-ex.Rout.save’ ... OK
Running ‘diana-boots.R’ [1s/2s]
Running ‘diana-ex.R’ [0s/1s]
Comparing ‘diana-ex.Rout’ to ‘diana-ex.Rout.save’ ... OK
Running ‘ellipsoid-ex.R’ [0s/1s]
Comparing ‘ellipsoid-ex.Rout’ to ‘ellipsoid-ex.Rout.save’ ... OK
Running ‘fanny-ex.R’ [1s/1s]
Comparing ‘fanny-ex.Rout’ to ‘fanny-ex.Rout.save’ ... OK
Running ‘mona.R’ [1s/1s]
Comparing ‘mona.Rout’ to ‘mona.Rout.save’ ...
213c213
< NAs :36 NAs :25
---
> NA's :36 NA's :25
218c218
< NAs :39 NAs :2 NAs :2
---
> NA's :39 NA's :2 NA's :2
222c222
< TRUE :73 TRUE :52 TRUE :97 NAs: 4
---
> TRUE :73 TRUE :52 TRUE :97 NA's: 4
227c227
< NAs: 4 NAs: 4 NAs: 4 NAs: 4
---
> NA's: 4 NA's: 4 NA's: 4 NA's: 4
229,232c229,232
< seasver everalw everparti elaio endozoo epizoo aquat windgl
< 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128 0 :108
< 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3 1 : 23
< NAs: 4 NAs: 4 NAs: 4 NAs: 5 NAs: 5 NAs: 5 NAs: 5 NAs: 5
---
> seasver everalw everparti elaio endozoo epizoo aquat
> 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128
> 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3
> NA's: 4 NA's: 4 NA's: 4 NA's: 5 NA's: 5 NA's: 5 NA's: 5
234,237c234,237
< unsp
< 0 :93
< 1 :38
< NAs: 5
---
> windgl unsp
> 0 :108 0 :93
> 1 : 23 1 :38
> NA's: 5 NA's: 5
Running ‘pam.R’ [43s/55s]
Comparing ‘pam.Rout’ to ‘pam.Rout.save’ ... OK
Running ‘silhouette-default.R’ [3s/4s]
Comparing ‘silhouette-default.Rout’ to ‘silhouette-default.Rout.save’ ... OK
Running ‘sweep-ex.R’ [0s/0s]
Flavor: r-devel-linux-x86_64-debian-clang
Version: 2.1.8.1
Check: tests
Result: NOTE
Running ‘agnes-ex.R’ [1s/2s]
Comparing ‘agnes-ex.Rout’ to ‘agnes-ex.Rout.save’ ... OK
Running ‘clara-NAs.R’ [0s/1s]
Comparing ‘clara-NAs.Rout’ to ‘clara-NAs.Rout.save’ ...
47c47
< Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
---
> Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
Running ‘clara-ex.R’ [2s/3s]
Comparing ‘clara-ex.Rout’ to ‘clara-ex.Rout.save’ ... OK
Running ‘clara-gower.R’ [0s/1s]
Running ‘clara.R’ [2s/4s]
Comparing ‘clara.Rout’ to ‘clara.Rout.save’ ... OK
Running ‘clusplot-out.R’ [0s/1s]
Comparing ‘clusplot-out.Rout’ to ‘clusplot-out.Rout.save’ ... OK
Running ‘daisy-ex.R’ [0s/1s]
Comparing ‘daisy-ex.Rout’ to ‘daisy-ex.Rout.save’ ... OK
Running ‘diana-boots.R’ [1s/1s]
Running ‘diana-ex.R’ [0s/1s]
Comparing ‘diana-ex.Rout’ to ‘diana-ex.Rout.save’ ... OK
Running ‘ellipsoid-ex.R’ [0s/0s]
Comparing ‘ellipsoid-ex.Rout’ to ‘ellipsoid-ex.Rout.save’ ... OK
Running ‘fanny-ex.R’ [0s/1s]
Comparing ‘fanny-ex.Rout’ to ‘fanny-ex.Rout.save’ ... OK
Running ‘mona.R’ [0s/1s]
Comparing ‘mona.Rout’ to ‘mona.Rout.save’ ...
213c213
< NAs :36 NAs :25
---
> NA's :36 NA's :25
218c218
< NAs :39 NAs :2 NAs :2
---
> NA's :39 NA's :2 NA's :2
222c222
< TRUE :73 TRUE :52 TRUE :97 NAs: 4
---
> TRUE :73 TRUE :52 TRUE :97 NA's: 4
227c227
< NAs: 4 NAs: 4 NAs: 4 NAs: 4
---
> NA's: 4 NA's: 4 NA's: 4 NA's: 4
229,232c229,232
< seasver everalw everparti elaio endozoo epizoo aquat windgl
< 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128 0 :108
< 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3 1 : 23
< NAs: 4 NAs: 4 NAs: 4 NAs: 5 NAs: 5 NAs: 5 NAs: 5 NAs: 5
---
> seasver everalw everparti elaio endozoo epizoo aquat
> 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128
> 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3
> NA's: 4 NA's: 4 NA's: 4 NA's: 5 NA's: 5 NA's: 5 NA's: 5
234,237c234,237
< unsp
< 0 :93
< 1 :38
< NAs: 5
---
> windgl unsp
> 0 :108 0 :93
> 1 : 23 1 :38
> NA's: 5 NA's: 5
Running ‘pam.R’ [26s/29s]
Comparing ‘pam.Rout’ to ‘pam.Rout.save’ ... OK
Running ‘silhouette-default.R’ [2s/2s]
Comparing ‘silhouette-default.Rout’ to ‘silhouette-default.Rout.save’ ... OK
Running ‘sweep-ex.R’ [0s/0s]
Flavor: r-devel-linux-x86_64-debian-gcc
Version: 2.1.8.1
Check: tests
Result: NOTE
Running ‘agnes-ex.R’
Comparing ‘agnes-ex.Rout’ to ‘agnes-ex.Rout.save’ ... OK
Running ‘clara-NAs.R’
Comparing ‘clara-NAs.Rout’ to ‘clara-NAs.Rout.save’ ...
47c47
< Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
---
> Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
Running ‘clara-ex.R’
Comparing ‘clara-ex.Rout’ to ‘clara-ex.Rout.save’ ... OK
Running ‘clara-gower.R’
Running ‘clara.R’
Comparing ‘clara.Rout’ to ‘clara.Rout.save’ ... OK
Running ‘clusplot-out.R’
Comparing ‘clusplot-out.Rout’ to ‘clusplot-out.Rout.save’ ... OK
Running ‘daisy-ex.R’
Comparing ‘daisy-ex.Rout’ to ‘daisy-ex.Rout.save’ ... OK
Running ‘diana-boots.R’
Running ‘diana-ex.R’
Comparing ‘diana-ex.Rout’ to ‘diana-ex.Rout.save’ ... OK
Running ‘ellipsoid-ex.R’
Comparing ‘ellipsoid-ex.Rout’ to ‘ellipsoid-ex.Rout.save’ ... OK
Running ‘fanny-ex.R’
Comparing ‘fanny-ex.Rout’ to ‘fanny-ex.Rout.save’ ... OK
Running ‘mona.R’
Comparing ‘mona.Rout’ to ‘mona.Rout.save’ ...
213c213
< NAs :36 NAs :25
---
> NA's :36 NA's :25
218c218
< NAs :39 NAs :2 NAs :2
---
> NA's :39 NA's :2 NA's :2
222c222
< TRUE :73 TRUE :52 TRUE :97 NAs: 4
---
> TRUE :73 TRUE :52 TRUE :97 NA's: 4
227c227
< NAs: 4 NAs: 4 NAs: 4 NAs: 4
---
> NA's: 4 NA's: 4 NA's: 4 NA's: 4
229,232c229,232
< seasver everalw everparti elaio endozoo epizoo aquat windgl
< 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128 0 :108
< 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3 1 : 23
< NAs: 4 NAs: 4 NAs: 4 NAs: 5 NAs: 5 NAs: 5 NAs: 5 NAs: 5
---
> seasver everalw everparti elaio endozoo epizoo aquat
> 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128
> 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3
> NA's: 4 NA's: 4 NA's: 4 NA's: 5 NA's: 5 NA's: 5 NA's: 5
234,237c234,237
< unsp
< 0 :93
< 1 :38
< NAs: 5
---
> windgl unsp
> 0 :108 0 :93
> 1 : 23 1 :38
> NA's: 5 NA's: 5
Running ‘pam.R’ [71s/88s]
Comparing ‘pam.Rout’ to ‘pam.Rout.save’ ... OK
Running ‘silhouette-default.R’
Comparing ‘silhouette-default.Rout’ to ‘silhouette-default.Rout.save’ ... OK
Running ‘sweep-ex.R’
Flavor: r-devel-linux-x86_64-fedora-clang
Version: 2.1.8.1
Check: tests
Result: NOTE
Running ‘agnes-ex.R’
Comparing ‘agnes-ex.Rout’ to ‘agnes-ex.Rout.save’ ... OK
Running ‘clara-NAs.R’
Comparing ‘clara-NAs.Rout’ to ‘clara-NAs.Rout.save’ ...
47c47
< Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
---
> Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
Running ‘clara-ex.R’
Comparing ‘clara-ex.Rout’ to ‘clara-ex.Rout.save’ ... OK
Running ‘clara-gower.R’
Running ‘clara.R’ [5s/11s]
Comparing ‘clara.Rout’ to ‘clara.Rout.save’ ... OK
Running ‘clusplot-out.R’
Comparing ‘clusplot-out.Rout’ to ‘clusplot-out.Rout.save’ ... OK
Running ‘daisy-ex.R’
Comparing ‘daisy-ex.Rout’ to ‘daisy-ex.Rout.save’ ... OK
Running ‘diana-boots.R’
Running ‘diana-ex.R’
Comparing ‘diana-ex.Rout’ to ‘diana-ex.Rout.save’ ... OK
Running ‘ellipsoid-ex.R’
Comparing ‘ellipsoid-ex.Rout’ to ‘ellipsoid-ex.Rout.save’ ... OK
Running ‘fanny-ex.R’
Comparing ‘fanny-ex.Rout’ to ‘fanny-ex.Rout.save’ ... OK
Running ‘mona.R’
Comparing ‘mona.Rout’ to ‘mona.Rout.save’ ...
213c213
< NAs :36 NAs :25
---
> NA's :36 NA's :25
218c218
< NAs :39 NAs :2 NAs :2
---
> NA's :39 NA's :2 NA's :2
222c222
< TRUE :73 TRUE :52 TRUE :97 NAs: 4
---
> TRUE :73 TRUE :52 TRUE :97 NA's: 4
227c227
< NAs: 4 NAs: 4 NAs: 4 NAs: 4
---
> NA's: 4 NA's: 4 NA's: 4 NA's: 4
229,232c229,232
< seasver everalw everparti elaio endozoo epizoo aquat windgl
< 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128 0 :108
< 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3 1 : 23
< NAs: 4 NAs: 4 NAs: 4 NAs: 5 NAs: 5 NAs: 5 NAs: 5 NAs: 5
---
> seasver everalw everparti elaio endozoo epizoo aquat
> 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128
> 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3
> NA's: 4 NA's: 4 NA's: 4 NA's: 5 NA's: 5 NA's: 5 NA's: 5
234,237c234,237
< unsp
< 0 :93
< 1 :38
< NAs: 5
---
> windgl unsp
> 0 :108 0 :93
> 1 : 23 1 :38
> NA's: 5 NA's: 5
Running ‘pam.R’ [85s/223s]
Comparing ‘pam.Rout’ to ‘pam.Rout.save’ ... OK
Running ‘silhouette-default.R’ [5s/13s]
Comparing ‘silhouette-default.Rout’ to ‘silhouette-default.Rout.save’ ... OK
Running ‘sweep-ex.R’
Flavor: r-devel-linux-x86_64-fedora-gcc
Version: 2.1.8.1
Check: tests
Result: NOTE
Running 'agnes-ex.R' [2s]
Comparing 'agnes-ex.Rout' to 'agnes-ex.Rout.save' ... OK
Running 'clara-NAs.R' [0s]
Comparing 'clara-NAs.Rout' to 'clara-NAs.Rout.save' ...
47c47
< Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
---
> Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
Running 'clara-ex.R' [2s]
Comparing 'clara-ex.Rout' to 'clara-ex.Rout.save' ... OK
Running 'clara-gower.R' [0s]
Running 'clara.R' [3s]
Comparing 'clara.Rout' to 'clara.Rout.save' ... OK
Running 'clusplot-out.R' [1s]
Comparing 'clusplot-out.Rout' to 'clusplot-out.Rout.save' ... OK
Running 'daisy-ex.R' [1s]
Comparing 'daisy-ex.Rout' to 'daisy-ex.Rout.save' ... OK
Running 'diana-boots.R' [2s]
Running 'diana-ex.R' [0s]
Comparing 'diana-ex.Rout' to 'diana-ex.Rout.save' ... OK
Running 'ellipsoid-ex.R' [0s]
Comparing 'ellipsoid-ex.Rout' to 'ellipsoid-ex.Rout.save' ... OK
Running 'fanny-ex.R' [1s]
Comparing 'fanny-ex.Rout' to 'fanny-ex.Rout.save' ... OK
Running 'mona.R' [1s]
Comparing 'mona.Rout' to 'mona.Rout.save' ...
213c213
< NAs :36 NAs :25
---
> NA's :36 NA's :25
218c218
< NAs :39 NAs :2 NAs :2
---
> NA's :39 NA's :2 NA's :2
222c222
< TRUE :73 TRUE :52 TRUE :97 NAs: 4
---
> TRUE :73 TRUE :52 TRUE :97 NA's: 4
227c227
< NAs: 4 NAs: 4 NAs: 4 NAs: 4
---
> NA's: 4 NA's: 4 NA's: 4 NA's: 4
229,232c229,232
< seasver everalw everparti elaio endozoo epizoo aquat windgl
< 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128 0 :108
< 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3 1 : 23
< NAs: 4 NAs: 4 NAs: 4 NAs: 5 NAs: 5 NAs: 5 NAs: 5 NAs: 5
---
> seasver everalw everparti elaio endozoo epizoo aquat
> 0 :126 0 :78 0 :112 0 :119 0 :93 0 :109 0 :128
> 1 : 6 1 :54 1 : 20 1 : 12 1 :38 1 : 22 1 : 3
> NA's: 4 NA's: 4 NA's: 4 NA's: 5 NA's: 5 NA's: 5 NA's: 5
234,237c234,237
< unsp
< 0 :93
< 1 :38
< NAs: 5
---
> windgl unsp
> 0 :108 0 :93
> 1 : 23 1 :38
> NA's: 5 NA's: 5
Running 'pam.R' [32s]
Comparing 'pam.Rout' to 'pam.Rout.save' ... OK
Running 'silhouette-default.R' [3s]
Comparing 'silhouette-default.Rout' to 'silhouette-default.Rout.save' ... OK
Running 'sweep-ex.R' [0s]
Flavor: r-devel-windows-x86_64
Current CRAN status: OK: 13
Current CRAN status: NOTE: 5, OK: 8
Version: 1.1-6
Check: relative paths in package URLs
Result: NOTE
Found the following (possibly) invalid URL:
URL: ../doc/Frank-Rmpfr
From: man/copula-package.Rd
Flavors: r-devel-linux-x86_64-debian-clang, r-devel-linux-x86_64-debian-gcc
Version: 1.1-6
Check: installed package size
Result: NOTE
installed size is 7.8Mb
sub-directories of 1Mb or more:
R 2.4Mb
doc 3.2Mb
Flavors: r-oldrel-macos-arm64, r-oldrel-macos-x86_64
Version: 1.1-6
Flags: --no-vignettes
Check: installed package size
Result: NOTE
installed size is 7.3Mb
sub-directories of 1Mb or more:
R 2.1Mb
doc 3.2Mb
Flavor: r-oldrel-windows-x86_64
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: NOTE: 2, OK: 11
Version: 0.8-5
Check: CRAN incoming feasibility
Result: NOTE
Maintainer: ‘Martin Maechler <maechler@stat.math.ethz.ch>’
No Authors@R field in DESCRIPTION.
Please add one, modifying
Authors@R: c(person(given = c("Eva", "Herrmann;", "R", "interface", "etc", "by", "Martin"),
family = "Maechler",
role = "aut",
email = "eherrmann@mathematik.tu-darmstadt.de",
comment = "C original"),
person(given = "Martin",
family = "Maechler",
role = "cre",
email = "maechler@stat.math.ethz.ch"))
as necessary.
Flavors: r-devel-linux-x86_64-debian-clang, r-devel-linux-x86_64-debian-gcc
Current CRAN status: OK: 13
Current CRAN status: ERROR: 1, OK: 12
Version: 0.99-6
Check: tests
Result: ERROR
Running ‘LTS-specials.R’ [1s/1s]
Running ‘MCD-specials.R’ [0s/1s]
Comparing ‘MCD-specials.Rout’ to ‘MCD-specials.Rout.save’ ... OK
Running ‘MT-tst.R’ [7s/9s]
Running ‘NAcoef.R’ [1s/1s]
Running ‘OGK-ex.R’ [1s/1s]
Comparing ‘OGK-ex.Rout’ to ‘OGK-ex.Rout.save’ ... OK
Running ‘Qn-Sn-plots.R’ [1s/1s]
Running ‘Rsquared.R’ [1s/1s]
Comparing ‘Rsquared.Rout’ to ‘Rsquared.Rout.save’ ... OK
Running ‘binom-ni-small.R’ [1s/1s]
Comparing ‘binom-ni-small.Rout’ to ‘binom-ni-small.Rout.save’ ... OK
Running ‘binom-no-x.R’ [0s/1s]
Running ‘comedian-tst.R’ [1s/1s]
Running ‘exact-fit-categorical.R’ [0s/1s]
Running ‘glmrob-1.R’ [8s/10s]
Running ‘glmrob-specials.R’ [0s/1s]
Running ‘huber-etc.R’ [0s/1s]
Comparing ‘huber-etc.Rout’ to ‘huber-etc.Rout.save’ ... OK
Running ‘large-values.R’ [1s/1s]
Running ‘lmrob-data.R’ [3s/3s]
Running ‘lmrob-ex12.R’ [2s/4s]
Running ‘lmrob-methods.R’ [1s/1s]
Comparing ‘lmrob-methods.Rout’ to ‘lmrob-methods.Rout.save’ ... OK
Running ‘lmrob-psifns.R’ [4s/5s]
Comparing ‘lmrob-psifns.Rout’ to ‘lmrob-psifns.Rout.save’ ... OK
Running ‘m-s-estimator.R’ [2s/2s]
Running ‘mc-etc.R’ [1s/1s]
Running ‘mc-strict.R’ [15s/20s]
Running ‘nlregrob-tst.R’ [16s/20s]
Running ‘nlrob-tst.R’ [3s/4s]
Running ‘poisson-ex.R’ [1s/1s]
Running ‘psi-rho-etc.R’ [1s/1s]
Comparing ‘psi-rho-etc.Rout’ to ‘psi-rho-etc.Rout.save’ ... OK
Running ‘small-sample.R’ [9s/12s]
Comparing ‘small-sample.Rout’ to ‘small-sample.Rout.save’ ... OK
Running ‘subsample.R’ [6s/9s]
Running ‘tlts.R’ [1s/1s]
Comparing ‘tlts.Rout’ to ‘tlts.Rout.save’ ... OK
Running ‘tmcd.R’ [6s/9s]
Running ‘weights.R’ [1s/2s]
Comparing ‘weights.Rout’ to ‘weights.Rout.save’ ... OK
Running ‘wgt-himed-xtra.R’ [3s/3s]
Running ‘wgt-himed.R’ [0s/1s]
Comparing ‘wgt-himed.Rout’ to ‘wgt-himed.Rout.save’ ... OK
Running the tests in ‘tests/NAcoef.R’ failed.
Complete output:
> ## test handing of NA coefficients / singular fits
> ## also check:
> ## -- what would have to be done if class "lm" was added.
> ## -- general compatibility to class lm.
> require(robustbase)
Loading required package: robustbase
> options(digits = 5)# -> higher chance of platform independence
>
> ## generate simple example data (almost as in ./weights.R )
> data <- expand.grid(x1=letters[1:3], x2=LETTERS[1:3], rep=1:3)
> set.seed(1)
> data$y <- rnorm(nrow(data))
> ## drop all combinations of one interaction:
> data <- subset(data, x1 != 'c' | (x2 != 'B' & x2 != 'C'))
> ## add collinear variables
> data$x3 <- rnorm(nrow(data))
> data$x4 <- rnorm(nrow(data))
> data$x5 <- data$x3 + data$x4
> ## add some NA terms
> data$y[1] <- NA
> data$x4[2:3] <- NA ## to test anova
>
> ## Classical models start with 'cm', robust just with 'rm' (or just 'm'):
> cm0 <- lm (y ~ x1*x2 + x3, data)
> cm1 <- lm (y ~ x1*x2 + x3 + x4 + x5, data)
> set.seed(2)
> rm1 <- lmrob(y ~ x1*x2 + x3 + x4 + x5, data)
> m3 <- lmrob(y ~ x1*x2 + x3 + x4, data) # same column space as rm1
> rm0 <- lmrob(y ~ x1*x2 + x3, data)
>
> ## clean version of rm1 (to check predict)
> data2 <- data.frame(y=data$y[-(1:3)], rm1$x[,!is.na(rm1$coef)])
> set.seed(2)
> rm1c <- lmrob(y ~ x1b + x1c + x2B + x2C + x3 + x4 + x1b:x2B + x1b:x2C, data2)
>
> ## add class lm to rm1 (for now)
> class(rm1) <- c(class(rm1), "lm")
> class(rm0) <- c(class(rm0), "lm")
>
> ## the full matrix (data) should be returned by model matrix (frame)
> stopifnot(all.equal(model.matrix(cm1), model.matrix(rm1)),
+ all.equal(model.frame (cm1), model.frame (rm1)))
> ## qr decomposition should be for the full data and pivots identical lm result
> qr.cm1 <- qr(cm1)$qr
> qr.rm1 <- rm1$qr$qr
> stopifnot(NCOL(qr.rm1) == NCOL(qr.cm1),
+ NROW(qr.rm1) == NROW(qr.cm1),
+ length(rm1$qr$qraux) == length(qr(cm1)$qraux),
+ all.equal(rm1$qr$pivot, qr(cm1)$pivot),
+ all.equal(dimnames(qr.rm1),dimnames(qr.cm1)))
> ## the alias function should return the same result
> stopifnot(all.equal(alias(cm1), alias(rm1)))
>
> ####
> ## these helper functions should print NAs for the dropped coefficients
> print(rm1)
Call:
lmrob(formula = y ~ x1 * x2 + x3 + x4 + x5, data = data)
\--> method = "MM"
Coefficients:
(Intercept) x1b x1c x2B x2C x3
0.4381 0.5968 0.0344 0.2012 0.1789 -0.1320
x4 x5 x1b:x2B x1c:x2B x1b:x2C x1c:x2C
-0.2155 NA -1.8763 NA -0.8651 NA
> summary(rm1) -> s1
> confint(rm1) -> ci1
> stopifnot(identical(is.na(coef(cm1)), apply(ci1, 1L, anyNA)),
+ identical(sigma(rm1), s1$ sigma),
+ identical(vcov(rm1, complete=FALSE), s1$ cov ),
+ TRUE)
>
> print(s1, showAlgo=FALSE)
Call:
lmrob(formula = y ~ x1 * x2 + x3 + x4 + x5, data = data)
\--> method = "MM"
Residuals:
Min 1Q Median 3Q Max
-1.4584 -0.3556 0.0246 0.3651 1.0296
Coefficients: (3 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.4381 0.5443 0.80 0.44
x1b 0.5968 0.6423 0.93 0.38
x1c 0.0344 0.6880 0.05 0.96
x2B 0.2012 0.7164 0.28 0.79
x2C 0.1789 0.6871 0.26 0.80
x3 -0.1320 0.4155 -0.32 0.76
x4 -0.2155 0.1694 -1.27 0.24
x5 NA NA NA NA
x1b:x2B -1.8763 1.2153 -1.54 0.16
x1c:x2B NA NA NA NA
x1b:x2C -0.8651 0.7466 -1.16 0.28
x1c:x2C NA NA NA NA
Robust residual standard error: 0.927
(3 observations deleted due to missingness)
Multiple R-squared: 0.338, Adjusted R-squared: -0.251
Convergence in 15 IRWLS iterations
Robustness weights:
2 weights are ~= 1. The remaining 16 ones are summarized as
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.787 0.937 0.985 0.952 0.988 0.994
> ci1
2.5 % 97.5 %
(Intercept) -0.79333 1.66946
x1b -0.85607 2.04973
x1c -1.52188 1.59076
x2B -1.41948 1.82189
x2C -1.37549 1.73320
x3 -1.07182 0.80783
x4 -0.59863 0.16756
x5 NA NA
x1b:x2B -4.62539 0.87283
x1c:x2B NA NA
x1b:x2C -2.55391 0.82381
x1c:x2C NA NA
> ## drop1 should return df = 0
> #drop1(rm1) ## drop.lm does not return valid results (yet)!
>
> ####
> ## methods that should just drop the NA coefficients
> ## m3 is actually the same as rm1, so anova should raise an error
> tools::assertError(anova(rm1, m3, test="Wald"))
> tools::assertError(anova(rm1, m3, test="Deviance"))
> ## but comparing rm1 and rm0 should be ok
> anova(rm1, rm0, test="Wald")
Robust Wald Test Table
Model 1: y ~ x1 * x2 + x3 + x4 + x5
Model 2: y ~ x1 * x2 + x3
Largest model fitted by lmrob(), i.e. SM
pseudoDf Test.Stat Df Pr(>chisq)
1 6
2 10 1.62 1 0.2
> anova(rm1, rm0, test="Deviance")
Robust Deviance Table
Model 1: y ~ x1 * x2 + x3 + x4 + x5
Model 2: y ~ x1 * x2 + x3
Largest model fitted by lmrob(), i.e. SM
pseudoDf Test.Stat Df Pr(>chisq)
1 6
2 10 1.4 1 0.24
> ## commands with single #:
> ## they do (or might) not return sensible results for robust fits
> ## and need to be checked again
> #cooks.distance(rm1)
> #deviance(rm1)
> #dfbeta(rm1)
> #dfbetas(rm1)
> #effects(rm1) ## fails
> #extractAIC(rm1)
> #influence(rm1)
> stopifnot(all.equal(hv1 <- hatvalues(rm1), .lmrob.hat(wqr=rm1$qr), tol=1e-15),
+ all.equal(hv1, stats:::hatvalues.lm(rm1), tol=1e-15),
+ all.equal(hat(cm1$qr), unname(hatvalues(cm1)), tol=1e-15),
+ all.equal(unname(hv1), hat(rm1$qr), tol=1e-15),
+ ## ditto :
+ all.equal(hv1c <- hatvalues(rm1c), stats:::hatvalues.lm(rm1c), tol=1e-15))
Error in match.arg(type) : 'arg' should be one of "prior", "robustness"
Calls: stopifnot ... weighted.residuals -> weights -> weights.lmrob -> match.arg
Execution halted
Flavor: r-devel-linux-x86_64-debian-clang
Current CRAN status: OK: 13
Current CRAN status: OK: 13
Current CRAN status: NOTE: 2, OK: 11
Version: 0.9-2
Check: CRAN incoming feasibility
Result: NOTE
Maintainer: ‘Martin Maechler <maechler@stat.math.ethz.ch>’
No Authors@R field in DESCRIPTION.
Please add one, modifying
Authors@R: c(person(given = "Valentin",
family = "Rousson",
role = "aut",
email = "rousson@ifspm.unizh.ch"),
person(given = "Martin",
family = "Maechler",
role = c("aut", "cre"),
email = "maechler@stat.math.ethz.ch"))
as necessary.
Flavors: r-devel-linux-x86_64-debian-clang, r-devel-linux-x86_64-debian-gcc
Current CRAN status: OK: 13
Current CRAN status: ERROR: 3, OK: 10
Version: 0.4-1-1
Check: tests
Result: ERROR
Running ‘ex_cvx.R’ [27s/27s]
Running ‘fastmerge-ex.R’ [1s/1s]
Running ‘smooth.pen_agcv.R’ [0s/1s]
Running the tests in ‘tests/smooth.pen_agcv.R’ failed.
Complete output:
> library(simest)
>
> ## Test smooth.pen.reg(*, agcv=TRUE) ---------------
>
> ## -- Using a subset of R's sunspots data ------
> osV <- abbreviate(gsub("[^[:alnum:]]", '', sub("\\(.*", '', osVersion)), 12)
>
> if(!dev.interactive(TRUE)) pdf(paste0("smooth.pen_sunsp__", osV, ".pdf"), width = 9, height=5)
>
> str(ssp <- window(sunspot.m2014, start = 1900))
Error: object 'sunspot.m2014' not found
Execution halted
Flavor: r-oldrel-macos-arm64
Version: 0.4-1-1
Check: tests
Result: ERROR
Running ‘ex_cvx.R’ [50s/71s]
Running ‘fastmerge-ex.R’ [4s/4s]
Running ‘smooth.pen_agcv.R’ [1s/1s]
Running the tests in ‘tests/smooth.pen_agcv.R’ failed.
Complete output:
> library(simest)
>
> ## Test smooth.pen.reg(*, agcv=TRUE) ---------------
>
> ## -- Using a subset of R's sunspots data ------
> osV <- abbreviate(gsub("[^[:alnum:]]", '', sub("\\(.*", '', osVersion)), 12)
>
> if(!dev.interactive(TRUE)) pdf(paste0("smooth.pen_sunsp__", osV, ".pdf"), width = 9, height=5)
>
> str(ssp <- window(sunspot.m2014, start = 1900))
Error: object 'sunspot.m2014' not found
Execution halted
Flavor: r-oldrel-macos-x86_64
Version: 0.4-1-1
Check: tests
Result: ERROR
Running 'ex_cvx.R' [87s]
Running 'fastmerge-ex.R' [7s]
Running 'smooth.pen_agcv.R' [3s]
Running the tests in 'tests/smooth.pen_agcv.R' failed.
Complete output:
> library(simest)
>
> ## Test smooth.pen.reg(*, agcv=TRUE) ---------------
>
> ## -- Using a subset of R's sunspots data ------
> osV <- abbreviate(gsub("[^[:alnum:]]", '', sub("\\(.*", '', osVersion)), 12)
>
> if(!dev.interactive(TRUE)) pdf(paste0("smooth.pen_sunsp__", osV, ".pdf"), width = 9, height=5)
>
> str(ssp <- window(sunspot.m2014, start = 1900))
Error: object 'sunspot.m2014' not found
Execution halted
Flavor: r-oldrel-windows-x86_64
Current CRAN status: OK: 13
Current CRAN status: NOTE: 2, OK: 11
Version: 1.1-1
Check: CRAN incoming feasibility
Result: NOTE
Maintainer: ‘Martin Maechler <maechler@stat.math.ethz.ch>’
No Authors@R field in DESCRIPTION.
Please add one, modifying
Authors@R: c(person(given = "Marcel",
family = "Dettling",
role = "aut",
email = "marcel.dettling@zhaw.ch"),
person(given = "Martin",
family = "Maechler",
role = c("aut", "cre"),
email = "maechler@stat.math.ethz.ch"))
as necessary.
Flavors: r-devel-linux-x86_64-debian-clang, r-devel-linux-x86_64-debian-gcc
Current CRAN status: OK: 13