## Lab 13 Code
diarrhea <- read.delim("http://myweb.uiowa.edu/pbreheny/data/diarrhea.txt")
attach(diarrhea)

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


logStool <- log(Stool)
require(lattice)

histogram(~logStool | Group)

var(logStool[Group == "Control"])
var(logStool[Group == "Treatment"])

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

## Point estimate
exp(5.212-4.871)

## 95% Lower Bound
exp(0.103)

## 95% Upper Bound
exp(0.581)


## Part 3

wilcox.test(Stool ~ Group)
t.test(Stool ~ Group, var.equal = TRUE)
t.test(logStool ~ Group, var.equal = TRUE)


## Part 4
detach(diarrhea)

oat <- read.delim("http://myweb.uiowa.edu/pbreheny/data/oatbran.txt")
attach(oat)

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