[Robast-commits] r133 - branches/robast-0.6/pkg/RobLox/R

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

Author: stamats
Date: 2008-07-31 18:23:52 +0200 (Thu, 31 Jul 2008)
New Revision: 133

Modified:
   branches/robast-0.6/pkg/RobLox/R/rowRoblox.R
Log:
minor change

Modified: branches/robast-0.6/pkg/RobLox/R/rowRoblox.R
===================================================================
--- branches/robast-0.6/pkg/RobLox/R/rowRoblox.R	2008-07-31 15:07:01 UTC (rev 132)
+++ branches/robast-0.6/pkg/RobLox/R/rowRoblox.R	2008-07-31 16:23:52 UTC (rev 133)
@@ -335,13 +335,14 @@
                     b <- sd*.getb.sc(r)
                 }
                 robEst <- .kstep.sc.matrix(x = x, initial.est = sd, A = A, a = a, b = b, mean = mean, k = k)
-                names(robEst$est) <- "sd"
+                robEst$est <- as.matrix(robEst$est)
+                colnames(robEst$est) <- "sd"
                 Info.matrix <- matrix(c("roblox", 
                                         paste("optimally robust estimates for contamination 'eps' =", round(eps, 3),
                                               "and 'asMSE'")),
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
                 return(new("kStepEstimate", name = "Optimally robust estimate",
-                           estimate = as.matrix(robEst$est), samplesize = ncol(x), steps = k, 
+                           estimate = robEst$est, samplesize = ncol(x), steps = k, 
                            pIC = NULL, Infos = Info.matrix))
 ## we need a class like "list of estimates" to set asvar and asbias consistently ...
 #                return(new("kStepEstimate", name = "Optimally robust estimate",
@@ -377,7 +378,8 @@
                     }
                 }
                 robEst <- .kstep.sc.matrix(x = x, initial.est = sd, A = A, a = a, b = b, mean = mean, k = k)
-                names(robEst$est) <- "sd"
+                robEst$est <- as.matrix(robEst$est)
+                colnames(robEst$est) <- "sd"
                 Info.matrix <- matrix(c(rep("roblox", 3), 
                                       paste("radius-minimax estimates for contamination interval [", 
                                         round(eps.lower, 3), ", ", round(eps.upper, 3), "]", sep = ""),
@@ -385,7 +387,7 @@
                                       paste("maximum MSE-inefficiency: ", round(ineff[1], 3), sep = "")), 
                                       ncol = 2, dimnames = list(NULL, c("method", "message")))
                 return(new("kStepEstimate", name = "Optimally robust estimate",
-                           estimate = as.matrix(robEst$est), samplesize = ncol(x), steps = k, 
+                           estimate = robEst$est, samplesize = ncol(x), steps = k, 
                            pIC = NULL, Infos = Info.matrix))
 ## we need a class like "list of estimates" to set asvar and asbias consistently ...
 #                return(new("kStepEstimate", name = "Optimally robust estimate",