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

tab <- table(lister)
tab

chisq.test(tab, correct = FALSE)
fisher.test(tab)

tab <- cbind(c(4475, 1597), c(65, 31))
fisher.test(tab)

tab[,2:1]
fisher.test(tab[,2:1])

# Power and sample size calculations
power.t.test(delta = 125, sd = 200, n = 14, type = "paired")

power.t.test(delta = 125, sd = 200, power = .8, type = "paired")
#power.t.test(delta = 125, sd = 200, n = 23, type = "paired")

n <- 10:50
t <- power.t.test(delta = 125, sd = 200, n = 10:50, type = "paired")

#setwd("M:/General Research Directory/171-161/Labs2014/Lab10")
#pdf("PowerPlotR.pdf")
plot(n, t$power, xlab = "Total Sample Size", ylab = "Power",
     main = "One Sample t Test for Mean", type = 'o')
#dev.off()