An Evolutionary Framework for Association Testing in Resequencing Studies
R code to estimate the combined rho model. This is absolutely not intended to be production code, and you should probably not try to use it. I'm writting an R package now that does it smarter and better.
SAS code implementing split rho calculation in a subset of the population. Variant data for non-synonymous snps. Full dataset not provided due to privacy restrictions. "mi" and "shat" are variant counts and estimated fitness burden per class of SNP.