R code for Ex 5.19 (3 factor factorial designs)

R code for Ex 5.19 (3 factor factorial designs)
# R code for 3 factor factorial design Ex 5.19
data=read.table("C:/Users/Mihinda/Desktop/ex519.txt", header=T) #the data file
data$CycleTime

Output
> # R code for 3 factor factorial design Ex 5.19
> data=read.table("C:/Users/Mihinda/Desktop/ex519.txt", header=T) #the
data file
> data$CycleTime data$Operator data$Temperature par(mfrow=c(2,2))
> plot.design(data)
> interaction.plot(data$CycleTime,data$Operator,data$Score)
> interaction.plot(data$CycleTime,data$Temperature,data$Score)
> interaction.plot(data$Operator,data$Temperature,data$Score)
mean of Score
29 30 31 32 33 34 35
50
2
3
40
60 1
CycleTime
350
300
Temperature
mean of data$Score
26 28 30 32 34 36
40 50 60
data$Operator
2
1
3
Factors
data$CycleTime
mean of data$Score
28 30 32 34
data$Temperature
350
300
mean of data$Score
30 32 34 36
data$Temperature
350
300
40 50 60
data$CycleTime
1 2 3
data$Operator

fit anova(fit)
Analysis of Variance Table
Response: Score
Df Sum Sq Mean Sq F value Pr(>F)
Operator 2 261.33 130.667 39.8644 7.439e-10 ***
CycleTime 2 436.00 218.000 66.5085 8.141e-13 ***
Temperature 1 50.07 50.074 15.2768 0.0003934 ***
Operator:CycleTime 4 355.67 88.917 27.1271 1.982e-10 ***
Operator:Temperature 2 11.26 5.630 1.7175 0.1938948
CycleTime:Temperature 2 78.81 39.407 12.0226 0.0001002 ***
Operator:CycleTime:Temperature 4 46.19 11.546 3.5226 0.0158701 *
Residuals 36 118.00 3.278
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1