chroma = read.table("chromatograph.txt",header=T)
chroma$means=apply(chroma[,4:6],1,mean)
chroma$vari=apply(chroma[,4:6],1,var)
chroma$stdv=apply(chroma[,4:6],1,sd)
attach(chroma)
round(means,2)
round(vari,2)
round(stdv,2)
chroma.sp=sqrt(mean(vari))
plot(means,stdv)
mean(means[CO==1])-mean(means[CO==-1])
chroma$CO.H=CO*H
chroma$CO.Fl=CO*Fl
chroma$H.Fl=H*Fl
chroma$CO.H.Fl=CO*H*Fl
2*coef(lm(means~CO*H*Fl))
par(mfrow=c(2,2))
boxplot(means~CO,ylab="Finish",xlab="CO")
boxplot(means~H,ylab="Finish",xlab="H")
boxplot(means~Fl,ylab="Finish",xlab="Fl")
par(mfrow=c(1,1))
chroma.se=2*chroma.sp/sqrt(24)
chroma.se
chroma.se*qt(.975,16)