Last updated on 2019-01-19 00:46:41 CET.
Flavor | Version | Tinstall | Tcheck | Ttotal | Status | Flags |
---|---|---|---|---|---|---|
r-devel-linux-x86_64-debian-clang | 1.2-2 | 2.49 | 34.42 | 36.91 | OK | |
r-devel-linux-x86_64-debian-gcc | 1.2-2 | 2.18 | 26.77 | 28.95 | OK | |
r-devel-linux-x86_64-fedora-clang | 1.2-2 | 42.92 | OK | |||
r-devel-linux-x86_64-fedora-gcc | 1.2-2 | 40.86 | OK | |||
r-devel-windows-ix86+x86_64 | 1.2-2 | 8.00 | 43.00 | 51.00 | OK | |
r-patched-linux-x86_64 | 1.2-2 | 2.56 | 31.74 | 34.30 | OK | |
r-patched-solaris-x86 | 1.2-2 | 62.80 | OK | |||
r-release-linux-x86_64 | 1.2-2 | 2.39 | 31.50 | 33.89 | OK | |
r-release-windows-ix86+x86_64 | 1.2-2 | 5.00 | 68.00 | 73.00 | OK | |
r-release-osx-x86_64 | 1.2-2 | ERROR | ||||
r-oldrel-windows-ix86+x86_64 | 1.2-2 | 3.00 | 46.00 | 49.00 | OK | |
r-oldrel-osx-x86_64 | 1.2-2 | ERROR |
Version: 1.2-2
Check: examples
Result: ERROR
Running examples in ‘robustX-Ex.R’ failed
The error most likely occurred in:
> ### Name: rbwheel
> ### Title: Multivariate Barrow Wheel Distribution Random Vectors
> ### Aliases: rbwheel
> ### Keywords: distribution robust
>
> ### ** Examples
>
> set.seed(17)
> rX8 <- rbwheel(1000,8, fullResult = TRUE, scaleAfter=FALSE)
> with(rX8, stopifnot(all.equal(X, X0 %*% A, tol = 1e-15),
+ all.equal(X0, X %*% t(A), tol = 1e-15)))
> ##--> here, don't need to keep X0 (nor A, since that is Qrot(p))
>
> ## for n = 100, you don't see "it", but may guess .. :
> n <- 100
> pairs(r <- rbwheel(n,6))
> n1 <- attr(r,"n1") ; pairs(r, col=1+((1:n) > n1))
>
> ## for n = 500, you *do* see it :
> n <- 500
> pairs(r <- rbwheel(n,6))
> ## show explicitly
> n1 <- attr(r,"n1") ; pairs(r, col=1+((1:n) > n1))
>
> ## but increasing sig2 does help:
> pairs(r <- rbwheel(n,6, sig2 = .2))
>
> ## show explicitly
> n1 <- attr(r,"n1") ; pairs(r, col=1+((1:n) > n1))
>
> set.seed(12)
> pairs(X <- rbwheel(n, 7, spherize=TRUE))
> colSums(X) # already centered
[1] -3.129441e-15 6.151329e-15 -8.555656e-15 2.352285e-15 2.428613e-16
[6] 3.505876e-15 -3.537101e-15
>
> if(require("ICS")) {
+ # ICS: Compare M-estimate [Max.Lik. of t_{df = 2}] with high-breakdown :
+ stopifnot(require("MASS"))
+ X.paM <- ics(X, S1 = cov, S2 = function(.) cov.trob(., nu=2)$cov, stdKurt = FALSE)
+ X.paM.<- ics(X, S1 = cov, S2 = function(.) tM(., df=2)$V, stdKurt = FALSE)
+ X.paR <- ics(X, S1 = cov, S2 = function(.) covMcd(.)$cov, stdKurt = FALSE)
+ plot(X.paM) # not at all clear
+ plot(X.paM.)# ditto
+ plot(X.paR)# very clear
+ }
Loading required package: ICS
Loading required package: mvtnorm
Loading required package: MASS
Error in covMcd(.) : could not find function "covMcd"
Calls: ics -> do.call -> S2
Execution halted
Flavors: r-release-osx-x86_64, r-oldrel-osx-x86_64