Hey everyone,

I had a quick conceptual question about admixMap() implemented in the GENESIS package.

The documentation clearly mentions that in case of three ancestries (eg: European, African, and Native American), only two of the three need to be sent into admixDataList option of admixMap() function. But what if my local ancestry dosages come from just two ancestries (eg: European and African). In this case should I just send one GDS object and not two?

I ran them both ways and the summary stats seem to differ (Joint pval when two GDS objects are sent and single pvals when only one GDS object is sent) which makes sense. But I'm wondering if there is an advantage of sending both the GDS objects since I'm getting different results.

#My dosage GDS objects. My local ancestry estimates are two-way *European* & *African*.
gdsList = lapply(list("Dosage.EUR.chrAll.gds",  "Dosage.AFR.chrAll.gds"), GdsGenotypeReader)

#Is this correct?
genoIterators <- lapply(genoDataList, GenotypeBlockIterator)

#Or is this correct?
genoIterators <- lapply(genoDataList[1], GenotypeBlockIterator)

#Test the model.
myassoc <- admixMap(genoIterators, null.model)

Thank you for all the help!

Best, Shivam