plot.hdmod: Diagnostic plot a Height-Diameter model residuals

In lmfor: Functions for Forest Biometrics

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 fithd.
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 (standd=TRUE) or against raw diameter (standd=FALSE)
cex	See par
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 <[e-mail]>

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)

lmfor documentation built on April 30, 2022, 1:08 a.m.