[Robast-commits] r170 - in pkg: ROptEst/R ROptEst/chm ROptEst/inst/scripts ROptEst/man RandVar/R RobAStBase/R RobAStBase/chm RobLox/chm
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Thu Aug 28 14:46:24 CEST 2008
Author: ruckdeschel
Date: 2008-08-28 14:46:23 +0200 (Thu, 28 Aug 2008)
New Revision: 170
Modified:
pkg/ROptEst/R/getModifyIC.R
pkg/ROptEst/R/getRiskIC.R
pkg/ROptEst/R/optIC.R
pkg/ROptEst/chm/00Index.html
pkg/ROptEst/chm/ROptEst.chm
pkg/ROptEst/chm/ROptEst.hhp
pkg/ROptEst/chm/ROptEst.toc
pkg/ROptEst/chm/getAsRisk.html
pkg/ROptEst/chm/getBiasIC.html
pkg/ROptEst/chm/getFixRobIC.html
pkg/ROptEst/chm/getIneffDiff.html
pkg/ROptEst/chm/getInfCent.html
pkg/ROptEst/chm/getInfClip.html
pkg/ROptEst/chm/getInfRobIC.html
pkg/ROptEst/chm/getInfStand.html
pkg/ROptEst/chm/getL1normL2deriv.html
pkg/ROptEst/chm/getRiskIC.html
pkg/ROptEst/chm/leastFavorableRadius.html
pkg/ROptEst/chm/minmaxBias.html
pkg/ROptEst/chm/optIC.html
pkg/ROptEst/chm/radiusMinimaxIC.html
pkg/ROptEst/inst/scripts/GammaModel.R
pkg/ROptEst/inst/scripts/GumbelLocationModel.R
pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R
pkg/ROptEst/inst/scripts/UnderOverShootRisk.R
pkg/ROptEst/man/getBiasIC.Rd
pkg/ROptEst/man/getModifyIC.Rd
pkg/ROptEst/man/getRiskIC.Rd
pkg/RandVar/R/util.R
pkg/RobAStBase/R/locMEstimator.R
pkg/RobAStBase/chm/00Index.html
pkg/RobAStBase/chm/IC-class.html
pkg/RobAStBase/chm/RobAStBase.chm
pkg/RobAStBase/chm/RobAStBase.toc
pkg/RobLox/chm/RobLox.chm
Log:
merged branch robast-0.6 into trunk
Modified: pkg/ROptEst/R/getModifyIC.R
===================================================================
--- pkg/ROptEst/R/getModifyIC.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/R/getModifyIC.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -27,6 +27,12 @@
return(modIC)
})
+setMethod("getModifyIC", signature(L2FamIC = "L2LocationFamily",
+ neighbor = "UncondNeighborhood", risk = "fiUnOvShoot"),
+ getMethod("getModifyIC",signature(L2FamIC = "L2LocationFamily",
+ neighbor = "UncondNeighborhood", risk = "asGRisk"))
+ )
+
setMethod("getModifyIC", signature(L2FamIC = "L2ScaleFamily",
neighbor = "ContNeighborhood", risk = "asGRisk"),
function(L2FamIC, neighbor, risk){
Modified: pkg/ROptEst/R/getRiskIC.R
===================================================================
--- pkg/ROptEst/R/getRiskIC.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/R/getRiskIC.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -19,6 +19,23 @@
return(list(asCov = list(distribution = .getDistr(L2Fam), value = Cov)))
})
+setMethod("getRiskIC", signature(IC = "TotalVarIC",
+ risk = "asCov",
+ neighbor = "missing",
+ L2Fam = "L2ParamFamily"),
+ function(IC, risk, L2Fam){
+ Cov <- IC at Risks[["asCov"]]
+ if (is.null(Cov)){
+ L2deriv <- L2Fam at L2derivDistr[[1]]
+ A <- IC at stand
+ c0 <- (IC at clipUp-IC@clipLo)/A
+ z <- IC at clipLo/A
+ neighbor <- TotalVarNeighborhood(1)
+ Cov <- getInfV(L2deriv = L2deriv, neighbor = neighbor,
+ biastype = biastype(IC), clip = c0, cent = z, stand = A)
+ }
+ return(list(asCov = list(distribution = .getDistr(L2Fam), value = Cov)))
+ })
###############################################################################
## asymptotic Bias for various types
@@ -27,10 +44,27 @@
neighbor = "UncondNeighborhood"),
function(IC, neighbor, L2Fam,...){
if(missing(L2Fam))
- L2Fam <- eval(IC at CallL2Fam)
+ L2Fam <- force(eval(IC at CallL2Fam))
+ Bias <- IC at Risks$asBias$value
+
return(list(asBias = list(distribution = .getDistr(L2Fam),
- neighborhood = neighbor at type, value = IC at Risks$asBias$value)))
+ neighborhood = neighbor at type, value = Bias)))
})
+setMethod("getBiasIC", signature(IC = "TotalVarIC",
+ neighbor = "UncondNeighborhood"),
+ function(IC, neighbor, L2Fam,...){
+ if(missing(L2Fam))
+ L2Fam <- force(eval(IC at CallL2Fam))
+ Bias <- IC at Risks$asBias$value
+ if (is.null(Bias)){
+ Bias <- if(is(neighbor,"ContNeighborhood"))
+ max(IC at clipUp,-IC at clipLo) else IC at clipUp-IC@clipLo
+ }
+
+ return(list(asBias = list(distribution = .getDistr(L2Fam),
+ neighborhood = neighbor at type, value = Bias)))
+ })
+
Modified: pkg/ROptEst/R/optIC.R
===================================================================
--- pkg/ROptEst/R/optIC.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/R/optIC.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -72,22 +72,26 @@
tol = .Machine$double.eps^0.4, warn = TRUE){
L2derivDistr <- model at center@L2derivDistr[[1]]
if((length(model at center@L2derivDistr) == 1) & is(L2derivDistr, "UnivariateDistribution")){
- ow <- options("warn")
- options(warn = -1)
- res <- getInfRobIC(L2deriv = L2derivDistr,
+ if(identical(all.equal(model at neighbor@radius, 0), TRUE)){
+ return(optIC(model at center, risk = asCov()))
+ }else{
+ ow <- options("warn")
+ options(warn = -1)
+ res <- getInfRobIC(L2deriv = L2derivDistr,
neighbor = model at neighbor, risk = risk,
symm = model at center@L2derivDistrSymm[[1]],
Finfo = model at center@FisherInfo, trafo = model at center@param at trafo,
upper = upper, maxiter = maxiter, tol = tol, warn = warn)
- options(ow)
- if(is(model at neighbor, "ContNeighborhood"))
- res$info <- c("optIC", "optIC", res$info, "Optimal IC for 'InfRobModel' with 'ContNeighborhood'!!!")
- else
- res$info <- c("optIC", res$info)
- res <- c(res, modifyIC = getModifyIC(L2FamIC = model at center,
- neighbor = model at neighbor,
- risk = risk))
- return(generateIC(TotalVarNeighborhood(radius = model at neighbor@radius), model at center, res))
+ options(ow)
+ if(is(model at neighbor, "ContNeighborhood"))
+ res$info <- c("optIC", "optIC", res$info, "Optimal IC for 'InfRobModel' with 'ContNeighborhood'!!!")
+ else
+ res$info <- c("optIC", res$info)
+ res <- c(res, modifyIC = getModifyIC(L2FamIC = model at center,
+ neighbor = model at neighbor,
+ risk = risk))
+ return(generateIC(TotalVarNeighborhood(radius = model at neighbor@radius), model at center, res))
+ }
}else{
stop("restricted to 1-dimensional parameteric models")
}
Modified: pkg/ROptEst/chm/00Index.html
===================================================================
--- pkg/ROptEst/chm/00Index.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/00Index.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -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>
@@ -198,12 +200,30 @@
<td>Calculation of L1 norm of L2derivative</td></tr>
<tr><td width="25%"><a href="getL2normL2deriv.html">getL2normL2deriv</a></td>
<td>Calculation of L2 norm of L2derivative</td></tr>
+<tr><td width="25%"><a href="getModifyIC.html">getModifyIC</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,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>
+<td>Generic Function for the Computation of Functions for Slot modifyIC</td></tr>
+<tr><td width="25%"><a href="getModifyIC.html">getModifyIC,L2ScaleFamily,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,L2ScaleFamily,TotalVarNeighborhood,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-methods</a></td>
+<td>Generic Function for the Computation of Functions for Slot modifyIC</td></tr>
<tr><td width="25%"><a href="getRiskIC.html">getRiskIC</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,HampIC,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,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: pkg/ROptEst/chm/ROptEst.chm
===================================================================
(Binary files differ)
Modified: pkg/ROptEst/chm/ROptEst.hhp
===================================================================
--- pkg/ROptEst/chm/ROptEst.hhp 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/ROptEst.hhp 2008-08-28 12:46:23 UTC (rev 170)
@@ -26,6 +26,7 @@
getInfV.html
getL1normL2deriv.html
getL2normL2deriv.html
+getModifyIC.html
getRiskIC.html
leastFavorableRadius.html
lowerCaseRadius.html
Modified: pkg/ROptEst/chm/ROptEst.toc
===================================================================
--- pkg/ROptEst/chm/ROptEst.toc 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/ROptEst.toc 2008-08-28 12:46:23 UTC (rev 170)
@@ -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>
@@ -354,6 +358,38 @@
<param name="Local" value="getL2normL2deriv.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC,L2LocationFamily,UncondNeighborhood,asGRisk-method">
+<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>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC,L2ParamFamily,Neighborhood,asRisk-method">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC,L2ScaleFamily,ContNeighborhood,asGRisk-method">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC,L2ScaleFamily,TotalVarNeighborhood,asGRisk-method">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="getModifyIC-methods">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="getRiskIC">
<param name="Local" value="getRiskIC.html">
</OBJECT>
@@ -366,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>
@@ -539,6 +579,10 @@
<param name="Local" value="minmaxBias.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Generic Function for the Computation of Functions for Slot modifyIC">
+<param name="Local" value="getModifyIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Generic Function for the Computation of Inefficiency Differences">
<param name="Local" value="getIneffDiff.html">
</OBJECT>
Modified: pkg/ROptEst/chm/getAsRisk.html
===================================================================
--- pkg/ROptEst/chm/getAsRisk.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getAsRisk.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -44,61 +44,75 @@
## S4 method for signature 'asMSE, UnivariateDistribution,
## Neighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, clip, cent, stand, trafo)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, clip, cent, stand, trafo)
## S4 method for signature 'asMSE, EuclRandVariable,
## Neighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, clip, cent, stand, trafo)
+getAsRisk(risk, L2deriv, neighbor,
+ biastype, clip, cent, stand, trafo)
## S4 method for signature 'asBias, UnivariateDistribution,
## ContNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, trafo)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, trafo)
## S4 method for signature 'asBias, UnivariateDistribution,
## ContNeighborhood, onesidedBias':
-getAsRisk(risk, L2deriv, neighbor, biastype, trafo)
+getAsRisk(risk,
+ L2deriv, neighbor, biastype, trafo)
## S4 method for signature 'asBias, UnivariateDistribution,
## ContNeighborhood, asymmetricBias':
-getAsRisk(risk, L2deriv, neighbor, biastype, trafo)
+getAsRisk(risk,
+ L2deriv, neighbor, biastype, trafo)
## S4 method for signature 'asBias, UnivariateDistribution,
## TotalVarNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, trafo)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, trafo)
## S4 method for signature 'asBias, RealRandVariable,
## ContNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, Distr, DistrSymm, L2derivSymm,
- L2derivDistrSymm, trafo, z.start, A.start, maxiter, tol, warn)
+getAsRisk(risk, L2deriv, neighbor,
+ biastype, Distr, DistrSymm, L2derivSymm, L2derivDistrSymm,
+ trafo, z.start, A.start, maxiter, tol, warn)
## S4 method for signature 'asCov, UnivariateDistribution,
## ContNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, clip, cent, stand)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, clip, cent, stand)
## S4 method for signature 'asCov, UnivariateDistribution,
## TotalVarNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, clip, cent, stand)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, clip, cent, stand)
## S4 method for signature 'asCov, RealRandVariable,
## ContNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, Distr, clip, cent, stand)
+getAsRisk(risk, L2deriv, neighbor,
+ biastype, Distr, cent, stand,
+ V.comp = matrix(TRUE, ncol = nrow(stand), nrow = nrow(stand)), w)
## S4 method for signature 'trAsCov,
## UnivariateDistribution, UncondNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, clip, cent, stand)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, clip, cent, stand)
## S4 method for signature 'trAsCov, RealRandVariable,
## ContNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, Distr, clip, cent, stand, normtype)
+getAsRisk(risk, L2deriv, neighbor,
+ biastype, Distr, clip, cent, stand, normtype)
## S4 method for signature 'asUnOvShoot,
## UnivariateDistribution, UncondNeighborhood, ANY':
-getAsRisk(risk, L2deriv, neighbor, biastype, clip, cent, stand, trafo)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, clip, cent, stand, trafo)
## S4 method for signature 'asSemivar,
## UnivariateDistribution, Neighborhood, onesidedBias':
-getAsRisk(risk, L2deriv, neighbor, biastype,
- clip, cent, stand, trafo)
+getAsRisk(risk, L2deriv,
+ neighbor, biastype, clip, cent, stand, trafo)
</pre>
@@ -163,6 +177,13 @@
<tr valign="top"><td><code>normtype</code></td>
<td>
object of class <code>"NormType"</code>. </td></tr>
+<tr valign="top"><td><code>V.comp</code></td>
+<td>
+matrix: indication which components of the standardizing
+matrix have to be computed. </td></tr>
+<tr valign="top"><td><code>w</code></td>
+<td>
+object of class <code>RobWeight</code>; current weight</td></tr>
</table>
<h3>Details</h3>
Modified: pkg/ROptEst/chm/getBiasIC.html
===================================================================
--- pkg/ROptEst/chm/getBiasIC.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getBiasIC.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -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>
@@ -28,7 +29,7 @@
getBiasIC(IC, neighbor, ...)
## S4 method for signature 'HampIC, UncondNeighborhood':
-getBiasIC(IC, neighbor, L2Fam)
+getBiasIC(IC, neighbor, L2Fam, ...)
</pre>
@@ -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: pkg/ROptEst/chm/getFixRobIC.html
===================================================================
--- pkg/ROptEst/chm/getFixRobIC.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getFixRobIC.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -74,6 +74,16 @@
"left" or "right". </td></tr>
</table>
+<h3>Details</h3>
+
+<p>
+Computation of the optimally robust IC in sense of Huber (1968) which
+is also treated in Kohl (2005). The Algorithm used to compute the exact
+finite sample risk is introduced and explained in Kohl (2005). It is
+based on FFT.
+</p>
+
+
<h3>Value</h3>
<p>
@@ -100,6 +110,9 @@
Verw. Geb. <B>10</B>:269–278.
</p>
<p>
+Rieder, H. (1980) Estimates derived from robust tests. Ann. Stats. <B>8</B>: 106-115.
+</p>
+<p>
Kohl, M. (2005) <EM>Numerical Contributions to the Asymptotic Theory of Robustness</EM>.
Bayreuth: Dissertation.
</p>
Modified: pkg/ROptEst/chm/getIneffDiff.html
===================================================================
--- pkg/ROptEst/chm/getIneffDiff.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getIneffDiff.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -34,7 +34,7 @@
getIneffDiff(
radius, L2Fam, neighbor, risk, loRad, upRad, loRisk, upRisk,
z.start = NULL, A.start = NULL, upper.b, MaxIter, eps, warn,
- loNorm = NULL, upNorm = NULL)
+ loNorm = NULL, upNorm = NULL, verbose = FALSE)
</pre>
@@ -96,6 +96,9 @@
object of class <code>"NormType"</code>; used in selfstandardization
to evaluate the bias of the current IC in the norm of the upper
bound</td></tr>
+<tr valign="top"><td><code>verbose</code></td>
+<td>
+logical: if <code>TRUE</code>, some messages are printed </td></tr>
</table>
<h3>Value</h3>
@@ -122,6 +125,10 @@
<h3>References</h3>
<p>
+Rieder, H., Kohl, M. and Ruckdeschel, P. (2008) The Costs of not Knowing
+the Radius. Statistical Methods and Applications, <EM>17</EM>(1) 13-40.
+</p>
+<p>
Rieder, H., Kohl, M. and Ruckdeschel, P. (2001) The Costs of not Knowing
the Radius. Submitted. Appeared as discussion paper Nr. 81.
SFB 373 (Quantification and Simulation of Economic Processes),
Modified: pkg/ROptEst/chm/getInfCent.html
===================================================================
--- pkg/ROptEst/chm/getInfCent.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getInfCent.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -48,7 +48,7 @@
## S4 method for signature 'RealRandVariable,
## ContNeighborhood, BiasType':
getInfCent(L2deriv,
- neighbor, biastype, z.comp, stand, cent, clip, w)
+ neighbor, biastype, Distr, z.comp, w)
## S4 method for signature 'UnivariateDistribution,
## ContNeighborhood, onesidedBias':
@@ -84,9 +84,6 @@
<tr valign="top"><td><code>cent</code></td>
<td>
optimal centering constant. </td></tr>
-<tr valign="top"><td><code>stand</code></td>
-<td>
-standardizing matrix. </td></tr>
<tr valign="top"><td><code>tol.z</code></td>
<td>
the desired accuracy (convergence tolerance). </td></tr>
@@ -96,6 +93,9 @@
<tr valign="top"><td><code>trafo</code></td>
<td>
matrix: transformation of the parameter. </td></tr>
+<tr valign="top"><td><code>Distr</code></td>
+<td>
+object of class <code>Distribution</code>. </td></tr>
<tr valign="top"><td><code>z.comp</code></td>
<td>
logical vector: indication which components of the
Modified: pkg/ROptEst/chm/getInfClip.html
===================================================================
--- pkg/ROptEst/chm/getInfClip.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getInfClip.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -35,24 +35,29 @@
## S4 method for signature 'numeric,
## UnivariateDistribution, asMSE, ContNeighborhood':
-getInfClip(clip, L2deriv, risk, neighbor, biastype, cent, symm, trafo)
+getInfClip(clip, L2deriv,
+ risk, neighbor, biastype, cent, symm, trafo)
## S4 method for signature 'numeric,
## UnivariateDistribution, asMSE, TotalVarNeighborhood':
-getInfClip(clip, L2deriv, risk, neighbor, biastype, cent, symm, trafo)
+getInfClip(clip, L2deriv,
+ risk, neighbor, biastype, cent, symm, trafo)
## S4 method for signature 'numeric, EuclRandVariable,
## asMSE, ContNeighborhood':
-getInfClip(clip, L2deriv, risk, neighbor, Distr, stand, biastype, cent, trafo)
+getInfClip(clip, L2deriv, risk,
+ neighbor, biastype, Distr, stand, cent, trafo)
## S4 method for signature 'numeric,
## UnivariateDistribution, asUnOvShoot,
## UncondNeighborhood':
-getInfClip(clip, L2deriv, risk, neighbor, biastype, cent, symm, trafo)
+getInfClip(clip, L2deriv,
+ risk, neighbor, biastype, cent, symm, trafo)
## S4 method for signature 'numeric,
## UnivariateDistribution, asSemivar, ContNeighborhood':
-getInfClip(clip, L2deriv, risk, neighbor, cent, symm, trafo)
+getInfClip(clip, L2deriv,
+ risk, neighbor, cent, symm, trafo)
</pre>
Modified: pkg/ROptEst/chm/getInfRobIC.html
===================================================================
--- pkg/ROptEst/chm/getInfRobIC.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getInfRobIC.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -40,46 +40,46 @@
## S4 method for signature 'UnivariateDistribution, asCov,
## ContNeighborhood':
-getInfRobIC(L2deriv, risk, neighbor, Finfo, trafo)
+getInfRobIC(L2deriv, risk, neighbor, Finfo, trafo, verbose = FALSE)
## S4 method for signature 'UnivariateDistribution, asCov,
## TotalVarNeighborhood':
-getInfRobIC(L2deriv, risk, neighbor, Finfo, trafo)
+getInfRobIC(L2deriv, risk, neighbor, Finfo, trafo, verbose = FALSE)
## S4 method for signature 'RealRandVariable, asCov,
## ContNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, Distr, Finfo, trafo,
- QuadForm = diag(nrow(trafo)))
+ QuadForm = diag(nrow(trafo)), verbose = FALSE)
## S4 method for signature 'UnivariateDistribution, asBias,
## UncondNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, symm, trafo,
- maxiter, tol, warn, Finfo)
+ maxiter, tol, warn, Finfo, verbose = FALSE, ...)
## S4 method for signature 'RealRandVariable, asBias,
## ContNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, Distr, DistrSymm, L2derivSymm,
- L2derivDistrSymm, z.start, A.start, Finfo, trafo, maxiter, tol, warn)
+ L2derivDistrSymm, z.start, A.start, Finfo, trafo, maxiter, tol, warn, verbose = FALSE, ...)
## S4 method for signature 'UnivariateDistribution,
## asHampel, UncondNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, symm, Finfo, trafo,
- upper, maxiter, tol, warn, noLow = FALSE)
+ upper, maxiter, tol, warn, noLow = FALSE, verbose = FALSE)
## S4 method for signature 'RealRandVariable, asHampel,
## ContNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, Distr, DistrSymm, L2derivSymm,
- L2derivDistrSymm, Finfo, z.start, A.start, trafo, upper, maxiter, tol, warn)
+ L2derivDistrSymm, Finfo, trafo, onesetLM = FALSE, z.start, A.start, upper, maxiter, tol, warn, verbose = FALSE)
## S4 method for signature 'UnivariateDistribution,
## asGRisk, UncondNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, symm, Finfo, trafo,
- upper, maxiter, tol, warn, noLow = FALSE)
+ upper, maxiter, tol, warn, noLow = FALSE, verbose = FALSE)
## S4 method for signature 'RealRandVariable, asGRisk,
## ContNeighborhood':
getInfRobIC(L2deriv, risk, neighbor, Distr, DistrSymm, L2derivSymm,
- L2derivDistrSymm, Finfo, z.start, A.start, trafo, upper, maxiter, tol, warn)
+ L2derivDistrSymm, Finfo, trafo, onesetLM = FALSE, z.start, A.start, upper, maxiter, tol, warn, verbose = FALSE)
## S4 method for signature 'UnivariateDistribution,
## asUnOvShoot, UncondNeighborhood':
@@ -146,11 +146,17 @@
<tr valign="top"><td><code>noLow</code></td>
<td>
logical: is lower case to be computed? </td></tr>
+<tr valign="top"><td><code>onesetLM</code></td>
+<td>
+logical: use one set of Lagrange multipliers? </td></tr>
<tr valign="top"><td><code>QuadForm</code></td>
<td>
matrix of (or which may coerced to) class
<code>PosSemDefSymmMatrix</code> for use of different
(standardizing) norm </td></tr>
+<tr valign="top"><td><code>verbose</code></td>
+<td>
+logical: if <code>TRUE</code>, some messages are printed </td></tr>
</table>
<h3>Value</h3>
Modified: pkg/ROptEst/chm/getInfStand.html
===================================================================
--- pkg/ROptEst/chm/getInfStand.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getInfStand.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -46,7 +46,7 @@
## S4 method for signature 'RealRandVariable,
## ContNeighborhood, BiasType':
getInfStand(L2deriv,
- neighbor, biastype, Distr, A.comp, stand, clip, cent, trafo, w)
+ neighbor, biastype, Distr, A.comp, cent, trafo, w)
## S4 method for signature 'UnivariateDistribution,
## ContNeighborhood, BiasType':
@@ -82,9 +82,6 @@
<tr valign="top"><td><code>cent</code></td>
<td>
optimal centering constant. </td></tr>
-<tr valign="top"><td><code>stand</code></td>
-<td>
-standardizing matrix. </td></tr>
<tr valign="top"><td><code>Distr</code></td>
<td>
object of class <code>"Distribution"</code>. </td></tr>
Modified: pkg/ROptEst/chm/getL1normL2deriv.html
===================================================================
--- pkg/ROptEst/chm/getL1normL2deriv.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getL1normL2deriv.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -30,10 +30,9 @@
getL1normL2deriv(L2deriv,
cent, ...)
-## S4 method for signature 'UnivariateDistribution':
+## S4 method for signature 'RealRandVariable':
getL1normL2deriv(L2deriv,
cent, stand, Distr, normtype, ...)
-
</pre>
Modified: pkg/ROptEst/chm/getRiskIC.html
===================================================================
--- pkg/ROptEst/chm/getRiskIC.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/getRiskIC.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -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: pkg/ROptEst/chm/leastFavorableRadius.html
===================================================================
--- pkg/ROptEst/chm/leastFavorableRadius.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/leastFavorableRadius.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -32,7 +32,7 @@
leastFavorableRadius(
L2Fam, neighbor, risk, rho, upRad = 1,
z.start = NULL, A.start = NULL, upper = 100, maxiter = 100,
- tol = .Machine$double.eps^0.4, warn = FALSE)
+ tol = .Machine$double.eps^0.4, warn = FALSE, verbose = FALSE)
</pre>
@@ -76,6 +76,9 @@
<tr valign="top"><td><code>warn</code></td>
<td>
logical: print warnings. </td></tr>
+<tr valign="top"><td><code>verbose</code></td>
+<td>
+logical: if <code>TRUE</code>, some messages are printed </td></tr>
</table>
<h3>Value</h3>
Modified: pkg/ROptEst/chm/minmaxBias.html
===================================================================
--- pkg/ROptEst/chm/minmaxBias.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/minmaxBias.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -36,7 +36,7 @@
## S4 method for signature 'UnivariateDistribution,
## ContNeighborhood, BiasType':
minmaxBias(L2deriv, neighbor, biastype, symm, trafo,
- maxiter, tol, Finfo)
+ maxiter, tol, warn, Finfo)
## S4 method for signature 'UnivariateDistribution,
## ContNeighborhood, asymmetricBias':
@@ -55,8 +55,8 @@
## S4 method for signature 'RealRandVariable,
## ContNeighborhood, BiasType':
-minmaxBias(L2deriv, neighbor, biastype, Distr,
- L2derivDistrSymm, z.start, A.start, z.comp, A.comp, trafo, maxiter, tol, warn)
+minmaxBias(L2deriv, neighbor, biastype, normtype, Distr,
+ z.start, A.start, z.comp, A.comp, trafo, maxiter, tol)
</pre>
@@ -74,6 +74,9 @@
<tr valign="top"><td><code>biastype</code></td>
<td>
object of class <code>"BiasType"</code>. </td></tr>
+<tr valign="top"><td><code>normtype</code></td>
+<td>
+object of class <code>"NormType"</code>. </td></tr>
<tr valign="top"><td><code>...</code></td>
<td>
additional parameters. </td></tr>
@@ -83,9 +86,6 @@
<tr valign="top"><td><code>symm</code></td>
<td>
logical: indicating symmetry of <code>L2deriv</code>. </td></tr>
-<tr valign="top"><td><code>L2derivDistrSymm</code></td>
-<td>
-object of class <code>"DistrSymmList"</code>. </td></tr>
<tr valign="top"><td><code>z.start</code></td>
<td>
initial value for the centering constant. </td></tr>
Modified: pkg/ROptEst/chm/optIC.html
===================================================================
--- pkg/ROptEst/chm/optIC.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/optIC.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -30,18 +30,20 @@
optIC(model, risk, ...)
## S4 method for signature 'InfRobModel, asRisk':
-optIC(model, risk,
- z.start = NULL, A.start = NULL, upper = 1e4,
- maxiter = 50, tol = .Machine$double.eps^0.4, warn = TRUE, noLow = FALSE)
+optIC(model, risk, z.start = NULL, A.start = NULL,
+ upper = 1e4, maxiter = 50,
+ tol = .Machine$double.eps^0.4, warn = TRUE,
+ noLow = FALSE, verbose = FALSE)
## S4 method for signature 'InfRobModel, asUnOvShoot':
-optIC(model, risk,
- upper = 1e4, maxiter = 50,
- tol = .Machine$double.eps^0.4, warn = TRUE)
+optIC(model, risk, upper = 1e4, maxiter = 50,
+ tol = .Machine$double.eps^0.4, warn = TRUE)
## S4 method for signature 'FixRobModel, fiUnOvShoot':
-optIC(model, risk, sampleSize, upper = 1e4, maxiter = 50,
- tol = .Machine$double.eps^0.4, warn = TRUE, Algo = "A", cont = "left")
+optIC(model, risk, sampleSize, upper = 1e4,
+ maxiter = 50, tol = .Machine$double.eps^0.4,
+ warn = TRUE, Algo = "A", cont = "left",
+ verbose = FALSE)
</pre>
@@ -87,6 +89,9 @@
<tr valign="top"><td><code>noLow</code></td>
<td>
logical: is lower case to be computed? </td></tr>
+<tr valign="top"><td><code>verbose</code></td>
+<td>
+logical: if <code>TRUE</code>, some messages are printed </td></tr>
</table>
<h3>Details</h3>
Modified: pkg/ROptEst/chm/radiusMinimaxIC.html
===================================================================
--- pkg/ROptEst/chm/radiusMinimaxIC.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/chm/radiusMinimaxIC.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -31,7 +31,8 @@
## UncondNeighborhood, asGRisk':
radiusMinimaxIC(
L2Fam, neighbor, risk, loRad, upRad, z.start = NULL, A.start = NULL,
- upper = 1e5, maxiter = 100, tol = .Machine$double.eps^0.4, warn = FALSE)
+ upper = 1e5, maxiter = 50, tol = .Machine$double.eps^0.4,
+ warn = FALSE, verbose = FALSE)
</pre>
@@ -74,8 +75,21 @@
<tr valign="top"><td><code>warn</code></td>
<td>
logical: print warnings. </td></tr>
+<tr valign="top"><td><code>verbose</code></td>
+<td>
+logical: if <code>TRUE</code>, some messages are printed </td></tr>
</table>
+<h3>Details</h3>
+
+<p>
+In case the neighborhood radius is unknown, Rieder et al. (2001, 2008)
+and Kohl (2005) show that there is nevertheless a way to compute an
+optimally robust IC - the so-called radius-minimax IC - which is
+optimal for some radius interval.
+</p>
+
+
<h3>Value</h3>
<p>
@@ -99,11 +113,11 @@
<p>
Rieder, H., Kohl, M. and Ruckdeschel, P. (2008) The Costs of not Knowing
-the Radius. Statistical Methods and Applications <EM>17</EM>(1) 13-40.
+the Radius. Statistical Methods and Applications, <EM>17</EM>(1) 13-40.
</p>
<p>
Rieder, H., Kohl, M. and Ruckdeschel, P. (2001) The Costs of not Knowing
-the Radius. Submitted. Appeared as discussion paper Nr. 81.
+the Radius. Appeared as discussion paper Nr. 81.
SFB 373 (Quantification and Simulation of Economic Processes),
Humboldt University, Berlin; also available under
<a href="www.uni-bayreuth.de/departments/math/org/mathe7/RIEDER/pubs/RR.pdf">www.uni-bayreuth.de/departments/math/org/mathe7/RIEDER/pubs/RR.pdf</a>
@@ -125,8 +139,9 @@
<pre>
N <- NormLocationFamily(mean=0, sd=1)
-radiusMinimaxIC(L2Fam=N, neighbor=ContNeighborhood(),
- risk=asMSE(), loRad=0.1, upRad=0.5)
+radIC <- radiusMinimaxIC(L2Fam=N, neighbor=ContNeighborhood(),
+ risk=asMSE(), loRad=0.1, upRad=0.5)
+checkIC(radIC)
</pre>
Modified: pkg/ROptEst/inst/scripts/GammaModel.R
===================================================================
--- pkg/ROptEst/inst/scripts/GammaModel.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/inst/scripts/GammaModel.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -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))
+(est1 <- oneStepEstimator(x, IC=IC9, start=estimate(est0)))
Modified: pkg/ROptEst/inst/scripts/GumbelLocationModel.R
===================================================================
--- pkg/ROptEst/inst/scripts/GumbelLocationModel.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/inst/scripts/GumbelLocationModel.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -131,4 +131,4 @@
neighbor=ContNeighborhood(), risk=asMSE(), loRad=0.5, upRad=1)
(G0.est21 <- locMEstimator(G0.x, IC=G0.IC101))
-distrExOptions(ElowerTruncQuantile, 0) # default
+distrExOptions(ElowerTruncQuantile=0) # default
Modified: pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R
===================================================================
--- pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/inst/scripts/NormalLocationScaleModel.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -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))
Modified: pkg/ROptEst/inst/scripts/UnderOverShootRisk.R
===================================================================
--- pkg/ROptEst/inst/scripts/UnderOverShootRisk.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/inst/scripts/UnderOverShootRisk.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -110,14 +110,14 @@
## 3. Kolmogorov(-Smirnov) minimum distance estimator
-(est0 <- MDEstimator(x=X, NormLocationFamily(), interval = c(-5, 5)))
+(est0 <- MDEstimator(x=X, NormLocationFamily()))
## 4. one-step estimation
-N0.Rob7 <- InfRobModel(center = NormLocationFamily(mean = est0$estimate),
+N0.Rob7 <- InfRobModel(center = NormLocationFamily(mean = estimate(est0)),
neighbor = ContNeighborhood(radius=0.5))
N0.IC9 <- optIC(model=N0.Rob7, risk=asUnOvShoot(width = 1.960))
-(est1 <- oneStepEstimator(X, IC = N0.IC9, start = est0$estimate))
-N0.Rob8 <- FixRobModel(center = NormLocationFamily(mean = est0$estimate),
+(est1 <- oneStepEstimator(X, IC = N0.IC9, start = estimate(est0)))
+N0.Rob8 <- FixRobModel(center = NormLocationFamily(mean = estimate(est0)),
neighbor = ContNeighborhood(radius=0.05))
N0.IC10 <- optIC(model=N0.Rob8, risk=fiUnOvShoot(width = 1.960/sqrt(n)), sampleSize = 1e2)
-(est2 <- oneStepEstimator(X, IC = N0.IC10, start = est0$estimate))
+(est2 <- oneStepEstimator(X, IC = N0.IC10, start = estimate(est0)))
Modified: pkg/ROptEst/man/getBiasIC.Rd
===================================================================
--- pkg/ROptEst/man/getBiasIC.Rd 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/man/getBiasIC.Rd 2008-08-28 12:46:23 UTC (rev 170)
@@ -3,6 +3,7 @@
\alias{getBiasIC}
\alias{getBiasIC-methods}
\alias{getBiasIC,HampIC,UncondNeighborhood-method}
+\alias{getBiasIC,TotalVarIC,UncondNeighborhood-method}
\title{Generic function for the computation of the asymptotic bias for an IC}
\description{
@@ -26,6 +27,10 @@
\describe{
\item{IC = "HampIC", neighbor = "UncondNeighborhood"}{
reads off the as. bias from the risks-slot of the IC. }
+ \item{IC = "TotalVarIC", neighbor = "UncondNeighborhood"}{
+ reads off the as. bias from the risks-slot of the IC,
+ resp. if this is \code{NULL} from the corresponding
+ Lagrange Multipliers. }
}}
\references{
Huber, P.J. (1968) Robust Confidence Limits. Z. Wahrscheinlichkeitstheor.
Modified: pkg/ROptEst/man/getModifyIC.Rd
===================================================================
--- pkg/ROptEst/man/getModifyIC.Rd 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/man/getModifyIC.Rd 2008-08-28 12:46:23 UTC (rev 170)
@@ -3,11 +3,12 @@
\alias{getModifyIC-methods}
\alias{getModifyIC,L2ParamFamily,Neighborhood,asRisk-method}
\alias{getModifyIC,L2LocationFamily,UncondNeighborhood,asGRisk-method}
+\alias{getModifyIC,L2LocationFamily,UncondNeighborhood,fiUnOvShoot-method}
\alias{getModifyIC,L2ScaleFamily,ContNeighborhood,asGRisk-method}
\alias{getModifyIC,L2ScaleFamily,TotalVarNeighborhood,asGRisk-method}
\alias{getModifyIC,L2LocationScaleFamily,ContNeighborhood,asGRisk-method}
+\title{Generic Function for the Computation of Functions for Slot modifyIC}
-\title{Generic Function for the Computation of Functions for Slot modifyIC}
\description{
This function is used by internal computations and is rarely called directly.
}
@@ -15,6 +16,7 @@
getModifyIC(L2FamIC, neighbor, risk)
\S4method{getModifyIC}{L2ParamFamily,Neighborhood,asRisk}(L2FamIC, neighbor, risk)
\S4method{getModifyIC}{L2LocationFamily,UncondNeighborhood,asGRisk}(L2FamIC, neighbor, risk)
+\S4method{getModifyIC}{L2LocationFamily,UncondNeighborhood,fiUnOvShoot}(L2FamIC, neighbor, risk)
\S4method{getModifyIC}{L2ScaleFamily,ContNeighborhood,asGRisk}(L2FamIC, neighbor, risk)
\S4method{getModifyIC}{L2ScaleFamily,TotalVarNeighborhood,asGRisk}(L2FamIC, neighbor, risk)
\S4method{getModifyIC}{L2LocationScaleFamily,ContNeighborhood,asGRisk}(L2FamIC, neighbor, risk)
Modified: pkg/ROptEst/man/getRiskIC.Rd
===================================================================
--- pkg/ROptEst/man/getRiskIC.Rd 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/ROptEst/man/getRiskIC.Rd 2008-08-28 12:46:23 UTC (rev 170)
@@ -4,6 +4,7 @@
\alias{getRiskIC-methods}
\alias{getRiskIC,HampIC,asCov,missing,missing-method}
\alias{getRiskIC,HampIC,asCov,missing,L2ParamFamily-method}
+\alias{getRiskIC,TotalVarIC,asCov,missing,L2ParamFamily-method}
\title{Generic function for the computation of a risk for an IC}
\description{
@@ -15,6 +16,7 @@
\S4method{getRiskIC}{HampIC,asCov,missing,missing}(IC, risk)
\S4method{getRiskIC}{HampIC,asCov,missing,L2ParamFamily}(IC, risk, L2Fam)
+\S4method{getRiskIC}{TotalVarIC,asCov,missing,L2ParamFamily}(IC, risk, L2Fam)
}
\arguments{
@@ -35,6 +37,9 @@
\item{IC = "HampIC", risk = "asCov", neighbor = "missing", L2Fam = "L2ParamFamily"}{
asymptotic covariance of \code{IC} under \code{L2Fam} read off from corresp. \code{Risks} slot. }
+ \item{IC = "TotalVarIC", risk = "asCov", neighbor = "missing", L2Fam = "L2ParamFamily"}{
+ asymptotic covariance of \code{IC} read off from corresp. \code{Risks} slot,
+ resp. if this is \code{NULL} calculates it via \code{\link{getInfV}}.}
}}
\references{
Huber, P.J. (1968) Robust Confidence Limits. Z. Wahrscheinlichkeitstheor.
Modified: pkg/RandVar/R/util.R
===================================================================
--- pkg/RandVar/R/util.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/RandVar/R/util.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -3,5 +3,5 @@
{ if (is(try(return(f(distr)), silent = TRUE),
"try-error")){
rl <- function(n) { xr <- r(distr)(n); f(xr) }
- return(new("AbscontDistribution", r = rl, .withArith = TRUE, .withSim = TRUE))}
+ return(AbscontDistribution( r = rl, .withArith = TRUE, .withSim = TRUE)))}
}
\ No newline at end of file
Modified: pkg/RobAStBase/R/locMEstimator.R
===================================================================
--- pkg/RobAStBase/R/locMEstimator.R 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/RobAStBase/R/locMEstimator.R 2008-08-28 12:46:23 UTC (rev 170)
@@ -21,7 +21,14 @@
ncol = 2)
colnames(Infos) <- c("method", "message")
asVar <- getRiskIC(IC, risk = asCov(), L2Fam = L2Fam)$asCov$value
- asBias <- getRiskIC(IC, risk = asBias(), L2Fam = L2Fam)$asBias$value
+ asBias <- getRiskIC(IC, risk = asBias(),
+ neighbor = ContNeighborhood(1),
+ L2Fam = L2Fam)$asBias$value
+
+ names(res$root) <- nms <- locscalename(L2Fam)
+ asVar <- PosDefSymmMatrix(asVar)
+ dimnames(asVar) <- list(nms, nms)
+ names(asBias) <- nms
}else{
Infos <- matrix(c("locMEstimator", "Location M estimate"), ncol = 2)
colnames(Infos) <- c("method", "message")
Modified: pkg/RobAStBase/chm/00Index.html
===================================================================
--- pkg/RobAStBase/chm/00Index.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/RobAStBase/chm/00Index.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -416,7 +416,7 @@
<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="IC-class.html">plot,IC-method</a></td>
+<tr><td width="25%"><a href="IC-class.html">plot,IC,ANY-method</a></td>
<td>Influence curve</td></tr>
</table>
@@ -448,6 +448,8 @@
<h2><a name="S">-- S --</a></h2>
<table width="100%">
+<tr><td width="25%"><a href="ALEstimate-class.html">show,ALEstimate-method</a></td>
+<td>ALEstimate-class.</td></tr>
<tr><td width="25%"><a href="ContIC-class.html">show,ContIC-method</a></td>
<td>Influence curve of contamination type</td></tr>
<tr><td width="25%"><a href="FixRobModel-class.html">show,FixRobModel-method</a></td>
@@ -458,6 +460,10 @@
<td>Influence curve</td></tr>
<tr><td width="25%"><a href="InfRobModel-class.html">show,InfRobModel-method</a></td>
<td>Robust model with infinitesimal (unconditional) neighborhood</td></tr>
+<tr><td width="25%"><a href="kStepEstimate-class.html">show,kStepEstimate-method</a></td>
+<td>kStepEstimate-class.</td></tr>
+<tr><td width="25%"><a href="MEstimate-class.html">show,MEstimate-method</a></td>
+<td>MEstimate-class.</td></tr>
<tr><td width="25%"><a href="Neighborhood-class.html">show,Neighborhood-method</a></td>
<td>Neighborhood</td></tr>
<tr><td width="25%"><a href="TotalVarIC-class.html">show,TotalVarIC-method</a></td>
Modified: pkg/RobAStBase/chm/IC-class.html
===================================================================
--- pkg/RobAStBase/chm/IC-class.html 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/RobAStBase/chm/IC-class.html 2008-08-28 12:46:23 UTC (rev 170)
@@ -17,7 +17,7 @@
<param name="keyword" value="R: evalIC,IC,numeric-method">
<param name="keyword" value="R: evalIC,IC,matrix-method">
<param name="keyword" value="R: infoPlot,IC-method">
-<param name="keyword" value="R: plot,IC-method">
+<param name="keyword" value="R: plot,IC,ANY-method">
<param name="keyword" value="R: show,IC-method">
<param name="keyword" value=" Influence curve">
</object>
Modified: pkg/RobAStBase/chm/RobAStBase.chm
===================================================================
(Binary files differ)
Modified: pkg/RobAStBase/chm/RobAStBase.toc
===================================================================
--- pkg/RobAStBase/chm/RobAStBase.toc 2008-08-28 08:23:52 UTC (rev 169)
+++ pkg/RobAStBase/chm/RobAStBase.toc 2008-08-28 12:46:23 UTC (rev 170)
@@ -682,7 +682,7 @@
<param name="Local" value="ALEstimate-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
-<param name="Name" value="plot,IC-method">
+<param name="Name" value="plot,IC,ANY-method">
<param name="Local" value="IC-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
@@ -730,6 +730,10 @@
<param name="Local" value="RobWeight-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="show,ALEstimate-method">
+<param name="Local" value="ALEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="show,ContIC-method">
<param name="Local" value="ContIC-class.html">
</OBJECT>
@@ -750,6 +754,14 @@
<param name="Local" value="InfRobModel-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="show,kStepEstimate-method">
+<param name="Local" value="kStepEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="show,MEstimate-method">
+<param name="Local" value="MEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="show,Neighborhood-method">
<param name="Local" value="Neighborhood-class.html">
</OBJECT>
Modified: pkg/RobLox/chm/RobLox.chm
===================================================================
(Binary files differ)
More information about the Robast-commits
mailing list