varcomp function

Estimates variance components for each allele of a locus