plot.MCMCBioGro. BioCro 0.93

Usage

"plot"(x, x2 = NULL, x3 = NULL, plot.kind = c("rss", "OF", "RF", "OFT", "trace", "density"), type = c("l", "p"), coef = 1, cols = c("blue", "green", "red", "magenta", "black", "purple"), ...)

Arguments

x: Object of class MCMCBioGro
x2: Optional object of class MCMCBioGro
x3: Optional object of class MCMCBioGro
plot.kind: Kind of plot. See details.
type: Point of line as in xyplot
coef: Whether to plot dry biomass partitioning coeficients (2) or Vmax and alpha (1).
cols: Colors. Modify if they don't suit you.
...: Additional arguments passed to the underlying xyplot function. Some can be really useful. See details.

Value

A lattice plot.

Description

Powerful plotting function to make a variety of plots regarding the MCMCBioGro class (output). It plots the residual sum of square progression, observed vs. fitted, residuals vs. fitted, trace of the coefficients and density.

Details

Kind of plots that can be produced

rss: Residual Sum of Squares progression. OF: Observed vs. Fitted RF: Residual vs. Fitted OFT: Observed and Fitted with thermal time as the x-axis. trace: trace for the parameters density: density for the parameters

biomass used. If coef=1 then Vmax instead. At this point it is highly and alpha are not optimized at this stage so shouldn't be plotted either (it will be a straight line if not

To choosing a subset of the 24 dry biomass partitioning coefficients use the subset option as you would in the xyplot using nams as the name. For example, plot(x,plot.kind="trace",subset=nams=="kLeaf_1") will select it for the leaf at the first phenological stage.

Examples

## See the MCMCBioGro function