根据回归模型估计参数(2.685)及其SE(0.319),我可以如下计算95%Wald置信区间:
> c("2.5%" = 2.685 - 1.96*0.319,
"97.5%" = 2.685 + 1.96*0.319)
2.5% 97.5%
2.05976 3.31024
CRAN上的现有软件包中是否有R功能可以执行此操作?
为简单起见,我们假设参数遵循正态分布。我正在寻找来自Rmisc
的等效CI()
,但它将估算值及其SE作为输入而不是矢量。
答案 0 :(得分:1)
假设正常,单线解决方案将是:
> sapply(c(0.025, 0.975), function(x) qnorm(x, 2.685, 0.319))
[1] 2.059771 3.310229
即使如评论中所指出的那样,它只是简单得多:
> qnorm(c(0.025, 0.975), 2.685, 0.319)
[1] 2.059771 3.310229