R语言rmcorr包CI Bootstrap参数不生效问题咨询

问题根源 #

你的参数语法没有错误，但rmcorr的默认输出只会展示analytic方法的置信区间，bootstrap的计算结果并不会直接显示在主输出里——它被存储在返回对象的boot子元素中。你只看了主输出的CI，所以误以为两种方法结果一致。

正确获取Bootstrap CI的步骤 #

1. 先将rmcorr的运行结果赋值给一个变量：

rmc_out <- rmcorr(Subject, PaCO2, pH, bland1995, CI.level=0.95, CIs=c("analytic","bootstrap"), nreps=100, bstrap.out=T)

2. 提取并查看bootstrap置信区间：

# 查看Bootstrap计算的相关结果
rmc_out$boot

这个子对象里会包含bootstrap得到的相关系数分布、置信区间等数据。

优化建议 #

• nreps=100的重复次数过低，实际分析中建议设置为1000次及以上（如nreps=1000），提升bootstrap结果的稳定性。
• 检查rmcorr包版本，旧版本可能存在参数逻辑问题，执行update.packages("rmcorr")更新到最新版。

内容的提问来源于stack exchange，提问作者cello