rlsOptIC.Ha3: Computation of the optimally robust IC for Ha3 estimators
Description
The function rlsOptIC.Ha3
computes the optimally robust IC for
Ha3 estimators in case of normal location with unknown scale and
(convex) contamination neighborhoods. The definition of
these estimators can be found in Subsection 8.5.2 of Kohl (2005).
Usage
rlsOptIC.Ha3(r, a.start = 0.25, b.start = 2.5, c.start = 5,
delta = 1e-06, MAX = 100)
Arguments
r
non-negative real: neighborhood radius.
a.start
positive real: starting value for a.
b.start
positive real: starting value for b.
c.start
positive real: starting value for c.
delta
the desired accuracy (convergence tolerance).
MAX
if a or b or c are beyond the admitted values,
MAX
is returned.
Value
Object of class "IC"
Details
The computation of the optimally robust IC for Ha3 estimators
is based on optim
where MAX
is used to
control the constraints on a, b and c. The optimal values of
the tuning constants a, b and c can be read off
from the slot Infos
of the resulting IC.
References
Kohl, M. (2005) Numerical Contributions to the Asymptotic Theory of Robustness.
Bayreuth: Dissertation.
Examples
Run this code# NOT RUN {
IC1 <- rlsOptIC.Ha3(r = 0.1)
checkIC(IC1)
Risks(IC1)
Infos(IC1)
## don't run to reduce check time on CRAN
# }
# NOT RUN {
plot(IC1)
infoPlot(IC1)
# }
Run the code above in your browser using DataLab