dta <- read.table("C:/Users/X510/Desktop/2020-04-13-In-Class-3-Trellis/beautyCourseEval.txt", header = T)[,c(1,2,7)] dta$courseID <- as.character(dta$courseID)
library(lattice) xyplot(eval ~ beauty | courseID, data=dta, xlab="Beauty judgement score", ylab="Average course evaluation score", type=c('p', 'g', 'r'), layout=c(6, 6), aspect=0.7, index.cond=list(c(9,31,5,17,4,8,12,29,7,11,24,15,14,18,1,25,16,20,27,13,10,23,6,3,26,22,28,2,21,30,19)))