See the introduction for an overview. Also see other analyses of this data. See the textbook for a full discussion.

Load the libraries:

library(ggplot2)
library(INLA)

Data

Load in and plot the data. We centre the Raven score and create unique labels for the classes within schools:

data(jsp, package="faraway")
jspr <- jsp[jsp$year==2,]
jspr$craven <- jspr$raven-mean(jspr$raven)
jspr$classch <- paste(jspr$school,jspr$class,sep=".")
ggplot(jspr, aes(x=raven, y=math))+xlab("Raven Score")+ylab("Math Score")+geom_point(position = position_jitter())

ggplot(jspr, aes(x=social, y=math))+xlab("Social Class")+ylab("Math Score")+geom_boxplot()