data <- matrix(scan(),ncol=6,byrow=TRUE)
1 1 1 1 1 14   1 1 1 0 0 11.5  1 1 0 1 0 12  1 1 0 0 1 12.5 
2 3 1 1 1 10   2 3 1 0 0 8     2 3 0 1 0 7   2 3 0 0 1 9 
1 2 1 1 0 13   1 2 1 0 1 14    1 2 0 1 1 13  1 2 0 0 0 10 
2 4 1 1 0 13   2 4 1 0 1 15    2 4 0 1 1 14  2 4 0 0 0 11 

data

Burgers1 <- data.frame(data)
colnames(Burgers1) <- c("rep","person","patty","cheese","mayo","taste")

Burgers1$rep       <- as.factor(Burgers1$rep)
Burgers1$person       <- as.factor(Burgers1$person)
Burgers1$patty       <- as.factor(Burgers1$patty)
Burgers1$mayo       <- as.factor(Burgers1$mayo)
Burgers1$cheese       <- as.factor(Burgers1$cheese)


aov(taste ~ rep + Error(rep/person) + patty*mayo*cheese, data = Burgers1)