confint | R Documentation |
The function confint
computes confindence interval for rmx estimates.
## S3 method for class 'rmx'
confint(object, parm, level = 0.95, method = "as", R = 9999, type = "all", ...)
object |
object of S3 class |
parm |
see |
level |
see |
method |
see details. |
R |
number of bootstrap replicates. |
type |
A vector of character strings representing the type of intervals required.
The value should be any subset of the values c("norm","basic", "stud", "perc", "bca")
or simply "all" which will compute all five types of intervals; see also
|
... |
further arguments passed through especially to |
The function is inspired by the respective function of the RobASt-family of packages.
In case of optimally-robust RMX estimators computed with function rmx
(S3 class rmx
), confint intervals are computed either only using
the asymptotic (co)variance (method = "as"
, default) or using the
asymptotic (co)variance and the maximum asymptotic bias (method = "as.bias"
)
or using bootstrap (method = "boot"
).
An object of class "rmxCI"
is returned. It contails at least the
following arguments:
method |
method selected for computing the confidence interval. |
conf.int |
numeric matrix with the confidence bounds. |
rmxEst |
point estimates. |
Matthias Kohl Matthias.Kohl@stamats.de
Kohl, M. (2005). Numerical Contributions to the Asymptotic Theory of Robustness. Bayreuth: Dissertation.
M. Kohl, P. Ruckdeschel, and H. Rieder (2010). Infinitesimally Robust Estimation in General Smoothly Parametrized Models. Statistical Methods and Application, 19(3):333-354.
Rieder, H. (1994) Robust Asymptotic Statistics. New York: Springer.
Rieder, H., Kohl, M. and Ruckdeschel, P. (2008) The Costs of not Knowing the Radius. Statistical Methods and Applications 17(1) 13-40. Extended version: http://r-kurs.de/RRlong.pdf.
rmx
, confint
ind <- rbinom(100, size=1, prob=0.05)
x <- rnorm(100, mean=ind*3, sd=(1-ind) + ind*9)
res <- rmx(x, eps.lower = 0.01, eps.upper = 0.1)
confint(res)
confint(res, method = "as.bias")
#confint(res, method = "boot", R = 999)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.