R code for AEM/AES algorithm

AES.R

Instruction:

Arguments of function aes (or aem) are

data: pooled genotype data, T rows (T=number of pools), q columns (q=number of loci);

rho: overdispersion factor, which can be estimated from data or from genome data;

N: pool size, number of individuals in each pool.