################################################### ### chunk number 1: ################################################### options("width"=70) set.seed(13131) load("students2010.rdata") attach(students) # load("students.RData") # attach(students) #library(rgdal) #nrw = readOGR("katha", "NRW", verbose=FALSE) # students$I.am.2 = as.character(students$I.am.) # students$I.am.2 = factor(students$I.am.2, levels = c("small", "intermediate", "tall")) ################################################### ### chunk number 2: ################################################### n1 = runif(1000) n2 = apply(replicate(1000, runif(2)), 2, mean) n3 = apply(replicate(1000, runif(3)), 2, mean) n5 = apply(replicate(1000, runif(5)), 2, mean) n10 = apply(replicate(1000, runif(10)), 2, mean) library(lattice) print(histogram(~values|ind, stack(data.frame('n=01'=n1,'n=02'=n2, 'n=03'=n3,'n=05'=n5,'n=10'=n10)), layout=c(1,5),breaks=seq(0,1,by=.05),as.table=TRUE)) ################################################### ### chunk number 3: ################################################### par(mfrow=c(2,1)) hist(rgamma(1000, .05), breaks=20, main = "gamma distribution, shape = 0.05") hist(apply(replicate(1000, rgamma(50, .05)), 2, mean), main = "means of random samples with size 50: still far from normal") ################################################### ### chunk number 4: ################################################### t.test(Length~Gender, var.equal = TRUE) ################################################### ### chunk number 5: ################################################### x1 = c(13.5, 15.3, 7.5, 10.3, 8.7) x2 = c(12.7, 15.1, 6.6, 8.5, 8.0) x1 - x2 ################################################### ### chunk number 6: ################################################### t.test(x1, x2, var.equal = TRUE) t.test(x1-x2) ################################################### ### chunk number 7: ################################################### qf(.95, 20, 30) v1 = var(Length[Gender == "male"]) v2 = var(Length[Gender == "female"]) v1 v2 v2/v1 qf(.95, length(Length[Gender == "female"]), length(Length[Gender == "male"])) ################################################### ### chunk number 8: ################################################### t.test(Length~Gender, var.equal = TRUE) t.test(Length~Gender)