library(earth) library(gam) data(airquality) airquality <- na.omit(airquality) # plotmo doesn't know how to deal with NAs plotmo(gam(Ozone^(1/3) ~ lo(Solar.R) + lo(Wind, Temp), data=airquality), caption="plotmo example 2: GAM model", pt.pch=20, pt.col="blue", type2="image", ngrid2=100, image.col=terrain.colors(100))