# Lab 10 supplemental code

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

t.test(TRG, mu = 120)
cf <- read.delim("http://myweb.uiowa.edu/pbreheny/data/cystic-fibrosis.txt")
attach(cf)

t.test(Drug,Placebo,paired=TRUE)
t.test(Drug - Placebo) ## Same thing

sd(Drug - Placebo)
power.t.test(delta = -136.5, power = .9, sd = 223.17,type = "paired")

power.t.test(delta = -136.5, n = 50, sd = 223.17,type = "paired")

lister <- read.delim("http://myweb.uiowa.edu/pbreheny/data/lister.txt")
attach(lister)
(tab <- table(lister))

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