Hipp metric

[willpearse]

require(ape)
w <- function(phy) {
  colSums(solve(vcv.phylo(phy)))
  }