plot.MCMCBioGro: Plotting function fo the MCMCBioGro class
In serbinsh/biocro: Simulation of Energycane and Sugarcane

Description Usage Arguments Details Value See Also Examples

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.

  plot.MCMCBioGro(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"),
    ...)

`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.

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 Vmax 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.

A lattice plot.

MCMCBioGro