data <- matrix(scan(),ncol=6,byrow=TRUE)
 0 0 0 0 1 10   0 0 0 0 1  15
 0 0 0 1 0 12   0 0 0 1 0  18
 0 0 1 0 0 18   0 0 1 0 0  13
 0 0 1 1 1 13   0 0 1 1 1  19
 0 1 0 0 0 15   0 1 0 0 0  11
 0 1 0 1 1  9   0 1 0 1 1  13
 0 1 1 0 1  9   0 1 1 0 1  16
 0 1 1 1 0 11   0 1 1 1 0  16
 1 0 0 0 0 13   1 0 0 0 0  18
 1 0 0 1 1 11   1 0 0 1 1  17
 1 0 1 0 1 14   1 0 1 0 1  20
 1 0 1 1 0 14   1 0 1 1 0  18
 1 1 0 0 1 10   1 1 0 0 1  14
 1 1 0 1 0 10   1 1 0 1 0  17
 1 1 1 0 0 13   1 1 1 0 0  17
 1 1 1 1 1 10   1 1 1 1 1  16

data
FracFact1  <- data.frame(data)
colnames(FracFact1)  <- c("a","b","c","d","e","y")
rm(data)

FracFact1$a  <- as.factor(FracFact1$a)
FracFact1$b  <- as.factor(FracFact1$b)
FracFact1$c  <- as.factor(FracFact1$c)
FracFact1$d  <- as.factor(FracFact1$d)
FracFact1$e  <- as.factor(FracFact1$e)

anova(lm(y ~ a*b*c*d, data = FracFact1))