[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