5.7 
> engineer<-read.table("C:/\--------------/prb0507.dat",header=T)
> engineer
  run x1 x2 x3  y
1   1 -1 -1 -1 60
2   2  1 -1 -1 61
3   3 -1  1 -1 54
4   4  1  1 -1 75
5   5 -1 -1  1 58
6   6  1 -1  1 61
7   7 -1  1  1 55
8   8  1  1  1 75

> interaction.plot(engineer$x1,engineer$x2,engineer$y)
> interaction.plot(engineer$x1,engineer$x3,engineer$y)
> interaction.plot(engineer$x2,engineer$x3,engineer$y)
> fit<-aov(y~factor(x1)*factor(x2)*factor(x3),data=engineer)
> summary(fit)
> model.tables(fit,type="effects",se=F)
> fit1<-lm(y~x1+x2 +x1:x2, data=engineer)
> summary(fit1)

5.10
> chemical<-read.table("C:/--==============/prb0510.dat",header=T)
> chemical
  run temp conc rate yield
1   1   50    6   60    54
2   2   60    6   60    57
3   3   50   10   60    69
4   4   60   10   60    70
5   5   50    6  100    55
6   6   60    6  100    54
7   7   50   10  100    80
8   8   60   10  100    81
> interaction.plot(chemical$temp,chemical$conc,chemical$yield)
> interaction.plot(chemical$temp,chemical$rate,chemical$yield)
> interaction.plot(chemical$conc,chemical$rate,chemical$yield)
> fit<-aov(yield~factor(temp)*factor(conc)*factor(rate),data=chemical)
> summary(fit)
> model.tables(fit,type="effects",se=T) 
 
> fit<-lm(yield~conc+rate+conc:rate, data=chemical)
> summary(fit)
> fit<-lm(yield~rate+conc:rate, data=chemical)
> summary(fit)


5.15 > lab<-read.table("C:/=====/515.txt",header=T)
> lab
   temp conc time  y
1    -1   -1   -1 80
2     1   -1   -1 65
3    -1    1   -1 69
4     1    1   -1 74
5    -1   -1    1 81
6     1   -1    1 84
7    -1    1    1 91
8     1    1    1 93
9    -1   -1   -1 62
10    1   -1   -1 63
11   -1    1   -1 73
12    1    1   -1 80
13   -1   -1    1 79
14    1   -1    1 86
15   -1    1    1 93
16    1    1    1 93

> interaction.plot(lab$temp,lab$conc,lab$y)
> interaction.plot(lab$temp,lab$time,lab$y)
> interaction.plot(lab$conc,lab$time,lab$y)
> fit<-aov(y~(factor(temp)*factor(conc)*factor(time)),data=lab)
> summary(fit)
> model.tables(fit,type="effects",se=T)
> fit<-lm(y~conc+time, data=lab)
> summary(fit)