[Robast-commits] r130 - in branches/robast-0.6/pkg/RobLox: . R inst/tests man

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Mon Jul 28 14:16:05 CEST 2008


Author: stamats
Date: 2008-07-28 14:16:05 +0200 (Mon, 28 Jul 2008)
New Revision: 130

Modified:
   branches/robast-0.6/pkg/RobLox/DESCRIPTION
   branches/robast-0.6/pkg/RobLox/NAMESPACE
   branches/robast-0.6/pkg/RobLox/R/roblox.R
   branches/robast-0.6/pkg/RobLox/R/rowRoblox.R
   branches/robast-0.6/pkg/RobLox/inst/tests/tests.R
   branches/robast-0.6/pkg/RobLox/man/rlOptIC.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AL.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An1.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An2.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AnMad.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.BM.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha3.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha4.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HaMad.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu1.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2a.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu3.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HuMad.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.M.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.MM2.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu1.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu2.Rd
   branches/robast-0.6/pkg/RobLox/man/rlsOptIC.TuMad.Rd
   branches/robast-0.6/pkg/RobLox/man/roblox.Rd
   branches/robast-0.6/pkg/RobLox/man/rowRoblox.Rd
   branches/robast-0.6/pkg/RobLox/man/rsOptIC.Rd
Log:
adaption to new implementation using kStepEstimate class; some minor corrections of Rd-files

Modified: branches/robast-0.6/pkg/RobLox/DESCRIPTION
===================================================================
--- branches/robast-0.6/pkg/RobLox/DESCRIPTION	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/DESCRIPTION	2008-07-28 12:16:05 UTC (rev 130)
@@ -1,6 +1,6 @@
 Package: RobLox
 Version: 0.6.0
-Date: 2008-07-21
+Date: 2008-07-28
 Title: Optimally robust influence curves for location and scale
 Description: functions for the determination of optimally 
     robust influence curves in case of normal
@@ -12,4 +12,3 @@
 LazyLoad: yes
 License: GPL version 2 or later
 URL: http://robast.r-forge.r-project.org/
-Packaged: Mon Aug  6 08:07:25 2007; kohl

Modified: branches/robast-0.6/pkg/RobLox/NAMESPACE
===================================================================
--- branches/robast-0.6/pkg/RobLox/NAMESPACE	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/NAMESPACE	2008-07-28 12:16:05 UTC (rev 130)
@@ -21,4 +21,3 @@
        roblox,
        rowRoblox,
        colRoblox)
-S3method(print, ALEstimate)

Modified: branches/robast-0.6/pkg/RobLox/R/roblox.R
===================================================================
--- branches/robast-0.6/pkg/RobLox/R/roblox.R	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/R/roblox.R	2008-07-28 12:16:05 UTC (rev 130)
@@ -255,11 +255,11 @@
                                       risk = list(asMSE = mse, asBias = b, asCov = mse - r^2*b^2), 
                                       info = c("roblox", "optimally robust IC for AL estimators and 'asMSE'"),
                                       w = w, biastype = symmetricBias(), normtype = NormType()))
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix, "optIC" = IC1), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = IC1, Infos = Info.matrix))
             }else
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
         }else{
             sqrtn <- sqrt(length(x))
             rlo <- sqrtn*eps.lower
@@ -323,11 +323,11 @@
                                       paste("least favorable contamination: ", round(r/sqrtn, 3), sep = ""),
                                       paste("maximum MSE-inefficiency: ", round(ineff, 3), sep = "")), 
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix, "optIC" = IC1),
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = IC1, Infos = Info.matrix))
             }else
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
         }
     }else{
         if(missing(mean)){
@@ -362,7 +362,7 @@
                     w <- new("HampelWeight")
                     clip(w) <- b
                     cent(w) <- 0
-                    stand(w) <- A
+                    stand(w) <- as.matrix(A)
                     weight(w) <- getweight(w, neighbor = ContNeighborhood(radius = r), 
                                        biastype = symmetricBias(), 
                                        normW = NormType())
@@ -372,11 +372,11 @@
                                           risk = list(asMSE = A, asBias = b, asCov = b^2), 
                                           info = c("roblox", "optimally robust IC for AL estimators and 'asMSE'"),
                                           w = w, biastype = symmetricBias(), normtype = NormType()))
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix, "optIC" = IC1),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = IC1, Infos = Info.matrix))
                 }else
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }else{
                 sqrtn <- sqrt(length(x))
                 rlo <- sqrtn*eps.lower
@@ -416,7 +416,7 @@
                     w <- new("HampelWeight")
                     clip(w) <- b
                     cent(w) <- 0
-                    stand(w) <- A
+                    stand(w) <- as.matrix(A)
                     weight(w) <- getweight(w, neighbor = ContNeighborhood(radius = r), 
                                        biastype = symmetricBias(), 
                                        normW = NormType())
@@ -432,11 +432,11 @@
                                  paste("least favorable contamination: ", round(r/sqrtn, 3), sep = ""),
                                  paste("maximum MSE-inefficiency: ", round(ineff, 3), sep = "")), 
                                  ncol = 2, dimnames = list(NULL, c("method", "message")))
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix, "optIC" = IC1),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = IC1, Infos = Info.matrix))
                 }else
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }
         }
         if(missing(sd)){
@@ -476,7 +476,7 @@
                     w <- new("HampelWeight")
                     clip(w) <- b
                     cent(w) <- a/A
-                    stand(w) <- A
+                    stand(w) <- as.matrix(A)
                     weight(w) <- getweight(w, neighbor = ContNeighborhood(radius = r), 
                                        biastype = symmetricBias(), 
                                        normW = NormType())
@@ -486,11 +486,11 @@
                                           risk = list(asMSE = A, asBias = b, asCov = b^2), 
                                           info = c("roblox", "optimally robust IC for AL estimators and 'asMSE'"),
                                           w = w, biastype = symmetricBias(), normtype = NormType()))
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix, "optIC" = IC1),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = IC1, Infos = Info.matrix))
                 }else
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }else{
                 sqrtn <- sqrt(length(x))
                 rlo <- sqrtn*eps.lower
@@ -532,7 +532,7 @@
                     w <- new("HampelWeight")
                     clip(w) <- b
                     cent(w) <- a/A
-                    stand(w) <- A
+                    stand(w) <- as.matrix(A)
                     weight(w) <- getweight(w, neighbor = ContNeighborhood(radius = r), 
                                        biastype = symmetricBias(), 
                                        normW = NormType())
@@ -548,18 +548,12 @@
                                  paste("least favorable contamination: ", round(r/sqrtn, 3), sep = ""),
                                  paste("maximum MSE-inefficiency: ", round(ineff, 3), sep = "")), 
                                  ncol = 2, dimnames = list(NULL, c("method", "message")))
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix, "optIC" = IC1),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = IC1, Infos = Info.matrix))
                 }else
-                    return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix),
-                                 class = c("ALEstimate", "Estimate")))
+                    return(new("kStepEstimate", name = "Optimally robust estimate",
+                               estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }
         }
     }
 }
-print.ALEstimate <- function(x, digits = getOption("digits"), ...){
-  print(x$estimate)
-  if(!is.null(x$Infos)){
-    print(x$Infos)
-  }
-}

Modified: branches/robast-0.6/pkg/RobLox/R/rowRoblox.R
===================================================================
--- branches/robast-0.6/pkg/RobLox/R/rowRoblox.R	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/R/rowRoblox.R	2008-07-28 12:16:05 UTC (rev 130)
@@ -153,8 +153,8 @@
                                     paste("optimally robust estimates for contamination 'eps' =", round(eps, 3),
                                           "and 'asMSE'")),
                                   ncol = 2, dimnames = list(NULL, c("method", "message")))
-            return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                             class = c("ALEstimate", "Estimate")))
+            return(new("kStepEstimate", name = "Optimally robust estimate",
+                       estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
         }else{
             sqrtn <- sqrt(ncol(x))
             rlo <- sqrtn*eps.lower
@@ -199,8 +199,8 @@
                                   paste("least favorable contamination: ", round(r/sqrtn, 3), sep = ""),
                                   paste("maximum MSE-inefficiency: ", round(ineff[1], 3), sep = "")), 
                                   ncol = 2, dimnames = list(NULL, c("method", "message")))
-            return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                             class = c("ALEstimate", "Estimate")))
+            return(new("kStepEstimate", name = "Optimally robust estimate",
+                       estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
         }
     }else{
         if(missing(mean)){
@@ -235,8 +235,8 @@
                                         paste("optimally robust estimates for contamination 'eps' =", round(eps, 3),
                                               "and 'asMSE'")),
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }else{
                 sqrtn <- sqrt(ncol(x))
                 rlo <- sqrtn*eps.lower
@@ -272,8 +272,8 @@
                                       paste("least favorable contamination: ", round(r/sqrtn, 3), sep = ""),
                                       paste("maximum MSE-inefficiency: ", round(ineff[1], 3), sep = "")), 
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }
         }
         if(missing(sd)){
@@ -314,8 +314,8 @@
                                         paste("optimally robust estimates for contamination 'eps' =", round(eps, 3),
                                               "and 'asMSE'")),
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }else{
                 sqrtn <- sqrt(ncol(x))
                 rlo <- sqrtn*eps.lower
@@ -353,8 +353,8 @@
                                       paste("least favorable contamination: ", round(r/sqrtn, 3), sep = ""),
                                       paste("maximum MSE-inefficiency: ", round(ineff[1], 3), sep = "")), 
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
-                return(structure(list("estimate" = robEst, "steps" = k, "Infos" = Info.matrix), 
-                                 class = c("ALEstimate", "Estimate")))
+                return(new("kStepEstimate", name = "Optimally robust estimate",
+                           estimate = robEst, steps = k, pIC = NULL, Infos = Info.matrix))
             }
         }
     }

Modified: branches/robast-0.6/pkg/RobLox/inst/tests/tests.R
===================================================================
--- branches/robast-0.6/pkg/RobLox/inst/tests/tests.R	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/inst/tests/tests.R	2008-07-28 12:16:05 UTC (rev 130)
@@ -9,13 +9,13 @@
 
 ## location and scale, radius unknown
 res1 <- roblox(x, returnIC = TRUE)
-checkIC(res1$optIC)
+checkIC(pIC(res1))
 
 res11 <- roblox(x, returnIC = TRUE, k = 2)
-checkIC(res11$optIC)
+checkIC(pIC(res11))
 
 res12 <- roblox(x, returnIC = TRUE, k = 5)
-checkIC(res12$optIC)
+checkIC(pIC(res12))
 
 roblox(x)
 roblox(x, k = 2)
@@ -24,13 +24,13 @@
 
 ## location and scale, radius interval
 res2 <- roblox(x, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE)
-checkIC(res2$optIC)
+checkIC(pIC(res2))
 
 res21 <- roblox(x, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE, k = 2)
-checkIC(res21$optIC)
+checkIC(pIC(res21))
 
 res22 <- roblox(x, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE, k = 4)
-checkIC(res22$optIC)
+checkIC(pIC(res22))
 
 roblox(x, eps.lower = 0.15, eps.upper = 0.3)
 roblox(x, eps.lower = 0.15, eps.upper = 0.3, k = 2)
@@ -39,13 +39,13 @@
 
 ## scale, radius interval
 res3 <- roblox(x, mean = -2, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE)
-checkIC(res3$optIC)
+checkIC(pIC(res3))
 
 res31 <- roblox(x, mean = -2, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE, k = 3)
-checkIC(res31$optIC)
+checkIC(pIC(res31))
 
 res32 <- roblox(x, mean = -2, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE, k = 6)
-checkIC(res32$optIC)
+checkIC(pIC(res32))
 
 roblox(x, mean = -2, eps.lower = 0.15, eps.upper = 0.3)
 roblox(x, mean = -2, eps.lower = 0.15, eps.upper = 0.3, k = 3)
@@ -54,13 +54,13 @@
 
 ## location, radius interval
 res4 <- roblox(x, sd = 3, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE)
-checkIC(res4$optIC)
+checkIC(pIC(res4))
 
 res41 <- roblox(x, sd = 3, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE, k = 2)
-checkIC(res41$optIC)
+checkIC(pIC(res41))
 
 res42 <- roblox(x, sd = 3, eps.lower = 0.15, eps.upper = 0.3, returnIC = TRUE, k = 5)
-checkIC(res42$optIC)
+checkIC(pIC(res42))
 
 roblox(x, sd = 3, eps.lower = 0.15, eps.upper = 0.3)
 roblox(x, sd = 3, eps.lower = 0.15, eps.upper = 0.3, k = 2)
@@ -215,7 +215,9 @@
 system.time(rowRoblox(X, eps = 0.02))
 
 system.time(apply(X, 1, roblox))
+## uses rowMedians of package Biobase if available
 system.time(rowRoblox(X))
+
 M <- apply(X, 1, median)
 S <- apply(X, 1, mad)
 init <- cbind(M, S)

Modified: branches/robast-0.6/pkg/RobLox/man/rlOptIC.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlOptIC.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlOptIC.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -37,7 +37,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{ContIC-class}}, \code{\link{roblox}}}
+\seealso{\code{\link[RobAStBase]{ContIC-class}}, \code{\link{roblox}}}
 \examples{
 IC1 <- rlOptIC(r = 0.1)
 distrExOptions("ErelativeTolerance" = 1e-12)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AL.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AL.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AL.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -50,7 +50,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{ContIC-class}}, \code{\link{roblox}}}
+\seealso{\code{\link[RobAStBase]{ContIC-class}}, \code{\link{roblox}}}
 \examples{
 IC1 <- rlsOptIC.AL(r = 0.1, check = TRUE)
 distrExOptions("ErelativeTolerance" = 1e-12)
@@ -78,14 +78,14 @@
 (est0 <- MDEstimator(x, ParamFamily = NormLocationScaleFamily(), distance = KolmogorovDist))
 
 ## 3. one-step estimation: radius known
-IC1 <- rlsOptIC.AL(r = 0.5, mean = est0$estimate[1], sd = est0$estimate[2])
-(est1 <- oneStepEstimator(x, IC1, est0$estimate))
+IC1 <- rlsOptIC.AL(r = 0.5, mean = estimate(est0)[1], sd = estimate(est0)[2])
+(est1 <- oneStepEstimator(x, IC1, est0))
 
 ## 4. one-step estimation: radius unknown
 ## take least favorable radius r = 0.579
 ## cf. Table 8.1 in Kohl(2005)
-IC2 <- rlsOptIC.AL(r = 0.579, mean = est0$estimate[1], sd = est0$estimate[2])
-(est2 <- oneStepEstimator(x, IC2, est0$estimate))
+IC2 <- rlsOptIC.AL(r = 0.579, mean = estimate(est0)[1], sd = estimate(est0)[2])
+(est2 <- oneStepEstimator(x, IC2, est0))
 }
 \concept{normal location and scale}
 \concept{influence curve}

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An1.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An1.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An1.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -29,7 +29,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.An1(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An2.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An2.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.An2.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -37,7 +37,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.An2(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AnMad.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AnMad.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.AnMad.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -30,7 +30,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.AnMad(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.BM.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.BM.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.BM.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -39,7 +39,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.BM(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha3.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha3.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha3.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -34,7 +34,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Ha3(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha4.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha4.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Ha4.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -38,7 +38,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Ha4(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HaMad.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HaMad.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HaMad.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -39,7 +39,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.HaMad(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu1.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu1.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu1.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -29,7 +29,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Hu1(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -35,7 +35,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Hu2(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2a.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2a.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu2a.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -37,7 +37,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Hu2a(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu3.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu3.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Hu3.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -36,7 +36,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Hu3(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HuMad.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HuMad.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.HuMad.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -30,7 +30,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.HuMad(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.M.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.M.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.M.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -37,7 +37,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.M(r = 0.1, check = TRUE)
 distrExOptions("ErelativeTolerance" = 1e-12)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.MM2.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.MM2.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.MM2.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -37,7 +37,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.MM2(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu1.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu1.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu1.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -29,7 +29,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Tu1(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu2.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu2.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.Tu2.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -37,7 +37,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.Tu2(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/rlsOptIC.TuMad.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rlsOptIC.TuMad.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rlsOptIC.TuMad.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -29,7 +29,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{IC-class}}}
+\seealso{\code{\link[RobAStBase]{IC-class}}}
 \examples{
 IC1 <- rlsOptIC.TuMad(r = 0.1)
 checkIC(IC1)

Modified: branches/robast-0.6/pkg/RobLox/man/roblox.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/roblox.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/roblox.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -55,16 +55,7 @@
   If \code{eps} is missing, the radius-minimax estimator in sense of 
   Rieder et al. (2001), respectively Section 2.2 of Kohl (2005) is returned.
 }
-\value{
-  An object of S3-class \code{"ALEstimate"} which inherits from 
-  class \code{"Estimate"}, a list with components
-  \item{estimate }{ location and/or scale estimate }
-  \item{steps }{ number of k-steps used to compute the estimate }
-  \item{Infos }{ additional information about the estimate }
-
-  if \code{returnIC} is \code{TRUE} the list also contains
-  \item{optIC}{ object of class \code{"ContIC"}; optimally robust IC }
-}
+\value{Object of class \code{"kStepEstimate"}. }
 \references{
   Kohl, M. (2005) \emph{Numerical Contributions to the Asymptotic Theory of Robustness}. 
   Bayreuth: Dissertation.
@@ -82,31 +73,32 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{ContIC-class}}, \code{\link{rlOptIC}}, 
-         \code{\link{rsOptIC}}, \code{\link{rlsOptIC.AL}}}
+\seealso{\code{\link[RobAStBase]{ContIC-class}}, \code{\link{rlOptIC}}, 
+         \code{\link{rsOptIC}}, \code{\link{rlsOptIC.AL}},
+         \code{\link[RobAStBase]{kStepEstimate-class}} }
 \examples{
 ind <- rbinom(100, size=1, prob=0.05) 
 x <- rnorm(100, mean=ind*3, sd=(1-ind) + ind*9)
 
 ## amount of gross errors known
 res1 <- roblox(x, eps = 0.05, returnIC = TRUE)
-res1$estimate
-res1$optIC
-checkIC(res1$optIC)
-Risks(res1$optIC)
-Infos(res1$optIC)
-plot(res1$optIC)
-infoPlot(res1$optIC)
+estimate(res1)
+pIC(res1)
+checkIC(pIC(res1))
+Risks(pIC(res1))
+Infos(pIC(res1))
+plot(pIC(res1))
+infoPlot(pIC(res1))
 
 ## amount of gross errors unknown
 res2 <- roblox(x, eps.lower = 0.01, eps.upper = 0.1, returnIC = TRUE)
-res2$estimate
-res2$optIC
-checkIC(res2$optIC)
-Risks(res2$optIC)
-Infos(res2$optIC)
-plot(res2$optIC)
-infoPlot(res2$optIC)
+estimate(res2)
+pIC(res2)
+checkIC(pIC(res2))
+Risks(pIC(res2))
+Infos(pIC(res2))
+plot(pIC(res2))
+infoPlot(pIC(res2))
 
 ## estimator comparison
 # classical optimal (non-robust)
@@ -116,22 +108,22 @@
 c(median(x), mad(x))
 
 # optimally robust (amount of gross errors known)
-res1$estimate
+estimate(res1)
 
 # optimally robust (amount of gross errors unknown)
-res2$estimate
+estimate(res2)
 
 # Kolmogorov(-Smirnov) minimum distance estimator (robust)
-(ks.est <- MDEstimator(x, ParamFamily = NormLocationScaleFamily(), distance = KolmogorovDist))
+(ks.est <- MDEstimator(x, ParamFamily = NormLocationScaleFamily()))
 
 # optimally robust (amount of gross errors known)
-roblox(x, eps = 0.05, initial.est = ks.est$estimate)
+roblox(x, eps = 0.05, initial.est = estimate(ks.est))
 
 # Cramer von Mises minimum distance estimator (robust)
 (CvM.est <- MDEstimator(x, ParamFamily = NormLocationScaleFamily(), distance = CvMDist))
 
 # optimally robust (amount of gross errors known)
-roblox(x, eps = 0.05, initial.est = CvM.est$estimate)
+roblox(x, eps = 0.05, initial.est = estimate(CvM.est))
 }
 \concept{normal location}
 \concept{normal scale}

Modified: branches/robast-0.6/pkg/RobLox/man/rowRoblox.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rowRoblox.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rowRoblox.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -1,10 +1,10 @@
 \name{rowRoblox and colRoblox}
 \alias{rowRoblox}
 \alias{colRoblox}
-\title{Optimally robust estimator for location and/or scale}
+\title{Optimally robust estimation for location and/or scale}
 \description{
   The functions \code{rowRoblox} and \code{colRoblox} compute 
-  the optimally robust estimator for normal location und/or scale and 
+  optimally robust estimates for normal location und/or scale and 
   (convex) contamination neighborhoods. The definition of 
   these estimators can be found in Rieder (1994) or Kohl (2005),
   respectively.
@@ -69,13 +69,7 @@
   or a vector with length identical to the number of rows, respectively 
   columns.
 }
-\value{
-  An object of S3-class \code{"ALEstimate"} which inherits from 
-  class \code{"Estimate"}, a list with components
-  \item{estimate }{ location and/or scale estimates }
-  \item{steps }{ number of k-steps used to compute the estimates }
-  \item{Infos }{ additional information about the estimates }
-}
+\value{Object of class \code{"kStepEstimate"}. }
 \references{
   Kohl, M. (2005) \emph{Numerical Contributions to the Asymptotic Theory of Robustness}. 
   Bayreuth: Dissertation.
@@ -93,7 +87,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link{roblox}}}
+\seealso{\code{\link{roblox}}, \code{\link[RobAStBase]{kStepEstimate-class}} }
 \examples{
 ind <- rbinom(200, size=1, prob=0.05) 
 X <- matrix(rnorm(200, mean=ind*3, sd=(1-ind) + ind*9), nrow = 2)

Modified: branches/robast-0.6/pkg/RobLox/man/rsOptIC.Rd
===================================================================
--- branches/robast-0.6/pkg/RobLox/man/rsOptIC.Rd	2008-07-28 11:33:52 UTC (rev 129)
+++ branches/robast-0.6/pkg/RobLox/man/rsOptIC.Rd	2008-07-28 12:16:05 UTC (rev 130)
@@ -39,7 +39,7 @@
 }
 \author{Matthias Kohl \email{Matthias.Kohl at stamats.de}}
 %\note{}
-\seealso{\code{\link[ROptEst]{ContIC-class}}, \code{\link{roblox}}}
+\seealso{\code{\link[RobAStBase]{ContIC-class}}, \code{\link{roblox}}}
 \examples{
 IC1 <- rsOptIC(r = 0.1)
 distrExOptions("ErelativeTolerance" = 1e-12)



More information about the Robast-commits mailing list