Image or Perspective plot of two-way model interactions. Ranges of two specified predictor variables are plotted on X and Y axis, and fitted model values are plotted on the Z axis. The remaining predictor variables are fixed at their mean (for continuous predictors) or their most common value (for categorical predictors).
model.interaction.plot(model.obj = NULL, x = NULL, y = NULL,
response.category=NULL, quantiles=NULL, all=FALSE, obs=1, qdata.trainfn = NULL,
folder = NULL, MODELfn = NULL, PLOTfn = NULL, pred.means = NULL,
xlab = NULL, ylab = NULL, x.range = NULL, y.range = NULL,
z.range = NULL, ticktype = "detailed", theta = 55, phi = 40,
smooth = "none", plot.type = NULL, device.type = NULL,
res=NULL, jpeg.res = 72, device.width = 7, device.height = 7,
units="in", pointsize=12, cex=par()$cex,
col = NULL, xlim = NULL, ylim = NULL, zlim = NULL, ...)
R model object. The model object to use for prediction. The model object must be of type "RF" (random forest), "QRF" (quantile random forest), or "CF" (conditional forest). The ModelMap package does not currently support SGB models.
String or Integer. Name of predictor variable to be plotted on the x axis. Alternativly, can be a number indicating a variable name from predList.
String or Integer. Name of predictor variable to be plotted on the y axis. Alternatively, can be a number indicating a variable name from predList.
String. Used for categorical response models. Specify which category of response variable to use. This category's probabilities will be plotted on the z axis.
Numeric. Used for QRF models. Specify which quantile of response variable to use. This quantile will be plotted on the z axis. Note: unlike other functions model.interaction.plot will only use a single quantile. If quantiles is a vector only the first value will be used.
Logical. Used for QRF models. A logical value. all=TRUE uses all observations for prediction. all=FALSE uses only a certain number of observations per node for prediction (set with argument obs). The default is all=FALSE.
Numeric. Used for QRF models. An integer number. Determines the maximal number of observations per node to use for prediction. The input is ignored for all=TRUE. The default is obs=1.
String. The name (full path or base name with path specified by folder) of the training data file used for building the model (file should include columns for both response and predictor variables). The file must be a comma-delimited file *.csv with column headings. qdata.trainfn can also be an R dataframe. If predictions will be made (predict = TRUE or map=TRUE) the predictor column headers must match the names of the raster layer files, or a rastLUT must be provided to match predictor columns to the appropriate raster and band. If qdata.trainfn = NULL (the default), a GUI interface prompts user to browse to the training data file.
String. The folder used for all output. Do not add ending slash to path string. If folder = NULL (default), a GUI interface prompts user to browse to a folder. To use the working directory, specify folder = getwd().
String. The file name used to save the generated model object, only used if PLOTfn = NULL. If MODELfn is supplied and If PLOTfn = NULL, a graphical file name is generated by pasting MODELfn_plot.type_x.name_y.name. If PLOTfn = NULL and MODELfn = NULL, a default name is generated by pasting model.type_response.type_response.name_plot.type_x.name_y.name. The filename can be the full path, or it can be the simple basename, in which case the output will be to the folder specified by folder.
String. The file name to use to save the generated graphical plots. The filename can be the full path, or it can be the simple basename, in which case the output will be to the folder specified by folder.
Vector. Allows specification of values for other predictor variables. If Null, other predictors are set to their mean value (for continuous predictors) or their most common value (for factored predictors).
String. Allows manual specification of the x label.
String. Allows manual specification of the y label.
Vector. Manual range specification for the x axis. Alternate argument name for xlim. Use one or the other. Do not provide both x.range and xlim.
Vector. Manual range specification for the y axis. Alternate argument name for ylim. Use one or the other. Do not provide both y.range and ylim.
Vector. Manual range specification for the z axis. Alternate argument name for zlim. Use one or the other. Do not provide both z.range and zlim.
Character: "simple" draws just an arrow parallel to the axis to indicate direction of increase; "detailed" (default) draws normal ticks as per 2D plots. If X or y is factored, ticks will be drawn on both axes.
Numeric. Angles defining the viewing direction. theta gives the azimuthal direction.
Numeric. Angles defining the viewing direction. phi gives the colatitude.
String. controls smoothing of the predicted surface. Options are "none" (default), "model" which uses a glm model to smooth the surface, and "average" which applies a 3x3 smoothing average. Note: smoothing is not appropriate if X or y is factored.
Character. "persp" gives a 3-D perspective plot. "image" gives an image plot.
String or vector of strings. Model validation. One or more device types for graphical output from model validation diagnostics.
Current choices:
"default" | default graphics device | |||
"jpeg" | *.jpg files | |||
"none" | no graphics device generated | |||
"pdf" | *.pdf files | |||
"postscript" | *.ps files | |||
"win.metafile" | *.emf files |
Integer. Model validation. Pixels per inch for jpeg, png, and tiff plots. The default is 72dpi, good for on screen viewing. For printing, suggested setting is 300dpi.
Integer. Model validation. Deprecated. Ignored unless res not provided.
Integer. Model validation. The device width for diagnostic plots in inches.
Integer. Model validation. The device height for diagnostic plots in inches.
Model validation. The units in which device.height and device.width are given. Can be "px" (pixels), "in" (inches, the default), "cm" or "mm".
Integer. Model validation. The default pointsize of plotted text, interpreted as big points (1/72 inch) at res ppi
Integer. Model validation. The cex for diagnostic plots.
Vector. Color table to use for image plots ( see help file on image for details).
Vector. X limits. Alternate argument name for x.range. Use one or the other. Do not provide both x.range and xlim.
Vector. Y limits. Alternate argument name for y.range. Use one or the other. Do not provide both y.range and ylim.
Vector. Z limits. Alternate argument name for z.range. Use one or the other. Do not provide both z.range and zlim.
additional graphical parameters (see par).
Elizabeth Freeman
This function provides a diagnostic plot useful in visualizing two-way interactions between predictor variables. Two of the predictor variables from the model are used to produce a grid of possible combinations of predictor values over the range of both variables. The remaining predictor variables from the model are fixed at either their means (for continuous predictors) or their most common value (for categorical predictors). Model predictions are generated over this grid and plotted as the z axis.
This function works with both continuous and categorical predictors, though the perspective plot should be interpreted with care for categorical predictors. In particular, the smooth option is not appropriate if either of the two selected predictor variables is categorical.
For categorical response models, a particular value must be specified for the response using the response.category argument.
This function is adapted from gbm.perspec version 2.9 April 2007, J Leathwick/J Elith. See appendix S3 from:
Elith, J., Leathwick, J. R. and Hastie, T. (2008). A working guide to boosted regression trees. Journal of Animal Ecology. 77:802-813.