我在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
有人能就这个问题给我一些建议吗? 。谢谢!