###################################################################### # Volleyball radius example # # # ###################################################################### ###################################################################### #Read in data volleyball <- read.csv(file = "volleyball.csv") head(volleyball) #Shows first 6 observations tail(volleyball) #Shows last 6 observations ###################################################################### #Longer way ybar <- mean(volleyball$radius) s <- sd(volleyball$radius) alpha <- 0.05 n <- length(volleyball$radius) #Confidence interval lower <- ybar - qt(p = 1 - alpha/2, df = n-1) * s / sqrt(n) upper <- ybar + qt(p = 1 - alpha/2, df = n-1) * s / sqrt(n) data.frame(lower, upper) #test statistic mu0 <- 11.6 t <- (ybar - mu0)/(s/sqrt(n)) crit.val <- qt(p = 1 - alpha/2, df = n-1) data.frame(t, crit.val) #P-value 2*(1 - pt(q = abs(t), df = n - 1)) ###################################################################### #Shorter way t.test(x = volleyball$radius, alternative = "two.sided", mu = 11.6, conf.level = 0.95) save.results <- t.test(x = volleyball$radius, alternative = "two.sided", mu = 11.6, conf.level = 0.95) save.results names(save.results) save.results$conf.int save.results$conf.int[1] save.results$conf.int[2] save.results$statistic save.results$null.value ###################################################################### #Ho: mu >= 11.6 vs. Ha: mu < 11.6 t.test(x = volleyball$radius, alternative = "less", mu = 11.6, conf.level = 0.95) # #