| plot.hdmod {lmfor} | R Documentation |
Diagnostic plot a Height-Diameter model residuals
Description
Plotting method for class hdmod
Usage
## S3 method for class 'hdmod'
plot(x, col.point = "blue", highlight = "red", standd = TRUE,
cex=1, corD=FALSE, ask=TRUE, ...)
Arguments
x |
A H-D model model fitted by |
col.point |
The color used for data points |
highlight |
The color used to higlight classes with mean sighnificantly different from zero |
standd |
Plot residuals against diameter standardized using the plot-specific mean and diameter ( |
cex |
See |
corD |
should predictions of random effects be plotted on mean diameter of the plot. |
ask |
ask before new plot. |
... |
Other arguments, currently ignored. |
Details
The function makes residual plots on a fitted H-D model, which can be used to explore whether the fixed part
satisfactorily models the shape of H-D models.
The residuals are plotted on diameters standardized at plot level (dsd) or on raw diameters (d) according to argument standd.
Here dsd=(d-D)/SDD, where d is tree diameter, D and SDD are the mean and standard deviation
of diameters on the plot in question. Using plot-specific standardized diameter ensures that e.g.,
the medium-sized trees of the plot are always in the middle of the plot, which provides
a better graph to explore the fit at the plot level in a dataset where the diameter range varies between plots.
Author(s)
Lauri Mehtatalo <lauri.mehtatalo@uef.fi>
See Also
The function plots model residuals on the required type of diameter and adds a whiskers plot
using mywhiskers with argument se=TRUE.
Examples
data(spati)
model<-fithd(spati$d,spati$h,spati$plot)
plot(model)