trplot                 package:VGAM                 R Documentation

_T_r_a_j_e_c_t_o_r_y _P_l_o_t

_D_e_s_c_r_i_p_t_i_o_n:

     Generic function for a trajectory plot.

_U_s_a_g_e:

     trplot(object, ...)

_A_r_g_u_m_e_n_t_s:

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

_D_e_t_a_i_l_s:

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

_V_a_l_u_e:

     The value returned depends specifically on the methods function
     invoked.

_A_u_t_h_o_r(_s):

     Thomas W. Yee

_R_e_f_e_r_e_n_c_e_s:

     Yee, T. W. (2005) On constrained and unconstrained quadratic
     ordination. _Manuscript in preparation_.

_S_e_e _A_l_s_o:

     'trplot.qrrvglm', 'persp.qrrvglm', 'lvplot'.

_E_x_a_m_p_l_e_s:

     ## Not run: 
     data(hspider)
     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")

     ## End(Not run)

