[Robast-commits] r359 - in branches/robast-0.7/pkg: ROptEst/R ROptEst/chm ROptEst/man RobAStBase RobAStBase/R RobAStBase/chm RobAStBase/man

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Mon Aug 31 19:04:47 CEST 2009


Author: ruckdeschel
Date: 2009-08-31 19:04:45 +0200 (Mon, 31 Aug 2009)
New Revision: 359

Added:
   branches/robast-0.7/pkg/RobAStBase/man/OptionalInfluenceCurve-Class.Rd
Modified:
   branches/robast-0.7/pkg/ROptEst/R/roptest.R
   branches/robast-0.7/pkg/ROptEst/chm/ROptEst.chm
   branches/robast-0.7/pkg/ROptEst/chm/roptest.html
   branches/robast-0.7/pkg/ROptEst/man/roptest.Rd
   branches/robast-0.7/pkg/RobAStBase/NAMESPACE
   branches/robast-0.7/pkg/RobAStBase/R/AllClass.R
   branches/robast-0.7/pkg/RobAStBase/R/AllGeneric.R
   branches/robast-0.7/pkg/RobAStBase/R/RobAStBaseOptions.R
   branches/robast-0.7/pkg/RobAStBase/R/kStepEstimator.R
   branches/robast-0.7/pkg/RobAStBase/R/oneStepEstimator.R
   branches/robast-0.7/pkg/RobAStBase/chm/00Index.html
   branches/robast-0.7/pkg/RobAStBase/chm/RobAStBase.chm
   branches/robast-0.7/pkg/RobAStBase/chm/RobAStBase.toc
   branches/robast-0.7/pkg/RobAStBase/chm/RobAStBaseOptions.html
   branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimate-class.html
   branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimator.html
   branches/robast-0.7/pkg/RobAStBase/man/RobAStBaseOptions.Rd
   branches/robast-0.7/pkg/RobAStBase/man/kStepEstimate-class.Rd
   branches/robast-0.7/pkg/RobAStBase/man/kStepEstimator.Rd
Log:
class kStepEstimator gains slots 
-start:   argument start of kStepEstimator
-startval:matrix of transformed starting values
-ustarval:matrix of untransformed starting values
-ksteps:  matrix of (intermediate) transformed estimates
-uksteps: matrix of (intermediate) untransformed estimates
-ICList:  optional: the (intermediate) (total) ICs
-pICList: optional: the (intermediate) [p]ICs
and corresponding accessors

Modified: branches/robast-0.7/pkg/ROptEst/R/roptest.R
===================================================================
--- branches/robast-0.7/pkg/ROptEst/R/roptest.R	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/ROptEst/R/roptest.R	2009-08-31 17:04:45 UTC (rev 359)
@@ -7,6 +7,8 @@
                     useLast = getRobAStBaseOption("kStepUseLast"),
                     withUpdateInKer = getRobAStBaseOption("withUpdateInKer"),
                     IC.UpdateInKer = getRobAStBaseOption("IC.UpdateInKer"),
+                    withICList = getRobAStBaseOption("withICList"),
+                    withPICList = getRobAStBaseOption("withPICList"),
                     na.rm = TRUE, initial.est.ArgList, ...){
     if(missing(verbose)|| is.null(verbose))
            verbose <- getRobAStBaseOption("all.verbose")
@@ -96,7 +98,7 @@
     }
     res <- kStepEstimator(x, IC = ICstart, start = initial.est, steps = steps, useLast = useLast,
                           withUpdateInKer = withUpdateInKer, IC.UpdateInKer = IC.UpdateInKer,
-                          na.rm = na.rm)
+                          withICList = withICList, withPICList = withPICList, na.rm = na.rm)
     res at estimate.call <- es.call
     Infos <- matrix(c("roptest", 
                       paste(steps, "-step estimate for ", name(L2Fam), sep = "")),

Modified: branches/robast-0.7/pkg/ROptEst/chm/ROptEst.chm
===================================================================
(Binary files differ)

Modified: branches/robast-0.7/pkg/ROptEst/chm/roptest.html
===================================================================
--- branches/robast-0.7/pkg/ROptEst/chm/roptest.html	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/ROptEst/chm/roptest.html	2009-08-31 17:04:45 UTC (rev 359)
@@ -30,6 +30,8 @@
         useLast = getRobAStBaseOption("kStepUseLast"),
         withUpdateInKer = getRobAStBaseOption("withUpdateInKer"),
         IC.UpdateInKer = getRobAStBaseOption("IC.UpdateInKer"),
+        withICList = getRobAStBaseOption("withICList"),
+        withPICList = getRobAStBaseOption("withPICList"),
         na.rm = TRUE, initial.est.ArgList, ...)
 </pre>
 
@@ -100,6 +102,14 @@
 if there is a non-trivial trafo in the model with matrix <i>D</i>,
 the IC to be used for this; if <code>NULL</code> the result of <code>getboundedIC(L2Fam,D)</code> is taken;
 this IC will then be projected onto <i>ker(D)</i>.</td></tr>
+<tr valign="top"><td><code>withPICList</code></td>
+<td>
+logical: shall slot <code>pICList</code> of return value
+be filled?</td></tr>
+<tr valign="top"><td><code>withICList</code></td>
+<td>
+logical: shall slot <code>ICList</code> of return value
+be filled?</td></tr>
 <tr valign="top"><td><code>na.rm</code></td>
 <td>
 logical: if  <code>TRUE</code>, the estimator is evaluated at <code>complete.cases(x)</code>.</td></tr>

Modified: branches/robast-0.7/pkg/ROptEst/man/roptest.Rd
===================================================================
--- branches/robast-0.7/pkg/ROptEst/man/roptest.Rd	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/ROptEst/man/roptest.Rd	2009-08-31 17:04:45 UTC (rev 359)
@@ -12,6 +12,8 @@
         useLast = getRobAStBaseOption("kStepUseLast"),
         withUpdateInKer = getRobAStBaseOption("withUpdateInKer"),
         IC.UpdateInKer = getRobAStBaseOption("IC.UpdateInKer"),
+        withICList = getRobAStBaseOption("withICList"),
+        withPICList = getRobAStBaseOption("withPICList"),
         na.rm = TRUE, initial.est.ArgList, ...)
 }
 \arguments{
@@ -46,6 +48,10 @@
   \item{IC.UpdateInKer}{if there is a non-trivial trafo in the model with matrix \eqn{D},
   the IC to be used for this; if \code{NULL} the result of \code{getboundedIC(L2Fam,D)} is taken;
   this IC will then be projected onto \eqn{{\rm ker}(D)}{ker(D)}.}
+  \item{withPICList}{logical: shall slot \code{pICList} of return value
+  be filled?}
+  \item{withICList}{logical: shall slot \code{ICList} of return value
+  be filled?}
   \item{na.rm}{logical: if  \code{TRUE}, the estimator is evaluated at \code{complete.cases(x)}.}
   \item{initial.est.ArgList}{a list of arguments to be given to argument \code{start} if the latter
    is a function; this list by default already starts with two unnamed items,

Modified: branches/robast-0.7/pkg/RobAStBase/NAMESPACE
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/NAMESPACE	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/NAMESPACE	2009-08-31 17:04:45 UTC (rev 359)
@@ -18,6 +18,7 @@
               "BdStWeight", "HampelWeight")
 exportClasses("ALEstimate", "kStepEstimate", "MEstimate")
 exportClasses("cutoff")
+exportClasses("StartClass", "pICList", "OptionalpICList")
 exportMethods("show", 
               "plot")
 exportMethods("type", "radius", "radius<-")
@@ -40,7 +41,7 @@
               "neighborRadius", "neighborRadius<-", 
               "clipLo", "clipLo<-", 
               "clipUp", "clipUp<-",
-              "optIC") 
+              "optIC", "start", "startval", "pICList") 
 exportMethods("locMEstimator")
 exportMethods("weight", "weight<-", 
               "getweight", 
@@ -53,6 +54,8 @@
 exportMethods("pIC", "asbias", 
               "steps", 
               "Mroot","kStepEstimator.start")
+exportMethods("pICList","ICList", "ksteps", "uksteps", 
+              "start", "startval", "ustartval")
 exportMethods("ddPlot")
 exportMethods("cutoff.quantile", "cutoff.quantile<-")
 export("oneStepEstimator", "kStepEstimator")

Modified: branches/robast-0.7/pkg/RobAStBase/R/AllClass.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/AllClass.R	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/R/AllClass.R	2009-08-31 17:04:45 UTC (rev 359)
@@ -198,7 +198,21 @@
 
 ## ALEstimate
 setClassUnion("OptionalInfluenceCurve", c("InfluenceCurve", "NULL"))
-setClass("ALEstimate", 
+setClassUnion("StartClass", c("numeric", "function", "Estimate"))
+setClass("pICList",
+          prototype = prototype(list()),
+            contains = "list",
+            validity = function(object){
+                nrvalues <- length(object)
+                if(nrvalues){
+                for(i in 1:nrvalues)
+                    if(!is(object[[i]], "OptionalInfluenceCurve"))
+                        stop("element ", i, " is no 'OptionalInfluenceCurve'")
+                }
+                return(TRUE)
+            })
+setClassUnion("OptionalpICList", c("pICList", "NULL"))
+setClass("ALEstimate",
          representation(pIC = "OptionalInfluenceCurve",
                         asbias = "OptionalNumeric"),
          prototype(name = "Asymptotically linear estimate",
@@ -219,7 +233,14 @@
                    untransformed.asvar = NULL),
          contains = "Estimate")
 setClass("kStepEstimate", 
-         representation(steps = "integer"),
+         representation(steps = "integer",
+                        pICList = "OptionalpICList",
+                        ICList = "OptionalpICList",
+                        start = "StartClass",
+                        startval = "matrix",
+                        ustartval = "matrix",
+                        ksteps = "OptionalMatrix",
+                        uksteps = "OptionalMatrix"),
          prototype(name = "Asymptotically linear estimate",
                    estimate = numeric(0),
                    samplesize = numeric(0),
@@ -229,6 +250,13 @@
                    asvar = NULL,
                    asbias = NULL,
                    pIC = NULL,
+                   pICList = NULL,
+                   ICList = NULL,
+                   ksteps = NULL,
+                   uksteps = NULL,
+                   start = NULL,
+                   startval = matrix(0),
+                   ustartval = matrix(0),
                    nuis.idx = NULL,
                    trafo = list(fct = function(x){
                                       list(fval = x, mat = matrix(1))},

Modified: branches/robast-0.7/pkg/RobAStBase/R/AllGeneric.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/AllGeneric.R	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/R/AllGeneric.R	2009-08-31 17:04:45 UTC (rev 359)
@@ -176,6 +176,27 @@
 if(!isGeneric("steps")){
     setGeneric("steps", function(object) standardGeneric("steps"))
 }
+if(!isGeneric("ksteps")){
+    setGeneric("ksteps", function(object,...) standardGeneric("ksteps"))
+}
+if(!isGeneric("uksteps")){
+    setGeneric("uksteps", function(object,...) standardGeneric("uksteps"))
+}
+if(!isGeneric("start")){
+    setGeneric("start", function(object) standardGeneric("start"))
+}
+if(!isGeneric("startval")){
+    setGeneric("startval", function(object) standardGeneric("startval"))
+}
+if(!isGeneric("ustartval")){
+    setGeneric("ustartval", function(object) standardGeneric("ustartval"))
+}
+if(!isGeneric("ICList")){
+    setGeneric("ICList", function(object) standardGeneric("ICList"))
+}
+if(!isGeneric("pICList")){
+    setGeneric("pICList", function(object) standardGeneric("pICList"))
+}
 if(!isGeneric("Mroot")){
     setGeneric("Mroot", function(object) standardGeneric("Mroot"))
 }

Modified: branches/robast-0.7/pkg/RobAStBase/R/RobAStBaseOptions.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/RobAStBaseOptions.R	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/R/RobAStBaseOptions.R	2009-08-31 17:04:45 UTC (rev 359)
@@ -2,7 +2,9 @@
     kStepUseLast = FALSE,
     withUpdateInKer = FALSE,
     IC.UpdateInKer = NULL,
-    all.verbose = FALSE
+    all.verbose = FALSE,
+    withICList = FALSE,
+    withPICList = FALSE
 )
 
 RobAStBaseOptions <- function(...) {

Modified: branches/robast-0.7/pkg/RobAStBase/R/kStepEstimator.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/kStepEstimator.R	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/R/kStepEstimator.R	2009-08-31 17:04:45 UTC (rev 359)
@@ -11,6 +11,8 @@
                            useLast = getRobAStBaseOption("kStepUseLast"),
                            withUpdateInKer = getRobAStBaseOption("withUpdateInKer"),
                            IC.UpdateInKer = getRobAStBaseOption("IC.UpdateInKer"),
+                           withICList = getRobAStBaseOption("withICList"),
+                           withPICList = getRobAStBaseOption("withPICList"),
                            na.rm = TRUE, startArgList = NULL, ...){
 ## save call
         es.call <- match.call()
@@ -74,7 +76,15 @@
         u.theta <- start.val
         theta <- if(is(start.val,"Estimate")) estimate(start.val)
                  else trafoF(u.theta)$fval
+        u.start.val <- matrix(start.val,nrow=1)
+        start.val <- matrix(theta,ncol=1)
+        rownames(u.start.val) <- u.est.names
+        rownames(start.val) <- est.names
 
+### shall intermediate IC's / pIC's be stored?
+        pICList <- if(withPICList) vector("list", steps) else NULL
+        ICList <- if(withICList) vector("list", steps) else NULL
+
         ### update - function
         updateStep <- function(u.theta, theta, IC, L2Fam, Param,
                                withModif = TRUE, with.u.var = FALSE){
@@ -141,7 +151,7 @@
 
                 return(list(IC = IC, Param = Param, L2Fam = L2Fam,
                             theta = theta, u.theta = u.theta, u.var = u.var,
-                            var = var0))
+                            var = var0, IC.tot = IC.tot, IC.c = IC.c))
         }
 
         Infos <- matrix(c("kStepEstimator",
@@ -155,6 +165,10 @@
 #        print(IC at Risks$asCov)
 #        print(Risks(IC)$asCov)
 
+        uksteps = matrix(0,ncol=steps, nrow = k)
+        ksteps = matrix(0,ncol=steps, nrow = p)
+        rownames(ksteps) <- est.names
+        rownames(uksteps) <- u.est.names
         if(!is(modifyIC(IC), "NULL") ){
            for(i in 1:steps){
                if(i>1){
@@ -166,11 +180,25 @@
                upd <- updateStep(u.theta,theta,IC, L2Fam, Param,
                                  withModif = (steps>1) | useLast,
                                  with.u.var = i==steps)
-               u.theta <- upd$u.theta
-               theta <- upd$theta
+               uksteps[,i] <- u.theta <- upd$u.theta
+               ksteps[,i] <- theta <- upd$theta
+               if(withICList)
+                  ICList[[i]] <- new("InfluenceCurve",
+                                      name = paste(gettext("(total) IC in step"),i),
+                                      Risks = list(),
+                                      Infos = matrix(c("",""),ncol=2),
+                                      Curve =  EuclRandVarList(upd$IC.tot))
+               if(withPICList)
+                  pICList[[i]] <- new("InfluenceCurve",
+                                      name = paste(gettext("pIC in step"),i),
+                                      Risks = list(),
+                                      Infos = matrix(c("",""),ncol=2),
+                                      Curve =  EuclRandVarList(upd$IC.c))
                u.var <- upd$u.var
                var0 <- upd$var
            }
+           if(withICList) ICList <- new("pICList",ICList)
+           if(withPICList) pICList <- new("pICList",pICList)
            if(useLast){
               IC <- upd$IC
               L2Fam <- upd$L2Fam
@@ -187,7 +215,9 @@
               stop("slot 'modifyIC' of 'IC' is 'NULL'!")
            upd <- updateStep(u.theta,theta,IC, L2Fam, Param, withModif = FALSE)
            u.theta <- upd$u.theta
+           uksteps <- NULL
            theta <- upd$theta
+           ksteps <- NULL
            u.var <- upd$u.var
            var0 <- upd$var
            if(useLast){
@@ -260,7 +290,10 @@
                        nuis.idx = nuis.idx, untransformed.estimate = u.theta,
                        completecases = completecases,
                        untransformed.asvar = u.var,
-                       asbias = asBias, pIC = IC, steps = steps, Infos = Infos))
+                       asbias = asBias, pIC = IC, steps = steps, Infos = Infos,
+                       start = start, startval = start.val, ustartval = u.start.val,
+                       ksteps = ksteps, uksteps = uksteps,
+                       pICList = pICList, ICList = ICList))
 }
 #  (est1.NS <- kStepEstimator(x, IC2.NS, est0, steps = 1))
 

Modified: branches/robast-0.7/pkg/RobAStBase/R/oneStepEstimator.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/oneStepEstimator.R	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/R/oneStepEstimator.R	2009-08-31 17:04:45 UTC (rev 359)
@@ -57,7 +57,8 @@
              fixed = fixed, nuis.idx = nuis.idx,
              completecases = completecases,
             estimate.call = es.call, samplesize = nrow(x0), asvar = asVar,
-            asbias = asBias, pIC = IC, steps = 1L, Infos = Infos)
+            asbias = asBias, pIC = IC, steps = 1L, Infos = Infos,
+            start = start, startval = start.val, ustartval = start.val)
     }
 
 

Modified: branches/robast-0.7/pkg/RobAStBase/chm/00Index.html
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/chm/00Index.html	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/chm/00Index.html	2009-08-31 17:04:45 UTC (rev 359)
@@ -309,6 +309,10 @@
 <td>Generating function for IC-class</td></tr>
 <tr><td width="25%"><a href="IC-class.html">IC-class</a></td>
 <td>Influence curve</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">ICList</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">ICList,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
 <tr><td width="25%"><a href="InfluenceCurve.html">InfluenceCurve</a></td>
 <td>Generating function for InfluenceCurve-class</td></tr>
 <tr><td width="25%"><a href="InfluenceCurve-class.html">InfluenceCurve-class</a></td>
@@ -350,6 +354,10 @@
 <td>Methods for function kStepEstimator.start in Package 'RobAStBase'</td></tr>
 <tr><td width="25%"><a href="kStepEstimatorStart-methods.html">kStepEstimator.start-methods</a></td>
 <td>Methods for function kStepEstimator.start in Package 'RobAStBase'</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">ksteps</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">ksteps,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
 <tr><td width="25%"><a href="RobAStBaseOptions.html">kStepUseLast</a></td>
 <td>Function to change the global variables of the package 'RobAStBase'</td></tr>
 </table>
@@ -472,6 +480,10 @@
 <td>Generic function for the computation of optimally robust ICs</td></tr>
 <tr><td width="25%"><a href="optIC.html">optIC-methods</a></td>
 <td>Generic function for the computation of optimally robust ICs</td></tr>
+<tr><td width="25%"><a href="OptionalInfluenceCurve-Class.html">OptionalInfluenceCurve-class</a></td>
+<td>Some helper Classes in package 'RobAStBase'</td></tr>
+<tr><td width="25%"><a href="OptionalInfluenceCurve-Class.html">OptionalpICList-class</a></td>
+<td>Some helper Classes in package 'RobAStBase'</td></tr>
 <tr><td width="25%"><a href="outlyingPlotIC.html">outlyingPlotIC</a></td>
 <td>Function outlyingPlotIC in Package 'RobAStBase'</td></tr>
 </table>
@@ -483,6 +495,12 @@
 <td>ALEstimate-class.</td></tr>
 <tr><td width="25%"><a href="ALEstimate-class.html">pIC,ALEstimate-method</a></td>
 <td>ALEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">pICList</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">pICList,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="OptionalInfluenceCurve-Class.html">pICList-class</a></td>
+<td>Some helper Classes in package 'RobAStBase'</td></tr>
 <tr><td width="25%"><a href="plot-methods.html">plot</a></td>
 <td>Methods for Function plot in Package 'RobAStBase'</td></tr>
 <tr><td width="25%"><a href="plot-methods.html">plot,IC,missing-method</a></td>
@@ -563,6 +581,16 @@
 <td>Influence curve of contamination type</td></tr>
 <tr><td width="25%"><a href="TotalVarIC-class.html">stand&lt;-,TotalVarIC-method</a></td>
 <td>Influence curve of total variation type</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">start</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">start,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="OptionalInfluenceCurve-Class.html">StartClass-class</a></td>
+<td>Some helper Classes in package 'RobAStBase'</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">startval</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">startval,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
 <tr><td width="25%"><a href="kStepEstimate-class.html">steps</a></td>
 <td>kStepEstimate-class.</td></tr>
 <tr><td width="25%"><a href="kStepEstimate-class.html">steps,kStepEstimate-method</a></td>
@@ -591,8 +619,16 @@
 <h2><a name="U">-- U --</a></h2>
 
 <table width="100%">
+<tr><td width="25%"><a href="kStepEstimate-class.html">uksteps</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">uksteps,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
 <tr><td width="25%"><a href="UncondNeighborhood-class.html">UncondNeighborhood-class</a></td>
 <td>Unconditional neighborhood</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">ustartval</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">ustartval,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
 </table>
 
 <h2><a name="W">-- W --</a></h2>

Modified: branches/robast-0.7/pkg/RobAStBase/chm/RobAStBase.chm
===================================================================
(Binary files differ)

Modified: branches/robast-0.7/pkg/RobAStBase/chm/RobAStBase.toc
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/chm/RobAStBase.toc	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/chm/RobAStBase.toc	2009-08-31 17:04:45 UTC (rev 359)
@@ -458,6 +458,14 @@
 <param name="Local" value="IC-class.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ICList">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ICList,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="InfluenceCurve">
 <param name="Local" value="InfluenceCurve.html">
 </OBJECT>
@@ -530,6 +538,14 @@
 <param name="Local" value="kStepEstimatorStart-methods.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ksteps">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ksteps,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="kStepUseLast">
 <param name="Local" value="RobAStBaseOptions.html">
 </OBJECT>
@@ -734,6 +750,14 @@
 <param name="Local" value="optIC.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="OptionalInfluenceCurve-class">
+<param name="Local" value="OptionalInfluenceCurve-Class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="OptionalpICList-class">
+<param name="Local" value="OptionalInfluenceCurve-Class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="outlyingPlotIC">
 <param name="Local" value="outlyingPlotIC.html">
 </OBJECT>
@@ -746,6 +770,18 @@
 <param name="Local" value="ALEstimate-class.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="pICList">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="pICList,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="pICList-class">
+<param name="Local" value="OptionalInfluenceCurve-Class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="plot">
 <param name="Local" value="plot-methods.html">
 </OBJECT>
@@ -886,6 +922,26 @@
 <param name="Local" value="TotalVarIC-class.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="start">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="start,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="StartClass-class">
+<param name="Local" value="OptionalInfluenceCurve-Class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="startval">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="startval,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="steps">
 <param name="Local" value="kStepEstimate-class.html">
 </OBJECT>
@@ -922,10 +978,26 @@
 <param name="Local" value="Neighborhood-class.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="uksteps">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="uksteps,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="UncondNeighborhood-class">
 <param name="Local" value="UncondNeighborhood-class.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ustartval">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ustartval,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="weight">
 <param name="Local" value="RobWeight-class.html">
 </OBJECT>
@@ -1143,6 +1215,10 @@
 <param name="Local" value="HampelWeight-class.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Some helper Classes in package 'RobAStBase'">
+<param name="Local" value="OptionalInfluenceCurve-Class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="Total variation neighborhood">
 <param name="Local" value="TotalVarNeighborhood-class.html">
 </OBJECT>

Modified: branches/robast-0.7/pkg/RobAStBase/chm/RobAStBaseOptions.html
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/chm/RobAStBaseOptions.html	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/chm/RobAStBaseOptions.html	2009-08-31 17:04:45 UTC (rev 359)
@@ -70,6 +70,10 @@
 <dt>all.verbose:</dt><dd>argument <code>verbose</code> passed on by default to many calls of
 <code>optIC</code>, <code>radiusminimaxIC</code>, <code>getinfRobIC</code> etc.; well suited
 for testing purposes.</dd>
+<dt>withPICList:</dt><dd>logical: shall slot <code>pICList</code> of return value of <code>kStepEstimator</code>
+be filled?</dd>
+<dt>withICList:</dt><dd>logical: shall slot <code>ICList</code> of return value of <code>kStepEstimator</code>
+be filled?</dd>
 </dl>
 
 </p>

Modified: branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimate-class.html
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimate-class.html	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimate-class.html	2009-08-31 17:04:45 UTC (rev 359)
@@ -8,6 +8,20 @@
 <param name="keyword" value="R:   kStepEstimate-class">
 <param name="keyword" value="R:   steps">
 <param name="keyword" value="R:   steps,kStepEstimate-method">
+<param name="keyword" value="R:   ksteps">
+<param name="keyword" value="R:   ksteps,kStepEstimate-method">
+<param name="keyword" value="R:   uksteps">
+<param name="keyword" value="R:   uksteps,kStepEstimate-method">
+<param name="keyword" value="R:   ICList">
+<param name="keyword" value="R:   ICList,kStepEstimate-method">
+<param name="keyword" value="R:   pICList">
+<param name="keyword" value="R:   pICList,kStepEstimate-method">
+<param name="keyword" value="R:   start">
+<param name="keyword" value="R:   start,kStepEstimate-method">
+<param name="keyword" value="R:   startval">
+<param name="keyword" value="R:   startval,kStepEstimate-method">
+<param name="keyword" value="R:   ustartval">
+<param name="keyword" value="R:   ustartval,kStepEstimate-method">
 <param name="keyword" value="R:   show,kStepEstimate-method">
 <param name="keyword" value=" kStepEstimate-class.">
 </object>
@@ -68,6 +82,28 @@
 <dt><code>untransformed.asvar</code>:</dt><dd> object of class <code>"OptionalNumericOrMatrix"</code>
 which may contain the asymptotic (co)variance of the untransformed
 estimator. </dd>
+<dt><code>pICList</code></dt><dd>Optional object of class <code>"OptionalpICList"</code>:
+the list of (intermediate) (partial) influence curves used;
+only filled when called from <code>kStepEstimator</code> with argument
+<code>withPICList==TRUE</code>. </dd>
+<dt><code>ICList</code></dt><dd>Optional object of class <code>"OptionalpICList"</code>:
+the list of (intermediate) (total) influence curves used;
+only filled when called from <code>kStepEstimator</code> with argument
+<code>withICList==TRUE</code>. </dd>
+<dt><code>start</code></dt><dd>The argument <code>start</code> &mdash; of class <code>"StartClass"</code>
+used in call to <code>kStepEstimator</code>. </dd>
+<dt><code>startval</code></dt><dd>Object of class <code>matrix</code>:
+the starting value with which the k-step Estimator was initialized
+(in <i>p</i>-space / transformed). </dd>
+<dt><code>ustartval</code></dt><dd>Object of class <code>matrix</code>:
+the starting value with which the k-step Estimator was initialized
+(in <i>k</i>-space / untransformed). </dd>
+<dt><code>ksteps</code></dt><dd>Object of class <code>"OptionalMatrix"</code>:
+the intermediate estimates (in <i>p</i>-space) for the parameter;
+only filled when called from <code>kStepEstimator</code>. </dd>
+<dt><code>uksteps</code></dt><dd>Object of class <code>"OptionalMatrix"</code>:
+the intermediate estimates (in <i>k</i>-space) for the parameter;
+only filled when called from <code>kStepEstimator</code>. </dd>
 </dl>
 
 </p>
@@ -90,6 +126,40 @@
 accessor function for slot <code>steps</code>. </dd>
 </p>
 <p>
+<dt>ksteps</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>ksteps</code>; has additional argument
+<code>diff</code>, defaulting to <code>FALSE</code>; if the latter is <code>TRUE</code>,
+the starting value from slot <code>startval</code> is prepended as first column;
+otherwise we return the corresponding increments in each step. </dd>
+</p>
+<p>
+<dt>uksteps</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>uksteps</code>; has additional argument
+<code>diff</code>, defaulting to <code>FALSE</code>; if the latter is <code>TRUE</code>,
+the starting value from slot <code>ustartval</code> is prepended as first column;
+otherwise we return the corresponding increments in each step. </dd>
+</p>
+<p>
+<dt>start</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>start</code>. </dd>
+</p>
+<p>
+<dt>startval</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>startval</code>. </dd>
+</p>
+<p>
+<dt>ustartval</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>startval</code>. </dd>
+</p>
+<p>
+<dt>ICList</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>ICList</code>. </dd>
+</p>
+<p>
+<dt>pICList</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>pICList</code>. </dd>
+</p>
+<p>
 <dt>show</dt><dd><code>signature(object = "kStepEstimate")</code>: a show method; </dd>
 </dl>
 

Modified: branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimator.html
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimator.html	2009-08-31 11:32:11 UTC (rev 358)
+++ branches/robast-0.7/pkg/RobAStBase/chm/kStepEstimator.html	2009-08-31 17:04:45 UTC (rev 359)
@@ -27,6 +27,8 @@
       useLast = getRobAStBaseOption("kStepUseLast"),
       withUpdateInKer = getRobAStBaseOption("withUpdateInKer"),
       IC.UpdateInKer = getRobAStBaseOption("IC.UpdateInKer"),
+      withICList = getRobAStBaseOption("withICList"),
+      withPICList = getRobAStBaseOption("withPICList"),
       na.rm = TRUE, startArgList = NULL, ...)
 </pre>
 
@@ -72,6 +74,14 @@
 a list of arguments to be given to argument <code>start</code> if the latter
 is a function; this list by default already starts with two unnamed items,
 the sample <code>x</code>, and the model <code>eval(CallL2Fam(IC))</code>.</td></tr>
+<tr valign="top"><td><code>withPICList</code></td>
+<td>
+logical: shall slot <code>pICList</code> of return value
+be filled?</td></tr>
+<tr valign="top"><td><code>withICList</code></td>
+<td>
+logical: shall slot <code>ICList</code> of return value
+be filled?</td></tr>
 <tr valign="top"><td><code>...</code></td>
 <td>
  additional parameters </td></tr>

Added: branches/robast-0.7/pkg/RobAStBase/man/OptionalInfluenceCurve-Class.Rd
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/man/OptionalInfluenceCurve-Class.Rd	                        (rev 0)
+++ branches/robast-0.7/pkg/RobAStBase/man/OptionalInfluenceCurve-Class.Rd	2009-08-31 17:04:45 UTC (rev 359)
@@ -0,0 +1,41 @@
+\name{OptionalInfluenceCurve-class}
+\docType{class}
+\alias{OptionalInfluenceCurve-class}
+\alias{OptionalpICList-class}
+\alias{StartClass-class}
+\alias{pICList-class}
+
+\title{Some helper Classes in package 'RobAStBase'}
+\description{Some helper Classes in package 'RobAStBase':
+Classes \code{OptionalInfluenceCurve}, \code{OptionalpICList},
+\code{StartClass}, \code{pICList}}
+\section{Class Unions}{
+  \code{OptionalInfluenceCurve} is a class union of classes
+  \code{InfluenceCurve} and \code{NULL} --- it is the slot
+  class of slot \code{pIC} in \code{ALEstimate};
+  \code{OptionalpICList} is a class union of classes
+  \code{pICList} and \code{NULL} --- it is the slot
+  class of slot \code{pICList} in \code{kStepEstimate};
+  \code{StartClass} is a class union of classes
+  \code{function}, \code{numeric} and \code{Estimate} --- it is the slot
[TRUNCATED]

To get the complete diff run:
    svnlook diff /svnroot/robast -r 359


More information about the Robast-commits mailing list