[Robast-commits] r450 - pkg/RobLox/tests/Examples

[noreply at r-forge.r-project.org]

Author: stamats
Date: 2011-01-06 11:41:06 +0100 (Thu, 06 Jan 2011)
New Revision: 450

Modified:
   pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save
Log:
forgot to update RobLox-Ex.Rout.save

Modified: pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save
===================================================================
--- pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save	2011-01-06 10:31:51 UTC (rev 449)
+++ pkg/RobLox/tests/Examples/RobLox-Ex.Rout.save	2011-01-06 10:41:06 UTC (rev 450)
@@ -22,14 +22,23 @@
 > source(file.path(R.home("share"), "R", "examples-header.R"))
 > options(warn = 1)
 > library('RobLox')
-Loading required package: distrMod
+Loading required package: lattice
+Loading required package: RColorBrewer
+Loading required package: Biobase
+
+Welcome to Bioconductor
+
+  Vignettes contain introductory material. To view, type
+  'openVignette()'. To cite Bioconductor, see
+  'citation("Biobase")' and for packages 'citation(pkgname)'.
+
+Loading required package: distr
 Loading required package: startupmsg
 :startupmsg>  Utilities for start-up messages (version 0.7.1)
 :startupmsg> 
 :startupmsg>  For more information see ?"startupmsg",
 :startupmsg>  NEWS("startupmsg")
 
-Loading required package: distr
 Loading required package: sfsmisc
 Loading required package: SweaveListingUtils
 :SweaveListingUtils>  Utilities for Sweave together with
@@ -84,8 +93,16 @@
 
     df, qqplot, sd
 
+Loading required package: distrMod
 Loading required package: distrEx
 Loading required package: evd
+
+Attaching package: 'evd'
+
+The following object(s) are masked from 'package:lattice':
+
+    qq
+
 Loading required package: actuar
 
 Attaching package: 'actuar'
@@ -325,14 +342,6 @@
 > ind <- rbinom(200, size=1, prob=0.05) 
 > X <- matrix(rnorm(200, mean=ind*3, sd=(1-ind) + ind*9), nrow = 2)
 > rowRoblox(X)
-Loading required package: Biobase
-
-Welcome to Bioconductor
-
-  Vignettes contain introductory material. To view, type
-  'openVignette()'. To cite Bioconductor, see
-  'citation("Biobase")' and for packages 'citation(pkgname)'.
-
 Evaluations of Optimally robust estimate:
 -----------------------------------------
 An object of class “Estimate” 
@@ -356,9 +365,6 @@
 > 
 > 
 > cleanEx()
-
-detaching ‘package:Biobase’
-
 > nameEx("finiteSampleCorrection")
 > ### * finiteSampleCorrection
 > 
@@ -1774,14 +1780,6 @@
 > ind <- rbinom(200, size=1, prob=0.05) 
 > X <- matrix(rnorm(200, mean=ind*3, sd=(1-ind) + ind*9), nrow = 2)
 > rowRoblox(X)
-Loading required package: Biobase
-
-Welcome to Bioconductor
-
-  Vignettes contain introductory material. To view, type
-  'openVignette()'. To cite Bioconductor, see
-  'citation("Biobase")' and for packages 'citation(pkgname)'.
-
 Evaluations of Optimally robust estimate:
 -----------------------------------------
 An object of class “Estimate” 
@@ -2111,9 +2109,6 @@
 > 
 > 
 > cleanEx()
-
-detaching ‘package:Biobase’
-
 > nameEx("rsOptIC")
 > ### * rsOptIC
 > 
@@ -2157,13 +2152,56 @@
 > 
 > 
 > 
+> cleanEx()
+> nameEx("showdown")
+> ### * showdown
+> 
+> flush(stderr()); flush(stdout())
+> 
+> ### Name: showdown
+> ### Title: Estimator Showdown by Monte-Carlo Study.
+> ### Aliases: showdown
+> ### Keywords: robust
+> 
+> ### ** Examples
+> 
+> library(MASS)
+> ## compare with Huber's Proposal 2
+> showdown(n = 20, M = 100, eps = 0.02, contD = Norm(mean = 3, sd = 3), 
++          estfun = function(x){ unlist(hubers(x)) },
++          plot1 = TRUE, plot2 = TRUE, plot3 = TRUE)
+dev.new(): using pdf(file="Rplots1.pdf")
+                               ML      Med Competitor       rmx
+n x empMSE (loc)         1.151729 1.571687  1.0217401 1.0412422
+relMSE (loc)             1.106111 1.509434  0.9812703 1.0000000
+n x empMSE (scale)       1.086260 1.493500  0.7204555 0.6113725
+relMSE (scale)           1.776757 2.442864  1.1784231 1.0000000
+n x empMSE (loc + scale) 2.237989 3.065187  1.7421956 1.6526147
+relMSE (loc + scale)     1.354211 1.854750  1.0542055 1.0000000
+> 
+> ## compare with Huber M estimator with MAD scale
+> showdown(n = 20, M = 100, eps = 0.02, contD = Norm(mean = 3, sd = 3), 
++          estfun = function(x){ unlist(huber(x)) },
++          plot1 = TRUE, plot2 = TRUE, plot3 = TRUE)
+dev.new(): using pdf(file="Rplots2.pdf")
+dev.new(): using pdf(file="Rplots3.pdf")
+                               ML      Med Competitor       rmx
+n x empMSE (loc)         1.151729 1.571687   1.091427 1.0412422
+relMSE (loc)             1.106111 1.509434   1.048197 1.0000000
+n x empMSE (scale)       1.086260 1.493500   1.493499 0.6113725
+relMSE (scale)           1.776757 2.442864   2.442862 1.0000000
+n x empMSE (loc + scale) 2.237989 3.065187   2.584925 1.6526147
+relMSE (loc + scale)     1.354211 1.854750   1.564143 1.0000000
+> 
+> 
+> 
 > ### * <FOOTER>
 > ###
 > cat("Time elapsed: ", proc.time() - get("ptime", pos = 'CheckExEnv'),"\n")
-Time elapsed:  63.3 0.31 63.862 0 0 
+Time elapsed:  64.36 0.06 64.627 0 0 
 > grDevices::dev.off()
-null device 
-          1 
+postscript 
+         2 
 > ###
 > ### Local variables: ***
 > ### mode: outline-minor ***