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.

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