[R] bootstrap for confidence intervals of the mean

[Robert A LaBudde]

See the help for boot().

The function in the 2nd argument has to be of a special form.

You need to define such a form, as in:

fmean<- function (x, i) mean(x[i]) #use data x[] and indices i

and then

boot.out<- boot(d, fmean, R=1000, sim='permutation')

At 12:59 PM 4/22/2008, stephen sefick wrote:
>d = c(0L, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0L, 0L, 7375L,
>NA, NA, 17092L, 0L, 0L, 32390L, 2326L, 22672L, 13550L, 18285L)
>
>boot.out <-boot(d, mean, R=1000, sim="permutation")
>
>Error in mean.default(data, original, ...) :
>   'trim' must be numeric of length one
>
>I know that I am missing something but I can't figure it out.
>thanks
>
>stephen
>
>--
>Let's not spend our time and resources thinking about things that are so
>little or so large that all they really do for us is puff us up and make us
>feel like gods. We are mammals, and have not exhausted the annoying little
>problems of being mammals.
>
>-K. Mullis
>
>         [[alternative HTML version deleted]]
>
>______________________________________________
>R-help at r-project.org mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.

================================================================
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS  e-mail: ral at lcfltd.com
Least Cost Formulations, Ltd.            URL: http://lcfltd.com/
824 Timberlake Drive                     Tel: 757-467-0954
Virginia Beach, VA 23464-3239            Fax: 757-467-2947

"Vere scire est per causas scire"