Sign Up to our social questions and Answers Engine to ask questions, answer people’s questions, and connect with other people.
Login to our social questions & Answers Engine to ask questions answer people’s questions & connect with other people.
Lost your password? Please enter your email address. You will receive a link and will create a new password via email.
Please briefly explain why you feel this question should be reported.
Please briefly explain why you feel this answer should be reported.
Please briefly explain why you feel this user should be reported.
I would like to use the boot() and boot.ci() functions from library("boot") for very large data sets(~50 000) with type="bca".
If R(number of bootstraps) is too small, it will give the following error:
Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o, :
estimated adjustment 'a' is NA
I wouldn’t want it to be too large either.
What is a good number for R? I know it would depend on the size of the data.
You need a greater number of bootstraps, in general, for confidence intervals than for estimates. I have heard suggestions of 1000 bootstraps for estimates and 10000 for confidence intervals since computing power has increased in the last decade.
The
boot()function nowadays has argumentparallelwhich can run the bootstraps on multiple cores, which if compute time is a concern for you, you might look to use to reduce significantly the time it takes to run the analysis on data sets of your size.