rmaia

require(rstan)
require(geiger)
require(MCMCglmm)

# load data
data(geospiza)
dat <- geospiza$geospiza.data

# create fake sample of trees
tr <- drop.tip(geospiza$geospiza.tree, 'olivacea')

rmaia

require(rstan)
require(geiger)
require(MCMCglmm)

# load data
data(geospiza)

# get the inverse of the phylogenetic variance-covariance matrix
tr <- drop.tip(geospiza$geospiza.tree, 'olivacea')
invA <- solve(vcv.phylo(tr))

rmaia

require(ggplot2)
require(plyr)

# make toy data
set.seed(1)
dataf <- data.frame(
  x = rnorm(100),
  y = c(rnorm(50,-1), rnorm(50,1)),
  cat = rep(c("A","B"), each=50)
  )