siteFrequency function

Get site frequency spectrum of a tree