gpt4 book ai didi

r - 从具有多个统计数据的 boot() 函数中获取绘图和 95% CI

转载 作者:行者123 更新时间:2023-12-04 09:23:39 26 4
gpt4 key购买 nike

我开始探索 boot() 的极限boot 中的函数包裹。我刚刚研究了如何从单个调用中提取多个统计信息,但是我不知道如何将它们编入索引。

bs <- function(formula, data, indices) {
d <- data[indices,] # allows boot to select sample
fit <- lm(formula, data=d)
return(c(coef(fit), summary(fit)$r.squared)) # four stats extracted
}

现在引导 1000 次复制。
(results <- boot(data = mtcars, statistic = bs, R = 1000, formula = mpg ~ wt + disp))

现在,如果我们绘制并尝试获得该对象的 95% CI,我们只会获得第一个统计数据,在这种情况下是模型的截距。
plot(results)
boot.ci(results, type="bca")

# BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
# Based on 1000 bootstrap replicates
#
# CALL :
# boot.ci(boot.out = results, type = "bca")
#
# Intervals :
# Level BCa
# 95% (30.27, 39.53 )
# Calculations and Intervals on Original Scale

我如何索引其他四个参数,一个一个和/或一起?

最佳答案

首先注意:

> head(results$t)
[,1] [,2] [,3] [,4]
[1,] 35.94765 -4.358146 -0.010605142 0.7471584
[2,] 33.54691 -3.319141 -0.011821276 0.8463980
[3,] 34.01732 -3.627647 -0.009797742 0.6727989
[4,] 32.26678 -1.717800 -0.030443778 0.7466273
[5,] 35.78895 -3.944054 -0.016471864 0.8277447
[6,] 33.99407 -3.560855 -0.014129072 0.8412922

所以你可以使用 result$t[,i]访问第 i 个元素。看起来第一列是截距,第二列是权重系数,第三列是位移系数,最后一列是 R 平方。

您也可以使用 boot.ci(results, type="bca", index=i)要访问第 i 列,情节 plot(results, index=i) 也是如此.

关于r - 从具有多个统计数据的 boot() 函数中获取绘图和 95% CI,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/51371307/

26 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com