data(BumpyGlossy)
plot(BumpyGlossy,
xlab = expression(paste("Surface ", S[ij],
"Gloss Level (i) and Bump level (j)")),
ylab = expression(paste("Surface ", S[kl],
"Gloss Level (k) and Bump level (l)"))
)
Run the code above in your browser using DataLab