data <- matrix(scan(),ncol=6,byrow=TRUE)
 0 0 0 0 1 41.5   1 0 0 0 0 39.6   0 1 0 0 0 43.9   1 1 0 0 1 38.8 
 0 0 1 0 0 48.7   1 0 1 0 1 52.0   0 1 1 0 1 55.8   1 1 1 0 0 43.2 
 0 0 0 1 0 39.5   1 0 0 1 1 42.6   0 1 0 1 1 44.0   1 1 0 1 0 33.8 
 0 0 1 1 1 53.6   1 0 1 1 0 48.1   0 1 1 1 0 51.3   1 1 1 1 1 48.7  

data


table124  <- data.frame(data)
colnames(table124)  <- c("a","s","m","c","t","y")
rm(data)


table124$a  <- as.factor(table124$a)
table124$s  <- as.factor(table124$s)
table124$m  <- as.factor(table124$m)
table124$c  <- as.factor(table124$c)
table124$t  <- as.factor(table124$t)


table124.aov  <-  aov(lm(y ~ (a+s+m+c+t)^2, data = table124))

library(gplots)
qqnorm(table124.aov)

 # with this command, you can click on and label points
 # to get out of this mode, you can select the "R Console" window and hit the Escape button
if (dev.interactive()) qqnorm(table124.aov, label=TRUE)


 # now fitting the model with only "active-appearing" effects

table124.lm2 <- lm(y ~ a +m +t +s +a:s, data = table124)

anova(table124.lm2)

par(mfrow=c(2,2))  
plot(table124.lm2)
par(mfrow=c(1,1))  


 # many other packages have features for FF designs, FrF2 is one of them

library(FrF2)

FrF2(nruns=4,nfactors=3,factor.names=list(Mayo=c("Regular","Spicy"),Patty=c("Veggie","Beef"),
                         Cheese=c("Absent","Present")),randomize=TRUE)