Displaying 1 result from an estimated 1 matches for "sampleskew".
2009 Jun 12
1
Studentized intervals
...d bootstrap intervals for the skewness of a data set. I have tried the following (nerve.dat is a set of data containing observations on one variable) (using Windows XP):
x <- scan("e:/Flashbackup2009/Nonparametrics/nerve.dat")
n <- length(x)
library(e1071)
skewness(x)
library(boot)
sampleskew <- function(x,d) {return(skewness(x[d]))}
bb <- boot(x,sampleskew,R=1000)
boot.ci(bb)
This gives four intervals (Normal, Basic, Percentile, BCa)
But gives the following warning message:
Warning message:
In boot.ci(bb, type = "all") :
bootstrap variances needed for studentized i...