library(simex)
load(url("http://192.38.117.59/~linearpredictors/datafiles/simex_model2.Rdata"))

par(mfrow = c(1,3))
plot(simex.model2,
     ## decide which of the plots to make
     show = c(T,T,T), 
     ## labels for each of the plots
     ylab = c("Intercept","Effect of BMI","Effect of log10(vitamin D intake)"),
     xlab = rep(expression(1 + lambda), 3)
     )