R codes for AEM/AES algorithm: AES.R

Instruction: Arguments of function aes (or aem):
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.


Example:

data: (genotypes of pools of N=2 individuals, 4 pools, 3 markers)
3 2 0
0 2 1
4 1 1
2 2 3


In R, run the following: (note: rho can be estimated from the data set using formula in the paper)

> aes(data, N=2, rho=estimated.rho)

An example