ca = read.csv("calcium.csv")
d = ca$oxalate - ca$flame
mean(d)
sd(d)
summary(d)
hist(d,10)

plot(calcium$oxalate,d)
t.test(d)
wilcox.test(d)

#ci for median
q025 = qbinom(.025,length(d),.5)
q975 = qbinom(.975,length(d),.5)
z = sort(d)
median(d)
cimed = c(z[q025], z[q975])