R script for grabbing sister pairs across a sample of trees

This R script grabs all sister pairs from across a sample of trees,computes the average branch length (dist), computes the posterior probability (pp = proportion of sampled trees with a given sister pair), and grabs the mating system of the sister pair.