Function to plot squared residuals and the error model for an mkin object
This function plots the squared residuals for the specified subset of the observed variables from an mkinfit object. In addition, one or more dashed line(s) show the fitted error model. A combined plot of the fitted model and this error model plot can be obtained with plot.mkinfit
using the argument show_errplot = TRUE
.
mkinerrplot( object, obs_vars = names(object$mkinmod$map), xlim = c(0, 1.1 * max(object$data$predicted)), xlab = "Predicted", ylab = "Squared residual", maxy = "auto", legend = TRUE, lpos = "topright", col_obs = "auto", pch_obs = "auto", frame = TRUE, ... )
object
: A fit represented in an mkinfit
object.obs_vars
: A character vector of names of the observed variables for which residuals should be plotted. Defaults to all observed variables in the modelxlim
: plot range in x direction.xlab
: Label for the x axis.ylab
: Label for the y axis.maxy
: Maximum value of the residuals. This is used for the scaling of the y axis and defaults to "auto".legend
: Should a legend be plotted?lpos
: Where should the legend be placed? Default is "topright". Will be passed on to legend
.col_obs
: Colors for the observed variables.pch_obs
: Symbols to be used for the observed variables.frame
: Should a frame be drawn around the plots?...
: further arguments passed to plot
.Nothing is returned by this function, as it is called for its side effect, namely to produce a plot.
## Not run: model <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO")) fit <- mkinfit(model, FOCUS_2006_D, error_model = "tc", quiet = TRUE) mkinerrplot(fit) ## End(Not run)
mkinplot
, for a way to plot the data and the fitted lines of the mkinfit object.
Johannes Ranke