#Partial dependence plot is a plot of the marginal effect of a variable on the class probability (classification) or response (regression).

imp <- importance(sahdrf)
impvar <- rownames(imp)[order(imp[, 1], decreasing=TRUE)]

par(mfrow=c(3, 3))

for (i in seq_along(impvar)) {
    partialPlot(sahdrf, sahd, impvar[i], xlab=impvar[i],
                main=paste("Partial Dependence on", impvar[i]),
        ylim=c(0, 1) )
}