[Robast-commits] r466 - branches/robast-0.9/pkg/RobAStBase/R

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

Author: ruckdeschel
Date: 2012-05-19 18:30:33 +0200 (Sat, 19 May 2012)
New Revision: 466

Modified:
   branches/robast-0.9/pkg/RobAStBase/R/kStepEstimator.R
Log:
RobAStbase: a small change in kStepEstimator -- now uses the (possibly nonlinear) trafo instead of the Jacobian to get the pIC-estimator from the untransformed estimator

Modified: branches/robast-0.9/pkg/RobAStBase/R/kStepEstimator.R
===================================================================
--- branches/robast-0.9/pkg/RobAStBase/R/kStepEstimator.R	2012-03-03 22:41:04 UTC (rev 465)
+++ branches/robast-0.9/pkg/RobAStBase/R/kStepEstimator.R	2012-05-19 16:30:33 UTC (rev 466)
@@ -95,8 +95,6 @@
                 IC.c <- as(diag(p) %*% IC at Curve, "EuclRandVariable")
 
 
-                theta <- theta + rowMeans(evalRandVar(IC.c, x0),
-                                          na.rm = na.rm )
 
 #                print(theta)
                 tf <- trafo(L2Fam, Param)
@@ -127,7 +125,10 @@
                      IC.tot <- IC.tot1 + IC.tot2
                      u.theta <- u.theta + rowMeans(evalRandVar(IC.tot, x0),
                                                    na.rm = na.rm)
+                     theta <- (tf$fct(u.theta))$fval
                 }else{
+                     theta <- theta + rowMeans(evalRandVar(IC.c, x0),
+                                               na.rm = na.rm )
                      IC.tot <- IC.c
                      u.theta <- theta
                 }