chickwts = read.csv("chickwts.csv") names(chickwts) boxplot(weight~feed) attach(chickwts) boxplot(weight~feed) a = lm(weight~feed) shapiro.test(a$resid) bartlett.test(weight~feed) data(ozone) data(Ozone) ??ozone ?kruskal.test data(airquality) allomet = read.csv("allomet.csv") iq = read.csv("IQdata.csv") names(allomet) attach(allomet) plot(length, mass) a = lm(mass~length) plot(a) shapiro.test(a$resid) names(iq) attach(iq) a = lm(offspringIQ ~ parentsIQ) plot(a$resid) a plot(a) attach(chickwts) boxplot(weight~feed) a = lm(weight~feed) shapiro.test(a$resid) plot(a) bartlett.test(weight~feed) attach(airquality) names(airquality) boxplot(Ozone~Month) a = lm(Ozone~Month) plot(a) shapiro.test(a$resid) bartlett.test(Ozone~Month) a = lm(log(Ozone) ~ Month) shapiro.test(a$resid) a = lm(Ozone^.2 ~ Month) shapiro.test(a$resid) bartlett.test(Ozone^.2 ~ Month) kruskal.test(Ozone ~ Month) anova(lm(Ozone^.2 ~ Month)) pairwise.wilcoxon.test(Ozone, Month) pairwise.wilcox.test(Ozone, Month) x = c(1, 2, 2, 3, 4, 4, 5) x rank(x) rO = rank(Ozone) a lm(rO~ Month) a lm(rO ~ Month) a = lm(rO ~ Month) anova(a)