load("bumpus.R") bumpus bumpus.R b attach(b) hist(X1) hist(X1, prob=T) hist(X1, prob=T) lines(density(X1)) lines(density(X1)) lines(density(X1[survive == "Y"]), col="green") lines(density(X1[survive == "N"]), col="red") t.test(X1~survive) e = read.csv("exam1R.csv") e save(e, file="exam1.R") load("exam1.R") e 434/176 pf(2.47, 27, 20) 1 - pf(2.47, 27, 20) 2*(1 - pf(2.47, 27, 20)) b X4 hist(X4) lines(density(X4[survive == "Y"]), col="green") hist(X4, prob=T) lines(density(X4[survive == "Y"]), col="green") hist(X4, prob=T, ylim=c(0,1)) lines(density(X4[survive == "Y"]), col="green") lines(density(X4[survive == "N"]), col="red") var.test(X4~survive) var.test(X1~survive) var.test(X2~survive) var.test(X3~survive) var.test(X3~survive) var.test(X4~survive) var.test(X5~survive) t.test(X1~survive) t.test(X2~survive) t.test(X3~survive) t.test(X4~survive) t.test(X5~survive) var.test(X4) var.test(X4~survive)