PDA

View Full Version : qpAdm: 50% Greek_Northern_Thrace (Lukasz') + 50% Pontic Greek Trebizond/Gumushane.



kingmob
02-07-2022, 06:49 PM
1) Turkey_N.SG 54.4% + Iran_GanjDareh_N 17.7% + Russia_Samara_EBA_Yamnaya 27.9%.



qpAdm: parameter file: parqpadm.txt
### THE INPUT PARAMETERS
##PARAMETER NAME: VALUE
fstatsname: fstatsa.txt
popleft: left.txt
popright: right.txt
details: YES
## qpAdm version: 1520
seed: 1182635581
*** recommended that inbreed be explicitly set ***
inbreed set NO

left pops:
eupator
Turkey_N.SG
Iran_GanjDareh_N
Russia_Samara_EBA_Yamnaya

right pops:
Mbuti.DG
Russia_Ust_Ishim.DG
China_Tianyuan
Goyet_Neanderthal.SG
Russia_Sunghir3.SG
Russia_Kostenki14.SG
Spain_ElMiron
Russia_MA1_HG.SG
Georgia_Satsurblia.SG
Russia_DevilsCave_N.SG
Papuan.DG
Switzerland_Bichon.SG
Israel_Natufian_published
Morocco_Iberomaurusian
Russia_Steppe_Eneolithic
ONG.SG

codimension 1
f4info:
f4rank: 2 dof: 13 chisq: 17.254 tail: 0.187941462 dofdiff: 15 chisqdiff: -17.254 taildiff: 1
B:
scale 1.000 1.000
Russia_Ust_Ishim.DG 0.469 -0.020
China_Tianyuan 0.527 -0.681
Goyet_Neanderthal.SG 0.135 0.117
Russia_Sunghir3.SG 1.082 0.663
Russia_Kostenki14.SG 0.868 0.288
Spain_ElMiron 1.476 0.751
Russia_MA1_HG.SG 1.756 -1.449
Georgia_Satsurblia.SG 0.406 -0.906
Russia_DevilsCave_N.SG 0.616 -0.451
Papuan.DG 0.421 -0.483
Switzerland_Bichon.SG 1.858 0.821
Israel_Natufian_published 0.482 2.372
Morocco_Iberomaurusian 0.484 1.176
Russia_Steppe_Eneolithic 1.545 -1.514
ONG.SG 0.521 -0.320
A:
scale 512.183 745.798
Turkey_N.SG -0.169 0.866
Iran_GanjDareh_N -1.286 -1.159
Russia_Samara_EBA_Yamnaya 1.148 -0.952


full rank
f4info:
f4rank: 3 dof: 0 chisq: 0.000 tail: 1 dofdiff: 13 chisqdiff: 17.254 taildiff: 0.187941462
B:
scale 836.495 355.469 350.294
Russia_Ust_Ishim.DG 0.253 -0.250 0.579
China_Tianyuan -0.377 0.052 0.911
Goyet_Neanderthal.SG -0.087 -0.257 -0.022
Russia_Sunghir3.SG 0.489 -1.260 0.632
Russia_Kostenki14.SG 0.646 -0.684 0.826
Spain_ElMiron 0.689 -1.587 1.004
Russia_MA1_HG.SG -1.546 -0.661 2.145
Georgia_Satsurblia.SG -1.167 0.067 0.683
Russia_DevilsCave_N.SG -0.142 -0.119 0.892
Papuan.DG -0.813 -0.212 0.478
Switzerland_Bichon.SG -0.110 -2.263 0.985
Israel_Natufian_published 2.116 -1.770 -0.680
Morocco_Iberomaurusian 0.955 -1.109 -0.153
Russia_Steppe_Eneolithic -1.952 -0.551 1.896
ONG.SG 0.092 -0.081 0.777
A:
scale 1.732 1.732 1.732
Turkey_N.SG 1.732 0.000 0.000
Iran_GanjDareh_N 0.000 1.732 0.000
Russia_Samara_EBA_Yamnaya 0.000 0.000 1.732


best coefficients: 0.544 0.177 0.279
totmean: 0.544 0.177 0.279
boot mean: 0.544 0.177 0.279
std. errors: 0.042 0.050 0.044

error covariance (* 1,000,000)
1725 -1164 -561
-1164 2506 -1342
-561 -1342 1902


summ: eupator 3 0.187941 0.544 0.177 0.279 1725 -1164 -561 2506 -1342 ...
1902

fixed pat wt dof chisq tail prob
000 0 13 17.254 0.187941 0.544 0.177 0.279
001 1 14 57.080 3.78963e-07 0.623 0.377 0.000
010 1 14 29.450 0.009077 0.631 0.000 0.369
100 1 14 178.311 1.4246e-30 0.000 0.580 0.420
011 2 15 124.155 2.9456e-19 1.000 0.000 0.000
101 2 15 246.997 0 0.000 1.000 0.000
110 2 15 314.468 0 0.000 0.000 1.000
best pat: 000 0.187941 - -
best pat: 010 0.009077 chi(nested): 12.196 p-value for nested model: 0.000478859
best pat: 011 2.9456e-19 chi(nested): 94.705 p-value for nested model: 2.21028e-22

coeffs: 0.544 0.177 0.279

## dscore:: f_4(Base, Fit, Rbase, right2)
## genstat:: f_4(Base, Fit, right1, right2)

details: Turkey_N.SG Russia_Ust_Ishim.DG 0.000302 0.618214
details: Iran_GanjDareh_N Russia_Ust_Ishim.DG -0.000704 -1.377963
details: Russia_Samara_EBA_Yamnaya Russia_Ust_Ishim.DG 0.001652 3.410899
dscore: Russia_Ust_Ishim.DG f4: 0.000500 Z: 1.113515

details: Turkey_N.SG China_Tianyuan -0.000451 -0.877752
details: Iran_GanjDareh_N China_Tianyuan 0.000147 0.274109
details: Russia_Samara_EBA_Yamnaya China_Tianyuan 0.002601 5.048198
dscore: China_Tianyuan f4: 0.000506 Z: 1.069118

details: Turkey_N.SG Goyet_Neanderthal.SG -0.000104 -0.307809
details: Iran_GanjDareh_N Goyet_Neanderthal.SG -0.000723 -2.045924
details: Russia_Samara_EBA_Yamnaya Goyet_Neanderthal.SG -0.000063 -0.192775
dscore: Goyet_Neanderthal.SG f4: -0.000202 Z: -0.675928

details: Turkey_N.SG Russia_Sunghir3.SG 0.000584 1.038869
details: Iran_GanjDareh_N Russia_Sunghir3.SG -0.003543 -6.018781
details: Russia_Samara_EBA_Yamnaya Russia_Sunghir3.SG 0.001803 3.353448
dscore: Russia_Sunghir3.SG f4: 0.000192 Z: 0.376621

details: Turkey_N.SG Russia_Kostenki14.SG 0.000772 1.421486
details: Iran_GanjDareh_N Russia_Kostenki14.SG -0.001923 -3.377255
details: Russia_Samara_EBA_Yamnaya Russia_Kostenki14.SG 0.002358 4.523081
dscore: Russia_Kostenki14.SG f4: 0.000736 Z: 1.482232

details: Turkey_N.SG Spain_ElMiron 0.000824 1.441107
details: Iran_GanjDareh_N Spain_ElMiron -0.004464 -7.371176
details: Russia_Samara_EBA_Yamnaya Spain_ElMiron 0.002865 5.169620
dscore: Spain_ElMiron f4: 0.000455 Z: 0.873834

details: Turkey_N.SG Russia_MA1_HG.SG -0.001849 -3.200808
details: Iran_GanjDareh_N Russia_MA1_HG.SG -0.001860 -3.124025
details: Russia_Samara_EBA_Yamnaya Russia_MA1_HG.SG 0.006124 11.118424
dscore: Russia_MA1_HG.SG f4: 0.000372 Z: 0.717092

details: Turkey_N.SG Georgia_Satsurblia.SG -0.001395 -2.444830
details: Iran_GanjDareh_N Georgia_Satsurblia.SG 0.000189 0.327158
details: Russia_Samara_EBA_Yamnaya Georgia_Satsurblia.SG 0.001950 3.676764
dscore: Georgia_Satsurblia.SG f4: -0.000182 Z: -0.355988

details: Turkey_N.SG Russia_DevilsCave_N.SG -0.000169 -0.400825
details: Iran_GanjDareh_N Russia_DevilsCave_N.SG -0.000334 -0.705857
details: Russia_Samara_EBA_Yamnaya Russia_DevilsCave_N.SG 0.002547 6.035887
dscore: Russia_DevilsCave_N.SG f4: 0.000558 Z: 1.431705

details: Turkey_N.SG Papuan.DG -0.000972 -2.348136
details: Iran_GanjDareh_N Papuan.DG -0.000595 -1.368727
details: Russia_Samara_EBA_Yamnaya Papuan.DG 0.001365 3.308717
dscore: Papuan.DG f4: -0.000253 Z: -0.663602

details: Turkey_N.SG Switzerland_Bichon.SG -0.000131 -0.235759
details: Iran_GanjDareh_N Switzerland_Bichon.SG -0.006368 -11.057025
details: Russia_Samara_EBA_Yamnaya Switzerland_Bichon.SG 0.002812 5.033356
dscore: Switzerland_Bichon.SG f4: -0.000417 Z: -0.815145

details: Turkey_N.SG Israel_Natufian_published 0.002530 4.137315
details: Iran_GanjDareh_N Israel_Natufian_published -0.004979 -8.172541
details: Russia_Samara_EBA_Yamnaya Israel_Natufian_published -0.001942 -3.430105
dscore: Israel_Natufian_published f4: -0.000048 Z: -0.089802

details: Turkey_N.SG Morocco_Iberomaurusian 0.001141 2.903000
details: Iran_GanjDareh_N Morocco_Iberomaurusian -0.003121 -7.724254
details: Russia_Samara_EBA_Yamnaya Morocco_Iberomaurusian -0.000437 -1.175811
dscore: Morocco_Iberomaurusian f4: -0.000055 Z: -0.155598

details: Turkey_N.SG Russia_Steppe_Eneolithic -0.002334 -5.431361
details: Iran_GanjDareh_N Russia_Steppe_Eneolithic -0.001550 -3.472956
details: Russia_Samara_EBA_Yamnaya Russia_Steppe_Eneolithic 0.005412 13.105774
dscore: Russia_Steppe_Eneolithic f4: -0.000036 Z: -0.092055

details: Turkey_N.SG ONG.SG 0.000110 0.256772
details: Iran_GanjDareh_N ONG.SG -0.000227 -0.505430
details: Russia_Samara_EBA_Yamnaya ONG.SG 0.002217 5.385728
dscore: ONG.SG f4: 0.000638 Z: 1.631663

gendstat: Mbuti.DG Russia_Ust_Ishim.DG 1.114
gendstat: Mbuti.DG China_Tianyuan 1.069
gendstat: Mbuti.DG Goyet_Neanderthal.SG -0.676
gendstat: Mbuti.DG Russia_Sunghir3.SG 0.377
gendstat: Mbuti.DG Russia_Kostenki14.SG 1.482
gendstat: Mbuti.DG Spain_ElMiron 0.874
gendstat: Mbuti.DG Russia_MA1_HG.SG 0.717
gendstat: Mbuti.DG Georgia_Satsurblia.SG -0.356
gendstat: Mbuti.DG Russia_DevilsCave_N.SG 1.432
gendstat: Mbuti.DG Papuan.DG -0.664
gendstat: Mbuti.DG Switzerland_Bichon.SG -0.815
gendstat: Mbuti.DG Israel_Natufian_published -0.090
gendstat: Mbuti.DG Morocco_Iberomaurusian -0.156
gendstat: Mbuti.DG Russia_Steppe_Eneolithic -0.092
gendstat: Mbuti.DG ONG.SG 1.632
gendstat: Russia_Ust_Ishim.DG China_Tianyuan 0.011
gendstat: Russia_Ust_Ishim.DG Goyet_Neanderthal.SG -1.348
gendstat: Russia_Ust_Ishim.DG Russia_Sunghir3.SG -0.549
gendstat: Russia_Ust_Ishim.DG Russia_Kostenki14.SG 0.419
gendstat: Russia_Ust_Ishim.DG Spain_ElMiron -0.075
gendstat: Russia_Ust_Ishim.DG Russia_MA1_HG.SG -0.220
gendstat: Russia_Ust_Ishim.DG Georgia_Satsurblia.SG -1.163
gendstat: Russia_Ust_Ishim.DG Russia_DevilsCave_N.SG 0.128
gendstat: Russia_Ust_Ishim.DG Papuan.DG -1.685
gendstat: Russia_Ust_Ishim.DG Switzerland_Bichon.SG -1.577
gendstat: Russia_Ust_Ishim.DG Israel_Natufian_published -0.879
gendstat: Russia_Ust_Ishim.DG Morocco_Iberomaurusian -1.104
gendstat: Russia_Ust_Ishim.DG Russia_Steppe_Eneolithic -1.132
gendstat: Russia_Ust_Ishim.DG ONG.SG 0.308
gendstat: China_Tianyuan Goyet_Neanderthal.SG -1.269
gendstat: China_Tianyuan Russia_Sunghir3.SG -0.558
gendstat: China_Tianyuan Russia_Kostenki14.SG 0.394
gendstat: China_Tianyuan Spain_ElMiron -0.084
gendstat: China_Tianyuan Russia_MA1_HG.SG -0.230
gendstat: China_Tianyuan Georgia_Satsurblia.SG -1.094
gendstat: China_Tianyuan Russia_DevilsCave_N.SG 0.115
gendstat: China_Tianyuan Papuan.DG -1.682
gendstat: China_Tianyuan Switzerland_Bichon.SG -1.545
gendstat: China_Tianyuan Israel_Natufian_published -0.833
gendstat: China_Tianyuan Morocco_Iberomaurusian -1.103
gendstat: China_Tianyuan Russia_Steppe_Eneolithic -1.071
gendstat: China_Tianyuan ONG.SG 0.293
gendstat: Goyet_Neanderthal.SG Russia_Sunghir3.SG 0.710
gendstat: Goyet_Neanderthal.SG Russia_Kostenki14.SG 1.655
gendstat: Goyet_Neanderthal.SG Spain_ElMiron 1.099
gendstat: Goyet_Neanderthal.SG Russia_MA1_HG.SG 0.981
gendstat: Goyet_Neanderthal.SG Georgia_Satsurblia.SG 0.035
gendstat: Goyet_Neanderthal.SG Russia_DevilsCave_N.SG 1.612
gendstat: Goyet_Neanderthal.SG Papuan.DG -0.108
gendstat: Goyet_Neanderthal.SG Switzerland_Bichon.SG -0.373
gendstat: Goyet_Neanderthal.SG Israel_Natufian_published 0.249
gendstat: Goyet_Neanderthal.SG Morocco_Iberomaurusian 0.323
gendstat: Goyet_Neanderthal.SG Russia_Steppe_Eneolithic 0.352
gendstat: Goyet_Neanderthal.SG ONG.SG 1.751
gendstat: Russia_Sunghir3.SG Russia_Kostenki14.SG 0.980
gendstat: Russia_Sunghir3.SG Spain_ElMiron 0.437
gendstat: Russia_Sunghir3.SG Russia_MA1_HG.SG 0.324
gendstat: Russia_Sunghir3.SG Georgia_Satsurblia.SG -0.578
gendstat: Russia_Sunghir3.SG Russia_DevilsCave_N.SG 0.714
gendstat: Russia_Sunghir3.SG Papuan.DG -0.878
gendstat: Russia_Sunghir3.SG Switzerland_Bichon.SG -1.077
gendstat: Russia_Sunghir3.SG Israel_Natufian_published -0.362
gendstat: Russia_Sunghir3.SG Morocco_Iberomaurusian -0.467
gendstat: Russia_Sunghir3.SG Russia_Steppe_Eneolithic -0.442
gendstat: Russia_Sunghir3.SG ONG.SG 0.872
gendstat: Russia_Kostenki14.SG Spain_ElMiron -0.475
gendstat: Russia_Kostenki14.SG Russia_MA1_HG.SG -0.638
gendstat: Russia_Kostenki14.SG Georgia_Satsurblia.SG -1.503
gendstat: Russia_Kostenki14.SG Russia_DevilsCave_N.SG -0.330
gendstat: Russia_Kostenki14.SG Papuan.DG -1.923
gendstat: Russia_Kostenki14.SG Switzerland_Bichon.SG -2.033
gendstat: Russia_Kostenki14.SG Israel_Natufian_published -1.173
gendstat: Russia_Kostenki14.SG Morocco_Iberomaurusian -1.526
gendstat: Russia_Kostenki14.SG Russia_Steppe_Eneolithic -1.556
gendstat: Russia_Kostenki14.SG ONG.SG -0.188
gendstat: Spain_ElMiron Russia_MA1_HG.SG -0.139
gendstat: Spain_ElMiron Georgia_Satsurblia.SG -1.007
gendstat: Spain_ElMiron Russia_DevilsCave_N.SG 0.189
gendstat: Spain_ElMiron Papuan.DG -1.343
gendstat: Spain_ElMiron Switzerland_Bichon.SG -1.603
gendstat: Spain_ElMiron Israel_Natufian_published -0.753
gendstat: Spain_ElMiron Morocco_Iberomaurusian -0.967
gendstat: Spain_ElMiron Russia_Steppe_Eneolithic -0.966
gendstat: Spain_ElMiron ONG.SG 0.335
gendstat: Russia_MA1_HG.SG Georgia_Satsurblia.SG -0.881
gendstat: Russia_MA1_HG.SG Russia_DevilsCave_N.SG 0.352
gendstat: Russia_MA1_HG.SG Papuan.DG -1.227
gendstat: Russia_MA1_HG.SG Switzerland_Bichon.SG -1.337
gendstat: Russia_MA1_HG.SG Israel_Natufian_published -0.618
gendstat: Russia_MA1_HG.SG Morocco_Iberomaurusian -0.808
gendstat: Russia_MA1_HG.SG Russia_Steppe_Eneolithic -0.840
gendstat: Russia_MA1_HG.SG ONG.SG 0.528
gendstat: Georgia_Satsurblia.SG Russia_DevilsCave_N.SG 1.375
gendstat: Georgia_Satsurblia.SG Papuan.DG -0.131
gendstat: Georgia_Satsurblia.SG Switzerland_Bichon.SG -0.391
gendstat: Georgia_Satsurblia.SG Israel_Natufian_published 0.203
gendstat: Georgia_Satsurblia.SG Morocco_Iberomaurusian 0.227
gendstat: Georgia_Satsurblia.SG Russia_Steppe_Eneolithic 0.301
gendstat: Georgia_Satsurblia.SG ONG.SG 1.581
gendstat: Russia_DevilsCave_N.SG Papuan.DG -2.172
gendstat: Russia_DevilsCave_N.SG Switzerland_Bichon.SG -1.820
gendstat: Russia_DevilsCave_N.SG Israel_Natufian_published -1.016
gendstat: Russia_DevilsCave_N.SG Morocco_Iberomaurusian -1.434
gendstat: Russia_DevilsCave_N.SG Russia_Steppe_Eneolithic -1.467
gendstat: Russia_DevilsCave_N.SG ONG.SG 0.217
gendstat: Papuan.DG Switzerland_Bichon.SG -0.315
gendstat: Papuan.DG Israel_Natufian_published 0.336
gendstat: Papuan.DG Morocco_Iberomaurusian 0.477
gendstat: Papuan.DG Russia_Steppe_Eneolithic 0.540
gendstat: Papuan.DG ONG.SG 2.614
gendstat: Switzerland_Bichon.SG Israel_Natufian_published 0.557
gendstat: Switzerland_Bichon.SG Morocco_Iberomaurusian 0.698
gendstat: Switzerland_Bichon.SG Russia_Steppe_Eneolithic 0.792
gendstat: Switzerland_Bichon.SG ONG.SG 1.963
gendstat: Israel_Natufian_published Morocco_Iberomaurusian -0.011
gendstat: Israel_Natufian_published Russia_Steppe_Eneolithic 0.023
gendstat: Israel_Natufian_published ONG.SG 1.136
gendstat: Morocco_Iberomaurusian Russia_Steppe_Eneolithic 0.045
gendstat: Morocco_Iberomaurusian ONG.SG 1.593
gendstat: Russia_Steppe_Eneolithic ONG.SG 1.682

##end of qpAdm: 9.267 seconds cpu 1.270 Mbytes in use


2) Iran_GanjDareh_N 26.2% + Greece_Peloponnese_N_all 53% + Russia_Samara_EBA_Yamnaya 20.9%.




qpAdm: parameter file: parqpadm.txt
### THE INPUT PARAMETERS
##PARAMETER NAME: VALUE
fstatsname: fstatsa.txt
popleft: left.txt
popright: right.txt
details: YES
## qpAdm version: 1520
seed: 1632377610
*** recommended that inbreed be explicitly set ***
inbreed set NO

left pops:
eupator
Iran_GanjDareh_N
Greece_Peloponnese_N_all
Russia_Samara_EBA_Yamnaya

right pops:
Mbuti.DG
Russia_Ust_Ishim.DG
China_Tianyuan
Goyet_Neanderthal.SG
Turkey_N_published
Russia_MA1_HG.SG
Georgia_Satsurblia.SG
Georgia_Kotias.SG
Papuan.DG
Russia_DevilsCave_N.SG
Israel_Natufian_published
Morocco_Iberomaurusian
Russia_Steppe_Eneolithic
ONG.SG

codimension 1
f4info:
f4rank: 2 dof: 11 chisq: 9.967 tail: 0.533388563 dofdiff: 13 chisqdiff: -9.967 taildiff: 1
B:
scale 1.000 1.000
Russia_Ust_Ishim.DG -0.089 0.615
China_Tianyuan -0.342 0.620
Goyet_Neanderthal.SG 0.105 0.198
Turkey_N_published 2.486 1.226
Russia_MA1_HG.SG -1.176 2.052
Georgia_Satsurblia.SG -0.472 0.407
Georgia_Kotias.SG -0.524 0.591
Papuan.DG -0.315 0.487
Russia_DevilsCave_N.SG -0.294 0.754
Israel_Natufian_published 1.652 0.970
Morocco_Iberomaurusian 0.901 0.802
Russia_Steppe_Eneolithic -0.988 1.787
ONG.SG -0.317 0.623
A:
scale 496.205 650.643
Iran_GanjDareh_N -0.956 -1.290
Greece_Peloponnese_N_all 0.913 0.189
Russia_Samara_EBA_Yamnaya -1.119 1.140


full rank
f4info:
f4rank: 3 dof: 0 chisq: 0.000 tail: 1 dofdiff: 11 chisqdiff: 9.967 taildiff: 0.533388563
B:
scale 375.218 517.856 336.564
Russia_Ust_Ishim.DG -0.287 0.178 0.539
China_Tianyuan 0.044 0.281 0.867
Goyet_Neanderthal.SG -0.275 0.057 -0.024
Turkey_N_published -2.651 2.605 -1.107
Russia_MA1_HG.SG -0.695 -0.881 2.062
Georgia_Satsurblia.SG 0.155 -0.206 0.696
Georgia_Kotias.SG 0.045 -0.298 0.817
Papuan.DG -0.215 -0.371 0.464
Russia_DevilsCave_N.SG -0.125 0.232 0.864
Israel_Natufian_published -1.866 1.850 -0.639
Morocco_Iberomaurusian -1.194 1.102 -0.163
Russia_Steppe_Eneolithic -0.569 -0.603 1.841
ONG.SG -0.095 0.020 0.736
A:
scale 1.732 1.732 1.732
Iran_GanjDareh_N 1.732 0.000 0.000
Greece_Peloponnese_N_all 0.000 1.732 0.000
Russia_Samara_EBA_Yamnaya 0.000 0.000 1.732


best coefficients: 0.262 0.530 0.209
totmean: 0.262 0.530 0.209
boot mean: 0.262 0.530 0.208
std. errors: 0.050 0.023 0.051

error covariance (* 1,000,000)
2546 -260 -2287
-260 535 -275
-2287 -275 2562


summ: eupator 3 0.533389 0.262 0.530 0.208 2546 -260 -2287 535 -275 ...
2562

fixed pat wt dof chisq tail prob
000 0 11 9.967 0.533389 0.262 0.530 0.209
001 1 12 25.946 0.0109239 0.447 0.553 0.000
010 1 12 536.144 0 1.724 0.000 -0.724 infeasible
100 1 12 36.058 0.000317153 0.000 0.565 0.435
011 2 13 531.843 0 1.000 0.000 0.000
101 2 13 261.162 0 0.000 1.000 0.000
110 2 13 677.525 0 0.000 0.000 1.000
best pat: 000 0.533389 - -
best pat: 100 0.000317153 chi(nested): 26.091 p-value for nested model: 3.25702e-07
best pat: 101 3.06274e-48 chi(nested): 225.104 p-value for nested model: 6.96805e-51

coeffs: 0.262 0.530 0.209

## dscore:: f_4(Base, Fit, Rbase, right2)
## genstat:: f_4(Base, Fit, right1, right2)

details: Iran_GanjDareh_N Russia_Ust_Ishim.DG -0.000765 -1.484520
details: Greece_Peloponnese_N_all Russia_Ust_Ishim.DG 0.000344 0.535334
details: Russia_Samara_EBA_Yamnaya Russia_Ust_Ishim.DG 0.001603 3.291657
dscore: Russia_Ust_Ishim.DG f4: 0.000316 Z: 0.618349

details: Iran_GanjDareh_N China_Tianyuan 0.000117 0.218656
details: Greece_Peloponnese_N_all China_Tianyuan 0.000543 0.832302
details: Russia_Samara_EBA_Yamnaya China_Tianyuan 0.002576 5.017138
dscore: China_Tianyuan f4: 0.000856 Z: 1.660690

details: Iran_GanjDareh_N Goyet_Neanderthal.SG -0.000732 -2.060796
details: Greece_Peloponnese_N_all Goyet_Neanderthal.SG 0.000111 0.256292
details: Russia_Samara_EBA_Yamnaya Goyet_Neanderthal.SG -0.000072 -0.217575
dscore: Goyet_Neanderthal.SG f4: -0.000148 Z: -0.444906

details: Iran_GanjDareh_N Turkey_N_published -0.007066 -17.580192
details: Greece_Peloponnese_N_all Turkey_N_published 0.005030 10.786634
details: Russia_Samara_EBA_Yamnaya Turkey_N_published -0.003290 -8.829826
dscore: Turkey_N_published f4: 0.000129 Z: 0.355120

details: Iran_GanjDareh_N Russia_MA1_HG.SG -0.001853 -3.104973
details: Greece_Peloponnese_N_all Russia_MA1_HG.SG -0.001701 -2.452482
details: Russia_Samara_EBA_Yamnaya Russia_MA1_HG.SG 0.006126 11.090729
dscore: Russia_MA1_HG.SG f4: -0.000108 Z: -0.195712

details: Iran_GanjDareh_N Georgia_Satsurblia.SG 0.000414 0.717062
details: Greece_Peloponnese_N_all Georgia_Satsurblia.SG -0.000398 -0.591303
details: Russia_Samara_EBA_Yamnaya Georgia_Satsurblia.SG 0.002069 3.868145
dscore: Georgia_Satsurblia.SG f4: 0.000329 Z: 0.615015

details: Iran_GanjDareh_N Georgia_Kotias.SG 0.000120 0.221597
details: Greece_Peloponnese_N_all Georgia_Kotias.SG -0.000575 -0.861967
details: Russia_Samara_EBA_Yamnaya Georgia_Kotias.SG 0.002427 4.775409
dscore: Georgia_Kotias.SG f4: 0.000233 Z: 0.441048

details: Iran_GanjDareh_N Papuan.DG -0.000573 -1.317589
details: Greece_Peloponnese_N_all Papuan.DG -0.000716 -1.397809
details: Russia_Samara_EBA_Yamnaya Papuan.DG 0.001379 3.361719
dscore: Papuan.DG f4: -0.000242 Z: -0.587946

details: Iran_GanjDareh_N Russia_DevilsCave_N.SG -0.000334 -0.700963
details: Greece_Peloponnese_N_all Russia_DevilsCave_N.SG 0.000447 0.843878
details: Russia_Samara_EBA_Yamnaya Russia_DevilsCave_N.SG 0.002567 6.061835
dscore: Russia_DevilsCave_N.SG f4: 0.000685 Z: 1.593348

details: Iran_GanjDareh_N Israel_Natufian_published -0.004972 -8.218806
details: Greece_Peloponnese_N_all Israel_Natufian_published 0.003573 4.809877
details: Russia_Samara_EBA_Yamnaya Israel_Natufian_published -0.001898 -3.378756
dscore: Israel_Natufian_published f4: 0.000196 Z: 0.339543

details: Iran_GanjDareh_N Morocco_Iberomaurusian -0.003181 -7.884584
details: Greece_Peloponnese_N_all Morocco_Iberomaurusian 0.002127 4.518811
details: Russia_Samara_EBA_Yamnaya Morocco_Iberomaurusian -0.000485 -1.306593
dscore: Morocco_Iberomaurusian f4: 0.000193 Z: 0.507130

details: Iran_GanjDareh_N Russia_Steppe_Eneolithic -0.001517 -3.400652
details: Greece_Peloponnese_N_all Russia_Steppe_Eneolithic -0.001165 -2.274589
details: Russia_Samara_EBA_Yamnaya Russia_Steppe_Eneolithic 0.005469 13.326350
dscore: Russia_Steppe_Eneolithic f4: 0.000126 Z: 0.309758

details: Iran_GanjDareh_N ONG.SG -0.000253 -0.563236
details: Greece_Peloponnese_N_all ONG.SG 0.000039 0.076001
details: Russia_Samara_EBA_Yamnaya ONG.SG 0.002188 5.296813
dscore: ONG.SG f4: 0.000411 Z: 0.986145

gendstat: Mbuti.DG Russia_Ust_Ishim.DG 0.618
gendstat: Mbuti.DG China_Tianyuan 1.661
gendstat: Mbuti.DG Goyet_Neanderthal.SG -0.445
gendstat: Mbuti.DG Turkey_N_published 0.355
gendstat: Mbuti.DG Russia_MA1_HG.SG -0.196
gendstat: Mbuti.DG Georgia_Satsurblia.SG 0.615
gendstat: Mbuti.DG Georgia_Kotias.SG 0.441
gendstat: Mbuti.DG Papuan.DG -0.588
gendstat: Mbuti.DG Russia_DevilsCave_N.SG 1.593
gendstat: Mbuti.DG Israel_Natufian_published 0.340
gendstat: Mbuti.DG Morocco_Iberomaurusian 0.507
gendstat: Mbuti.DG Russia_Steppe_Eneolithic 0.310
gendstat: Mbuti.DG ONG.SG 0.986
gendstat: Russia_Ust_Ishim.DG China_Tianyuan 0.896
gendstat: Russia_Ust_Ishim.DG Goyet_Neanderthal.SG -0.790
gendstat: Russia_Ust_Ishim.DG Turkey_N_published -0.360
gendstat: Russia_Ust_Ishim.DG Russia_MA1_HG.SG -0.641
gendstat: Russia_Ust_Ishim.DG Georgia_Satsurblia.SG 0.020
gendstat: Russia_Ust_Ishim.DG Georgia_Kotias.SG -0.128
gendstat: Russia_Ust_Ishim.DG Papuan.DG -1.134
gendstat: Russia_Ust_Ishim.DG Russia_DevilsCave_N.SG 0.703
gendstat: Russia_Ust_Ishim.DG Israel_Natufian_published -0.174
gendstat: Russia_Ust_Ishim.DG Morocco_Iberomaurusian -0.219
gendstat: Russia_Ust_Ishim.DG Russia_Steppe_Eneolithic -0.361
gendstat: Russia_Ust_Ishim.DG ONG.SG 0.186
gendstat: China_Tianyuan Goyet_Neanderthal.SG -1.704
gendstat: China_Tianyuan Turkey_N_published -1.385
gendstat: China_Tianyuan Russia_MA1_HG.SG -1.492
gendstat: China_Tianyuan Georgia_Satsurblia.SG -0.760
gendstat: China_Tianyuan Georgia_Kotias.SG -0.943
gendstat: China_Tianyuan Papuan.DG -2.263
gendstat: China_Tianyuan Russia_DevilsCave_N.SG -0.332
gendstat: China_Tianyuan Israel_Natufian_published -0.915
gendstat: China_Tianyuan Morocco_Iberomaurusian -1.181
gendstat: China_Tianyuan Russia_Steppe_Eneolithic -1.323
gendstat: China_Tianyuan ONG.SG -0.872
gendstat: Goyet_Neanderthal.SG Turkey_N_published 0.598
gendstat: Goyet_Neanderthal.SG Russia_MA1_HG.SG 0.063
gendstat: Goyet_Neanderthal.SG Georgia_Satsurblia.SG 0.789
gendstat: Goyet_Neanderthal.SG Georgia_Kotias.SG 0.628
gendstat: Goyet_Neanderthal.SG Papuan.DG -0.185
gendstat: Goyet_Neanderthal.SG Russia_DevilsCave_N.SG 1.631
gendstat: Goyet_Neanderthal.SG Israel_Natufian_published 0.518
gendstat: Goyet_Neanderthal.SG Morocco_Iberomaurusian 0.691
gendstat: Goyet_Neanderthal.SG Russia_Steppe_Eneolithic 0.548
gendstat: Goyet_Neanderthal.SG ONG.SG 1.083
gendstat: Turkey_N_published Russia_MA1_HG.SG -0.432
gendstat: Turkey_N_published Georgia_Satsurblia.SG 0.389
gendstat: Turkey_N_published Georgia_Kotias.SG 0.209
gendstat: Turkey_N_published Papuan.DG -0.885
gendstat: Turkey_N_published Russia_DevilsCave_N.SG 1.276
gendstat: Turkey_N_published Israel_Natufian_published 0.116
gendstat: Turkey_N_published Morocco_Iberomaurusian 0.151
gendstat: Turkey_N_published Russia_Steppe_Eneolithic -0.007
gendstat: Turkey_N_published ONG.SG 0.652
gendstat: Russia_MA1_HG.SG Georgia_Satsurblia.SG 0.632
gendstat: Russia_MA1_HG.SG Georgia_Kotias.SG 0.527
gendstat: Russia_MA1_HG.SG Papuan.DG -0.238
gendstat: Russia_MA1_HG.SG Russia_DevilsCave_N.SG 1.371
gendstat: Russia_MA1_HG.SG Israel_Natufian_published 0.414
gendstat: Russia_MA1_HG.SG Morocco_Iberomaurusian 0.523
gendstat: Russia_MA1_HG.SG Russia_Steppe_Eneolithic 0.452
gendstat: Russia_MA1_HG.SG ONG.SG 0.921
gendstat: Georgia_Satsurblia.SG Georgia_Kotias.SG -0.174
gendstat: Georgia_Satsurblia.SG Papuan.DG -0.974
gendstat: Georgia_Satsurblia.SG Russia_DevilsCave_N.SG 0.607
gendstat: Georgia_Satsurblia.SG Israel_Natufian_published -0.193
gendstat: Georgia_Satsurblia.SG Morocco_Iberomaurusian -0.229
gendstat: Georgia_Satsurblia.SG Russia_Steppe_Eneolithic -0.388
gendstat: Georgia_Satsurblia.SG ONG.SG 0.145
gendstat: Georgia_Kotias.SG Papuan.DG -0.830
gendstat: Georgia_Kotias.SG Russia_DevilsCave_N.SG 0.795
gendstat: Georgia_Kotias.SG Israel_Natufian_published -0.054
gendstat: Georgia_Kotias.SG Morocco_Iberomaurusian -0.070
gendstat: Georgia_Kotias.SG Russia_Steppe_Eneolithic -0.218
gendstat: Georgia_Kotias.SG ONG.SG 0.310
gendstat: Papuan.DG Russia_DevilsCave_N.SG 2.224
gendstat: Papuan.DG Israel_Natufian_published 0.671
gendstat: Papuan.DG Morocco_Iberomaurusian 0.938
gendstat: Papuan.DG Russia_Steppe_Eneolithic 0.847
gendstat: Papuan.DG ONG.SG 1.713
gendstat: Russia_DevilsCave_N.SG Israel_Natufian_published -0.752
gendstat: Russia_DevilsCave_N.SG Morocco_Iberomaurusian -1.054
gendstat: Russia_DevilsCave_N.SG Russia_Steppe_Eneolithic -1.290
gendstat: Russia_DevilsCave_N.SG ONG.SG -0.660
gendstat: Israel_Natufian_published Morocco_Iberomaurusian -0.004
gendstat: Israel_Natufian_published Russia_Steppe_Eneolithic -0.113
gendstat: Israel_Natufian_published ONG.SG 0.324
gendstat: Morocco_Iberomaurusian Russia_Steppe_Eneolithic -0.145
gendstat: Morocco_Iberomaurusian ONG.SG 0.445
gendstat: Russia_Steppe_Eneolithic ONG.SG 0.669

##end of qpAdm: 7.179 seconds cpu 1.114 Mbytes in use



3) Bulgaria_IA 79.1% + 20.9% Ukraine_IA_WesternScythian_Cimmerians_o2.SG. This model puzzles me as it rejects South Caucasian/Anatolian references.


qpAdm: parameter file: parqpadm.txt
### THE INPUT PARAMETERS
##PARAMETER NAME: VALUE
fstatsname: fstatsa.txt
popleft: left.txt
popright: right.txt
details: YES
## qpAdm version: 1520
seed: 1774793406
*** recommended that inbreed be explicitly set ***
inbreed set NO

left pops:
eupator
Bulgaria_IA
Ukraine_IA_WesternScythian_Cimmerians_o2.SG

right pops:
Mbuti.DG
Russia_Ust_Ishim.DG
China_Tianyuan
Goyet_Neanderthal.SG
Turkey_N_published
Russia_MA1_HG.SG
Georgia_Satsurblia.SG
Papuan.DG
Israel_Natufian_published
Russia_DevilsCave_N.SG
Russia_Steppe_Eneolithic
Russia_HG_Samara
ONG.SG

codimension 1
f4info:
f4rank: 1 dof: 11 chisq: 15.490 tail: 0.161129019 dofdiff: 13 chisqdiff: -15.490 taildiff: 1
B:
scale 1.000
Russia_Ust_Ishim.DG 0.407
China_Tianyuan 0.465
Goyet_Neanderthal.SG -0.073
Turkey_N_published -1.609
Russia_MA1_HG.SG 1.441
Georgia_Satsurblia.SG -0.031
Papuan.DG 0.607
Israel_Natufian_published -1.599
Russia_DevilsCave_N.SG 1.291
Russia_Steppe_Eneolithic 0.398
Russia_HG_Samara 1.245
ONG.SG 0.804
A:
scale 434.269
Bulgaria_IA -0.360
Ukraine_IA_WesternScythian_Cimmerians_o2.SG 1.368


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 11 chisqdiff: 15.490 taildiff: 0.161129019
B:
scale 1001.280 283.461
Russia_Ust_Ishim.DG 0.330 0.490
China_Tianyuan 0.713 0.698
Goyet_Neanderthal.SG 0.129 -0.037
Turkey_N_published 2.506 -1.192
Russia_MA1_HG.SG -0.453 1.476
Georgia_Satsurblia.SG 0.505 0.097
Papuan.DG 0.062 0.667
Israel_Natufian_published 1.963 -1.361
Russia_DevilsCave_N.SG 0.420 1.505
Russia_Steppe_Eneolithic 0.246 0.450
Russia_HG_Samara 0.706 1.486
ONG.SG -0.170 0.831
A:
scale 1.414 1.414
Bulgaria_IA 1.414 0.000
Ukraine_IA_WesternScythian_Cimmerians_o2.SG 0.000 1.414


best coefficients: 0.791 0.209
totmean: 0.791 0.209
boot mean: 0.792 0.208
std. errors: 0.038 0.038

error covariance (* 1,000,000)
1443 -1443
-1443 1443


summ: eupator 2 0.161129 0.792 0.208 1443 -1443 1443

fixed pat wt dof chisq tail prob
00 0 11 15.490 0.161129 0.791 0.209
01 1 12 42.261 3.00953e-05 1.000 0.000
10 1 12 303.078 0 0.000 1.000
best pat: 00 0.161129 - -
best pat: 01 3.00953e-05 chi(nested): 26.771 p-value for nested model: 2.29057e-07

coeffs: 0.791 0.209

## dscore:: f_4(Base, Fit, Rbase, right2)
## genstat:: f_4(Base, Fit, right1, right2)

details: Bulgaria_IA Russia_Ust_Ishim.DG 0.000330 0.525065
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Russia_Ust_Ishim.DG 0.001729 2.377501
dscore: Russia_Ust_Ishim.DG f4: 0.000622 Z: 1.088903

details: Bulgaria_IA China_Tianyuan 0.000712 1.070427
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG China_Tianyuan 0.002462 3.169412
dscore: China_Tianyuan f4: 0.001077 Z: 1.770830

details: Bulgaria_IA Goyet_Neanderthal.SG 0.000129 0.277287
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Goyet_Neanderthal.SG -0.000132 -0.228959
dscore: Goyet_Neanderthal.SG f4: 0.000075 Z: 0.181094

details: Bulgaria_IA Turkey_N_published 0.002503 5.308251
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Turkey_N_published -0.004204 -7.283983
dscore: Turkey_N_published f4: 0.001105 Z: 2.562839

details: Bulgaria_IA Russia_MA1_HG.SG -0.000453 -0.615556
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Russia_MA1_HG.SG 0.005209 5.905472
dscore: Russia_MA1_HG.SG f4: 0.000728 Z: 1.074979

details: Bulgaria_IA Georgia_Satsurblia.SG 0.000505 0.727765
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Georgia_Satsurblia.SG 0.000343 0.412815
dscore: Georgia_Satsurblia.SG f4: 0.000471 Z: 0.734669

details: Bulgaria_IA Papuan.DG 0.000062 0.118459
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Papuan.DG 0.002352 4.031770
dscore: Papuan.DG f4: 0.000539 Z: 1.150661

details: Bulgaria_IA Israel_Natufian_published 0.001961 2.525743
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Israel_Natufian_published -0.004800 -4.258779
dscore: Israel_Natufian_published f4: 0.000551 Z: 0.779963

details: Bulgaria_IA Russia_DevilsCave_N.SG 0.000419 0.769164
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Russia_DevilsCave_N.SG 0.005310 7.985637
dscore: Russia_DevilsCave_N.SG f4: 0.001439 Z: 2.888103

details: Bulgaria_IA Russia_Steppe_Eneolithic 0.000246 0.468313
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Russia_Steppe_Eneolithic 0.001589 2.535702
dscore: Russia_Steppe_Eneolithic f4: 0.000526 Z: 1.098264

details: Bulgaria_IA Russia_HG_Samara 0.000705 0.979550
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG Russia_HG_Samara 0.005242 5.795423
dscore: Russia_HG_Samara f4: 0.001651 Z: 2.510476

details: Bulgaria_IA ONG.SG -0.000170 -0.323666
details: Ukraine_IA_WesternScythian_Cimmerians_o2.SG ONG.SG 0.002931 4.964356
dscore: ONG.SG f4: 0.000477 Z: 0.994976

gendstat: Mbuti.DG Russia_Ust_Ishim.DG 1.089
gendstat: Mbuti.DG China_Tianyuan 1.771
gendstat: Mbuti.DG Goyet_Neanderthal.SG 0.181
gendstat: Mbuti.DG Turkey_N_published 2.563
gendstat: Mbuti.DG Russia_MA1_HG.SG 1.075
gendstat: Mbuti.DG Georgia_Satsurblia.SG 0.735
gendstat: Mbuti.DG Papuan.DG 1.151
gendstat: Mbuti.DG Israel_Natufian_published 0.780
gendstat: Mbuti.DG Russia_DevilsCave_N.SG 2.888
gendstat: Mbuti.DG Russia_Steppe_Eneolithic 1.098
gendstat: Mbuti.DG Russia_HG_Samara 2.510
gendstat: Mbuti.DG ONG.SG 0.995
gendstat: Russia_Ust_Ishim.DG China_Tianyuan 0.633
gendstat: Russia_Ust_Ishim.DG Goyet_Neanderthal.SG -0.808
gendstat: Russia_Ust_Ishim.DG Turkey_N_published 0.821
gendstat: Russia_Ust_Ishim.DG Russia_MA1_HG.SG 0.139
gendstat: Russia_Ust_Ishim.DG Georgia_Satsurblia.SG -0.204
gendstat: Russia_Ust_Ishim.DG Papuan.DG -0.136
gendstat: Russia_Ust_Ishim.DG Israel_Natufian_published -0.084
gendstat: Russia_Ust_Ishim.DG Russia_DevilsCave_N.SG 1.292
gendstat: Russia_Ust_Ishim.DG Russia_Steppe_Eneolithic -0.156
gendstat: Russia_Ust_Ishim.DG Russia_HG_Samara 1.313
gendstat: Russia_Ust_Ishim.DG ONG.SG -0.237
gendstat: China_Tianyuan Goyet_Neanderthal.SG -1.359
gendstat: China_Tianyuan Turkey_N_published 0.043
gendstat: China_Tianyuan Russia_MA1_HG.SG -0.447
gendstat: China_Tianyuan Georgia_Satsurblia.SG -0.761
gendstat: China_Tianyuan Papuan.DG -0.890
gendstat: China_Tianyuan Israel_Natufian_published -0.606
gendstat: China_Tianyuan Russia_DevilsCave_N.SG 0.586
gendstat: China_Tianyuan Russia_Steppe_Eneolithic -0.878
gendstat: China_Tianyuan Russia_HG_Samara 0.734
gendstat: China_Tianyuan ONG.SG -1.009
gendstat: Goyet_Neanderthal.SG Turkey_N_published 1.808
gendstat: Goyet_Neanderthal.SG Russia_MA1_HG.SG 0.851
gendstat: Goyet_Neanderthal.SG Georgia_Satsurblia.SG 0.532
gendstat: Goyet_Neanderthal.SG Papuan.DG 0.765
gendstat: Goyet_Neanderthal.SG Israel_Natufian_published 0.600
gendstat: Goyet_Neanderthal.SG Russia_DevilsCave_N.SG 2.156
gendstat: Goyet_Neanderthal.SG Russia_Steppe_Eneolithic 0.735
gendstat: Goyet_Neanderthal.SG Russia_HG_Samara 2.082
gendstat: Goyet_Neanderthal.SG ONG.SG 0.632
gendstat: Turkey_N_published Russia_MA1_HG.SG -0.555
gendstat: Turkey_N_published Georgia_Satsurblia.SG -1.010
gendstat: Turkey_N_published Papuan.DG -1.149
gendstat: Turkey_N_published Israel_Natufian_published -0.784
gendstat: Turkey_N_published Russia_DevilsCave_N.SG 0.655
gendstat: Turkey_N_published Russia_Steppe_Eneolithic -1.286
gendstat: Turkey_N_published Russia_HG_Samara 0.874
gendstat: Turkey_N_published ONG.SG -1.272
gendstat: Russia_MA1_HG.SG Georgia_Satsurblia.SG -0.314
gendstat: Russia_MA1_HG.SG Papuan.DG -0.285
gendstat: Russia_MA1_HG.SG Israel_Natufian_published -0.202
gendstat: Russia_MA1_HG.SG Russia_DevilsCave_N.SG 1.055
gendstat: Russia_MA1_HG.SG Russia_Steppe_Eneolithic -0.326
gendstat: Russia_MA1_HG.SG Russia_HG_Samara 1.212
gendstat: Russia_MA1_HG.SG ONG.SG -0.387
gendstat: Georgia_Satsurblia.SG Papuan.DG 0.103
gendstat: Georgia_Satsurblia.SG Israel_Natufian_published 0.093
gendstat: Georgia_Satsurblia.SG Russia_DevilsCave_N.SG 1.458
gendstat: Georgia_Satsurblia.SG Russia_Steppe_Eneolithic 0.088
gendstat: Georgia_Satsurblia.SG Russia_HG_Samara 1.495
gendstat: Georgia_Satsurblia.SG ONG.SG 0.009
gendstat: Papuan.DG Israel_Natufian_published 0.015
gendstat: Papuan.DG Russia_DevilsCave_N.SG 1.915
gendstat: Papuan.DG Russia_Steppe_Eneolithic -0.028
gendstat: Papuan.DG Russia_HG_Samara 1.669
gendstat: Papuan.DG ONG.SG -0.141
gendstat: Israel_Natufian_published Russia_DevilsCave_N.SG 1.128
gendstat: Israel_Natufian_published Russia_Steppe_Eneolithic -0.034
gendstat: Israel_Natufian_published Russia_HG_Samara 1.324
gendstat: Israel_Natufian_published ONG.SG -0.096
gendstat: Russia_DevilsCave_N.SG Russia_Steppe_Eneolithic -1.794
gendstat: Russia_DevilsCave_N.SG Russia_HG_Samara 0.326
gendstat: Russia_DevilsCave_N.SG ONG.SG -2.090
gendstat: Russia_Steppe_Eneolithic Russia_HG_Samara 1.784
gendstat: Russia_Steppe_Eneolithic ONG.SG -0.101
gendstat: Russia_HG_Samara ONG.SG -1.764

##end of qpAdm: 1.925 seconds cpu 0.725 Mbytes in use



4) Greece_BA_Mycenaean 77.4% + Russia_Samara_EBA_Yamnaya 5% + Iran_GanjDareh_N 17.7%.



qpAdm: parameter file: parqpadm.txt
### THE INPUT PARAMETERS
##PARAMETER NAME: VALUE
fstatsname: fstatsa.txt
popleft: left.txt
popright: right.txt
details: YES
## qpAdm version: 1520
seed: 1860613956
*** recommended that inbreed be explicitly set ***
inbreed set NO

left pops:
eupator
Greece_BA_Mycenaean
Russia_Samara_EBA_Yamnaya
Iran_GanjDareh_N

right pops:
Mbuti.DG
Russia_Ust_Ishim.DG
China_Tianyuan
Goyet_Neanderthal.SG
Turkey_N_published
Russia_MA1_HG.SG
Georgia_Satsurblia.SG
Papuan.DG
Israel_Natufian_published
Russia_DevilsCave_N.SG
Georgia_Kotias.SG
Russia_Steppe_Eneolithic
Greece_Minoan_Lassithi
ONG.SG

codimension 1
f4info:
f4rank: 2 dof: 11 chisq: 10.305 tail: 0.503151124 dofdiff: 13 chisqdiff: -10.305 taildiff: 1
B:
scale 1.000 1.000
Russia_Ust_Ishim.DG -0.097 0.625
China_Tianyuan -0.323 0.679
Goyet_Neanderthal.SG 0.083 0.165
Turkey_N_published 2.147 0.748
Russia_MA1_HG.SG -0.874 2.195
Georgia_Satsurblia.SG -0.474 0.485
Papuan.DG -0.099 0.523
Israel_Natufian_published 1.684 0.640
Russia_DevilsCave_N.SG -0.204 0.806
Georgia_Kotias.SG -0.538 0.660
Russia_Steppe_Eneolithic -0.736 1.919
Greece_Minoan_Lassithi 1.871 0.781
ONG.SG -0.243 0.666
A:
scale 479.711 639.859
Greece_BA_Mycenaean 0.366 0.192
Russia_Samara_EBA_Yamnaya -1.092 1.245
Iran_GanjDareh_N -1.294 -1.189


full rank
f4info:
f4rank: 3 dof: 0 chisq: 0.000 tail: 1 dofdiff: 11 chisqdiff: 10.305 taildiff: 0.503151124
B:
scale 864.297 322.019 321.857
Russia_Ust_Ishim.DG 0.300 0.533 -0.234
China_Tianyuan 0.533 0.810 0.029
Goyet_Neanderthal.SG -0.097 -0.017 -0.236
Turkey_N_published 1.712 -1.076 -2.287
Russia_MA1_HG.SG -0.116 2.015 -0.571
Georgia_Satsurblia.SG -0.210 0.662 0.124
Papuan.DG 0.227 0.457 -0.169
Israel_Natufian_published 2.587 -0.622 -1.606
Russia_DevilsCave_N.SG 0.838 0.843 -0.088
Georgia_Kotias.SG -0.273 0.777 0.050
Russia_Steppe_Eneolithic 0.103 1.789 -0.468
Greece_Minoan_Lassithi 1.393 -0.895 -2.114
ONG.SG 0.397 0.701 -0.085
A:
scale 1.732 1.732 1.732
Greece_BA_Mycenaean 1.732 0.000 0.000
Russia_Samara_EBA_Yamnaya 0.000 1.732 0.000
Iran_GanjDareh_N 0.000 0.000 1.732


best coefficients: 0.774 0.050 0.177
totmean: 0.774 0.050 0.177
boot mean: 0.774 0.049 0.177
std. errors: 0.033 0.053 0.053

error covariance (* 1,000,000)
1101 -511 -590
-511 2759 -2248
-590 -2248 2839


summ: eupator 3 0.503151 0.774 0.049 0.177 1101 -511 -590 2759 -2248 ...
2839

fixed pat wt dof chisq tail prob
000 0 11 10.305 0.503151 0.774 0.050 0.177
001 1 12 20.956 0.0510355 0.814 0.186 0.000
010 1 12 13.788 0.314478 0.785 0.000 0.215
100 1 12 571.114 0 0.000 -3.927 4.927 infeasible
011 2 13 48.693 4.98778e-06 1.000 0.000 0.000
101 2 13 741.780 0 0.000 1.000 0.000
110 2 13 611.970 0 0.000 0.000 1.000
best pat: 000 0.503151 - -
best pat: 010 0.314478 chi(nested): 3.482 p-value for nested model: 0.0620295
best pat: 011 4.98778e-06 chi(nested): 34.905 p-value for nested model: 3.4613e-09

coeffs: 0.774 0.050 0.177

## dscore:: f_4(Base, Fit, Rbase, right2)
## genstat:: f_4(Base, Fit, right1, right2)

details: Greece_BA_Mycenaean Russia_Ust_Ishim.DG 0.000347 0.629351
details: Russia_Samara_EBA_Yamnaya Russia_Ust_Ishim.DG 0.001654 3.383814
details: Iran_GanjDareh_N Russia_Ust_Ishim.DG -0.000727 -1.407089
dscore: Russia_Ust_Ishim.DG f4: 0.000222 Z: 0.435703

details: Greece_BA_Mycenaean China_Tianyuan 0.000616 1.013922
details: Russia_Samara_EBA_Yamnaya China_Tianyuan 0.002514 4.821634
details: Iran_GanjDareh_N China_Tianyuan 0.000090 0.165611
dscore: China_Tianyuan f4: 0.000617 Z: 1.106874

details: Greece_BA_Mycenaean Goyet_Neanderthal.SG -0.000112 -0.276273
details: Russia_Samara_EBA_Yamnaya Goyet_Neanderthal.SG -0.000054 -0.163010
details: Iran_GanjDareh_N Goyet_Neanderthal.SG -0.000734 -2.049430
dscore: Goyet_Neanderthal.SG f4: -0.000219 Z: -0.601257

details: Greece_BA_Mycenaean Turkey_N_published 0.001980 4.550947
details: Russia_Samara_EBA_Yamnaya Turkey_N_published -0.003342 -8.904749
details: Iran_GanjDareh_N Turkey_N_published -0.007107 -17.497838
dscore: Turkey_N_published f4: 0.000109 Z: 0.273756

details: Greece_BA_Mycenaean Russia_MA1_HG.SG -0.000134 -0.213026
details: Russia_Samara_EBA_Yamnaya Russia_MA1_HG.SG 0.006258 11.273361
details: Iran_GanjDareh_N Russia_MA1_HG.SG -0.001775 -2.925709
dscore: Russia_MA1_HG.SG f4: -0.000107 Z: -0.183042

details: Greece_BA_Mycenaean Georgia_Satsurblia.SG -0.000243 -0.403118
details: Russia_Samara_EBA_Yamnaya Georgia_Satsurblia.SG 0.002055 3.820292
details: Iran_GanjDareh_N Georgia_Satsurblia.SG 0.000384 0.659024
dscore: Georgia_Satsurblia.SG f4: -0.000018 Z: -0.032347

details: Greece_BA_Mycenaean Papuan.DG 0.000262 0.567161
details: Russia_Samara_EBA_Yamnaya Papuan.DG 0.001419 3.423867
details: Iran_GanjDareh_N Papuan.DG -0.000524 -1.191074
dscore: Papuan.DG f4: 0.000181 Z: 0.421911

details: Greece_BA_Mycenaean Israel_Natufian_published 0.002993 3.872332
details: Russia_Samara_EBA_Yamnaya Israel_Natufian_published -0.001931 -3.417930
details: Iran_GanjDareh_N Israel_Natufian_published -0.004991 -8.164318
dscore: Israel_Natufian_published f4: 0.001337 Z: 1.970652

details: Greece_BA_Mycenaean Russia_DevilsCave_N.SG 0.000970 1.960136
details: Russia_Samara_EBA_Yamnaya Russia_DevilsCave_N.SG 0.002618 6.121602
details: Iran_GanjDareh_N Russia_DevilsCave_N.SG -0.000274 -0.566581
dscore: Russia_DevilsCave_N.SG f4: 0.000831 Z: 1.815686

details: Greece_BA_Mycenaean Georgia_Kotias.SG -0.000315 -0.533495
details: Russia_Samara_EBA_Yamnaya Georgia_Kotias.SG 0.002413 4.730769
details: Iran_GanjDareh_N Georgia_Kotias.SG 0.000156 0.284895
dscore: Georgia_Kotias.SG f4: -0.000097 Z: -0.177323

details: Greece_BA_Mycenaean Russia_Steppe_Eneolithic 0.000119 0.256099
details: Russia_Samara_EBA_Yamnaya Russia_Steppe_Eneolithic 0.005556 13.362992
details: Iran_GanjDareh_N Russia_Steppe_Eneolithic -0.001454 -3.216285
dscore: Russia_Steppe_Eneolithic f4: 0.000111 Z: 0.256668

details: Greece_BA_Mycenaean Greece_Minoan_Lassithi 0.001612 3.555173
details: Russia_Samara_EBA_Yamnaya Greece_Minoan_Lassithi -0.002779 -7.177125
details: Iran_GanjDareh_N Greece_Minoan_Lassithi -0.006568 -16.139615
dscore: Greece_Minoan_Lassithi f4: -0.000053 Z: -0.128028

details: Greece_BA_Mycenaean ONG.SG 0.000459 0.997032
details: Russia_Samara_EBA_Yamnaya ONG.SG 0.002176 5.227393
details: Iran_GanjDareh_N ONG.SG -0.000264 -0.581036
dscore: ONG.SG f4: 0.000417 Z: 0.968407

gendstat: Mbuti.DG Russia_Ust_Ishim.DG 0.436
gendstat: Mbuti.DG China_Tianyuan 1.107
gendstat: Mbuti.DG Goyet_Neanderthal.SG -0.601
gendstat: Mbuti.DG Turkey_N_published 0.274
gendstat: Mbuti.DG Russia_MA1_HG.SG -0.183
gendstat: Mbuti.DG Georgia_Satsurblia.SG -0.032
gendstat: Mbuti.DG Papuan.DG 0.422
gendstat: Mbuti.DG Israel_Natufian_published 1.971
gendstat: Mbuti.DG Russia_DevilsCave_N.SG 1.816
gendstat: Mbuti.DG Georgia_Kotias.SG -0.177
gendstat: Mbuti.DG Russia_Steppe_Eneolithic 0.257
gendstat: Mbuti.DG Greece_Minoan_Lassithi -0.128
gendstat: Mbuti.DG ONG.SG 0.968
gendstat: Russia_Ust_Ishim.DG China_Tianyuan 0.649
gendstat: Russia_Ust_Ishim.DG Goyet_Neanderthal.SG -0.746
gendstat: Russia_Ust_Ishim.DG Turkey_N_published -0.226
gendstat: Russia_Ust_Ishim.DG Russia_MA1_HG.SG -0.507
gendstat: Russia_Ust_Ishim.DG Georgia_Satsurblia.SG -0.372
gendstat: Russia_Ust_Ishim.DG Papuan.DG -0.083
gendstat: Russia_Ust_Ishim.DG Israel_Natufian_published 1.474
gendstat: Russia_Ust_Ishim.DG Russia_DevilsCave_N.SG 1.186
gendstat: Russia_Ust_Ishim.DG Georgia_Kotias.SG -0.515
gendstat: Russia_Ust_Ishim.DG Russia_Steppe_Eneolithic -0.212
gendstat: Russia_Ust_Ishim.DG Greece_Minoan_Lassithi -0.536
gendstat: Russia_Ust_Ishim.DG ONG.SG 0.399
gendstat: China_Tianyuan Goyet_Neanderthal.SG -1.289
gendstat: China_Tianyuan Turkey_N_published -0.943
gendstat: China_Tianyuan Russia_MA1_HG.SG -1.050
gendstat: China_Tianyuan Georgia_Satsurblia.SG -0.903
gendstat: China_Tianyuan Papuan.DG -0.827
gendstat: China_Tianyuan Israel_Natufian_published 0.898
gendstat: China_Tianyuan Russia_DevilsCave_N.SG 0.390
gendstat: China_Tianyuan Georgia_Kotias.SG -1.074
gendstat: China_Tianyuan Russia_Steppe_Eneolithic -0.885
gendstat: China_Tianyuan Greece_Minoan_Lassithi -1.170
gendstat: China_Tianyuan ONG.SG -0.377
gendstat: Goyet_Neanderthal.SG Turkey_N_published 0.645
gendstat: Goyet_Neanderthal.SG Russia_MA1_HG.SG 0.171
gendstat: Goyet_Neanderthal.SG Georgia_Satsurblia.SG 0.313
gendstat: Goyet_Neanderthal.SG Papuan.DG 0.761
gendstat: Goyet_Neanderthal.SG Israel_Natufian_published 2.056
gendstat: Goyet_Neanderthal.SG Russia_DevilsCave_N.SG 1.938
gendstat: Goyet_Neanderthal.SG Georgia_Kotias.SG 0.192
gendstat: Goyet_Neanderthal.SG Russia_Steppe_Eneolithic 0.616
gendstat: Goyet_Neanderthal.SG Greece_Minoan_Lassithi 0.318
gendstat: Goyet_Neanderthal.SG ONG.SG 1.192
gendstat: Turkey_N_published Russia_MA1_HG.SG -0.364
gendstat: Turkey_N_published Georgia_Satsurblia.SG -0.233
gendstat: Turkey_N_published Papuan.DG 0.162
gendstat: Turkey_N_published Israel_Natufian_published 1.883
gendstat: Turkey_N_published Russia_DevilsCave_N.SG 1.564
gendstat: Turkey_N_published Georgia_Kotias.SG -0.410
gendstat: Turkey_N_published Russia_Steppe_Eneolithic 0.004
gendstat: Turkey_N_published Greece_Minoan_Lassithi -0.523
gendstat: Turkey_N_published ONG.SG 0.715
gendstat: Russia_MA1_HG.SG Georgia_Satsurblia.SG 0.124
gendstat: Russia_MA1_HG.SG Papuan.DG 0.496
gendstat: Russia_MA1_HG.SG Israel_Natufian_published 1.777
gendstat: Russia_MA1_HG.SG Russia_DevilsCave_N.SG 1.534
gendstat: Russia_MA1_HG.SG Georgia_Kotias.SG 0.015
gendstat: Russia_MA1_HG.SG Russia_Steppe_Eneolithic 0.397
gendstat: Russia_MA1_HG.SG Greece_Minoan_Lassithi 0.092
gendstat: Russia_MA1_HG.SG ONG.SG 0.913
gendstat: Georgia_Satsurblia.SG Papuan.DG 0.319
gendstat: Georgia_Satsurblia.SG Israel_Natufian_published 1.739
gendstat: Georgia_Satsurblia.SG Russia_DevilsCave_N.SG 1.392
gendstat: Georgia_Satsurblia.SG Georgia_Kotias.SG -0.138
gendstat: Georgia_Satsurblia.SG Russia_Steppe_Eneolithic 0.236
gendstat: Georgia_Satsurblia.SG Greece_Minoan_Lassithi -0.064
gendstat: Georgia_Satsurblia.SG ONG.SG 0.755
gendstat: Papuan.DG Israel_Natufian_published 1.587
gendstat: Papuan.DG Russia_DevilsCave_N.SG 1.550
gendstat: Papuan.DG Georgia_Kotias.SG -0.477
gendstat: Papuan.DG Russia_Steppe_Eneolithic -0.159
gendstat: Papuan.DG Greece_Minoan_Lassithi -0.518
gendstat: Papuan.DG ONG.SG 0.593
gendstat: Israel_Natufian_published Russia_DevilsCave_N.SG -0.711
gendstat: Israel_Natufian_published Georgia_Kotias.SG -1.879
gendstat: Israel_Natufian_published Russia_Steppe_Eneolithic -1.765
gendstat: Israel_Natufian_published Greece_Minoan_Lassithi -2.075
gendstat: Israel_Natufian_published ONG.SG -1.275
gendstat: Russia_DevilsCave_N.SG Georgia_Kotias.SG -1.643
gendstat: Russia_DevilsCave_N.SG Russia_Steppe_Eneolithic -1.569
gendstat: Russia_DevilsCave_N.SG Greece_Minoan_Lassithi -1.913
gendstat: Russia_DevilsCave_N.SG ONG.SG -1.018
gendstat: Georgia_Kotias.SG Russia_Steppe_Eneolithic 0.407
gendstat: Georgia_Kotias.SG Greece_Minoan_Lassithi 0.085
gendstat: Georgia_Kotias.SG ONG.SG 0.891
gendstat: Russia_Steppe_Eneolithic Greece_Minoan_Lassithi -0.410
gendstat: Russia_Steppe_Eneolithic ONG.SG 0.724
gendstat: Greece_Minoan_Lassithi ONG.SG 1.095

##end of qpAdm: 7.184 seconds cpu 1.114 Mbytes in use

kingmob
07-05-2022, 05:53 PM
FST distances run, using my own sample, in order to check how it fares against the Greek samples in the Reich dataset (the ones used in Lazaridis' paper, GREEKALPOP, etc).

Compared to Myceneans:



fst(prefix, pop1 = "Greece_BA_Mycenaean", pop2 = c("eptr", "Greek"))
ℹ Reading allele frequencies from packedancestrymap files...
ℹ eptrfamilyHO.geno has 14317 samples and 597573 SNPs
ℹ Calculating allele frequencies from 24 samples in 3 populations
ℹ Expected size of allele frequency data: 81 MB
597k SNPs read...
✔ 597573 SNPs read in total
! 169155 SNPs remain after filtering. 156076 are polymorphic.
ℹ Allele frequency matrix for 169155 SNPs and 3 populations is 16 MB
ℹ Computing pairwise f2 for all SNPs and population pairs requires 97 MB RAM without splitting
ℹ Computing without splitting since 97 < 8000 (maxmem)...
ℹ Returning fst blocks
# A tibble: 2 × 4
pop1 pop2 est se
<chr> <chr> <dbl> <dbl>
1 Greece_BA_Mycenaean eptr 0.00227 0.00325
2 Greece_BA_Mycenaean Greek 0.00792 0.00136

Compared to Imperial_Roman:



fst(prefix, pop1 = "Italy_Imperial.SG", pop2 = c("eptr", "Greek"), adjust_pseudohaploid = FALSE)
ℹ Reading allele frequencies from packedancestrymap files...
ℹ eptrfamilyHO.geno has 14317 samples and 597573 SNPs
ℹ Calculating allele frequencies from 56 samples in 3 populations
ℹ Expected size of allele frequency data: 81 MB
597k SNPs read...
✔ 597573 SNPs read in total
! 245195 SNPs remain after filtering. 231896 are polymorphic.
ℹ Allele frequency matrix for 245195 SNPs and 3 populations is 24 MB
ℹ Computing pairwise f2 for all SNPs and population pairs requires 141 MB RAM without splitting
ℹ Computing without splitting since 141 < 8000 (maxmem)...
ℹ Returning fst blocks
# A tibble: 2 × 4
pop1 pop2 est se
<chr> <chr> <dbl> <dbl>
1 Italy_Imperial.SG eptr 0.00957 0.00232
2 Italy_Imperial.SG Greek 0.0141 0.000197

Looks like I win against the average, in both counts.

Voskos
07-05-2022, 06:03 PM
What does the first post mean? I get you're closer than average to Myc. and Rome, that's post 2. I don't understand the goal/context of post 1.

kingmob
07-05-2022, 06:08 PM
What does the first post mean? I get you're closer than average to Myc. and Rome, that's post 2. I don't understand the goal/context of post 1.

It is different models using several samples as sources with varying levels of success (p values), similar to what people usually do with g25 and admixture calculators.

If you have any requests (not for me necessarily) I can give it a try.

Roy
07-05-2022, 08:24 PM
What is it supposed to be - that simulation? About you? I am out of the loop.

kingmob
07-06-2022, 08:46 AM
What is it supposed to be - that simulation? About you? I am out of the loop.


Yes, it's modelling of myself using admixtools/qpAdm, the packages that are used for the academic papers produced by Reich Labs and the likes.