library(lattice)

mtcars$transmission <- factor(mtcars$am, levels=c(0,1),

labels=c("Automatic", "Manual"))

panel.smoother <- function(x, y) {

panel.grid(h=-1, v=-1)

panel.xyplot(x, y)

panel.loess(x, y)

panel.abline(h=mean(y), lwd=2, lty=2, col="darkgreen")

}

xyplot(mpg~disp|transmission,data=mtcars,

scales=list(cex=.8, col="red"),

panel=panel.smoother,

xlab="Displacement", ylab="Miles per Gallon",

main="MPG vs Displacement by Transmission Type",

sub = "Dotted lines are Group Means", aspect=1)

在上面的代碼中有幾個地方需要指出。 panel.xyplot()函數(shù)畫出了個別點， panel.

loess()函數(shù)在每個面板圖中畫出了非參數(shù)擬合曲線。 panel.abline()函數(shù)在調(diào)解變量的每個水平中添加了水平參考線（mpg的均值）。（如果你用h=mean(mtcars$mpg)代替h=mean(y)，在

整個樣本中將產(chǎn)生以mpg均值為基礎(chǔ)的單個參考線。）scales=選項呈現(xiàn)大小為默認字體80%的紅

色刻度注釋（坐標(biāo)軸數(shù)字和刻度線）。