prob5_15<-read.table("E:/FSU-FALL07-School/Fall07-school/STA5166/hw5/prb0515.dat",header=T)
prob5_15


#prob5_15= data.frame(apply(prob5_15_1[,2:4],2, as.factor), y= prob5_15_1[,5])

par(mfrow=c(2,2))
interaction.plot(prob5_15$temp,prob5_15$conc,prob5_15$y)
interaction.plot(prob5_15$temp,prob5_15$time,prob5_15$y)
interaction.plot(prob5_15$conc,prob5_15$time,prob5_15$y)


fit1 = aov(y~factor(temp)*factor(conc)*factor(time), data=prob5_15)
summary(fit1)
par(mfrow=c(2,2))
model.tables(fit1, type="effects", se=T)
qqnorm(fit1$fitted)





fit2<-lm(y~factor(conc)+factor(time)+ factor(conc):factor(time),data=prob5_15)

summary(fit2)



fit3<-lm(y~factor(conc)+factor(time),data=prob5_15)

summary(fit3)

par(mfrow=c(1,1))
plot(prob5_15$y,fit1$resid) 

prob5_15_2<-read.table("E:/FSU-FALL07-School/Fall07-school/STA5166/hw5/prb0515_2.dat",header=T)
prob5_15_2

fit1 = aov(y~factor(temp)*factor(conc)*factor(time), data=prob5_15_2)
summary(fit1)
plot(prob5_15_2$y,fit1$resid)
qqnorm(fit1$fitted)
source("C:\\Documents and Settings\\Jaime\\Desktop\\FALL07\\STA5166\\hw1eda1-2.txt")
par(mfrow=c(2,2))
eda1(fit1$fitted)