Right now, when I was performing some preliminary data analysis (data obtained during my Ms that I was so late to publish the paper), I came across with something interesting. As we have already known GENEPOP and FREENA are brother softwares developed for the analysis of microsatellite and some co-dominante markers (as Allozymes and RFLPs). Also, as we have already known, as a first step analysis we must have to calculate the null allele frequency on each loci through each population to confirm the presence of null alleles. The most powerful way to estimate the frequency of the null allele is using maximum likelihood estimator with EM algorithm. Both softwares have this approach (sure!!), and in both you can assign the null alleles (if you previously know) as 999 (or “99”) and the genotyping failures as 000 (or “00”). As a good approach if you really don´t know how occurred you must assign 000. Of course we expect to have the same null allele’s frequency estimative (or approximately). For the majority of the loci I had the same estimative (the MLE with EM is sufficiently robust) but for some loci I had a super estimative of null alleles. It occurs with the GENEPOP. I have seem that for some locus the highest allele (>pb) have taken as a null. For that reason at the output file the null allele is the bigger one (considered as null). I have been using both, for some analysis as null allele estimative and of course for FST with ENA correction. But I prefer GENEPOP to other analysis as HW test specially.
Nenhum comentário:
Postar um comentário