R confint 模型参数的置信区间

说明

计算拟合模型中一个或多个参数的置信区间。从类 "lm" 继承的对象有一个默认值和一个方法。

用法

confint(object, parm, level = 0.95, ...)
## Default S3 method:
confint(object, parm, level = 0.95, ...)
## S3 method for class 'lm'
confint(object, parm, level = 0.95, ...)
## S3 method for class 'glm'
confint(object, parm, level = 0.95, trace = FALSE, test=c("LRT", "Rao"), ...)
## S3 method for class 'nls'
confint(object, parm, level = 0.95, ...)

参数

object	拟合的模型对象。
parm	指定哪些参数被赋予置信区间，可以是数字向量或名称向量。如果缺失，则考虑所有参数。
level	所需的置信度。
trace	合乎逻辑的。是否应该追踪分析？
test	在分析中使用似然比或 Rao 分数检验。
...	方法的附加参数。

细节

confint 是一个通用函数。默认方法假定正态性，并且需要合适的 coef 和 vcov 方法可用。可以直接调用默认方法与其他方法进行比较。

对于"lm" 类的对象，使用基于t 值的直接公式。

类 "glm" 和 "nls" 的方法调用适当的配置文件方法，然后通过在配置文件轨迹中插值来查找置信区间。如果轮廓对象已经可用，则可以将其用作主要参数，而不是拟合模型对象本身。

值

矩阵(或向量)，其中的列给出每个参数的置信下限和上限。这些将以百分比标记为 (1-level)/2 和 1 - (1-level)/2(默认为 2.5% 和 97.5%)。

例子

fit <- lm(100/mpg ~ disp + hp + wt + am, data = mtcars)
confint(fit)
confint(fit, "wt")

## from example(glm)
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3, 1, 9); treatment <- gl(3, 3)
glm.D93 <- glm(counts ~ outcome + treatment, family = poisson())
confint(glm.D93) # needs MASS to be installed
confint.default(glm.D93)  # based on asymptotic normality

参考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.

也可以看看

原始版本：包 MASS 中的confint.glm和confint.nls。