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.
______________________________________________
R-help@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.
