[Robast-commits] r820 - branches/robast-1.0/pkg/RobLox/tests/Examples
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Sun May 3 15:34:28 CEST 2015
Author: stamats
Date: 2015-05-03 15:34:27 +0200 (Sun, 03 May 2015)
New Revision: 820
Modified:
branches/robast-1.0/pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save
Log:
updated Rout.save file
Modified: branches/robast-1.0/pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save
===================================================================
--- branches/robast-1.0/pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save 2015-05-03 13:29:47 UTC (rev 819)
+++ branches/robast-1.0/pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save 2015-05-03 13:34:27 UTC (rev 820)
@@ -1,7 +1,6 @@
-R Under development (unstable) (2013-02-06 r61845) -- "Unsuffered Consequences"
-Copyright (C) 2013 The R Foundation for Statistical Computing
-ISBN 3-900051-07-0
+R Under development (unstable) (2015-05-02 r68310) -- "Unsuffered Consequences"
+Copyright (C) 2015 The R Foundation for Statistical Computing
Platform: x86_64-unknown-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
@@ -22,49 +21,19 @@
> source(file.path(R.home("share"), "R", "examples-header.R"))
> options(warn = 1)
> library('RobLox')
-Loading required package: lattice
-Loading required package: RColorBrewer
-Loading required package: Biobase
-Loading required package: BiocGenerics
-Loading required package: parallel
-
-Attaching package: ‘BiocGenerics’
-
-The following objects are masked from ‘package:parallel’:
-
- clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
- clusterExport, clusterMap, parApply, parCapply, parLapply,
- parLapplyLB, parRapply, parSapply, parSapplyLB
-
-The following object is masked from ‘package:stats’:
-
- xtabs
-
-The following objects are masked from ‘package:base’:
-
- Filter, Find, Map, Position, Reduce, anyDuplicated, as.data.frame,
- cbind, colnames, duplicated, eval, get, intersect, lapply, mapply,
- match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, rank,
- rbind, rep.int, rownames, sapply, setdiff, sort, table, tapply,
- union, unique, unlist
-
-Welcome to Bioconductor
-
- Vignettes contain introductory material; view with
- 'browseVignettes()'. To cite Bioconductor, see
- 'citation("Biobase")', and for packages 'citation("pkgname")'.
-
+Loading required package: distrMod
Loading required package: distr
Loading required package: startupmsg
-:startupmsg> Utilities for start-up messages (version 0.8)
+:startupmsg> Utilities for Start-Up Messages (version 0.9.1)
:startupmsg>
:startupmsg> For more information see ?"startupmsg",
:startupmsg> NEWS("startupmsg")
Loading required package: sfsmisc
Loading required package: SweaveListingUtils
-:SweaveListingUtils> Utilities for Sweave together with
-:SweaveListingUtils> TeX listings package (version 0.6)
+:SweaveListingUtils> Utilities for Sweave Together with
+:SweaveListingUtils> TeX 'listings' Package (version
+:SweaveListingUtils> 0.7)
:SweaveListingUtils>
:SweaveListingUtils> NOTE: Support for this package
:SweaveListingUtils> will stop soon.
@@ -92,12 +61,12 @@
Attaching package: ‘SweaveListingUtils’
-The following object is masked from ‘package:base’:
+The following objects are masked from ‘package:base’:
library, require
-:distr> Object oriented implementation of distributions (version
-:distr> 2.4)
+:distr> Object Oriented Implementation of Distributions (version
+:distr> 2.6)
:distr>
:distr> Attention: Arithmetics on distribution objects are
:distr> understood as operations on corresponding random variables
@@ -118,13 +87,12 @@
Attaching package: ‘distr’
-The following object is masked from ‘package:stats’:
+The following objects are masked from ‘package:stats’:
df, qqplot, sd
-Loading required package: distrMod
Loading required package: distrEx
-:distrEx> Extensions of package distr (version 2.4)
+:distrEx> Extensions of Package 'distr' (version 2.6)
:distrEx>
:distrEx> Note: Packages "e1071", "moments", "fBasics" should be
:distrEx> attached /before/ package "distrEx". See
@@ -143,12 +111,12 @@
Attaching package: ‘distrEx’
-The following object is masked from ‘package:stats’:
+The following objects are masked from ‘package:stats’:
IQR, mad, median, var
Loading required package: RandVar
-:RandVar> Implementation of random variables (version 0.9)
+:RandVar> Implementation of Random Variables (version 1.0)
:RandVar>
:RandVar> For more information see ?"RandVar", NEWS("RandVar"), as
:RandVar> well as
@@ -158,8 +126,8 @@
Loading required package: MASS
Loading required package: stats4
-:distrMod> Object oriented implementation of probability models
-:distrMod> (version 2.4)
+:distrMod> Object Oriented Implementation of Probability Models
+:distrMod> (version 2.6)
:distrMod>
:distrMod> Some functions from pkg's 'base' and 'stats' are
:distrMod> intentionally masked ---see distrModMASK().
@@ -195,11 +163,9 @@
Loading required package: RobAStBase
Loading required package: rrcov
Loading required package: robustbase
-Loading required package: pcaPP
-Loading required package: mvtnorm
-Scalable Robust Estimators with High Breakdown Point (version 1.3-02)
+Scalable Robust Estimators with High Breakdown Point (version 1.3-8)
-:RobAStBase> Robust Asymptotic Statistics (version 0.9)
+:RobAStBase> Robust Asymptotic Statistics (version 1.0)
:RobAStBase>
:RobAStBase> Some functions from pkg's 'stats' and 'graphics'
:RobAStBase> are intentionally masked ---see RobAStBaseMASK().
@@ -220,23 +186,12 @@
>
> base::assign(".oldSearch", base::search(), pos = 'CheckExEnv')
-> base::assign(".ExTimings", "RobLox-Ex.timings", pos = 'CheckExEnv')
-> base::cat("name\tuser\tsystem\telapsed\n", file=base::get(".ExTimings", pos = 'CheckExEnv'))
-> base::assign(".format_ptime",
-+ function(x) {
-+ if(!is.na(x[4L])) x[1L] <- x[1L] + x[4L]
-+ if(!is.na(x[5L])) x[2L] <- x[2L] + x[5L]
-+ format(x[1L:3L])
-+ },
-+ pos = 'CheckExEnv')
->
> cleanEx()
> nameEx("0RobLox-package")
> ### * 0RobLox-package
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: RobLox-package
> ### Title: Optimally robust influence curves and estimators for location
> ### and scale
@@ -369,16 +324,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("0RobLox-package", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("finiteSampleCorrection")
> ### * finiteSampleCorrection
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: finiteSampleCorrection
> ### Title: Function to compute finite-sample corrected radii
> ### Aliases: finiteSampleCorrection
@@ -395,16 +346,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("finiteSampleCorrection", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlOptIC")
> ### * rlOptIC
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlOptIC
> ### Title: Computation of the optimally robust IC for AL estimators
> ### Aliases: rlOptIC
@@ -419,6 +366,9 @@
precision of Fisher consistency:
mean
mean -7.580384e-06
+precision of Fisher consistency - relativ error [%]:
+ mean
+mean -0.0007580384
maximum deviation
7.580384e-06
> distrExOptions("ErelativeTolerance" = .Machine$double.eps^0.25) # default
@@ -444,16 +394,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlOptIC", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.AL")
> ### * rlsOptIC.AL
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.AL
> ### Title: Computation of the optimally robust IC for AL estimators
> ### Aliases: rlsOptIC.AL
@@ -462,17 +408,21 @@
> ### ** Examples
>
> IC1 <- rlsOptIC.AL(r = 0.1, check = TRUE)
-Fisher consistency of eta.loc: -1.743783e-10
-centering of eta.sc: -3.904033e-10
-Fisher consistency of eta.sc: 2.926179e-09
+Fisher consistency of eta.loc: -1.743714e-10
+centering of eta.sc: -3.903789e-10
+Fisher consistency of eta.sc: 2.926104e-09
MSE equation: 1.207368e-14
> distrExOptions("ErelativeTolerance" = 1e-12)
> checkIC(IC1)
-precision of centering: 0 -6.039278e-07
+precision of centering: 0 -6.039298e-07
precision of Fisher consistency:
- mean sd
-mean -1.10248e-06 0.000000e+00
-sd 0.00000e+00 -1.685676e-05
+ mean sd
+mean -1.102483e-06 0.000000e+00
+sd 0.000000e+00 -1.685676e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0001102483 NaN
+sd NaN -0.001685676
maximum deviation
1.685676e-05
> distrExOptions("ErelativeTolerance" = .Machine$double.eps^0.25) # default
@@ -542,16 +492,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.AL", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.An1")
> ### * rlsOptIC.An1
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.An1
> ### Title: Computation of the optimally robust IC for An1 estimators
> ### Aliases: rlsOptIC.An1
@@ -566,6 +512,10 @@
mean sd
mean 2.530605e-08 0.000000e+00
sd 0.000000e+00 3.875467e-06
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean 2.530605e-06 NaN
+sd NaN 0.0003875467
maximum deviation
3.875467e-06
> Risks(IC1)
@@ -590,16 +540,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.An1", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.An2")
> ### * rlsOptIC.An2
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.An2
> ### Title: Computation of the optimally robust IC for An2 estimators
> ### Aliases: rlsOptIC.An2
@@ -614,6 +560,10 @@
mean sd
mean -4.006394e-07 0.000000e+00
sd 0.000000e+00 -1.496408e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -4.006394e-05 NaN
+sd NaN -0.001496408
maximum deviation
0.0001066371
> Risks(IC1)
@@ -637,16 +587,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.An2", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.AnMad")
> ### * rlsOptIC.AnMad
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.AnMad
> ### Title: Computation of the optimally robust IC for AnMad estimators
> ### Aliases: rlsOptIC.AnMad
@@ -661,6 +607,10 @@
mean sd
mean -4.044557e-07 0.000000e+00
sd 0.000000e+00 -6.821323e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -4.044557e-05 NaN
+sd NaN -0.006821323
maximum deviation
6.821323e-05
> Risks(IC1)
@@ -684,16 +634,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.AnMad", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.BM")
> ### * rlsOptIC.BM
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.BM
> ### Title: Computation of the optimally robust IC for BM estimators
> ### Aliases: rlsOptIC.BM
@@ -707,9 +653,13 @@
precision of Fisher consistency:
mean sd
mean -2.355426e-06 0.000000e+00
-sd 0.000000e+00 -1.737383e-05
+sd 0.000000e+00 -1.732518e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0002355426 NaN
+sd NaN -0.001732518
maximum deviation
- 1.737383e-05
+ 1.732518e-05
> Risks(IC1)
$asMSE
[1] 1.659506
@@ -734,16 +684,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.BM", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Ha3")
> ### * rlsOptIC.Ha3
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Ha3
> ### Title: Computation of the optimally robust IC for Ha3 estimators
> ### Aliases: rlsOptIC.Ha3
@@ -753,11 +699,15 @@
>
> IC1 <- rlsOptIC.Ha3(r = 0.1)
> checkIC(IC1)
-precision of centering: 0 -2.659214e-07
+precision of centering: 0 -3.191376e-07
precision of Fisher consistency:
mean sd
mean -8.538336e-07 0.000000e+00
sd 0.000000e+00 -9.691769e-06
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -8.538336e-05 NaN
+sd NaN -0.0009691769
maximum deviation
9.691769e-06
> Risks(IC1)
@@ -782,16 +732,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Ha3", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Ha4")
> ### * rlsOptIC.Ha4
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Ha4
> ### Title: Computation of the optimally robust IC for Ha4 estimators
> ### Aliases: rlsOptIC.Ha4
@@ -806,6 +752,10 @@
mean sd
mean -2.265319e-06 0.000000e+00
sd 0.000000e+00 -1.498731e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0002265319 NaN
+sd NaN -0.001498731
maximum deviation
0.0001069318
> Risks(IC1)
@@ -829,16 +779,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Ha4", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.HaMad")
> ### * rlsOptIC.HaMad
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.HaMad
> ### Title: Computation of the optimally robust IC for HuMad estimators
> ### Aliases: rlsOptIC.HaMad
@@ -853,6 +799,10 @@
mean sd
mean -2.041724e-06 0.000000e+00
sd 0.000000e+00 -6.821323e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0002041724 NaN
+sd NaN -0.006821323
maximum deviation
6.821323e-05
> Risks(IC1)
@@ -876,16 +826,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.HaMad", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Hu1")
> ### * rlsOptIC.Hu1
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Hu1
> ### Title: Computation of the optimally robust IC for Hu1 estimators
> ### Aliases: rlsOptIC.Hu1
@@ -900,6 +846,10 @@
mean sd
mean -3.746843e-06 0.000000e+00
sd 0.000000e+00 -1.350353e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0003746843 NaN
+sd NaN -0.001350353
maximum deviation
1.350353e-05
> Risks(IC1)
@@ -923,16 +873,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Hu1", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Hu2")
> ### * rlsOptIC.Hu2
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Hu2
> ### Title: Computation of the optimally robust IC for Hu2 estimators
> ### Aliases: rlsOptIC.Hu2
@@ -947,6 +893,10 @@
mean sd
mean -2.563384e-06 0.000000e+00
sd 0.000000e+00 -1.321969e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0002563384 NaN
+sd NaN -0.001321969
maximum deviation
1.321969e-05
> Risks(IC1)
@@ -970,16 +920,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Hu2", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Hu2a")
> ### * rlsOptIC.Hu2a
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Hu2a
> ### Title: Computation of the optimally robust IC for Hu2a estimators
> ### Aliases: rlsOptIC.Hu2a
@@ -994,6 +940,10 @@
mean sd
mean -3.15714e-06 0.000000e+00
sd 0.00000e+00 -1.321476e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.000315714 NaN
+sd NaN -0.001321476
maximum deviation
1.321476e-05
> Risks(IC1)
@@ -1017,16 +967,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Hu2a", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Hu3")
> ### * rlsOptIC.Hu3
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Hu3
> ### Title: Computation of the optimally robust IC for Hu3 estimators
> ### Aliases: rlsOptIC.Hu3
@@ -1041,6 +987,10 @@
mean sd
mean -2.132286e-06 0.000000e+00
sd 0.000000e+00 -1.574855e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0002132286 NaN
+sd NaN -0.001574855
maximum deviation
8.525774e-05
> Risks(IC1)
@@ -1064,16 +1014,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Hu3", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.HuMad")
> ### * rlsOptIC.HuMad
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.HuMad
> ### Title: Computation of the optimally robust IC for HuMad estimators
> ### Aliases: rlsOptIC.HuMad
@@ -1088,6 +1034,10 @@
mean sd
mean -1.987572e-06 0.000000e+00
sd 0.000000e+00 -6.821323e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0001987572 NaN
+sd NaN -0.006821323
maximum deviation
6.821323e-05
> Risks(IC1)
@@ -1111,16 +1061,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.HuMad", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.M")
> ### * rlsOptIC.M
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.M
> ### Title: Computation of the optimally robust IC for M estimators
> ### Aliases: rlsOptIC.M
@@ -1131,14 +1077,18 @@
> IC1 <- rlsOptIC.M(r = 0.1, check = TRUE)
constraint a1: 4.682285e-09
constraint a3: -1.033529e-08
-MSE equation: -3.383516e-12
+MSE equation: -3.383432e-12
> distrExOptions("ErelativeTolerance" = 1e-12)
> checkIC(IC1, NormLocationScaleFamily())
-precision of centering: -1.418367e-20 9.579943e-06
+precision of centering: 4.510461e-21 9.579943e-06
precision of Fisher consistency:
mean sd
-mean 1.646253e-05 -8.358660e-22
-sd 1.313976e-21 2.094876e-05
+mean 1.646253e-05 -5.609939e-21
+sd 5.208550e-21 2.094876e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean 0.001646253 -Inf
+sd Inf 0.002094876
maximum deviation
2.094876e-05
> distrExOptions("ErelativeTolerance" = .Machine$double.eps^0.25)
@@ -1163,16 +1113,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.M", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.MM2")
> ### * rlsOptIC.MM2
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.MM2
> ### Title: Computation of the optimally robust IC for MM2 estimators
> ### Aliases: rlsOptIC.MM2
@@ -1187,6 +1133,10 @@
mean sd
mean -2.225204e-06 0.000000e+00
sd 0.000000e+00 -1.566028e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -0.0002225204 NaN
+sd NaN -0.001566028
maximum deviation
1.566028e-05
> Risks(IC1)
@@ -1210,16 +1160,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.MM2", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Tu1")
> ### * rlsOptIC.Tu1
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Tu1
> ### Title: Computation of the optimally robust IC for Tu1 estimators
> ### Aliases: rlsOptIC.Tu1
@@ -1234,6 +1180,10 @@
mean sd
mean -3.649066e-08 0.000000e+00
sd 0.000000e+00 2.888236e-06
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -3.649066e-06 NaN
+sd NaN 0.0002888236
maximum deviation
2.888236e-06
> Risks(IC1)
@@ -1257,16 +1207,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Tu1", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.Tu2")
> ### * rlsOptIC.Tu2
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.Tu2
> ### Title: Computation of the optimally robust IC for Tu2 estimators
> ### Aliases: rlsOptIC.Tu2
@@ -1276,11 +1222,15 @@
>
> IC1 <- rlsOptIC.Tu2(r = 0.1)
> checkIC(IC1)
-precision of centering: 0 -3.919236e-07
+precision of centering: 0 -1.429842e-07
precision of Fisher consistency:
mean sd
mean -5.447399e-07 0.000000e+00
sd 0.000000e+00 -1.437669e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -5.447399e-05 NaN
+sd NaN -0.001437669
maximum deviation
1.437669e-05
> Risks(IC1)
@@ -1304,16 +1254,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.Tu2", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rlsOptIC.TuMad")
> ### * rlsOptIC.TuMad
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rlsOptIC.TuMad
> ### Title: Computation of the optimally robust IC for TuMad estimators
> ### Aliases: rlsOptIC.TuMad
@@ -1328,6 +1274,10 @@
mean sd
mean -5.410012e-07 0.000000e+00
sd 0.000000e+00 -6.821323e-05
+precision of Fisher consistency - relativ error [%]:
+ mean sd
+mean -5.410012e-05 NaN
+sd NaN -0.006821323
maximum deviation
6.821323e-05
> Risks(IC1)
@@ -1351,16 +1301,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rlsOptIC.TuMad", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("roblox")
> ### * roblox
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: roblox
> ### Title: Optimally robust estimator for location and/or scale
> ### Aliases: roblox
@@ -1434,10 +1380,10 @@
samplesize: 100
estimate:
mean sd
--0.03980472 0.90817692
+-0.08687757 0.83485767
Criterion:
Kolmogorov distance
- 0.0568421
+ 1e+20
>
> # optimally robust (amount of gross errors known)
> roblox(x, eps = 0.05, initial.est = estimate(ks.est))
@@ -1449,12 +1395,12 @@
samplesize: 100
estimate:
mean sd
- -0.11586692 0.90365737
- ( 0.10307092) ( 0.07651628)
+ -0.11229870 0.88801808
+ ( 0.10128710) ( 0.07519204)
asymptotic (co)variance (multiplied with samplesize):
[,1] [,2]
-[1,] 1.062361 0.0000000
-[2,] 0.000000 0.5854741
+[1,] 1.025908 0.0000000
+[2,] 0.000000 0.5653843
Infos:
method
[1,] "roblox"
@@ -1462,7 +1408,7 @@
[1,] "finite-sample corrected optimally robust estimate for contamination 'eps' = 0.05 and 'asMSE'"
asymptotic bias:
sd
-1.035576
+1.017654
steps:
[1] 1
>
@@ -1476,10 +1422,10 @@
samplesize: 100
estimate:
mean sd
--0.06995363 0.89265661
+-0.08687757 0.83485767
Criterion:
CvM distance
- 0.02436665
+ 1e+20
>
> # optimally robust (amount of gross errors known)
> roblox(x, eps = 0.05, initial.est = estimate(CvM.est))
@@ -1491,12 +1437,12 @@
samplesize: 100
estimate:
mean sd
- -0.11421421 0.89779415
- ( 0.10240216) ( 0.07601982)
+ -0.11229870 0.88801808
+ ( 0.10128710) ( 0.07519204)
asymptotic (co)variance (multiplied with samplesize):
- [,1] [,2]
-[1,] 1.04862 0.0000000
-[2,] 0.00000 0.5779013
+ [,1] [,2]
+[1,] 1.025908 0.0000000
+[2,] 0.000000 0.5653843
Infos:
method
[1,] "roblox"
@@ -1504,22 +1450,18 @@
[1,] "finite-sample corrected optimally robust estimate for contamination 'eps' = 0.05 and 'asMSE'"
asymptotic bias:
sd
-1.028857
+1.017654
steps:
[1] 1
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("roblox", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rowRoblox")
> ### * rowRoblox
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rowRoblox and colRoblox
> ### Title: Optimally robust estimation for location and/or scale
> ### Aliases: rowRoblox colRoblox
@@ -1858,16 +1800,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rowRoblox", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("rsOptIC")
> ### * rsOptIC
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: rsOptIC
> ### Title: Computation of the optimally robust IC for AL estimators
> ### Aliases: rsOptIC
@@ -1882,6 +1820,9 @@
precision of Fisher consistency:
scale
scale -1.394372e-05
+precision of Fisher consistency - relativ error [%]:
+ scale
+scale -0.001394372
maximum deviation
1.394372e-05
> distrExOptions("ErelativeTolerance" = .Machine$double.eps^0.25) # default
@@ -1907,16 +1848,12 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("rsOptIC", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> cleanEx()
> nameEx("showdown")
> ### * showdown
>
> flush(stderr()); flush(stdout())
>
-> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: showdown
> ### Title: Estimator Showdown by Monte-Carlo Study.
> ### Aliases: showdown
@@ -1954,13 +1891,11 @@
>
>
>
->
-> base::assign(".dptime", (proc.time() - get(".ptime", pos = "CheckExEnv")), pos = "CheckExEnv")
-> base::cat("showdown", base::get(".format_ptime", pos = 'CheckExEnv')(get(".dptime", pos = "CheckExEnv")), "\n", file=base::get(".ExTimings", pos = 'CheckExEnv'), append=TRUE, sep="\t")
> ### * <FOOTER>
> ###
+> options(digits = 7L)
> base::cat("Time elapsed: ", proc.time() - base::get("ptime", pos = 'CheckExEnv'),"\n")
-Time elapsed: 72.904 0.34 73.52 0 0
+Time elapsed: 19.658 0.08 19.723 0 0
> grDevices::dev.off()
pdf
2
More information about the Robast-commits
mailing list