Generic function for the computation of optimally robust ICs.
optIC(model, risk, ...)# S4 method for L2ParamFamily,asCov
optIC(model, risk)
# S4 method for InfRobModel,asRisk
optIC(model, risk, z.start = NULL, A.start = NULL, upper = 1e4,
maxiter = 50, tol = .Machine$double.eps^0.4, warn = TRUE)
# S4 method for InfRobModel,asUnOvShoot
optIC(model, risk, upper = 1e4, maxiter = 50,
tol = .Machine$double.eps^0.4, warn = TRUE)
# S4 method for FixRobModel,fiUnOvShoot
optIC(model, risk, sampleSize, upper = 1e4, maxiter = 50,
tol = .Machine$double.eps^0.4, warn = TRUE, Algo = "A", cont = "left")
probability model.
object of class "RiskType"
.
additional parameters.
initial value for the centering constant.
initial value for the standardizing matrix.
upper bound for the optimal clipping bound.
the maximum number of iterations.
the desired accuracy (convergence tolerance).
logical: print warnings.
integer: sample size.
"A" or "B".
"left" or "right".
Some optimally robust IC is computed.
computes classical optimal influence curve for L2 differentiable parametric families.
computes optimally robust influence curve for robust models with infinitesimal neighborhoods and various asymptotic risks.
computes optimally robust influence curve for robust models with infinitesimal neighborhoods and asymptotic under-/overshoot risk.
computes optimally robust influence curve for robust models with fixed neighborhoods and finite-sample under-/overshoot risk.
In case of the finite-sample risk "fiUnOvShoot"
one can choose
between two algorithms for the computation of this risk where the least favorable
contamination is assumed to be left or right of some bound. For more details
we refer to Section 11.3 of Kohl (2005).
Huber, P.J. (1968) Robust Confidence Limits. Z. Wahrscheinlichkeitstheor. Verw. Geb. 10:269--278.
Rieder, H. (1980) Estimates derived from robust tests. Ann. Stats. 8: 106--115.
Rieder, H. (1994) Robust Asymptotic Statistics. New York: Springer.
Kohl, M. (2005) Numerical Contributions to the Asymptotic Theory of Robustness. Bayreuth: Dissertation.
# NOT RUN {
B <- BinomFamily(size = 25, prob = 0.25)
## classical optimal IC
IC0 <- optIC(model = B, risk = asCov())
plot(IC0) # plot IC
checkIC(IC0, B)
# }
Run the code above in your browser using DataLab