### Example
#data (3 loci/columns, 4 pools/rows, Each pool has N=2 individuals)

3 2 0 
0 2 1 
4 1 1 
2 2 3

## R commands:

genodata=matrix(c(3,2,0,  0,2,1,  4,1,1,  2,2,3), nrow=4, byrow=T)

N=2 # number of individuals in each pool

V=apply(genodata, 2, var) # variances of genotype data 
mu=apply(genodata,2,mean)/(2*N) #  allele frequencies

estimated.rho = sum(V)/sum(2*N*mu*(1-mu))-1 # estimated rho

aem(genodata, N=2) ## run the aem program (no variance inflation factor/inbreeding coefficient)

aes(genodata, N=2, rho= estimated.rho) ## run the aes program, with rho