hw1prob3 = read.table("e:\\FSU-school\\statapp2\\HW1prob3.txt", header = TRUE)
hw1prob3
fit_1 = aov(tensile ~ factor(alloy)/factor(casting)/factor(bar),data = hw1prob3 )

summary(fit_1)

fit_1$resid

eda1(fit_1$resid)