Duarte
Regular Member
- Messages
- 2,771
- Reaction score
- 2,053
- Points
- 113
- Ethnic group
- Brazilian
- Y-DNA haplogroup
- R-FGC35133
- mtDNA haplogroup
- B2hh1
I believe there is an incompatibility regarding SNP IDs (different formats) between my raw data—which uses standard SNP IDs (e.g., rs12345)—and the archaeogenetic reference dataset (HO - Human Origins), which employs a proprietary SNP ID format. During the file merging process, the `plink --bmerge` command likely excluded SNPs due to strand mismatches (requiring a "flip") or other incompatibilities, resulting in a very low number of remaining SNPs in the merged file. When using the standard qpAdm input script (Target, Left Sources, Right Sources) with a merged dataset—specifically "MADA — Brazil_XXXXXXXXX + AADR v66 HO (PLINK)," which combines my raw data (submitted to MADAI) with the v66 HO (PLINK) dataset—execution halts with a series of errors:Although this qpAdm model may pass the statistical test, is not a good model, this estimated admixture is implausible.
Iberians, from whom you inherited most of your genome, don’t have ~50% direct (or indirect) Iran_GanjDareh_N admixture.
The use of a modern Iranian population (Iranian.DG) as an outgroup is incorrect, since modern Iranians are themselves admixed and are not ideal for distinguishing ancient Near Eastern ancestries.
This model is not historically or genetically correct. qpAdm tests only whether the proposed sources can statistically explain the target relative to the chosen outgroups. It does not test if the sources are genetically appropriate and if the model is historically realistic.
A more appropriate model would use populations directly ancestral to the Target. A statistically acceptable model is not necessarily a biologically meaningful one, and careful selection of sources and outgroups is just as important as obtaining a good p-value.
Calculating 0 f4 statistics for block 706 of 710...
Calculating 0 f4 statistics for block 707 of 710...
Calculating 0 f4 statistics for block 708 of 710...
Calculating 0 f4 statistics for block 709 of 710...
Calculating 0 f4 statistics for block 710 of 710...
Error in
Consequently, I can only obtain models that are statistically reliable but contextually (historically) unreliable when running the tool in "Rotation Analysis" mode using all non-missing SNPs on a pairwise population basis.
In short, I end up with randomly generated models that are statistically sound but reveal nothing about historical admixture.
For now, I'm giving up on this little experiment, lol.