Someone might say, "But do these models work in qpAdm?".
Well, here we are presenting a model for Meena using Ror and Bhil as proxies (p-value > 0.95, very strong model).
left pops:
Meena 5
Ror 6
Bhil 8
right pops:
Ethiopia_4500BP.SG 2
Ust_Ishim.DG 1
Kostenki14 1
Mongolia_North_N 5
EHG 8
ONG.SG 6
PPNB 1
Turkey_N 28
Xinjiang_Xiaohe_BA 4
jackknife block size: 0.050
snps: 748920 indivs: 75
number of blocks for block jackknife: 714
f4info:
f4rank: 1 dof: 7 chisq: 0.998 tail: 0.994855599 dofdiff: 9 chisqdiff: -0.998 taildiff: 1
B:
scale 1.000
Ust_Ishim.DG -0.413
Kostenki14 0.432
Mongolia_North_N -1.010
EHG 1.248
ONG.SG -1.369
PPNB 0.993
Turkey_N 1.224
Xinjiang_Xiaohe_BA 0.841
A:
scale 602.595
Ror 1.153
Bhil -0.819
full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 7 chisqdiff: 0.998 taildiff: 0.994855599
B:
scale 519.184 741.162
Ust_Ishim.DG -0.375 0.474
Kostenki14 0.434 -0.420
Mongolia_North_N -0.944 1.102
EHG 1.241 -1.252
ONG.SG -1.348 1.394
PPNB 1.080 -0.863
Turkey_N 1.276 -1.142
Xinjiang_Xiaohe_BA 0.791 -0.909
A:
scale 1.414 1.414
Ror 1.414 0.000
Bhil 0.000 1.414
best coefficients: 0.415 0.585
Jackknife mean: 0.415149944 0.584850056
std. errors: 0.025 0.025
error covariance (* 1,000,000)
633 -633
-633 633
summ: Meena 2 0.994856 0.415 0.585 633 -633 633
fixed pat wt dof chisq tail prob
00 0 7 0.998 0.994856 0.415 0.585
01 1 8 360.891 0 1.000 0.000
10 1 8 206.778 0 0.000 1.000
best pat: 00 0.994856 - -
best pat: 10 2.3799e-40 chi(nested): 205.780 p-value for nested model: 1.14435e-46
coeffs: 0.415 0.585
## dscore:: f_4(Base, Fit, Rbase, right2)
## genstat:: f_4(Base, Fit, right1, right2)
details: Ror Ust_Ishim.DG -0.000721 -1.645461
details: Bhil Ust_Ishim.DG 0.000640 1.544434
dscore: Ust_Ishim.DG f4: 0.000075 Z: 0.202023
details: Ror Kostenki14 0.000835 1.854928
details: Bhil Kostenki14 -0.000567 -1.386642
dscore: Kostenki14 f4: 0.000015 Z: 0.041166
details: Ror Mongolia_North_N -0.001819 -5.104447
details: Bhil Mongolia_North_N 0.001486 4.536255
dscore: Mongolia_North_N f4: 0.000114 Z: 0.390186
details: Ror EHG 0.002391 7.020298
details: Bhil EHG -0.001689 -5.039449
dscore: EHG f4: 0.000004 Z: 0.015385
details: Ror
ONG.SG -0.002597 -7.219302
details: Bhil
ONG.SG 0.001881 5.487526
dscore:
ONG.SG f4: 0.000022 Z: 0.072929
details: Ror PPNB 0.002080 4.899661
details: Bhil PPNB -0.001164 -2.802844
dscore: PPNB f4: 0.000183 Z: 0.498940
details: Ror Turkey_N 0.002458 7.792816
details: Bhil Turkey_N -0.001540 -5.035728
dscore: Turkey_N f4: 0.000120 Z: 0.444917
details: Ror Xinjiang_Xiaohe_BA 0.001523 3.676915
details: Bhil Xinjiang_Xiaohe_BA -0.001227 -3.099144
dscore: Xinjiang_Xiaohe_BA f4: -0.000085 Z: -0.242810
gendstat: Ethiopia_4500BP.SG Ust_Ishim.DG 0.202
gendstat: Ethiopia_4500BP.SG Kostenki14 0.041
gendstat: Ethiopia_4500BP.SG Mongolia_North_N 0.390
gendstat: Ethiopia_4500BP.SG EHG 0.015
gendstat: Ethiopia_4500BP.SG
ONG.SG 0.073
gendstat: Ethiopia_4500BP.SG PPNB 0.499
gendstat: Ethiopia_4500BP.SG Turkey_N 0.445
gendstat: Ethiopia_4500BP.SG Xinjiang_Xiaohe_BA -0.243
gendstat: Ust_Ishim.DG Kostenki14 -0.145
gendstat: Ust_Ishim.DG Mongolia_North_N 0.116
gendstat: Ust_Ishim.DG EHG -0.207
gendstat: Ust_Ishim.DG
ONG.SG -0.154
gendstat: Ust_Ishim.DG PPNB 0.274
gendstat: Ust_Ishim.DG Turkey_N 0.138
gendstat: Ust_Ishim.DG Xinjiang_Xiaohe_BA -0.412
gendstat: Kostenki14 Mongolia_North_N 0.280
gendstat: Kostenki14 EHG -0.033
gendstat: Kostenki14
ONG.SG 0.019
gendstat: Kostenki14 PPNB 0.420
gendstat: Kostenki14 Turkey_N 0.343
gendstat: Kostenki14 Xinjiang_Xiaohe_BA -0.262
gendstat: Mongolia_North_N EHG -0.439
gendstat: Mongolia_North_N
ONG.SG -0.378
gendstat: Mongolia_North_N PPNB 0.201
gendstat: Mongolia_North_N Turkey_N 0.023
gendstat: Mongolia_North_N Xinjiang_Xiaohe_BA -0.680
gendstat: EHG
ONG.SG 0.063
gendstat: EHG PPNB 0.554
gendstat: EHG Turkey_N 0.599
gendstat: EHG Xinjiang_Xiaohe_BA -0.362
gendstat:
ONG.SG PPNB 0.453
gendstat:
ONG.SG Turkey_N 0.375
gendstat:
ONG.SG Xinjiang_Xiaohe_BA -0.328
gendstat: PPNB Turkey_N -0.242
gendstat: PPNB Xinjiang_Xiaohe_BA -0.716
gendstat: Turkey_N Xinjiang_Xiaohe_BA -0.709
worst Z-score with right hand mix
f4(Target, Fit, Base, mix of Right pops; Z: 1.000 sum: 1.000
Ust_Ishim.DG 0.099
Kostenki14 -0.517
Mongolia_North_N 2.086
EHG -1.190
ONG.SG -1.061
PPNB 0.938
Turkey_N 2.428
Xinjiang_Xiaohe_BA -1.783