DF = read.csv('../../Data/ASCII_Comma/Chapter_12/waferlayers.txt', quote='\'', header=TRUE)
DF$Furnace = as.factor(DF$Furnace)

fit = aov(Thickness ~ Furnace * Wafer.Type, data=DF)
print(summary(fit))

fit = aov(Thickness ~ Furnace + Wafer.Type, data=DF)
print(summary(fit))

#postscript("../../WriteUp/Graphics/Chapter12/prob_29_interactionplot.eps", onefile=FALSE, horizontal=FALSE)
interaction.plot(DF$Furnace, DF$Wafer.Type, DF$Thickness, type='b', col=c('red', 'blue', 'green'), pch=c(16, 18, 20), xlab='Furnace number', ylab='Thickness' , main='Interaction between Furnace and Wafer Type')
#dev.off()