使用R语言bootstrap方法时出现维度数量错误的问题
复现《Eur. Phys. J. C (2022) 82:1106》中MFV统计量Bootstrap置信区间的问题及解决
- 本人是R编程新手,正在复现《Eur. Phys. J. C (2022) 82:1106》(DOI:10.1140/epjc/s10052-022-11071-9)的研究数值
- 该论文应用了最大频率值(MFV)统计量,作者使用Bootstrap方法估计MFV的置信区间(详见论文第5页)
- 最初在应用Bootstrap方法时遇到问题,已准备好最小工作示例复现该问题
- 后续已找到解决方案,核心实现步骤如下:
- 编写自定义函数,用于计算数据集的MFV统计量
- 安装并加载
boot包,使用boot()函数生成Bootstrap样本,同时批量计算每个样本的MFV值 - 调用
boot.ci()函数,基于Bootstrap结果计算置信区间,可选择百分位法、BCa法等合适的类型
内容的提问来源于stack exchange,提问作者Victor Golovko
相关产品推荐
相关产品推荐

