Blocking in a 2^k Factorial Design
getwd()
#Question:3
#Reading the data:
A<-c(-1,1,-1,1,-1,1,-1,1)
B<-c(-1,-1,1,1,-1,-1,1,1)
AB<-c(1,-1,-1,1,1,-1,-1,1)
C<-c(-1,-1,-1,-1,1,1,1,1)
AC<-c(1,-1,1,-1,-1,1,-1,1)
BC<-c(1,1,-1,-1,-1,-1,1,1)
ABC<-c(-1,1,1,-1,1,-1,-1,1)
obs<-c(22,32,35,55,44,40,60,39)
Data<-data.frame(A,B,AB,C,AC,BC,ABC)
#Half Normal Plot:
library(DoE.base)
?halfnormal
Model <- lm(obs~A*B*C,data = Data)
coef(Model)
halfnormal(Model)