control = scan("ozonecontrol.csv")
treat = scan("ozonetreat.csv")
summary(treat)
summary(control)

t.test(control, treat)
t.test(control, treat,var.equal=T)

boxplot(control, treat, names=c("control", "treatment"))
stripchart(list(control,treat), group.names=c("control","treatment"))
qqplot(control,treat)

m=length(control)
n = length(treat)
z=rep(0,2000)
for (i in 1:2000){
	z[i]=sum(sample(m+n,m))
	}
hist(z)
qqnorm(z)

wilcox.test(control,treat, conf.int=T)
pihat=385/(length(control)*length(treat))


x = read.csv("calcium.csv")
oxalate=x[,1]
flame = x[,2]
mean(oxalate)
mean(flame)
sd(oxalate)
sd(flame)
t.test(oxalate,flame)

plot(oxalate, flame)
abline(0,1)
boxplot(oxalate-flame)
hist(oxalate - flame)
plot(oxalate, oxalate - flame)
t.test(oxalate,flame, paired=T)
wilcox.test(oxalate, flame, paired=T)

#bootstrap
diff=oxalate-flame
med(diff)
mean(diff)
sd(diff)
sd(diff)/sqrt(n)
n=length(diff)

#bootstrap mean
mn = rep(0,1000)
for (i in 1:1000){
	mn[i]=mean(sample(diff,n,replace=T))
	}
sd(mn)



#bootstrap median
bmed = rep(0,1000)
for (i in 1:1000){
	bmed[i] = median(sample(diff,n, replace=T))
	}
sd(bmed)