data <- matrix(scan(),ncol=2,byrow=TRUE)
      1         8.2     1         8.5     1         8.4
      2         7.1     2         7.8     2         6.8
      3         8.5     3         8.8     3         8.5 
      4         8.2     4         8.3     4         9.1  
      5         7.4     5         6.7     5         7.2 
      6         7.9     6         8.2     6         8.0
      7         8.4     7         7.5     7         7.7

data
chicken <- data.frame(data)
colnames(chicken) <- c("restaurant","taste")
rm(data)

chicken$restaurant <- as.factor(chicken$restaurant)

 # to get variance component estimates

library(nlme)

chicken.lme  <- lme(taste ~ 1, data = chicken, random = ~ 1 | restaurant )

summary(chicken.lme)



 # to get the ANOVA table

anova(lm(taste ~ restaurant, data = chicken))