我想使用boot()和boot.ci()函数从library("boot")非常大的数据集(~50 000)type="bca".
如果R(bootstraps的数量)太小,则会出现以下错误:
Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o, :
estimated adjustment 'a' is NA
我也不希望它太大.
R的好数字是多少？我知道这将取决于数据的大小.
Rei*_*son 10
一般而言,您需要更多数量的自举,以获得置信区间而不是估计值.自从计算能力在过去十年中有所增加以来,我已经听到了1000个用于估计的自举和10000个用于置信区间的建议.
boot()现在这个函数有parallel一个可以在多个内核上运行bootstrap的参数,如果计算时间对你来说很重要,你可能会使用它来大大减少在你的大小的数据集上运行分析所需的时间.
| 归档时间： |
|
| 查看次数： |
3205 次 |
| 最近记录： |