Originally Posted by
Zoro
Totally agree. I wouldn't call him an idiot. Clearly he's not a scientist if he's relying exclusively on G25. I would just call him a misinformed amateur who doesn't understand the issues with Davidski's G25.
The Hurrian he's referring to I think is Haji-Firuz-IA ! and the the Archemenid I believe is Hasanlu-IA. Clearly the G25 is giving totally off results. Davidski's G25 is a mish mosh of modern and ancient samples. Obviously the moderns he has in the run to generate coordinates are having all sorts of effects on the test samples.
Anyways if he was a scientist he would use qpAdm and if he had he would have seen the following for Haji-Firuz-IA (courtesy EurasianDNA)
numthreads: 8
## qpAdm version: 634
seed: 1295334595
left pops:
Iran_IA_HajjiFiruz
Iran_C_HajjiFiruz
Kazakhstan_Sarmatian.SG
right pops:
Mbuti.DG
Russia_Ust_Ishim.DG
Onge_1000G
Russia_MA1_HG.SG
Chinese_Han_1000G
Karitiana.DG
Anatolia_N
0 Iran_IA_HajjiFiruz 1
1 Iran_C_HajjiFiruz 4
2 Kazakhstan_Sarmatian.SG 2
3 Mbuti.DG 4
4 Russia_Ust_Ishim.DG 1
5 Onge_1000G 6
6 Russia_MA1_HG.SG 1
7 Chinese_Han_1000G 211
8 Karitiana.DG 3
9 Anatolia_N 29
jackknife block size: 0.050
snps: 709309 indivs: 262
number of blocks for block jackknife: 711
## ncols: 709309
coverage: Iran_IA_HajjiFiruz 292807
coverage: Iran_C_HajjiFiruz 624648
coverage: Kazakhstan_Sarmatian.SG 474638
coverage: Mbuti.DG 690580
coverage: Russia_Ust_Ishim.DG 654138
coverage: Onge_1000G 704593
coverage: Russia_MA1_HG.SG 493659
coverage: Chinese_Han_1000G 709306
coverage: Karitiana.DG 690550
coverage: Anatolia_N 709269
dof (jackknife): 639.251
numsnps used: 155390
codimension 1
f4info:
f4rank: 1 dof: 5 chisq: 2.444 tail: 0.784873416 dofdiff: 7 chisqdiff: -2.444 taildiff: 1
B:
scale 1.000
Russia_Ust_Ishim.DG 0.201
Onge_1000G 0.337
Russia_MA1_HG.SG 1.217
Chinese_Han_1000G 0.919
Karitiana.DG 1.512
Anatolia_N -1.111
A:
scale 723.320
Iran_C_HajjiFiruz -0.613
Kazakhstan_Sarmatian.SG 1.274
full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 5 chisqdiff: 2.444 taildiff: 0.784873416
B:
scale 1.000 1.000
Russia_Ust_Ishim.DG 0.200 0.200
Onge_1000G 0.330 0.790
Russia_MA1_HG.SG 1.203 -0.685
Chinese_Han_1000G 0.901 1.264
Karitiana.DG 1.520 -0.226
Anatolia_N -1.132 1.794
A:
scale 746.734 3511.876
Iran_C_HajjiFiruz -0.755 -1.065
Kazakhstan_Sarmatian.SG 1.196 -0.931
best coefficients: 0.675 0.325
Jackknife mean: 0.675276889 0.324723111
std. errors: 0.069 0.069
error covariance (* 1000000)
4774 -4774
-4774 4774
summ: Iran_IA_HajjiFiruz 2 0.784873 0.675 0.325 4774 -4774 4774
fixed pat wt dof chisq tail prob
00 0 5 2.444 0.784873 0.675 0.325
01 1 6 22.856 0.000846255 1.000 0.000
10 1 6 69.650 4.82473e-13 0.000 1.000
best pat: 00 0.784873 - -
best pat: 01 0.000846255 chi(nested): 20.412 p-value for nested model: 6.24515e-06
coeffs: 0.675 0.325
## dscore:: f_4(Base, Fit, Rbase, right2)
## genstat:: f_4(Base, Fit, right1, right2)
details: Iran_C_HajjiFiruz Russia_Ust_Ishim.DG -0.000306 -0.488565
details: Kazakhstan_Sarmatian.SG Russia_Ust_Ishim.DG 0.000229 0.320063
dscore: Russia_Ust_Ishim.DG f4: -0.000132 Z: -0.221248
details: Iran_C_HajjiFiruz Onge_1000G -0.000616 -1.068608
details: Kazakhstan_Sarmatian.SG Onge_1000G 0.000281 0.433744
dscore: Onge_1000G f4: -0.000325 Z: -0.594148
details: Iran_C_HajjiFiruz Russia_MA1_HG.SG -0.001029 -1.356903
details: Kazakhstan_Sarmatian.SG Russia_MA1_HG.SG 0.002090 2.335904
dscore: Russia_MA1_HG.SG f4: -0.000015 Z: -0.020790
details: Iran_C_HajjiFiruz Chinese_Han_1000G -0.001319 -2.756043
details: Kazakhstan_Sarmatian.SG Chinese_Han_1000G 0.001087 1.983415
dscore: Chinese_Han_1000G f4: -0.000538 Z: -1.181505
details: Iran_C_HajjiFiruz Karitiana.DG -0.001515 -2.368054
details: Kazakhstan_Sarmatian.SG Karitiana.DG 0.002453 3.290047
dscore: Karitiana.DG f4: -0.000225 Z: -0.364006
details: Iran_C_HajjiFiruz Anatolia_N 0.000570 1.157302
details: Kazakhstan_Sarmatian.SG Anatolia_N -0.002317 -4.105960
dscore: Anatolia_N f4: -0.000368 Z: -0.780912
gendstat: Mbuti.DG Russia_Ust_Ishim.DG -0.221
gendstat: Mbuti.DG Onge_1000G -0.594
gendstat: Mbuti.DG Russia_MA1_HG.SG -0.021
gendstat: Mbuti.DG Chinese_Han_1000G -1.182
gendstat: Mbuti.DG Karitiana.DG -0.364
gendstat: Mbuti.DG Anatolia_N -0.781
gendstat: Russia_Ust_Ishim.DG Onge_1000G -0.326
gendstat: Russia_Ust_Ishim.DG Russia_MA1_HG.SG 0.154
gendstat: Russia_Ust_Ishim.DG Chinese_Han_1000G -0.792
gendstat: Russia_Ust_Ishim.DG Karitiana.DG -0.147
gendstat: Russia_Ust_Ishim.DG Anatolia_N -0.431
gendstat: Onge_1000G Russia_MA1_HG.SG 0.409
gendstat: Onge_1000G Chinese_Han_1000G -0.514
gendstat: Onge_1000G Karitiana.DG 0.180
gendstat: Onge_1000G Anatolia_N -0.087
gendstat: Russia_MA1_HG.SG Chinese_Han_1000G -0.785
gendstat: Russia_MA1_HG.SG Karitiana.DG -0.294
gendstat: Russia_MA1_HG.SG Anatolia_N -0.507
gendstat: Chinese_Han_1000G Karitiana.DG 0.748
gendstat: Chinese_Han_1000G Anatolia_N 0.421
gendstat: Karitiana.DG Anatolia_N -0.266
## end of run
So what the G25 shows is irrelevant. Haji-Firuz-IA is actually 67.5% Iran-Chl + 32.5% Central Asian Sarmatian related with a p-value of 0.78 (excellent p-value, model passing p-values used in scientific papers is 0.05)