[Robast-commits] r168 - in branches/robast-0.6/pkg/ROptEst: chm inst/scripts

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Thu Aug 28 10:23:17 CEST 2008


Author: ruckdeschel
Date: 2008-08-28 10:23:16 +0200 (Thu, 28 Aug 2008)
New Revision: 168

Modified:
   branches/robast-0.6/pkg/ROptEst/chm/00Index.html
   branches/robast-0.6/pkg/ROptEst/chm/ROptEst.chm
   branches/robast-0.6/pkg/ROptEst/chm/ROptEst.toc
   branches/robast-0.6/pkg/ROptEst/chm/getBiasIC.html
   branches/robast-0.6/pkg/ROptEst/chm/getRiskIC.html
   branches/robast-0.6/pkg/ROptEst/inst/scripts/GammaModel.R
   branches/robast-0.6/pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R
Log:
yet another small changes in GammaModel -- all the script took about 20min on my Laptop
and also similar change in NormalLocationScaleModel.R


Modified: branches/robast-0.6/pkg/ROptEst/chm/00Index.html
===================================================================
--- branches/robast-0.6/pkg/ROptEst/chm/00Index.html	2008-08-27 19:59:47 UTC (rev 167)
+++ branches/robast-0.6/pkg/ROptEst/chm/00Index.html	2008-08-28 08:23:16 UTC (rev 168)
@@ -62,6 +62,8 @@
 <td>Generic function for the computation of the asymptotic bias for an IC</td></tr>
 <tr><td width="25%"><a href="getBiasIC.html">getBiasIC,HampIC,UncondNeighborhood-method</a></td>
 <td>Generic function for the computation of the asymptotic bias for an IC</td></tr>
+<tr><td width="25%"><a href="getBiasIC.html">getBiasIC,TotalVarIC,UncondNeighborhood-method</a></td>
+<td>Generic function for the computation of the asymptotic bias for an IC</td></tr>
 <tr><td width="25%"><a href="getBiasIC.html">getBiasIC-methods</a></td>
 <td>Generic function for the computation of the asymptotic bias for an IC</td></tr>
 <tr><td width="25%"><a href="getFiRisk.html">getFiRisk</a></td>
@@ -202,6 +204,8 @@
 <td>Generic Function for the Computation of Functions for Slot modifyIC</td></tr>
 <tr><td width="25%"><a href="getModifyIC.html">getModifyIC,L2LocationFamily,UncondNeighborhood,asGRisk-method</a></td>
 <td>Generic Function for the Computation of Functions for Slot modifyIC</td></tr>
+<tr><td width="25%"><a href="getModifyIC.html">getModifyIC,L2LocationFamily,UncondNeighborhood,fiUnOvShoot-method</a></td>
+<td>Generic Function for the Computation of Functions for Slot modifyIC</td></tr>
 <tr><td width="25%"><a href="getModifyIC.html">getModifyIC,L2LocationScaleFamily,ContNeighborhood,asGRisk-method</a></td>
 <td>Generic Function for the Computation of Functions for Slot modifyIC</td></tr>
 <tr><td width="25%"><a href="getModifyIC.html">getModifyIC,L2ParamFamily,Neighborhood,asRisk-method</a></td>
@@ -218,6 +222,8 @@
 <td>Generic function for the computation of a risk for an IC</td></tr>
 <tr><td width="25%"><a href="getRiskIC.html">getRiskIC,HampIC,asCov,missing,missing-method</a></td>
 <td>Generic function for the computation of a risk for an IC</td></tr>
+<tr><td width="25%"><a href="getRiskIC.html">getRiskIC,TotalVarIC,asCov,missing,L2ParamFamily-method</a></td>
+<td>Generic function for the computation of a risk for an IC</td></tr>
 <tr><td width="25%"><a href="getRiskIC.html">getRiskIC-methods</a></td>
 <td>Generic function for the computation of a risk for an IC</td></tr>
 </table>

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

Modified: branches/robast-0.6/pkg/ROptEst/chm/ROptEst.toc
===================================================================
--- branches/robast-0.6/pkg/ROptEst/chm/ROptEst.toc	2008-08-27 19:59:47 UTC (rev 167)
+++ branches/robast-0.6/pkg/ROptEst/chm/ROptEst.toc	2008-08-28 08:23:16 UTC (rev 168)
@@ -82,6 +82,10 @@
 <param name="Local" value="getBiasIC.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getBiasIC,TotalVarIC,UncondNeighborhood-method">
+<param name="Local" value="getBiasIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="getBiasIC-methods">
 <param name="Local" value="getBiasIC.html">
 </OBJECT>
@@ -362,6 +366,10 @@
 <param name="Local" value="getModifyIC.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC,L2LocationFamily,UncondNeighborhood,fiUnOvShoot-method">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="getModifyIC,L2LocationScaleFamily,ContNeighborhood,asGRisk-method">
 <param name="Local" value="getModifyIC.html">
 </OBJECT>
@@ -394,6 +402,10 @@
 <param name="Local" value="getRiskIC.html">
 </OBJECT>
 <LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getRiskIC,TotalVarIC,asCov,missing,L2ParamFamily-method">
+<param name="Local" value="getRiskIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
 <param name="Name" value="getRiskIC-methods">
 <param name="Local" value="getRiskIC.html">
 </OBJECT>

Modified: branches/robast-0.6/pkg/ROptEst/chm/getBiasIC.html
===================================================================
--- branches/robast-0.6/pkg/ROptEst/chm/getBiasIC.html	2008-08-27 19:59:47 UTC (rev 167)
+++ branches/robast-0.6/pkg/ROptEst/chm/getBiasIC.html	2008-08-28 08:23:16 UTC (rev 168)
@@ -8,6 +8,7 @@
 <param name="keyword" value="R:   getBiasIC">
 <param name="keyword" value="R:   getBiasIC-methods">
 <param name="keyword" value="R:   getBiasIC,HampIC,UncondNeighborhood-method">
+<param name="keyword" value="R:   getBiasIC,TotalVarIC,UncondNeighborhood-method">
 <param name="keyword" value=" Generic function for the computation of the asymptotic bias for an IC">
 </object>
 
@@ -66,6 +67,9 @@
 
 <dl>
 <dt>IC = "HampIC", neighbor = "UncondNeighborhood"</dt><dd>reads off the as. bias from the risks-slot of the IC. </dd>
+<dt>IC = "TotalVarIC", neighbor = "UncondNeighborhood"</dt><dd>reads off the as. bias from the risks-slot of the IC,
+resp. if this is <code>NULL</code> from the corresponding
+Lagrange Multipliers. </dd>
 </dl>
 
 <h3>Note</h3>

Modified: branches/robast-0.6/pkg/ROptEst/chm/getRiskIC.html
===================================================================
--- branches/robast-0.6/pkg/ROptEst/chm/getRiskIC.html	2008-08-27 19:59:47 UTC (rev 167)
+++ branches/robast-0.6/pkg/ROptEst/chm/getRiskIC.html	2008-08-28 08:23:16 UTC (rev 168)
@@ -9,6 +9,7 @@
 <param name="keyword" value="R:   getRiskIC-methods">
 <param name="keyword" value="R:   getRiskIC,HampIC,asCov,missing,missing-method">
 <param name="keyword" value="R:   getRiskIC,HampIC,asCov,missing,L2ParamFamily-method">
+<param name="keyword" value="R:   getRiskIC,TotalVarIC,asCov,missing,L2ParamFamily-method">
 <param name="keyword" value=" Generic function for the computation of a risk for an IC">
 </object>
 
@@ -35,6 +36,9 @@
 ## S4 method for signature 'HampIC, asCov, missing,
 ##   L2ParamFamily':
 getRiskIC(IC, risk, L2Fam)
+## S4 method for signature 'TotalVarIC, asCov, missing,
+##   L2ParamFamily':
+getRiskIC(IC, risk, L2Fam)
 
 </pre>
 
@@ -80,6 +84,8 @@
 
 
 <dt>IC = "HampIC", risk = "asCov", neighbor = "missing", L2Fam = "L2ParamFamily"</dt><dd>asymptotic covariance of <code>IC</code> under <code>L2Fam</code> read off from corresp. <code>Risks</code> slot. </dd>
+<dt>IC = "TotalVarIC", risk = "asCov", neighbor = "missing", L2Fam = "L2ParamFamily"</dt><dd>asymptotic covariance of <code>IC</code> read off from corresp. <code>Risks</code> slot,
+resp. if this is <code>NULL</code> calculates it via <code><a href="getInfV.html">getInfV</a></code>.</dd>
 </dl>
 
 <h3>Note</h3>

Modified: branches/robast-0.6/pkg/ROptEst/inst/scripts/GammaModel.R
===================================================================
--- branches/robast-0.6/pkg/ROptEst/inst/scripts/GammaModel.R	2008-08-27 19:59:47 UTC (rev 167)
+++ branches/robast-0.6/pkg/ROptEst/inst/scripts/GammaModel.R	2008-08-28 08:23:16 UTC (rev 168)
@@ -71,7 +71,8 @@
 MLEstimator(x=x, GammaFamily())
 
 ## 3. one-step estimation: radius known
-RobG3 <- InfRobModel(center=GammaFamily(scale = est0$estimate[1], shape = est0$estimate[2]), 
+RobG3 <- InfRobModel(center=GammaFamily(scale = estimate(est0)["scale"], 
+                     shape = estimate(est0)["shape"]), 
                 neighbor=ContNeighborhood(radius=0.5))
 IC9 <- optIC(model=RobG3, risk=asMSE())
 (est1 <- oneStepEstimator(x, IC=IC9, start=est0$estimate))

Modified: branches/robast-0.6/pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R
===================================================================
--- branches/robast-0.6/pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R	2008-08-27 19:59:47 UTC (rev 167)
+++ branches/robast-0.6/pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R	2008-08-28 08:23:16 UTC (rev 168)
@@ -112,7 +112,7 @@
 (est0 <- MDEstimator(x=x, NormLocationScaleFamily()))
 
 ## 3. k-step estimation: radius known
-N1 <- NormLocationScaleFamily(mean=estimate(est0)[1], sd=estimate(est0)[2])
+N1 <- NormLocationScaleFamily(mean=estimate(est0)["mean"], sd=estimate(est0)["sd"])
 N1.Rob <- InfRobModel(center = N1, neighbor = ContNeighborhood(radius = 0.5))
 IC1 <- optIC(model = N1.Rob, risk = asMSE())
 (est1 <- kStepEstimator(x, IC1, est0, steps = 3))



More information about the Robast-commits mailing list