对于R平方 - 自举的R - 循环为95%置信区间

时间:2018-04-05 11:25:46

标签: r loops linear-regression bootstrapping confidence-interval

我在R没有多少经验,所以我可能要求一些明显的东西。我有一个数据框(INs3),每列都有一个变量。

head(INs3)

Dependant <- (1280.75, 1317.70, 1166.35, 1663.10, ...)  
X1 <- (0.029, 0.033, ...)  
X2 <- (0.076, 0.78, ...)  
Xn

我正在尝试使用95% confident interval估算R squared的{​​{1}}。我设法编写了用两个变量进行计算的代码。

bootstrap
boot.test <- boot(INs3,function(data,indices)                          
summary(lm(Dependant~X1,data[indices,]))$r.squared,R=100)

ci <- boot.ci(boot.test, type="bca")   
ci

我想估计将变量“Dependent”与每个变量(BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS Based on 100 bootstrap replicates CALL : boot.ci(boot.out = boot.test, type = "bca") Intervals : Level BCa 95% ( 0.0070, 0.6189 ) X1等)配对的置信区间。我的数据集有超过X2个变量,所以我想实现一个循环来进行所有这些计算并将它们存储在数据框或列表中。我设法生成了一个列表,其中存储了与100000相关的数据摘要,但我只对bootstrap的置信区间感兴趣。另外,我无法在下面的代码中使用合适的函数来估计那些置信区间。

R squared

有人能就这个问题给我一些建议吗? 。谢谢!

0 个答案:

没有答案