我使用混合模型(带有交互作用的多个变量和一个随机变量)进行了引导。我得到了这个结果(只是部分):
> boot_out
ORDINARY NONPARAMETRIC BOOTSTRAP
Call:
boot(data = a001a1, statistic = bootReg, R = 1000)
Bootstrap Statistics :
original bias std. error
t1* 4.887383e+01 -1.677061e+00 4.362948e-01
t2* 3.066825e+01 1.264024e+00 5.328387e-01
t3* 8.105422e+01 2.368599e+00 6.789091e-01
t4* 1.620562e+02 4.908711e+00 1.779522e+00
......
现在,我想获取截距的置信区间:
> boot.ci(boot_out,type=c("norm","basic","perc"), index=1)
BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
Based on 1000 bootstrap replicates
CALL :
boot.ci(boot.out = boot_out, type = c("norm", "basic",
"perc"), index = 1)
Intervals :
Level Normal Basic Percentile
95% (49.70, 51.41 ) (49.70, 51.41 ) (46.34, 48.05 )
Calculations and Intervals on Original Scale
校正后的偏差估计为:
48.873 -1.677
1 47.196
我的问题是正常CI和基本CI在估计值之外(原始值和校正值)。我只是想知道如何应对。
更新1:
这是一个类似的问题,有很多答复。
2
只需评论一下:经典著作中的Efron&Tibshirani(1993)主张反对偏差校正,称这是“危险”和“问题”做法,可能会导致标准误差增加。
—
蒂姆
@Tim感谢您的评论。我会看这本书。也许一种解决方案是使用估计值和引导程序se *来计算置信区间。在我的情况下,偏差只会对估算产生轻微影响。
—
giordano