man/corClasses.Rd

   1 \name{corClasses}
   2 \alias{corClasses}
   3 \alias{corPhyl}
   4 \title{Phylogenetic Correlation Structures}
   5 \description{
   6   Classes of phylogenetic correlation structures (\code{"corPhyl"})
   7   available in \pkg{ape}.
   8
   9   \itemize{
  10     \item{corBrownian}{Brownian motion model (Felsenstein 1985)}
  11     \item{corMartins}{The covariance matrix defined in Martins and Hansen
  12       (1997)}
  13     \item{corGrafen}{The covariance matrix defined in Grafen (1989)}
  14     \item{corPagel}{The covariance matrix defined in Freckelton et al. (2002)}
  15     \item{corBlomberg}{The covariance matrix defined in Blomberg et al. (2003)}
  16   }
  17
  18   See the help page of each class for references and detailed
  19   description.
  20 }
  21 \seealso{
  22   \code{\link[nlme]{corClasses}} and \code{\link[nlme]{gls}} in the
  23   \pkg{nlme} librarie, \code{\link{corBrownian}},
  24   \code{\link{corMartins}}, \code{\link{corGrafen}},
  25   \code{\link{corPagel}}, \code{\link{corBlomberg}}
  26 }
  27 \author{Julien Dutheil \email{julien.dutheil@univ-montp2.fr}, Emmanuel
  28   Paradis}
  29 \examples{
  30 library(nlme)
  31 cat("((((Homo:0.21,Pongo:0.21):0.28,",
  32 "Macaca:0.49):0.13,Ateles:0.62):0.38,Galago:1.00);",
  33 file = "ex.tre", sep = "\n")
  34 tree.primates <- read.tree("ex.tre")
  35 X <- c(4.09434, 3.61092, 2.37024, 2.02815, -1.46968)
  36 Y <- c(4.74493, 3.33220, 3.36730, 2.89037, 2.30259)
  37 unlink("ex.tre") # delete the file "ex.tre"
  38 m1 <- gls(Y~X, correlation=corBrownian(1, tree.primates))
  39 summary(m1)
  40 m2 <- gls(Y~X, correlation=corMartins(1, tree.primates))
  41 summary(m2)
  42 corMatrix(m2$modelStruct$corStruct)
  43 m3 <- gls(Y~X, correlation=corGrafen(1, tree.primates))
  44 summary(m3)
  45 corMatrix(m3$modelStruct$corStruct)
  46 }
  47 \keyword{models}
  48