8.2 Bootstrap in hypothesis testing
skew <- function(x){
m3 <- sum((x-mean(x))ˆ3)/length(x)
s3 <- sqrt(var(x))ˆ3
m3/s3
}
data <- read.table("c:\\temp\\skewdata.txt",header=T)
attach(data)
names(data)
skew(values)
skew(values)/sqrt(6/length(values))
1-t(2.94928)
skew(sqrt(values))/sqrt(6/length(values))
skew(log(values))/sqrt(6/length(values))