[R] Bootstrap Confidence Intervals

[_Fede_]

Hi all.

This is my first post in this forum. Finally I find a forum in the web about
R, although is not in my language. 

Now I'm working with Bootstrap CI. I'd like to know how I can calculate a
Bootstrap CI for any statistic, in particular, for Kurtosis Coeficient. I
have done the following code lines:

> library(boot)
> x=rnorm(20)
> kurtosis=function(x) (mean((x-mean(x))^4))/(sd(x)^4)
> z <- numeric(10000)
> for(i in 1:10000)
> z[i]=kurtosis(sample(x, replace=TRUE))
> boot.ci(z, conf = 0.95,type = c("norm","basic","perc","bca"))

But the output shows the next error:

Error en if (ncol(boot.out$t) < max(index)) { : 
        argumento tiene longitud cero

I don't know what is wrong.

I hope that somebody can help me. Sorry for my english.

All have a nice new year.

_Fede_

-- 
View this message in context: http://www.nabble.com/Bootstrap-Confidence-Intervals-tp14550471p14550471.html
Sent from the R help mailing list archive at Nabble.com.