lmfor (version 1.6)

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

Description

Plotting method for class hdmod

Usage

# S3 method for 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.

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

# NOT RUN {
data(spati)

model<-fithd(spati$d,spati$h,spati$plot)

plot(model)
# }