plotqrrvglm               package:VGAM               R Documentation

_M_o_d_e_l _D_i_a_g_n_o_s_t_i_c _P_l_o_t_s _f_o_r _Q_R_R-_V_G_L_M_s

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

     The residuals of a QRR-VGLM are plotted for model diagnostic
     purposes.

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

     plotqrrvglm(object, 
                 rtype = c("pearson", "response", "deviance", "working"), 
                 ask = FALSE, 
                 main = paste(Rtype, "residuals vs latent variable(s)"), 
                 xlab = "Latent Variable", 
                 ITolerances = object@control$EqualTolerances, ...)

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

  object: An object of class '"qrrvglm"'. 

   rtype: Character string giving residual type. By default, the first
          one is chosen. 

     ask: Logical. If 'TRUE', the user is asked to hit the return key
          for the next plot. 

    main: Character string giving the title of the plot. 

    xlab: Character string giving the x-axis caption. 

ITolerances: Logical. This argument is fed into 'Coef(object,
          ITolerances=ITolerances)'. 

     ...: Other plotting arguments (see 'par'). 

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

     Plotting the residuals can be potentially very useful for checking
     that the model fit is adequate.

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

     The original object.

_N_o_t_e:

     An ordination plot of a QRR-VGLM can be obtained by
     'lvplot.qrrvglm'.

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

     Thomas W. Yee

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

     Yee, T. W. (2004) A new technique for maximum-likelihood canonical
     Gaussian ordination. _Ecological Monographs_, *74*, 685-701.

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

     'lvplot.qrrvglm', 'cqo'.

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

     ## Not run: 
     # QRR-VGLM on the hunting spiders data
     # This is computationally expensive
     data(hspider)
     set.seed(111)  # This leads to the global solution
     # hspider[,1:6]=scale(hspider[,1:6]) # Standardize the environmental variables
     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)
     par(mfrow=c(3,4)) 
     plot(p1, rtype="d", col="blue", pch=4, las=1)
     ## End(Not run)

