Data > require(DoE.wrapper) > Data Data.lm <- lm(y~.^4, data=Data)
Data
> require(DoE.wrapper)
> Data <- data.frame(
+ qua.design(
+ fac.design(nlevels = 2, nfactors = 4, randomize = FALSE, factor.names = list(A=c(-1,1), B=c(-1,1), C=c(-1,1), D=c(-1,1))), quantitative=”all”), y = c(1.68,1.98,3.28,3.44,4.98,5.70,9.97,9.07,2.07,2.44,4.09,4.53,7.77,9.43,11.75,16.30))
creating full factorial with 16 runs …
> Data.lm <- lm(y~.^4, data=Data)