diarrhea <- read.table("http://myweb.uiowa.edu/pbreheny/161/data/diarrhea.txt", header = T)
attach(diarrhea)

boxplot(Stool~Group,col="gray")

t.test(Stool~Group, diarrhea[Stool < 750,], var.equal=TRUE)

logStool <- log(Stool)

require(lattice)
histogram(~logStool|Group)

t.test(logStool~Group, var.equal=FALSE)
t.test(logStool~Group)

wilcox.test(Stool~Group)


oatbran <- read.table("http://myweb.uiowa.edu/pbreheny/161/data/oatbran.txt", header = T)
attach(oatbran)


wilcox.test(CornFlakes, OatBran, paired = TRUE)