VGAM (version 0.8-1)

trplot: Trajectory Plot

Description

Generic function for a trajectory plot.

Usage

trplot(object, ...)

Arguments

object

An object for which a trajectory plot is meaningful.

...

Other arguments fed into the specific methods function of the model. They usually are graphical parameters, and sometimes they are fed into the methods function for Coef.

Value

• The value returned depends specifically on the methods function invoked.

Details

Trajectory plots can be defined in different ways for different models. Many models have no such notion or definition.

For quadratic and additive ordination models they plot the fitted values of two species against each other (more than two is theoretically possible, but not implemented in this software yet).

References

Yee, T. W. (2005) On constrained and unconstrained quadratic ordination. Manuscript in preparation.

See Also

trplot.qrrvglm, persp.qrrvglm, lvplot.

Examples

hspider[,1:6] = scale(hspider[,1:6]) # Standardized environmental vars
set.seed(123)

p1 = cqo(cbind(Alopacce, Alopcune, Alopfabr, Arctlute, Arctperi,
               Auloalbi, Pardlugu, Pardmont, Pardnigr, Pardpull,
               Trocterr, Zoraspin) ~
         WaterCon + BareSand + FallTwig + CoveMoss + CoveHerb + ReflLux,
         fam = quasipoissonff, data = hspider, Crow1positive=FALSE)

nos = ncol(p1@y)
clr = 1:nos   # OR (1:(nos+1))[-7]  to omit yellow

trplot(p1, whichSpecies=1:3, log="xy",
       col=c("blue","red","green"), lwd=2, label=TRUE) -> ii
legend(0.00005, 0.3, paste(ii$species[,1], ii$species[,2], sep=" and "),
       lwd=2, lty=1, col=c("blue","red","green"))
abline(a=0, b=1, lty="dashed")