ungulates                package:ade4                R Documentation

_P_h_y_l_o_g_e_n_y _a_n_d _q_u_a_n_t_i_t_a_t_i_v_e _t_r_a_i_t_s _o_f _u_n_g_u_l_a_t_e_s.

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

     This data set describes the phylogeny of 18 ungulates as reported
     by Plabon et al. (1995). It also gives 4 traits corresponding to
     these 18 species.

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

     data(ungulates)

_F_o_r_m_a_t:

     'fission' is a list containing the 2 following objects : 

     _t_r_e is a character string giving the phylogenetic tree in Newick
          format. 

     _t_a_b is a data frame with 18 species and 4 traits    

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

     Variables of 'ungulates$tab' are the following ones : 
      afbw: is a numeric vector that describes the adult female body
     weight (g) 
      mnw: is a numeric vector that describes the male neonatal weight
     (g) 
      fnw: is a numeric vector that describes the female neonatal
     weight (g) 
      ls: is a numeric vector that describes the litter size 

_S_o_u_r_c_e:

     Data were obtained from Plabon, C., Gaillard, J.M., Loison, A.
     and Portier, A. (1995)  Is sex-biased maternal care limited by
     total maternal expenditure in polygynous ungulates?  _Behavioral
     Ecology and Sociobiology_, *37*, 311-319.

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

     data(ungulates)
     ung.phy <- newick2phylog(ungulates$tre)
     plot.phylog(ung.phy,clabel.l=1.25,clabel.n=0.75)
     ung.x <- log(ungulates$tab[,1])
     ung.y <- log((ungulates$tab[,2]+ungulates$tab[,3])/2)
     names(ung.x) <- names(ung.phy$leaves)
     names(ung.y) <- names(ung.x)
     plot(ung.x,ung.y)
     abline(lm(ung.y~ung.x))
     symbols.phylog(ung.phy,ung.x-mean(ung.x))
     dotchart.phylog(ung.phy,ung.x,cle=1.5,cno=1.5,cdot=1)
     orthogram(ung.x,ung.phy$Bscores,nrep=9999)
     ung.z <- residuals(lm(ung.y~ung.x))
     names(ung.z) <- names(ung.phy$leaves)
     dotchart.phylog(ung.phy,ung.z,cle=1.5,cno=1.5,cdot=1,ceti=0.75)
     orthogram(ung.z,ung.phy$Bscores,nrep=9999)

