PDA

View Full Version : Modern model in 2018 natufians. How accurate is this



Gentica277282
02-13-2022, 09:22 PM
A 2018 analysis of autosomal DNA using modern populations as a reference, found The Natufian sample consisted of 61.2% Arabian, 21.2% Northern African, 10.9% Western Asian, and 6.8% Omotic-related ancestry (related to the Omotic peoples of southern Ethiopia).

Hitchcock
02-14-2022, 01:56 PM
This is modern admixture estimator , is not meaningful with deep-lineages like the Natufians , it would just show the ancestry of modern as ancient.
The correct composition of Natufians is Basal Eurasian + West Eurasian.
The Natufians cannot be modelled with modern , but there west eurasian ancestry is related to the common Ancestry found in Mesolithic Europe.
they don't have Iberomaurusian admixture as the Natufians don't have the common drift with West african lineages that the Iberomaurusians do.

Ruderico
02-14-2022, 02:04 PM
Chad made a blog post which included Natufians in qpGraph a few years ago
https://populationgenomics.blog/2019/05/23/the-big-picture-for-west-asia-and-specifically-anatolia-after-pinarbasi/


To begin this simple graph, I started with the older samples from each region, with Natufians being a mix between Iberomaurusians and ancient Anatolians, which does fit fine and first demonstrated by Lazaridis et al (2018)(pre-print). Since I do not have the Dzudzuana genomes, I used Anatolia_HG (Pinarbasi) as the stand-in for that admixture into the Iberomaurusians (IBM).

Gentica277282
02-14-2022, 02:47 PM
Chad made a blog post which included Natufians in qpGraph a few years ago
https://populationgenomics.blog/2019/05/23/the-big-picture-for-west-asia-and-specifically-anatolia-after-pinarbasi/

Target: Levant_Natufian:I0861
Distance: 3.0161% / 0.03016148
80.4 Dzudzuana
12.0 Ancestral_North-African(DZU-55%-45%)
7.6 Iran_Neolithic


Target: Levant_Natufian_contam:I1072
Distance: 4.1042% / 0.04104191
86.4 Dzudzuana
13.6 Ancestral_North-African(DZU-55%-45%)

Mnemonics
02-14-2022, 03:02 PM
A 2018 analysis of autosomal DNA using modern populations as a reference, found The Natufian sample consisted of 61.2% Arabian, 21.2% Northern African, 10.9% Western Asian, and 6.8% Omotic-related ancestry (related to the Omotic peoples of southern Ethiopia).


It seems roughly in line with what I get for them in qpAdm despite all those modern samples having some degree of Natufian admixture.







left pops:
Natufian
Iberomaurusian
Anatolia_Epipaleolithic
Iran_N
Tanzania_Kisese

right pops:
Cameroon_SMA.SG
South_Africa_400BP.SG
South_Africa_HG
Eastern_Euro_UP.SG
BachoKiro
ZlatyKun
DevilsGate_N.SG
Yana_UP
CHG
Iran_N.SG
British_Mesolithic.SG
Mota
Morocco_EN.SG

0 Natufian 6
1 Iberomaurusian 7
2 Anatolia_Epipaleolithic 1
3 Iran_N 4
4 Tanzania_Kisese 1
5 Cameroon_SMA.SG 3
6 South_Africa_400BP.SG 4
7 South_Africa_HG 4
8 Eastern_Euro_UP.SG 3
9 BachoKiro 3
10 ZlatyKun 1
11 DevilsGate_N.SG 4
12 Yana_UP 2
13 CHG 2
14 Iran_N.SG 4
15 British_Mesolithic.SG 3
16 Mota 1
17 Morocco_EN.SG 3
jackknife block size: 0.050
snps: 1150461 indivs: 56
number of blocks for block jackknife: 714
## ncols: 1150461
coverage: Natufian 529354
coverage: Iberomaurusian 1088914
coverage: Anatolia_Epipaleolithic 838782
coverage: Iran_N 1002874
coverage: Tanzania_Kisese 938012
coverage: Cameroon_SMA.SG 1149654
coverage: South_Africa_400BP.SG 1150180
coverage: South_Africa_HG 1135129
coverage: Eastern_Euro_UP.SG 1150149
coverage: BachoKiro 990561
coverage: ZlatyKun 1120502
coverage: DevilsGate_N.SG 1134715
coverage: Yana_UP 1135189
coverage: CHG 1134595
coverage: Iran_N.SG 1149571
coverage: British_Mesolithic.SG 1149936
coverage: Mota 1134612
coverage: Morocco_EN.SG 653097
Effective number of blocks: 616.513
numsnps used: 319495
codimension 1
f4info:
f4rank: 3 dof: 9 chisq: 6.409 tail: 0.698383677 dofdiff: 11 chisqdiff: -6.409 taildiff: 1
B:
scale 1.000 1.000 1.000
South_Africa_400BP.SG 0.088 0.136 -0.078
South_Africa_HG -0.208 -0.165 -0.436
Eastern_Euro_UP.SG 1.221 -0.260 -0.913
BachoKiro 0.969 -0.184 0.055
ZlatyKun 1.011 -0.273 0.008
DevilsGate_N.SG 1.073 -0.234 0.201
Yana_UP 1.151 -0.244 -0.492
CHG 1.325 -0.629 0.802
Iran_N.SG 1.336 -0.834 2.737
British_Mesolithic.SG 1.334 -0.278 -1.591
Mota 0.103 0.076 -0.124
Morocco_EN.SG 0.832 3.239 0.065
A:
scale 105.341 314.313 879.994
Iberomaurusian -0.364 1.828 0.731
Anatolia_Epipaleolithic 0.221 -0.344 -0.507
Iran_N 0.080 -0.587 1.723
Tanzania_Kisese -1.952 -0.444 0.489


full rank
f4info:
f4rank: 4 dof: 0 chisq: 0.000 tail: 1 dofdiff: 9 chisqdiff: 6.409 taildiff: 0.698383677
B:
scale 142.536 460.665 326.101 52.741
South_Africa_400BP.SG 0.028 -0.142 -0.168 -0.120
South_Africa_HG -0.114 -0.169 -0.279 0.181
Eastern_Euro_UP.SG -0.980 1.326 -0.234 -1.237
BachoKiro -0.665 0.833 0.312 -0.963
ZlatyKun -0.793 0.803 0.309 -1.013
DevilsGate_N.SG -0.736 0.963 0.453 -1.055
Yana_UP -0.879 1.184 -0.010 -1.160
CHG -1.094 1.302 1.231 -1.249
Iran_N.SG -1.064 0.823 2.523 -1.194
British_Mesolithic.SG -1.103 1.784 -0.579 -1.358
Mota -0.012 0.030 -0.107 -0.125
Morocco_EN.SG 2.257 -0.984 -1.793 -1.079
A:
scale 2.000 2.000 2.000 2.000
Iberomaurusian 2.000 0.000 0.000 0.000
Anatolia_Epipaleolithic 0.000 2.000 0.000 0.000
Iran_N 0.000 0.000 2.000 0.000
Tanzania_Kisese 0.000 0.000 0.000 2.000


best coefficients: 0.173 0.669 0.110 0.048
Jackknife mean: 0.172923656 0.668521805 0.110520948 0.048033592
std. errors: 0.015 0.032 0.028 0.010

error covariance (* 1,000,000)
225 -225 69 -69
-225 1056 -841 9
69 -841 812 -40
-69 9 -40 99


summ: Natufian 4 0.698384 0.173 0.669 0.111 0.048 225 -225 69 -69 1056 ...
-841 9 812 -40 99

fixed pat wt dof chisq tail prob
0000 0 9 6.409 0.698384 0.173 0.669 0.110 0.048
0001 1 10 24.432 0.00653235 0.203 0.676 0.122 0.000
0010 1 10 19.974 0.0295034 0.162 0.787 0.000 0.051
0100 1 10 490.746 0 0.302 0.000 0.661 0.036
1000 1 10 128.264 1.05533e-22 0.000 0.947 -0.037 0.090 infeasible
0011 2 11 39.104 5.08558e-05 0.195 0.805 0.000 0.000
0101 2 11 511.532 0 0.328 0.000 0.672 0.000
0110 2 11 2621.128 0 0.884 0.000 0.000 0.116
1001 2 11 186.717 0 0.000 1.061 -0.061 0.000 infeasible
1010 2 11 116.243 1.03094e-19 0.000 0.904 0.000 0.096
1100 2 11 910.480 0 0.000 0.000 0.847 0.153
0111 3 12 2629.130 0 1.000 0.000 0.000 0.000
1011 3 12 190.436 0 0.000 1.000 0.000 0.000
1101 3 12 1166.674 0 0.000 0.000 1.000 0.000
1110 3 12 4833.535 0 0.000 0.000 0.000 1.000
best pat: 0000 0.698384 - -
best pat: 0010 0.0295034 chi(nested): 13.565 p-value for nested model: 0.000230496
best pat: 0011 5.08558e-05 chi(nested): 19.131 p-value for nested model: 1.2207e-05
best pat: 1011 3.05494e-34 chi(nested): 151.331 p-value for nested model: 8.87069e-35

coeffs: 0.173 0.669 0.110 0.048

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

details: Iberomaurusian South_Africa_400BP.SG 0.000199 0.738494
details: Anatolia_Epipaleolithic South_Africa_400BP.SG -0.000308 -0.998514
details: Iran_N South_Africa_400BP.SG -0.000515 -1.896781
details: Tanzania_Kisese South_Africa_400BP.SG -0.002281 -6.450911
dscore: South_Africa_400BP.SG f4: -0.000338 Z: -1.315428

details: Iberomaurusian South_Africa_HG -0.000797 -2.738949
details: Anatolia_Epipaleolithic South_Africa_HG -0.000366 -1.144672
details: Iran_N South_Africa_HG -0.000856 -2.982723
details: Tanzania_Kisese South_Africa_HG 0.003428 9.316969
dscore: South_Africa_HG f4: -0.000312 Z: -1.150735

details: Iberomaurusian Eastern_Euro_UP.SG -0.006876 -17.155451
details: Anatolia_Epipaleolithic Eastern_Euro_UP.SG 0.002878 6.122408
details: Iran_N Eastern_Euro_UP.SG -0.000719 -1.736251
details: Tanzania_Kisese Eastern_Euro_UP.SG -0.023455 -46.269567
dscore: Eastern_Euro_UP.SG f4: -0.000471 Z: -1.205933

details: Iberomaurusian BachoKiro -0.004663 -11.914505
details: Anatolia_Epipaleolithic BachoKiro 0.001808 4.245605
details: Iran_N BachoKiro 0.000957 2.340372
details: Tanzania_Kisese BachoKiro -0.018266 -35.609375
dscore: BachoKiro f4: -0.000370 Z: -1.033257

details: Iberomaurusian ZlatyKun -0.005564 -12.774158
details: Anatolia_Epipaleolithic ZlatyKun 0.001744 3.159573
details: Iran_N ZlatyKun 0.000946 2.021091
details: Tanzania_Kisese ZlatyKun -0.019203 -32.439416
dscore: ZlatyKun f4: -0.000615 Z: -1.369185

details: Iberomaurusian DevilsGate_N.SG -0.005161 -12.755288
details: Anatolia_Epipaleolithic DevilsGate_N.SG 0.002090 4.722907
details: Iran_N DevilsGate_N.SG 0.001390 3.195445
details: Tanzania_Kisese DevilsGate_N.SG -0.020012 -38.637022
dscore: DevilsGate_N.SG f4: -0.000305 Z: -0.819733

details: Iberomaurusian Yana_UP -0.006163 -15.511576
details: Anatolia_Epipaleolithic Yana_UP 0.002570 5.751333
details: Iran_N Yana_UP -0.000030 -0.069252
details: Tanzania_Kisese Yana_UP -0.021989 -42.350806
dscore: Yana_UP f4: -0.000408 Z: -1.097100

details: Iberomaurusian CHG -0.007678 -18.043397
details: Anatolia_Epipaleolithic CHG 0.002827 5.996761
details: Iran_N CHG 0.003774 8.754750
details: Tanzania_Kisese CHG -0.023673 -43.313328
dscore: CHG f4: -0.000160 Z: -0.404856

details: Iberomaurusian Iran_N.SG -0.007464 -20.059089
details: Anatolia_Epipaleolithic Iran_N.SG 0.001786 4.393056
details: Iran_N Iran_N.SG 0.007736 20.830514
details: Tanzania_Kisese Iran_N.SG -0.022635 -46.852167
dscore: Iran_N.SG f4: -0.000334 Z: -0.971357

details: Iberomaurusian British_Mesolithic.SG -0.007741 -17.833344
details: Anatolia_Epipaleolithic British_Mesolithic.SG 0.003873 7.991567
details: Iran_N British_Mesolithic.SG -0.001774 -4.212199
details: Tanzania_Kisese British_Mesolithic.SG -0.025756 -49.085063
dscore: British_Mesolithic.SG f4: -0.000182 Z: -0.450950

details: Iberomaurusian Mota -0.000087 -0.234863
details: Anatolia_Epipaleolithic Mota 0.000065 0.153641
details: Iran_N Mota -0.000329 -0.894946
details: Tanzania_Kisese Mota -0.002374 -4.912219
dscore: Mota f4: -0.000122 Z: -0.344200

details: Iberomaurusian Morocco_EN.SG 0.015836 32.126795
details: Anatolia_Epipaleolithic Morocco_EN.SG -0.002136 -4.293190
details: Iran_N Morocco_EN.SG -0.005499 -11.643657
details: Tanzania_Kisese Morocco_EN.SG -0.020464 -36.906489
dscore: Morocco_EN.SG f4: -0.000282 Z: -0.671605

gendstat: Cameroon_SMA.SG South_Africa_400BP.SG -1.315
gendstat: Cameroon_SMA.SG South_Africa_HG -1.151
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG -1.206
gendstat: Cameroon_SMA.SG BachoKiro -1.033
gendstat: Cameroon_SMA.SG ZlatyKun -1.369
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.820
gendstat: Cameroon_SMA.SG Yana_UP -1.097
gendstat: Cameroon_SMA.SG CHG -0.405
gendstat: Cameroon_SMA.SG Iran_N.SG -0.971
gendstat: Cameroon_SMA.SG British_Mesolithic.SG -0.451
gendstat: Cameroon_SMA.SG Mota -0.344
gendstat: Cameroon_SMA.SG Morocco_EN.SG -0.672
gendstat: South_Africa_400BP.SG South_Africa_HG 0.106
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG -0.344
gendstat: South_Africa_400BP.SG BachoKiro -0.091
gendstat: South_Africa_400BP.SG ZlatyKun -0.644
gendstat: South_Africa_400BP.SG DevilsGate_N.SG 0.092
gendstat: South_Africa_400BP.SG Yana_UP -0.196
gendstat: South_Africa_400BP.SG CHG 0.476
gendstat: South_Africa_400BP.SG Iran_N.SG 0.013
gendstat: South_Africa_400BP.SG British_Mesolithic.SG 0.400
gendstat: South_Africa_400BP.SG Mota 0.609
gendstat: South_Africa_400BP.SG Morocco_EN.SG 0.138
gendstat: South_Africa_HG Eastern_Euro_UP.SG -0.408
gendstat: South_Africa_HG BachoKiro -0.161
gendstat: South_Africa_HG ZlatyKun -0.675
gendstat: South_Africa_HG DevilsGate_N.SG 0.020
gendstat: South_Africa_HG Yana_UP -0.258
gendstat: South_Africa_HG CHG 0.381
gendstat: South_Africa_HG Iran_N.SG -0.063
gendstat: South_Africa_HG British_Mesolithic.SG 0.308
gendstat: South_Africa_HG Mota 0.513
gendstat: South_Africa_HG Morocco_EN.SG 0.067
gendstat: Eastern_Euro_UP.SG BachoKiro 0.255
gendstat: Eastern_Euro_UP.SG ZlatyKun -0.302
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG 0.448
gendstat: Eastern_Euro_UP.SG Yana_UP 0.172
gendstat: Eastern_Euro_UP.SG CHG 0.760
gendstat: Eastern_Euro_UP.SG Iran_N.SG 0.375
gendstat: Eastern_Euro_UP.SG British_Mesolithic.SG 0.754
gendstat: Eastern_Euro_UP.SG Mota 0.813
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG 0.420
gendstat: BachoKiro ZlatyKun -0.513
gendstat: BachoKiro DevilsGate_N.SG 0.184
gendstat: BachoKiro Yana_UP -0.102
gendstat: BachoKiro CHG 0.509
gendstat: BachoKiro Iran_N.SG 0.101
gendstat: BachoKiro British_Mesolithic.SG 0.447
gendstat: BachoKiro Mota 0.585
gendstat: BachoKiro Morocco_EN.SG 0.192
gendstat: ZlatyKun DevilsGate_N.SG 0.662
gendstat: ZlatyKun Yana_UP 0.455
gendstat: ZlatyKun CHG 0.937
gendstat: ZlatyKun Iran_N.SG 0.640
gendstat: ZlatyKun British_Mesolithic.SG 0.928
gendstat: ZlatyKun Mota 1.029
gendstat: ZlatyKun Morocco_EN.SG 0.647
gendstat: DevilsGate_N.SG Yana_UP -0.293
gendstat: DevilsGate_N.SG CHG 0.356
gendstat: DevilsGate_N.SG Iran_N.SG -0.084
gendstat: DevilsGate_N.SG British_Mesolithic.SG 0.310
gendstat: DevilsGate_N.SG Mota 0.416
gendstat: DevilsGate_N.SG Morocco_EN.SG 0.048
gendstat: Yana_UP CHG 0.614
gendstat: Yana_UP Iran_N.SG 0.212
gendstat: Yana_UP British_Mesolithic.SG 0.585
gendstat: Yana_UP Mota 0.659
gendstat: Yana_UP Morocco_EN.SG 0.275
gendstat: CHG Iran_N.SG -0.499
gendstat: CHG British_Mesolithic.SG -0.052
gendstat: CHG Mota 0.085
gendstat: CHG Morocco_EN.SG -0.258
gendstat: Iran_N.SG British_Mesolithic.SG 0.407
gendstat: Iran_N.SG Mota 0.508
gendstat: Iran_N.SG Morocco_EN.SG 0.117
gendstat: British_Mesolithic.SG Mota 0.133
gendstat: British_Mesolithic.SG Morocco_EN.SG -0.212
gendstat: Mota Morocco_EN.SG -0.343

Ahuwarhd
02-14-2022, 05:44 PM
Target: Levant_Natufian:I0861
Distance: 3.0161% / 0.03016148
80.4 Dzudzuana
12.0 Ancestral_North-African(DZU-55%-45%)
7.6 Iran_Neolithic


Target: Levant_Natufian_contam:I1072
Distance: 4.1042% / 0.04104191
86.4 Dzudzuana
13.6 Ancestral_North-African(DZU-55%-45%)

This is similar to what lazaridis estimated natufian composition to be in the dzudzuana study as well. Majority dzudzuana/anatolian HG like with a 20-25% minority of north african HG ancestry to distinguish it.

Gentica277282
02-15-2022, 01:18 PM
48443
Facial reconstruction of natufian

Typic
02-16-2022, 02:42 PM
It seems roughly in line with what I get for them in qpAdm despite all those modern samples having some degree of Natufian admixture.







left pops:
Natufian
Iberomaurusian
Anatolia_Epipaleolithic
Iran_N
Tanzania_Kisese

right pops:
Cameroon_SMA.SG
South_Africa_400BP.SG
South_Africa_HG
Eastern_Euro_UP.SG
BachoKiro
ZlatyKun
DevilsGate_N.SG
Yana_UP
CHG
Iran_N.SG
British_Mesolithic.SG
Mota
Morocco_EN.SG

0 Natufian 6
1 Iberomaurusian 7
2 Anatolia_Epipaleolithic 1
3 Iran_N 4
4 Tanzania_Kisese 1
5 Cameroon_SMA.SG 3
6 South_Africa_400BP.SG 4
7 South_Africa_HG 4
8 Eastern_Euro_UP.SG 3
9 BachoKiro 3
10 ZlatyKun 1
11 DevilsGate_N.SG 4
12 Yana_UP 2
13 CHG 2
14 Iran_N.SG 4
15 British_Mesolithic.SG 3
16 Mota 1
17 Morocco_EN.SG 3
jackknife block size: 0.050
snps: 1150461 indivs: 56
number of blocks for block jackknife: 714
## ncols: 1150461
coverage: Natufian 529354
coverage: Iberomaurusian 1088914
coverage: Anatolia_Epipaleolithic 838782
coverage: Iran_N 1002874
coverage: Tanzania_Kisese 938012
coverage: Cameroon_SMA.SG 1149654
coverage: South_Africa_400BP.SG 1150180
coverage: South_Africa_HG 1135129
coverage: Eastern_Euro_UP.SG 1150149
coverage: BachoKiro 990561
coverage: ZlatyKun 1120502
coverage: DevilsGate_N.SG 1134715
coverage: Yana_UP 1135189
coverage: CHG 1134595
coverage: Iran_N.SG 1149571
coverage: British_Mesolithic.SG 1149936
coverage: Mota 1134612
coverage: Morocco_EN.SG 653097
Effective number of blocks: 616.513
numsnps used: 319495
codimension 1
f4info:
f4rank: 3 dof: 9 chisq: 6.409 tail: 0.698383677 dofdiff: 11 chisqdiff: -6.409 taildiff: 1
B:
scale 1.000 1.000 1.000
South_Africa_400BP.SG 0.088 0.136 -0.078
South_Africa_HG -0.208 -0.165 -0.436
Eastern_Euro_UP.SG 1.221 -0.260 -0.913
BachoKiro 0.969 -0.184 0.055
ZlatyKun 1.011 -0.273 0.008
DevilsGate_N.SG 1.073 -0.234 0.201
Yana_UP 1.151 -0.244 -0.492
CHG 1.325 -0.629 0.802
Iran_N.SG 1.336 -0.834 2.737
British_Mesolithic.SG 1.334 -0.278 -1.591
Mota 0.103 0.076 -0.124
Morocco_EN.SG 0.832 3.239 0.065
A:
scale 105.341 314.313 879.994
Iberomaurusian -0.364 1.828 0.731
Anatolia_Epipaleolithic 0.221 -0.344 -0.507
Iran_N 0.080 -0.587 1.723
Tanzania_Kisese -1.952 -0.444 0.489


full rank
f4info:
f4rank: 4 dof: 0 chisq: 0.000 tail: 1 dofdiff: 9 chisqdiff: 6.409 taildiff: 0.698383677
B:
scale 142.536 460.665 326.101 52.741
South_Africa_400BP.SG 0.028 -0.142 -0.168 -0.120
South_Africa_HG -0.114 -0.169 -0.279 0.181
Eastern_Euro_UP.SG -0.980 1.326 -0.234 -1.237
BachoKiro -0.665 0.833 0.312 -0.963
ZlatyKun -0.793 0.803 0.309 -1.013
DevilsGate_N.SG -0.736 0.963 0.453 -1.055
Yana_UP -0.879 1.184 -0.010 -1.160
CHG -1.094 1.302 1.231 -1.249
Iran_N.SG -1.064 0.823 2.523 -1.194
British_Mesolithic.SG -1.103 1.784 -0.579 -1.358
Mota -0.012 0.030 -0.107 -0.125
Morocco_EN.SG 2.257 -0.984 -1.793 -1.079
A:
scale 2.000 2.000 2.000 2.000
Iberomaurusian 2.000 0.000 0.000 0.000
Anatolia_Epipaleolithic 0.000 2.000 0.000 0.000
Iran_N 0.000 0.000 2.000 0.000
Tanzania_Kisese 0.000 0.000 0.000 2.000


best coefficients: 0.173 0.669 0.110 0.048
Jackknife mean: 0.172923656 0.668521805 0.110520948 0.048033592
std. errors: 0.015 0.032 0.028 0.010

error covariance (* 1,000,000)
225 -225 69 -69
-225 1056 -841 9
69 -841 812 -40
-69 9 -40 99


summ: Natufian 4 0.698384 0.173 0.669 0.111 0.048 225 -225 69 -69 1056 ...
-841 9 812 -40 99

fixed pat wt dof chisq tail prob
0000 0 9 6.409 0.698384 0.173 0.669 0.110 0.048
0001 1 10 24.432 0.00653235 0.203 0.676 0.122 0.000
0010 1 10 19.974 0.0295034 0.162 0.787 0.000 0.051
0100 1 10 490.746 0 0.302 0.000 0.661 0.036
1000 1 10 128.264 1.05533e-22 0.000 0.947 -0.037 0.090 infeasible
0011 2 11 39.104 5.08558e-05 0.195 0.805 0.000 0.000
0101 2 11 511.532 0 0.328 0.000 0.672 0.000
0110 2 11 2621.128 0 0.884 0.000 0.000 0.116
1001 2 11 186.717 0 0.000 1.061 -0.061 0.000 infeasible
1010 2 11 116.243 1.03094e-19 0.000 0.904 0.000 0.096
1100 2 11 910.480 0 0.000 0.000 0.847 0.153
0111 3 12 2629.130 0 1.000 0.000 0.000 0.000
1011 3 12 190.436 0 0.000 1.000 0.000 0.000
1101 3 12 1166.674 0 0.000 0.000 1.000 0.000
1110 3 12 4833.535 0 0.000 0.000 0.000 1.000
best pat: 0000 0.698384 - -
best pat: 0010 0.0295034 chi(nested): 13.565 p-value for nested model: 0.000230496
best pat: 0011 5.08558e-05 chi(nested): 19.131 p-value for nested model: 1.2207e-05
best pat: 1011 3.05494e-34 chi(nested): 151.331 p-value for nested model: 8.87069e-35

coeffs: 0.173 0.669 0.110 0.048

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

details: Iberomaurusian South_Africa_400BP.SG 0.000199 0.738494
details: Anatolia_Epipaleolithic South_Africa_400BP.SG -0.000308 -0.998514
details: Iran_N South_Africa_400BP.SG -0.000515 -1.896781
details: Tanzania_Kisese South_Africa_400BP.SG -0.002281 -6.450911
dscore: South_Africa_400BP.SG f4: -0.000338 Z: -1.315428

details: Iberomaurusian South_Africa_HG -0.000797 -2.738949
details: Anatolia_Epipaleolithic South_Africa_HG -0.000366 -1.144672
details: Iran_N South_Africa_HG -0.000856 -2.982723
details: Tanzania_Kisese South_Africa_HG 0.003428 9.316969
dscore: South_Africa_HG f4: -0.000312 Z: -1.150735

details: Iberomaurusian Eastern_Euro_UP.SG -0.006876 -17.155451
details: Anatolia_Epipaleolithic Eastern_Euro_UP.SG 0.002878 6.122408
details: Iran_N Eastern_Euro_UP.SG -0.000719 -1.736251
details: Tanzania_Kisese Eastern_Euro_UP.SG -0.023455 -46.269567
dscore: Eastern_Euro_UP.SG f4: -0.000471 Z: -1.205933

details: Iberomaurusian BachoKiro -0.004663 -11.914505
details: Anatolia_Epipaleolithic BachoKiro 0.001808 4.245605
details: Iran_N BachoKiro 0.000957 2.340372
details: Tanzania_Kisese BachoKiro -0.018266 -35.609375
dscore: BachoKiro f4: -0.000370 Z: -1.033257

details: Iberomaurusian ZlatyKun -0.005564 -12.774158
details: Anatolia_Epipaleolithic ZlatyKun 0.001744 3.159573
details: Iran_N ZlatyKun 0.000946 2.021091
details: Tanzania_Kisese ZlatyKun -0.019203 -32.439416
dscore: ZlatyKun f4: -0.000615 Z: -1.369185

details: Iberomaurusian DevilsGate_N.SG -0.005161 -12.755288
details: Anatolia_Epipaleolithic DevilsGate_N.SG 0.002090 4.722907
details: Iran_N DevilsGate_N.SG 0.001390 3.195445
details: Tanzania_Kisese DevilsGate_N.SG -0.020012 -38.637022
dscore: DevilsGate_N.SG f4: -0.000305 Z: -0.819733

details: Iberomaurusian Yana_UP -0.006163 -15.511576
details: Anatolia_Epipaleolithic Yana_UP 0.002570 5.751333
details: Iran_N Yana_UP -0.000030 -0.069252
details: Tanzania_Kisese Yana_UP -0.021989 -42.350806
dscore: Yana_UP f4: -0.000408 Z: -1.097100

details: Iberomaurusian CHG -0.007678 -18.043397
details: Anatolia_Epipaleolithic CHG 0.002827 5.996761
details: Iran_N CHG 0.003774 8.754750
details: Tanzania_Kisese CHG -0.023673 -43.313328
dscore: CHG f4: -0.000160 Z: -0.404856

details: Iberomaurusian Iran_N.SG -0.007464 -20.059089
details: Anatolia_Epipaleolithic Iran_N.SG 0.001786 4.393056
details: Iran_N Iran_N.SG 0.007736 20.830514
details: Tanzania_Kisese Iran_N.SG -0.022635 -46.852167
dscore: Iran_N.SG f4: -0.000334 Z: -0.971357

details: Iberomaurusian British_Mesolithic.SG -0.007741 -17.833344
details: Anatolia_Epipaleolithic British_Mesolithic.SG 0.003873 7.991567
details: Iran_N British_Mesolithic.SG -0.001774 -4.212199
details: Tanzania_Kisese British_Mesolithic.SG -0.025756 -49.085063
dscore: British_Mesolithic.SG f4: -0.000182 Z: -0.450950

details: Iberomaurusian Mota -0.000087 -0.234863
details: Anatolia_Epipaleolithic Mota 0.000065 0.153641
details: Iran_N Mota -0.000329 -0.894946
details: Tanzania_Kisese Mota -0.002374 -4.912219
dscore: Mota f4: -0.000122 Z: -0.344200

details: Iberomaurusian Morocco_EN.SG 0.015836 32.126795
details: Anatolia_Epipaleolithic Morocco_EN.SG -0.002136 -4.293190
details: Iran_N Morocco_EN.SG -0.005499 -11.643657
details: Tanzania_Kisese Morocco_EN.SG -0.020464 -36.906489
dscore: Morocco_EN.SG f4: -0.000282 Z: -0.671605

gendstat: Cameroon_SMA.SG South_Africa_400BP.SG -1.315
gendstat: Cameroon_SMA.SG South_Africa_HG -1.151
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG -1.206
gendstat: Cameroon_SMA.SG BachoKiro -1.033
gendstat: Cameroon_SMA.SG ZlatyKun -1.369
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.820
gendstat: Cameroon_SMA.SG Yana_UP -1.097
gendstat: Cameroon_SMA.SG CHG -0.405
gendstat: Cameroon_SMA.SG Iran_N.SG -0.971
gendstat: Cameroon_SMA.SG British_Mesolithic.SG -0.451
gendstat: Cameroon_SMA.SG Mota -0.344
gendstat: Cameroon_SMA.SG Morocco_EN.SG -0.672
gendstat: South_Africa_400BP.SG South_Africa_HG 0.106
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG -0.344
gendstat: South_Africa_400BP.SG BachoKiro -0.091
gendstat: South_Africa_400BP.SG ZlatyKun -0.644
gendstat: South_Africa_400BP.SG DevilsGate_N.SG 0.092
gendstat: South_Africa_400BP.SG Yana_UP -0.196
gendstat: South_Africa_400BP.SG CHG 0.476
gendstat: South_Africa_400BP.SG Iran_N.SG 0.013
gendstat: South_Africa_400BP.SG British_Mesolithic.SG 0.400
gendstat: South_Africa_400BP.SG Mota 0.609
gendstat: South_Africa_400BP.SG Morocco_EN.SG 0.138
gendstat: South_Africa_HG Eastern_Euro_UP.SG -0.408
gendstat: South_Africa_HG BachoKiro -0.161
gendstat: South_Africa_HG ZlatyKun -0.675
gendstat: South_Africa_HG DevilsGate_N.SG 0.020
gendstat: South_Africa_HG Yana_UP -0.258
gendstat: South_Africa_HG CHG 0.381
gendstat: South_Africa_HG Iran_N.SG -0.063
gendstat: South_Africa_HG British_Mesolithic.SG 0.308
gendstat: South_Africa_HG Mota 0.513
gendstat: South_Africa_HG Morocco_EN.SG 0.067
gendstat: Eastern_Euro_UP.SG BachoKiro 0.255
gendstat: Eastern_Euro_UP.SG ZlatyKun -0.302
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG 0.448
gendstat: Eastern_Euro_UP.SG Yana_UP 0.172
gendstat: Eastern_Euro_UP.SG CHG 0.760
gendstat: Eastern_Euro_UP.SG Iran_N.SG 0.375
gendstat: Eastern_Euro_UP.SG British_Mesolithic.SG 0.754
gendstat: Eastern_Euro_UP.SG Mota 0.813
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG 0.420
gendstat: BachoKiro ZlatyKun -0.513
gendstat: BachoKiro DevilsGate_N.SG 0.184
gendstat: BachoKiro Yana_UP -0.102
gendstat: BachoKiro CHG 0.509
gendstat: BachoKiro Iran_N.SG 0.101
gendstat: BachoKiro British_Mesolithic.SG 0.447
gendstat: BachoKiro Mota 0.585
gendstat: BachoKiro Morocco_EN.SG 0.192
gendstat: ZlatyKun DevilsGate_N.SG 0.662
gendstat: ZlatyKun Yana_UP 0.455
gendstat: ZlatyKun CHG 0.937
gendstat: ZlatyKun Iran_N.SG 0.640
gendstat: ZlatyKun British_Mesolithic.SG 0.928
gendstat: ZlatyKun Mota 1.029
gendstat: ZlatyKun Morocco_EN.SG 0.647
gendstat: DevilsGate_N.SG Yana_UP -0.293
gendstat: DevilsGate_N.SG CHG 0.356
gendstat: DevilsGate_N.SG Iran_N.SG -0.084
gendstat: DevilsGate_N.SG British_Mesolithic.SG 0.310
gendstat: DevilsGate_N.SG Mota 0.416
gendstat: DevilsGate_N.SG Morocco_EN.SG 0.048
gendstat: Yana_UP CHG 0.614
gendstat: Yana_UP Iran_N.SG 0.212
gendstat: Yana_UP British_Mesolithic.SG 0.585
gendstat: Yana_UP Mota 0.659
gendstat: Yana_UP Morocco_EN.SG 0.275
gendstat: CHG Iran_N.SG -0.499
gendstat: CHG British_Mesolithic.SG -0.052
gendstat: CHG Mota 0.085
gendstat: CHG Morocco_EN.SG -0.258
gendstat: Iran_N.SG British_Mesolithic.SG 0.407
gendstat: Iran_N.SG Mota 0.508
gendstat: Iran_N.SG Morocco_EN.SG 0.117
gendstat: British_Mesolithic.SG Mota 0.133
gendstat: British_Mesolithic.SG Morocco_EN.SG -0.212
gendstat: Mota Morocco_EN.SG -0.343






If it is not asking too much, could you please help me understand this?

Mnemonics
02-16-2022, 04:32 PM
If it is not asking too much, could you please help me understand this?

They get roughly 17.3% Iberomaurusian, 66.9% Anatolia_Epipaleolithic, 11.0% Iran_N related ancestry, and 4.8% Tanzania_Kisese ( a more Southern shifted Mota-like HG ).

Gentica277282
02-16-2022, 06:23 PM
They get roughly 17.3% Iberomaurusian, 66.9 Anatolia_Epipaleolithic, 11.0% Iran_N related ancestry, and 4.8% Tanzania_Kisese ( a more Southern shifted Mota-like HG ).

Makes sense as to why they plot away from mahra

Gentica277282
02-16-2022, 10:50 PM
They get roughly 17.3% Iberomaurusian, 66.9% Anatolia_Epipaleolithic, 11.0% Iran_N related ancestry, and 4.8% Tanzania_Kisese ( a more Southern shifted Mota-like HG ).

Could you run Arabs and see if they score any ana please I was reading a paper the other day and saw something interesting about Arabs.

The ADMIXTURE analysis identified K = 12 ancestral populations as having the lowest cross-validation error (Supplemental Fig. 7A). At this level of resolution, the Q1 (Bedouin) had a high average (84%) proportion of ancestry that was also present in the Human Origins Bedouin B population at a high average proportion (93%) (Supplemental Fig. 7B,C), in which this same ancestry was also shared with Saudis, and at lower levels among other Middle Eastern populations. This ancestry therefore appears to be the signal of an indigenous Arab ancestral population.

Also I noticed when modelling a Saudi who's eef around 20 percent and comparing distances rich he ended up plotting further away from us then normal Arabs so I know the lack of eef ancestry is not an issue.

Mnemonics
02-17-2022, 12:10 AM
Could you run Arabs and see if they score any ana please I was reading a paper the other day and saw something interesting about Arabs.

The ADMIXTURE analysis identified K = 12 ancestral populations as having the lowest cross-validation error (Supplemental Fig. 7A). At this level of resolution, the Q1 (Bedouin) had a high average (84%) proportion of ancestry that was also present in the Human Origins Bedouin B population at a high average proportion (93%) (Supplemental Fig. 7B,C), in which this same ancestry was also shared with Saudis, and at lower levels among other Middle Eastern populations. This ancestry therefore appears to be the signal of an indigenous Arab ancestral population.

Also I noticed when modelling a Saudi who's eef rich around 20 percent and comparing distances rich he ended up plotting further away from us then normal Arabs so I know the lack of eef ancestry is not an issue.



The BedouinB diploid samples seem to mostly Jordan_EBA-like with extra Kenya_PastoralN. They also seem to have some extra Yana_UP affinity I can't model.





left pops:
BedouinB.DG
Levant_BA
Kenya_PastoralN

right pops:
South_Africa_400BP.SG
Cameroon_SMA.SG
Mota
South_Africa_HG
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 BedouinB.DG 2
1 Levant_BA 3
2 Kenya_PastoralN 8
3 South_Africa_400BP.SG 4
4 Cameroon_SMA.SG 3
5 Mota 1
6 South_Africa_HG 4
7 ZlatyKun 1
8 BachoKiro 3
9 Eastern_Euro_UP.SG 3
10 DevilsGate_N.SG 4
11 Kazakhstan_Botai_Eneolithic.SG 3
12 Barcin_N.SG 3
13 Iran_N.SG 4
14 MA1 1
15 Morocco_EN.SG 3
16 CHG 2
17 Italy_Mesolithic.SG 3
18 Romania_IronGates_Mesolithic.SG 3
19 Russia_BA_Okunevo.SG 14
20 Germany_EarlyMedieval.SG 21
21 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150471 indivs: 106
number of blocks for block jackknife: 714
## ncols: 1150471
coverage: BedouinB.DG 1121254
coverage: Levant_BA 870086
coverage: Kenya_PastoralN 1061629
coverage: South_Africa_400BP.SG 1150180
coverage: Cameroon_SMA.SG 1149654
coverage: Mota 1134612
coverage: South_Africa_HG 1135129
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 623.267
numsnps used: 472322
codimension 1
f4info:
f4rank: 1 dof: 17 chisq: 17.995 tail: 0.389136695 dofdiff: 19 chisqdiff: -17.995 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG -0.021
Mota 0.039
South_Africa_HG -0.077
ZlatyKun 0.869
BachoKiro 0.840
Eastern_Euro_UP.SG 1.084
DevilsGate_N.SG 0.953
Kazakhstan_Botai_Eneolithic.SG 1.106
Barcin_N.SG 1.289
Iran_N.SG 1.159
MA1 1.106
Morocco_EN.SG 0.631
CHG 1.220
Italy_Mesolithic.SG 1.177
Romania_IronGates_Mesolithic.SG 1.184
Russia_BA_Okunevo.SG 1.090
Germany_EarlyMedieval.SG 1.271
Poland_Southeast_CordedWare.SG 1.232
A:
scale 132.284
Levant_BA 0.130
Kenya_PastoralN -1.408


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 17 chisqdiff: 17.995 taildiff: 0.389136695
B:
scale 1005.459 92.843
Cameroon_SMA.SG -0.003 0.011
Mota 0.161 -0.042
South_Africa_HG 0.097 0.070
ZlatyKun 0.671 -0.882
BachoKiro 0.967 -0.838
Eastern_Euro_UP.SG 1.320 -1.078
DevilsGate_N.SG 0.768 -0.965
Kazakhstan_Botai_Eneolithic.SG 0.905 -1.120
Barcin_N.SG 1.445 -1.274
Iran_N.SG 1.380 -1.143
MA1 0.953 -1.118
Morocco_EN.SG 0.758 -0.633
CHG 1.487 -1.202
Italy_Mesolithic.SG 1.216 -1.175
Romania_IronGates_Mesolithic.SG 1.047 -1.190
Russia_BA_Okunevo.SG 0.874 -1.101
Germany_EarlyMedieval.SG 1.164 -1.267
Poland_Southeast_CordedWare.SG 1.021 -1.237
A:
scale 1.414 1.414
Levant_BA 1.414 0.000
Kenya_PastoralN 0.000 1.414


best coefficients: 0.915 0.085
Jackknife mean: 0.915340507 0.084659493
std. errors: 0.012 0.012

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


summ: BedouinB.DG 2 0.389137 0.915 0.085 152 -152 152

fixed pat wt dof chisq tail prob
00 0 17 17.995 0.389137 0.915 0.085
01 1 18 59.603 2.37107e-06 1.000 0.000
10 1 18 4315.882 0 0.000 1.000
best pat: 00 0.389137 - -
best pat: 01 2.37107e-06 chi(nested): 41.607 p-value for nested model: 1.11567e-10

coeffs: 0.915 0.085

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

details: Levant_BA Cameroon_SMA.SG -0.000003 -0.019341
details: Kenya_PastoralN Cameroon_SMA.SG 0.000119 0.732389
dscore: Cameroon_SMA.SG f4: 0.000007 Z: 0.041879

details: Levant_BA Mota 0.000160 0.715608
details: Kenya_PastoralN Mota -0.000454 -2.095219
dscore: Mota f4: 0.000108 Z: 0.505451

details: Levant_BA South_Africa_HG 0.000097 0.566199
details: Kenya_PastoralN South_Africa_HG 0.000756 4.758071
dscore: South_Africa_HG f4: 0.000153 Z: 0.935890

details: Levant_BA ZlatyKun 0.000668 2.221524
details: Kenya_PastoralN ZlatyKun -0.009497 -32.148828
dscore: ZlatyKun f4: -0.000192 Z: -0.665088

details: Levant_BA BachoKiro 0.000962 4.047038
details: Kenya_PastoralN BachoKiro -0.009024 -36.335295
dscore: BachoKiro f4: 0.000117 Z: 0.516532

details: Levant_BA Eastern_Euro_UP.SG 0.001313 5.106083
details: Kenya_PastoralN Eastern_Euro_UP.SG -0.011606 -42.370973
dscore: Eastern_Euro_UP.SG f4: 0.000220 Z: 0.891249

details: Levant_BA DevilsGate_N.SG 0.000764 3.024543
details: Kenya_PastoralN DevilsGate_N.SG -0.010391 -39.774104
dscore: DevilsGate_N.SG f4: -0.000179 Z: -0.739968

details: Levant_BA Kazakhstan_Botai_Eneolithic.SG 0.000900 3.636812
details: Kenya_PastoralN Kazakhstan_Botai_Eneolithic.SG -0.012063 -46.286225
dscore: Kazakhstan_Botai_Eneolithic.SG f4: -0.000196 Z: -0.824513

details: Levant_BA Barcin_N.SG 0.001437 5.586478
details: Kenya_PastoralN Barcin_N.SG -0.013724 -51.492202
dscore: Barcin_N.SG f4: 0.000155 Z: 0.626290

details: Levant_BA Iran_N.SG 0.001372 5.838001
details: Kenya_PastoralN Iran_N.SG -0.012308 -50.079052
dscore: Iran_N.SG f4: 0.000215 Z: 0.952358

details: Levant_BA MA1 0.000948 2.979787
details: Kenya_PastoralN MA1 -0.012040 -38.431862
dscore: MA1 f4: -0.000151 Z: -0.494533

details: Levant_BA Morocco_EN.SG 0.000754 2.933128
details: Kenya_PastoralN Morocco_EN.SG -0.006819 -25.142524
dscore: Morocco_EN.SG f4: 0.000113 Z: 0.458658

details: Levant_BA CHG 0.001479 5.303404
details: Kenya_PastoralN CHG -0.012943 -47.054098
dscore: CHG f4: 0.000259 Z: 0.967891

details: Levant_BA Italy_Mesolithic.SG 0.001210 4.474158
details: Kenya_PastoralN Italy_Mesolithic.SG -0.012655 -45.198473
dscore: Italy_Mesolithic.SG f4: 0.000037 Z: 0.142190

details: Levant_BA Romania_IronGates_Mesolithic.SG 0.001041 3.944235
details: Kenya_PastoralN Romania_IronGates_Mesolithic.SG -0.012821 -47.617020
dscore: Romania_IronGates_Mesolithic.SG f4: -0.000131 Z: -0.519257

details: Levant_BA Russia_BA_Okunevo.SG 0.000869 3.896733
details: Kenya_PastoralN Russia_BA_Okunevo.SG -0.011853 -49.602181
dscore: Russia_BA_Okunevo.SG f4: -0.000207 Z: -0.965882

details: Levant_BA Germany_EarlyMedieval.SG 0.001158 5.105662
details: Kenya_PastoralN Germany_EarlyMedieval.SG -0.013643 -57.669211
dscore: Germany_EarlyMedieval.SG f4: -0.000094 Z: -0.432724

details: Levant_BA Poland_Southeast_CordedWare.SG 0.001016 4.398323
details: Kenya_PastoralN Poland_Southeast_CordedWare.SG -0.013326 -54.627733
dscore: Poland_Southeast_CordedWare.SG f4: -0.000197 Z: -0.888079

gendstat: South_Africa_400BP.SG Cameroon_SMA.SG 0.042
gendstat: South_Africa_400BP.SG Mota 0.505
gendstat: South_Africa_400BP.SG South_Africa_HG 0.936
gendstat: South_Africa_400BP.SG ZlatyKun -0.665
gendstat: South_Africa_400BP.SG BachoKiro 0.517
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG 0.891
gendstat: South_Africa_400BP.SG DevilsGate_N.SG -0.740
gendstat: South_Africa_400BP.SG Kazakhstan_Botai_Eneolithic.SG -0.825
gendstat: South_Africa_400BP.SG Barcin_N.SG 0.626
gendstat: South_Africa_400BP.SG Iran_N.SG 0.952
gendstat: South_Africa_400BP.SG MA1 -0.495
gendstat: South_Africa_400BP.SG Morocco_EN.SG 0.459
gendstat: South_Africa_400BP.SG CHG 0.968
gendstat: South_Africa_400BP.SG Italy_Mesolithic.SG 0.142
gendstat: South_Africa_400BP.SG Romania_IronGates_Mesolithic.SG -0.519
gendstat: South_Africa_400BP.SG Russia_BA_Okunevo.SG -0.966
gendstat: South_Africa_400BP.SG Germany_EarlyMedieval.SG -0.433
gendstat: South_Africa_400BP.SG Poland_Southeast_CordedWare.SG -0.888
gendstat: Cameroon_SMA.SG Mota 0.430
gendstat: Cameroon_SMA.SG South_Africa_HG 0.841
gendstat: Cameroon_SMA.SG ZlatyKun -0.649
gendstat: Cameroon_SMA.SG BachoKiro 0.427
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG 0.797
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.700
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.772
gendstat: Cameroon_SMA.SG Barcin_N.SG 0.558
gendstat: Cameroon_SMA.SG Iran_N.SG 0.826
gendstat: Cameroon_SMA.SG MA1 -0.503
gendstat: Cameroon_SMA.SG Morocco_EN.SG 0.388
gendstat: Cameroon_SMA.SG CHG 0.870
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG 0.108
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -0.513
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -0.899
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -0.420
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -0.830
gendstat: Mota South_Africa_HG 0.194
gendstat: Mota ZlatyKun -0.898
gendstat: Mota BachoKiro 0.034
gendstat: Mota Eastern_Euro_UP.SG 0.371
gendstat: Mota DevilsGate_N.SG -0.946
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG -1.051
gendstat: Mota Barcin_N.SG 0.157
gendstat: Mota Iran_N.SG 0.389
gendstat: Mota MA1 -0.752
gendstat: Mota Morocco_EN.SG 0.018
gendstat: Mota CHG 0.482
gendstat: Mota Italy_Mesolithic.SG -0.226
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.769
gendstat: Mota Russia_BA_Okunevo.SG -1.193
gendstat: Mota Germany_EarlyMedieval.SG -0.747
gendstat: Mota Poland_Southeast_CordedWare.SG -1.105
gendstat: South_Africa_HG ZlatyKun -1.125
gendstat: South_Africa_HG BachoKiro -0.143
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.254
gendstat: South_Africa_HG DevilsGate_N.SG -1.271
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG -1.335
gendstat: South_Africa_HG Barcin_N.SG 0.009
gendstat: South_Africa_HG Iran_N.SG 0.257
gendstat: South_Africa_HG MA1 -0.961
gendstat: South_Africa_HG Morocco_EN.SG -0.144
gendstat: South_Africa_HG CHG 0.371
gendstat: South_Africa_HG Italy_Mesolithic.SG -0.423
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -1.054
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -1.559
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -1.042
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG -1.455
gendstat: ZlatyKun BachoKiro 1.055
gendstat: ZlatyKun Eastern_Euro_UP.SG 1.364
gendstat: ZlatyKun DevilsGate_N.SG 0.042
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG -0.014
gendstat: ZlatyKun Barcin_N.SG 1.130
gendstat: ZlatyKun Iran_N.SG 1.437
gendstat: ZlatyKun MA1 0.118
gendstat: ZlatyKun Morocco_EN.SG 0.896
gendstat: ZlatyKun CHG 1.427
gendstat: ZlatyKun Italy_Mesolithic.SG 0.703
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG 0.195
gendstat: ZlatyKun Russia_BA_Okunevo.SG -0.056
gendstat: ZlatyKun Germany_EarlyMedieval.SG 0.361
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.017
gendstat: BachoKiro Eastern_Euro_UP.SG 0.419
gendstat: BachoKiro DevilsGate_N.SG -1.219
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -1.290
gendstat: BachoKiro Barcin_N.SG 0.145
gendstat: BachoKiro Iran_N.SG 0.421
gendstat: BachoKiro MA1 -0.867
gendstat: BachoKiro Morocco_EN.SG -0.014
gendstat: BachoKiro CHG 0.532
gendstat: BachoKiro Italy_Mesolithic.SG -0.301
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -0.940
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.536
gendstat: BachoKiro Germany_EarlyMedieval.SG -0.933
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.349
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -1.572
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -1.788
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -0.272
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.020
gendstat: Eastern_Euro_UP.SG MA1 -1.302
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -0.390
gendstat: Eastern_Euro_UP.SG CHG 0.145
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG -0.763
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -1.468
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -2.088
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -1.469
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -1.862
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG -0.070
gendstat: DevilsGate_N.SG Barcin_N.SG 1.282
gendstat: DevilsGate_N.SG Iran_N.SG 1.627
gendstat: DevilsGate_N.SG MA1 0.092
gendstat: DevilsGate_N.SG Morocco_EN.SG 1.012
gendstat: DevilsGate_N.SG CHG 1.618
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.792
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG 0.183
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.136
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG 0.382
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG -0.075
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG 1.457
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 1.832
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 0.174
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG 1.107
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG 1.740
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.951
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG 0.298
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -0.066
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG 0.529
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG -0.004
gendstat: Barcin_N.SG Iran_N.SG 0.261
gendstat: Barcin_N.SG MA1 -1.017
gendstat: Barcin_N.SG Morocco_EN.SG -0.151
gendstat: Barcin_N.SG CHG 0.416
gendstat: Barcin_N.SG Italy_Mesolithic.SG -0.491
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG -1.248
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG -1.694
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG -1.338
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG -1.884
gendstat: Iran_N.SG MA1 -1.295
gendstat: Iran_N.SG Morocco_EN.SG -0.358
gendstat: Iran_N.SG CHG 0.180
gendstat: Iran_N.SG Italy_Mesolithic.SG -0.701
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -1.428
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -2.142
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -1.588
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -2.038
gendstat: MA1 Morocco_EN.SG 0.783
gendstat: MA1 CHG 1.311
gendstat: MA1 Italy_Mesolithic.SG 0.635
gendstat: MA1 Romania_IronGates_Mesolithic.SG 0.070
gendstat: MA1 Russia_BA_Okunevo.SG -0.223
gendstat: MA1 Germany_EarlyMedieval.SG 0.215
gendstat: MA1 Poland_Southeast_CordedWare.SG -0.175
gendstat: Morocco_EN.SG CHG 0.465
gendstat: Morocco_EN.SG Italy_Mesolithic.SG -0.267
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG -0.858
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG -1.242
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG -0.817
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG -1.159
gendstat: CHG Italy_Mesolithic.SG -0.805
gendstat: CHG Romania_IronGates_Mesolithic.SG -1.492
gendstat: CHG Russia_BA_Okunevo.SG -1.981
gendstat: CHG Germany_EarlyMedieval.SG -1.590
gendstat: CHG Poland_Southeast_CordedWare.SG -1.955
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -0.817
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -1.124
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -0.643
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -1.088
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.376
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG 0.193
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.340
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG 0.713
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG 0.061
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG -0.666

venustas
02-17-2022, 03:21 AM
The BedouinB diploid samples seem to mostly Jordan_EBA-like with extra Kenya_PastoralN. They also seem to have some extra Yana_UP affinity I can't model.





left pops:
BedouinB.DG
Levant_BA
Kenya_PastoralN

right pops:
South_Africa_400BP.SG
Cameroon_SMA.SG
Mota
South_Africa_HG
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 BedouinB.DG 2
1 Levant_BA 3
2 Kenya_PastoralN 8
3 South_Africa_400BP.SG 4
4 Cameroon_SMA.SG 3
5 Mota 1
6 South_Africa_HG 4
7 ZlatyKun 1
8 BachoKiro 3
9 Eastern_Euro_UP.SG 3
10 DevilsGate_N.SG 4
11 Kazakhstan_Botai_Eneolithic.SG 3
12 Barcin_N.SG 3
13 Iran_N.SG 4
14 MA1 1
15 Morocco_EN.SG 3
16 CHG 2
17 Italy_Mesolithic.SG 3
18 Romania_IronGates_Mesolithic.SG 3
19 Russia_BA_Okunevo.SG 14
20 Germany_EarlyMedieval.SG 21
21 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150471 indivs: 106
number of blocks for block jackknife: 714
## ncols: 1150471
coverage: BedouinB.DG 1121254
coverage: Levant_BA 870086
coverage: Kenya_PastoralN 1061629
coverage: South_Africa_400BP.SG 1150180
coverage: Cameroon_SMA.SG 1149654
coverage: Mota 1134612
coverage: South_Africa_HG 1135129
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 623.267
numsnps used: 472322
codimension 1
f4info:
f4rank: 1 dof: 17 chisq: 17.995 tail: 0.389136695 dofdiff: 19 chisqdiff: -17.995 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG -0.021
Mota 0.039
South_Africa_HG -0.077
ZlatyKun 0.869
BachoKiro 0.840
Eastern_Euro_UP.SG 1.084
DevilsGate_N.SG 0.953
Kazakhstan_Botai_Eneolithic.SG 1.106
Barcin_N.SG 1.289
Iran_N.SG 1.159
MA1 1.106
Morocco_EN.SG 0.631
CHG 1.220
Italy_Mesolithic.SG 1.177
Romania_IronGates_Mesolithic.SG 1.184
Russia_BA_Okunevo.SG 1.090
Germany_EarlyMedieval.SG 1.271
Poland_Southeast_CordedWare.SG 1.232
A:
scale 132.284
Levant_BA 0.130
Kenya_PastoralN -1.408


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 17 chisqdiff: 17.995 taildiff: 0.389136695
B:
scale 1005.459 92.843
Cameroon_SMA.SG -0.003 0.011
Mota 0.161 -0.042
South_Africa_HG 0.097 0.070
ZlatyKun 0.671 -0.882
BachoKiro 0.967 -0.838
Eastern_Euro_UP.SG 1.320 -1.078
DevilsGate_N.SG 0.768 -0.965
Kazakhstan_Botai_Eneolithic.SG 0.905 -1.120
Barcin_N.SG 1.445 -1.274
Iran_N.SG 1.380 -1.143
MA1 0.953 -1.118
Morocco_EN.SG 0.758 -0.633
CHG 1.487 -1.202
Italy_Mesolithic.SG 1.216 -1.175
Romania_IronGates_Mesolithic.SG 1.047 -1.190
Russia_BA_Okunevo.SG 0.874 -1.101
Germany_EarlyMedieval.SG 1.164 -1.267
Poland_Southeast_CordedWare.SG 1.021 -1.237
A:
scale 1.414 1.414
Levant_BA 1.414 0.000
Kenya_PastoralN 0.000 1.414


best coefficients: 0.915 0.085
Jackknife mean: 0.915340507 0.084659493
std. errors: 0.012 0.012

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


summ: BedouinB.DG 2 0.389137 0.915 0.085 152 -152 152

fixed pat wt dof chisq tail prob
00 0 17 17.995 0.389137 0.915 0.085
01 1 18 59.603 2.37107e-06 1.000 0.000
10 1 18 4315.882 0 0.000 1.000
best pat: 00 0.389137 - -
best pat: 01 2.37107e-06 chi(nested): 41.607 p-value for nested model: 1.11567e-10

coeffs: 0.915 0.085

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

details: Levant_BA Cameroon_SMA.SG -0.000003 -0.019341
details: Kenya_PastoralN Cameroon_SMA.SG 0.000119 0.732389
dscore: Cameroon_SMA.SG f4: 0.000007 Z: 0.041879

details: Levant_BA Mota 0.000160 0.715608
details: Kenya_PastoralN Mota -0.000454 -2.095219
dscore: Mota f4: 0.000108 Z: 0.505451

details: Levant_BA South_Africa_HG 0.000097 0.566199
details: Kenya_PastoralN South_Africa_HG 0.000756 4.758071
dscore: South_Africa_HG f4: 0.000153 Z: 0.935890

details: Levant_BA ZlatyKun 0.000668 2.221524
details: Kenya_PastoralN ZlatyKun -0.009497 -32.148828
dscore: ZlatyKun f4: -0.000192 Z: -0.665088

details: Levant_BA BachoKiro 0.000962 4.047038
details: Kenya_PastoralN BachoKiro -0.009024 -36.335295
dscore: BachoKiro f4: 0.000117 Z: 0.516532

details: Levant_BA Eastern_Euro_UP.SG 0.001313 5.106083
details: Kenya_PastoralN Eastern_Euro_UP.SG -0.011606 -42.370973
dscore: Eastern_Euro_UP.SG f4: 0.000220 Z: 0.891249

details: Levant_BA DevilsGate_N.SG 0.000764 3.024543
details: Kenya_PastoralN DevilsGate_N.SG -0.010391 -39.774104
dscore: DevilsGate_N.SG f4: -0.000179 Z: -0.739968

details: Levant_BA Kazakhstan_Botai_Eneolithic.SG 0.000900 3.636812
details: Kenya_PastoralN Kazakhstan_Botai_Eneolithic.SG -0.012063 -46.286225
dscore: Kazakhstan_Botai_Eneolithic.SG f4: -0.000196 Z: -0.824513

details: Levant_BA Barcin_N.SG 0.001437 5.586478
details: Kenya_PastoralN Barcin_N.SG -0.013724 -51.492202
dscore: Barcin_N.SG f4: 0.000155 Z: 0.626290

details: Levant_BA Iran_N.SG 0.001372 5.838001
details: Kenya_PastoralN Iran_N.SG -0.012308 -50.079052
dscore: Iran_N.SG f4: 0.000215 Z: 0.952358

details: Levant_BA MA1 0.000948 2.979787
details: Kenya_PastoralN MA1 -0.012040 -38.431862
dscore: MA1 f4: -0.000151 Z: -0.494533

details: Levant_BA Morocco_EN.SG 0.000754 2.933128
details: Kenya_PastoralN Morocco_EN.SG -0.006819 -25.142524
dscore: Morocco_EN.SG f4: 0.000113 Z: 0.458658

details: Levant_BA CHG 0.001479 5.303404
details: Kenya_PastoralN CHG -0.012943 -47.054098
dscore: CHG f4: 0.000259 Z: 0.967891

details: Levant_BA Italy_Mesolithic.SG 0.001210 4.474158
details: Kenya_PastoralN Italy_Mesolithic.SG -0.012655 -45.198473
dscore: Italy_Mesolithic.SG f4: 0.000037 Z: 0.142190

details: Levant_BA Romania_IronGates_Mesolithic.SG 0.001041 3.944235
details: Kenya_PastoralN Romania_IronGates_Mesolithic.SG -0.012821 -47.617020
dscore: Romania_IronGates_Mesolithic.SG f4: -0.000131 Z: -0.519257

details: Levant_BA Russia_BA_Okunevo.SG 0.000869 3.896733
details: Kenya_PastoralN Russia_BA_Okunevo.SG -0.011853 -49.602181
dscore: Russia_BA_Okunevo.SG f4: -0.000207 Z: -0.965882

details: Levant_BA Germany_EarlyMedieval.SG 0.001158 5.105662
details: Kenya_PastoralN Germany_EarlyMedieval.SG -0.013643 -57.669211
dscore: Germany_EarlyMedieval.SG f4: -0.000094 Z: -0.432724

details: Levant_BA Poland_Southeast_CordedWare.SG 0.001016 4.398323
details: Kenya_PastoralN Poland_Southeast_CordedWare.SG -0.013326 -54.627733
dscore: Poland_Southeast_CordedWare.SG f4: -0.000197 Z: -0.888079

gendstat: South_Africa_400BP.SG Cameroon_SMA.SG 0.042
gendstat: South_Africa_400BP.SG Mota 0.505
gendstat: South_Africa_400BP.SG South_Africa_HG 0.936
gendstat: South_Africa_400BP.SG ZlatyKun -0.665
gendstat: South_Africa_400BP.SG BachoKiro 0.517
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG 0.891
gendstat: South_Africa_400BP.SG DevilsGate_N.SG -0.740
gendstat: South_Africa_400BP.SG Kazakhstan_Botai_Eneolithic.SG -0.825
gendstat: South_Africa_400BP.SG Barcin_N.SG 0.626
gendstat: South_Africa_400BP.SG Iran_N.SG 0.952
gendstat: South_Africa_400BP.SG MA1 -0.495
gendstat: South_Africa_400BP.SG Morocco_EN.SG 0.459
gendstat: South_Africa_400BP.SG CHG 0.968
gendstat: South_Africa_400BP.SG Italy_Mesolithic.SG 0.142
gendstat: South_Africa_400BP.SG Romania_IronGates_Mesolithic.SG -0.519
gendstat: South_Africa_400BP.SG Russia_BA_Okunevo.SG -0.966
gendstat: South_Africa_400BP.SG Germany_EarlyMedieval.SG -0.433
gendstat: South_Africa_400BP.SG Poland_Southeast_CordedWare.SG -0.888
gendstat: Cameroon_SMA.SG Mota 0.430
gendstat: Cameroon_SMA.SG South_Africa_HG 0.841
gendstat: Cameroon_SMA.SG ZlatyKun -0.649
gendstat: Cameroon_SMA.SG BachoKiro 0.427
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG 0.797
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.700
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.772
gendstat: Cameroon_SMA.SG Barcin_N.SG 0.558
gendstat: Cameroon_SMA.SG Iran_N.SG 0.826
gendstat: Cameroon_SMA.SG MA1 -0.503
gendstat: Cameroon_SMA.SG Morocco_EN.SG 0.388
gendstat: Cameroon_SMA.SG CHG 0.870
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG 0.108
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -0.513
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -0.899
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -0.420
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -0.830
gendstat: Mota South_Africa_HG 0.194
gendstat: Mota ZlatyKun -0.898
gendstat: Mota BachoKiro 0.034
gendstat: Mota Eastern_Euro_UP.SG 0.371
gendstat: Mota DevilsGate_N.SG -0.946
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG -1.051
gendstat: Mota Barcin_N.SG 0.157
gendstat: Mota Iran_N.SG 0.389
gendstat: Mota MA1 -0.752
gendstat: Mota Morocco_EN.SG 0.018
gendstat: Mota CHG 0.482
gendstat: Mota Italy_Mesolithic.SG -0.226
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.769
gendstat: Mota Russia_BA_Okunevo.SG -1.193
gendstat: Mota Germany_EarlyMedieval.SG -0.747
gendstat: Mota Poland_Southeast_CordedWare.SG -1.105
gendstat: South_Africa_HG ZlatyKun -1.125
gendstat: South_Africa_HG BachoKiro -0.143
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.254
gendstat: South_Africa_HG DevilsGate_N.SG -1.271
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG -1.335
gendstat: South_Africa_HG Barcin_N.SG 0.009
gendstat: South_Africa_HG Iran_N.SG 0.257
gendstat: South_Africa_HG MA1 -0.961
gendstat: South_Africa_HG Morocco_EN.SG -0.144
gendstat: South_Africa_HG CHG 0.371
gendstat: South_Africa_HG Italy_Mesolithic.SG -0.423
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -1.054
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -1.559
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -1.042
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG -1.455
gendstat: ZlatyKun BachoKiro 1.055
gendstat: ZlatyKun Eastern_Euro_UP.SG 1.364
gendstat: ZlatyKun DevilsGate_N.SG 0.042
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG -0.014
gendstat: ZlatyKun Barcin_N.SG 1.130
gendstat: ZlatyKun Iran_N.SG 1.437
gendstat: ZlatyKun MA1 0.118
gendstat: ZlatyKun Morocco_EN.SG 0.896
gendstat: ZlatyKun CHG 1.427
gendstat: ZlatyKun Italy_Mesolithic.SG 0.703
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG 0.195
gendstat: ZlatyKun Russia_BA_Okunevo.SG -0.056
gendstat: ZlatyKun Germany_EarlyMedieval.SG 0.361
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.017
gendstat: BachoKiro Eastern_Euro_UP.SG 0.419
gendstat: BachoKiro DevilsGate_N.SG -1.219
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -1.290
gendstat: BachoKiro Barcin_N.SG 0.145
gendstat: BachoKiro Iran_N.SG 0.421
gendstat: BachoKiro MA1 -0.867
gendstat: BachoKiro Morocco_EN.SG -0.014
gendstat: BachoKiro CHG 0.532
gendstat: BachoKiro Italy_Mesolithic.SG -0.301
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -0.940
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.536
gendstat: BachoKiro Germany_EarlyMedieval.SG -0.933
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.349
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -1.572
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -1.788
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -0.272
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.020
gendstat: Eastern_Euro_UP.SG MA1 -1.302
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -0.390
gendstat: Eastern_Euro_UP.SG CHG 0.145
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG -0.763
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -1.468
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -2.088
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -1.469
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -1.862
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG -0.070
gendstat: DevilsGate_N.SG Barcin_N.SG 1.282
gendstat: DevilsGate_N.SG Iran_N.SG 1.627
gendstat: DevilsGate_N.SG MA1 0.092
gendstat: DevilsGate_N.SG Morocco_EN.SG 1.012
gendstat: DevilsGate_N.SG CHG 1.618
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.792
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG 0.183
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.136
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG 0.382
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG -0.075
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG 1.457
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 1.832
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 0.174
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG 1.107
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG 1.740
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.951
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG 0.298
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -0.066
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG 0.529
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG -0.004
gendstat: Barcin_N.SG Iran_N.SG 0.261
gendstat: Barcin_N.SG MA1 -1.017
gendstat: Barcin_N.SG Morocco_EN.SG -0.151
gendstat: Barcin_N.SG CHG 0.416
gendstat: Barcin_N.SG Italy_Mesolithic.SG -0.491
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG -1.248
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG -1.694
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG -1.338
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG -1.884
gendstat: Iran_N.SG MA1 -1.295
gendstat: Iran_N.SG Morocco_EN.SG -0.358
gendstat: Iran_N.SG CHG 0.180
gendstat: Iran_N.SG Italy_Mesolithic.SG -0.701
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -1.428
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -2.142
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -1.588
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -2.038
gendstat: MA1 Morocco_EN.SG 0.783
gendstat: MA1 CHG 1.311
gendstat: MA1 Italy_Mesolithic.SG 0.635
gendstat: MA1 Romania_IronGates_Mesolithic.SG 0.070
gendstat: MA1 Russia_BA_Okunevo.SG -0.223
gendstat: MA1 Germany_EarlyMedieval.SG 0.215
gendstat: MA1 Poland_Southeast_CordedWare.SG -0.175
gendstat: Morocco_EN.SG CHG 0.465
gendstat: Morocco_EN.SG Italy_Mesolithic.SG -0.267
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG -0.858
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG -1.242
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG -0.817
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG -1.159
gendstat: CHG Italy_Mesolithic.SG -0.805
gendstat: CHG Romania_IronGates_Mesolithic.SG -1.492
gendstat: CHG Russia_BA_Okunevo.SG -1.981
gendstat: CHG Germany_EarlyMedieval.SG -1.590
gendstat: CHG Poland_Southeast_CordedWare.SG -1.955
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -0.817
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -1.124
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -0.643
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -1.088
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.376
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG 0.193
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.340
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG 0.713
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG 0.061
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG -0.666



you guys frequently give out interesting tables but omit the haplotype data this is annoying what is the reason for this. Thanks

altvred
02-17-2022, 09:47 AM
you guys frequently give out interesting tables but omit the haplotype data this is annoying what is the reason for this. Thanks
What? Haplotype as in Y-DNA/mtDNA?

It's qpAdm. It is a tool used to model autosomal ancestry; it has nothing to do with uniparental markers.

So nobody is omitting anything because qpAdm isn't even looking at the Y or mtDNA.

If you're interested in a particular sample or set of samples being modeled, look it up on the annotation file of the Reich Dataset. They usually have some basic Y-haplogroup info on Male samples.

https://reichdata.hms.harvard.edu/pub/datasets/amh_repo/curated_releases/V50/V50.0/SHARE/public.dir/v50.0_1240K_public.anno

@Mnemonics

Do you have the HO-dataset? Please run the same model for Yemeni Jews if you can. Too lazy to do it myself lol...

Btw, if you're interested, I've got samples from the "Genomic History of the Middle-East" in eigenstrat format with around 1 million coverage for the 1240k dataset.

Take a look at the Saudi samples with less SSA-admixture (SaudiA).

https://drive.google.com/file/d/1_cwsHtWFkCr-XM8EfgiArphlMNgpdY5k/view?usp=sharing

Gentica277282
02-17-2022, 01:00 PM
The BedouinB diploid samples seem to mostly Jordan_EBA-like with extra Kenya_PastoralN. They also seem to have some extra Yana_UP affinity I can't model.





left pops:
BedouinB.DG
Levant_BA
Kenya_PastoralN

right pops:
South_Africa_400BP.SG
Cameroon_SMA.SG
Mota
South_Africa_HG
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 BedouinB.DG 2
1 Levant_BA 3
2 Kenya_PastoralN 8
3 South_Africa_400BP.SG 4
4 Cameroon_SMA.SG 3
5 Mota 1
6 South_Africa_HG 4
7 ZlatyKun 1
8 BachoKiro 3
9 Eastern_Euro_UP.SG 3
10 DevilsGate_N.SG 4
11 Kazakhstan_Botai_Eneolithic.SG 3
12 Barcin_N.SG 3
13 Iran_N.SG 4
14 MA1 1
15 Morocco_EN.SG 3
16 CHG 2
17 Italy_Mesolithic.SG 3
18 Romania_IronGates_Mesolithic.SG 3
19 Russia_BA_Okunevo.SG 14
20 Germany_EarlyMedieval.SG 21
21 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150471 indivs: 106
number of blocks for block jackknife: 714
## ncols: 1150471
coverage: BedouinB.DG 1121254
coverage: Levant_BA 870086
coverage: Kenya_PastoralN 1061629
coverage: South_Africa_400BP.SG 1150180
coverage: Cameroon_SMA.SG 1149654
coverage: Mota 1134612
coverage: South_Africa_HG 1135129
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 623.267
numsnps used: 472322
codimension 1
f4info:
f4rank: 1 dof: 17 chisq: 17.995 tail: 0.389136695 dofdiff: 19 chisqdiff: -17.995 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG -0.021
Mota 0.039
South_Africa_HG -0.077
ZlatyKun 0.869
BachoKiro 0.840
Eastern_Euro_UP.SG 1.084
DevilsGate_N.SG 0.953
Kazakhstan_Botai_Eneolithic.SG 1.106
Barcin_N.SG 1.289
Iran_N.SG 1.159
MA1 1.106
Morocco_EN.SG 0.631
CHG 1.220
Italy_Mesolithic.SG 1.177
Romania_IronGates_Mesolithic.SG 1.184
Russia_BA_Okunevo.SG 1.090
Germany_EarlyMedieval.SG 1.271
Poland_Southeast_CordedWare.SG 1.232
A:
scale 132.284
Levant_BA 0.130
Kenya_PastoralN -1.408


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 17 chisqdiff: 17.995 taildiff: 0.389136695
B:
scale 1005.459 92.843
Cameroon_SMA.SG -0.003 0.011
Mota 0.161 -0.042
South_Africa_HG 0.097 0.070
ZlatyKun 0.671 -0.882
BachoKiro 0.967 -0.838
Eastern_Euro_UP.SG 1.320 -1.078
DevilsGate_N.SG 0.768 -0.965
Kazakhstan_Botai_Eneolithic.SG 0.905 -1.120
Barcin_N.SG 1.445 -1.274
Iran_N.SG 1.380 -1.143
MA1 0.953 -1.118
Morocco_EN.SG 0.758 -0.633
CHG 1.487 -1.202
Italy_Mesolithic.SG 1.216 -1.175
Romania_IronGates_Mesolithic.SG 1.047 -1.190
Russia_BA_Okunevo.SG 0.874 -1.101
Germany_EarlyMedieval.SG 1.164 -1.267
Poland_Southeast_CordedWare.SG 1.021 -1.237
A:
scale 1.414 1.414
Levant_BA 1.414 0.000
Kenya_PastoralN 0.000 1.414


best coefficients: 0.915 0.085
Jackknife mean: 0.915340507 0.084659493
std. errors: 0.012 0.012

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


summ: BedouinB.DG 2 0.389137 0.915 0.085 152 -152 152

fixed pat wt dof chisq tail prob
00 0 17 17.995 0.389137 0.915 0.085
01 1 18 59.603 2.37107e-06 1.000 0.000
10 1 18 4315.882 0 0.000 1.000
best pat: 00 0.389137 - -
best pat: 01 2.37107e-06 chi(nested): 41.607 p-value for nested model: 1.11567e-10

coeffs: 0.915 0.085

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

details: Levant_BA Cameroon_SMA.SG -0.000003 -0.019341
details: Kenya_PastoralN Cameroon_SMA.SG 0.000119 0.732389
dscore: Cameroon_SMA.SG f4: 0.000007 Z: 0.041879

details: Levant_BA Mota 0.000160 0.715608
details: Kenya_PastoralN Mota -0.000454 -2.095219
dscore: Mota f4: 0.000108 Z: 0.505451

details: Levant_BA South_Africa_HG 0.000097 0.566199
details: Kenya_PastoralN South_Africa_HG 0.000756 4.758071
dscore: South_Africa_HG f4: 0.000153 Z: 0.935890

details: Levant_BA ZlatyKun 0.000668 2.221524
details: Kenya_PastoralN ZlatyKun -0.009497 -32.148828
dscore: ZlatyKun f4: -0.000192 Z: -0.665088

details: Levant_BA BachoKiro 0.000962 4.047038
details: Kenya_PastoralN BachoKiro -0.009024 -36.335295
dscore: BachoKiro f4: 0.000117 Z: 0.516532

details: Levant_BA Eastern_Euro_UP.SG 0.001313 5.106083
details: Kenya_PastoralN Eastern_Euro_UP.SG -0.011606 -42.370973
dscore: Eastern_Euro_UP.SG f4: 0.000220 Z: 0.891249

details: Levant_BA DevilsGate_N.SG 0.000764 3.024543
details: Kenya_PastoralN DevilsGate_N.SG -0.010391 -39.774104
dscore: DevilsGate_N.SG f4: -0.000179 Z: -0.739968

details: Levant_BA Kazakhstan_Botai_Eneolithic.SG 0.000900 3.636812
details: Kenya_PastoralN Kazakhstan_Botai_Eneolithic.SG -0.012063 -46.286225
dscore: Kazakhstan_Botai_Eneolithic.SG f4: -0.000196 Z: -0.824513

details: Levant_BA Barcin_N.SG 0.001437 5.586478
details: Kenya_PastoralN Barcin_N.SG -0.013724 -51.492202
dscore: Barcin_N.SG f4: 0.000155 Z: 0.626290

details: Levant_BA Iran_N.SG 0.001372 5.838001
details: Kenya_PastoralN Iran_N.SG -0.012308 -50.079052
dscore: Iran_N.SG f4: 0.000215 Z: 0.952358

details: Levant_BA MA1 0.000948 2.979787
details: Kenya_PastoralN MA1 -0.012040 -38.431862
dscore: MA1 f4: -0.000151 Z: -0.494533

details: Levant_BA Morocco_EN.SG 0.000754 2.933128
details: Kenya_PastoralN Morocco_EN.SG -0.006819 -25.142524
dscore: Morocco_EN.SG f4: 0.000113 Z: 0.458658

details: Levant_BA CHG 0.001479 5.303404
details: Kenya_PastoralN CHG -0.012943 -47.054098
dscore: CHG f4: 0.000259 Z: 0.967891

details: Levant_BA Italy_Mesolithic.SG 0.001210 4.474158
details: Kenya_PastoralN Italy_Mesolithic.SG -0.012655 -45.198473
dscore: Italy_Mesolithic.SG f4: 0.000037 Z: 0.142190

details: Levant_BA Romania_IronGates_Mesolithic.SG 0.001041 3.944235
details: Kenya_PastoralN Romania_IronGates_Mesolithic.SG -0.012821 -47.617020
dscore: Romania_IronGates_Mesolithic.SG f4: -0.000131 Z: -0.519257

details: Levant_BA Russia_BA_Okunevo.SG 0.000869 3.896733
details: Kenya_PastoralN Russia_BA_Okunevo.SG -0.011853 -49.602181
dscore: Russia_BA_Okunevo.SG f4: -0.000207 Z: -0.965882

details: Levant_BA Germany_EarlyMedieval.SG 0.001158 5.105662
details: Kenya_PastoralN Germany_EarlyMedieval.SG -0.013643 -57.669211
dscore: Germany_EarlyMedieval.SG f4: -0.000094 Z: -0.432724

details: Levant_BA Poland_Southeast_CordedWare.SG 0.001016 4.398323
details: Kenya_PastoralN Poland_Southeast_CordedWare.SG -0.013326 -54.627733
dscore: Poland_Southeast_CordedWare.SG f4: -0.000197 Z: -0.888079

gendstat: South_Africa_400BP.SG Cameroon_SMA.SG 0.042
gendstat: South_Africa_400BP.SG Mota 0.505
gendstat: South_Africa_400BP.SG South_Africa_HG 0.936
gendstat: South_Africa_400BP.SG ZlatyKun -0.665
gendstat: South_Africa_400BP.SG BachoKiro 0.517
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG 0.891
gendstat: South_Africa_400BP.SG DevilsGate_N.SG -0.740
gendstat: South_Africa_400BP.SG Kazakhstan_Botai_Eneolithic.SG -0.825
gendstat: South_Africa_400BP.SG Barcin_N.SG 0.626
gendstat: South_Africa_400BP.SG Iran_N.SG 0.952
gendstat: South_Africa_400BP.SG MA1 -0.495
gendstat: South_Africa_400BP.SG Morocco_EN.SG 0.459
gendstat: South_Africa_400BP.SG CHG 0.968
gendstat: South_Africa_400BP.SG Italy_Mesolithic.SG 0.142
gendstat: South_Africa_400BP.SG Romania_IronGates_Mesolithic.SG -0.519
gendstat: South_Africa_400BP.SG Russia_BA_Okunevo.SG -0.966
gendstat: South_Africa_400BP.SG Germany_EarlyMedieval.SG -0.433
gendstat: South_Africa_400BP.SG Poland_Southeast_CordedWare.SG -0.888
gendstat: Cameroon_SMA.SG Mota 0.430
gendstat: Cameroon_SMA.SG South_Africa_HG 0.841
gendstat: Cameroon_SMA.SG ZlatyKun -0.649
gendstat: Cameroon_SMA.SG BachoKiro 0.427
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG 0.797
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.700
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.772
gendstat: Cameroon_SMA.SG Barcin_N.SG 0.558
gendstat: Cameroon_SMA.SG Iran_N.SG 0.826
gendstat: Cameroon_SMA.SG MA1 -0.503
gendstat: Cameroon_SMA.SG Morocco_EN.SG 0.388
gendstat: Cameroon_SMA.SG CHG 0.870
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG 0.108
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -0.513
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -0.899
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -0.420
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -0.830
gendstat: Mota South_Africa_HG 0.194
gendstat: Mota ZlatyKun -0.898
gendstat: Mota BachoKiro 0.034
gendstat: Mota Eastern_Euro_UP.SG 0.371
gendstat: Mota DevilsGate_N.SG -0.946
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG -1.051
gendstat: Mota Barcin_N.SG 0.157
gendstat: Mota Iran_N.SG 0.389
gendstat: Mota MA1 -0.752
gendstat: Mota Morocco_EN.SG 0.018
gendstat: Mota CHG 0.482
gendstat: Mota Italy_Mesolithic.SG -0.226
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.769
gendstat: Mota Russia_BA_Okunevo.SG -1.193
gendstat: Mota Germany_EarlyMedieval.SG -0.747
gendstat: Mota Poland_Southeast_CordedWare.SG -1.105
gendstat: South_Africa_HG ZlatyKun -1.125
gendstat: South_Africa_HG BachoKiro -0.143
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.254
gendstat: South_Africa_HG DevilsGate_N.SG -1.271
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG -1.335
gendstat: South_Africa_HG Barcin_N.SG 0.009
gendstat: South_Africa_HG Iran_N.SG 0.257
gendstat: South_Africa_HG MA1 -0.961
gendstat: South_Africa_HG Morocco_EN.SG -0.144
gendstat: South_Africa_HG CHG 0.371
gendstat: South_Africa_HG Italy_Mesolithic.SG -0.423
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -1.054
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -1.559
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -1.042
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG -1.455
gendstat: ZlatyKun BachoKiro 1.055
gendstat: ZlatyKun Eastern_Euro_UP.SG 1.364
gendstat: ZlatyKun DevilsGate_N.SG 0.042
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG -0.014
gendstat: ZlatyKun Barcin_N.SG 1.130
gendstat: ZlatyKun Iran_N.SG 1.437
gendstat: ZlatyKun MA1 0.118
gendstat: ZlatyKun Morocco_EN.SG 0.896
gendstat: ZlatyKun CHG 1.427
gendstat: ZlatyKun Italy_Mesolithic.SG 0.703
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG 0.195
gendstat: ZlatyKun Russia_BA_Okunevo.SG -0.056
gendstat: ZlatyKun Germany_EarlyMedieval.SG 0.361
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.017
gendstat: BachoKiro Eastern_Euro_UP.SG 0.419
gendstat: BachoKiro DevilsGate_N.SG -1.219
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -1.290
gendstat: BachoKiro Barcin_N.SG 0.145
gendstat: BachoKiro Iran_N.SG 0.421
gendstat: BachoKiro MA1 -0.867
gendstat: BachoKiro Morocco_EN.SG -0.014
gendstat: BachoKiro CHG 0.532
gendstat: BachoKiro Italy_Mesolithic.SG -0.301
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -0.940
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.536
gendstat: BachoKiro Germany_EarlyMedieval.SG -0.933
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.349
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -1.572
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -1.788
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -0.272
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.020
gendstat: Eastern_Euro_UP.SG MA1 -1.302
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -0.390
gendstat: Eastern_Euro_UP.SG CHG 0.145
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG -0.763
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -1.468
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -2.088
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -1.469
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -1.862
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG -0.070
gendstat: DevilsGate_N.SG Barcin_N.SG 1.282
gendstat: DevilsGate_N.SG Iran_N.SG 1.627
gendstat: DevilsGate_N.SG MA1 0.092
gendstat: DevilsGate_N.SG Morocco_EN.SG 1.012
gendstat: DevilsGate_N.SG CHG 1.618
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.792
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG 0.183
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.136
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG 0.382
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG -0.075
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG 1.457
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 1.832
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 0.174
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG 1.107
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG 1.740
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.951
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG 0.298
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -0.066
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG 0.529
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG -0.004
gendstat: Barcin_N.SG Iran_N.SG 0.261
gendstat: Barcin_N.SG MA1 -1.017
gendstat: Barcin_N.SG Morocco_EN.SG -0.151
gendstat: Barcin_N.SG CHG 0.416
gendstat: Barcin_N.SG Italy_Mesolithic.SG -0.491
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG -1.248
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG -1.694
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG -1.338
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG -1.884
gendstat: Iran_N.SG MA1 -1.295
gendstat: Iran_N.SG Morocco_EN.SG -0.358
gendstat: Iran_N.SG CHG 0.180
gendstat: Iran_N.SG Italy_Mesolithic.SG -0.701
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -1.428
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -2.142
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -1.588
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -2.038
gendstat: MA1 Morocco_EN.SG 0.783
gendstat: MA1 CHG 1.311
gendstat: MA1 Italy_Mesolithic.SG 0.635
gendstat: MA1 Romania_IronGates_Mesolithic.SG 0.070
gendstat: MA1 Russia_BA_Okunevo.SG -0.223
gendstat: MA1 Germany_EarlyMedieval.SG 0.215
gendstat: MA1 Poland_Southeast_CordedWare.SG -0.175
gendstat: Morocco_EN.SG CHG 0.465
gendstat: Morocco_EN.SG Italy_Mesolithic.SG -0.267
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG -0.858
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG -1.242
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG -0.817
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG -1.159
gendstat: CHG Italy_Mesolithic.SG -0.805
gendstat: CHG Romania_IronGates_Mesolithic.SG -1.492
gendstat: CHG Russia_BA_Okunevo.SG -1.981
gendstat: CHG Germany_EarlyMedieval.SG -1.590
gendstat: CHG Poland_Southeast_CordedWare.SG -1.955
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -0.817
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -1.124
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -0.643
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -1.088
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.376
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG 0.193
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.340
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG 0.713
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG 0.061
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG -0.666





Thanks. Makes sense even though I don’t like mixing moderns with ancients I notice an increase in my ibm when using a modern Arabian sample as opposed to natufian. I guess we will have to wait to see if they find this indigenous Arabian sample.


Distance: 1.4533% / 0.01453315
30.8 Early_European_Farmer
25.6 Saudi
23.4 Iberomaurusian
12.6 ZAF_400BP
5.6 Steppe_Pastoralist
1.4 Iran_Neolithic
0.6 Western_Hunter-Gatherer

Using Levant Jor
Distance: 1.4802% / 0.01480177
35.4 Levant_JOR_EBA
23.2 Iberomaurusian
22.2 Early_European_Farmer
13.2 ZAF_400BP
4.2 Steppe_Pastoralist
1.8 Western_Hunter-Gatherer

Mnemonics
02-17-2022, 03:10 PM
What? Haplotype as in Y-DNA/mtDNA?

It's qpAdm. It is a tool used to model autosomal ancestry; it has nothing to do with uniparental markers.

So nobody is omitting anything because qpAdm isn't even looking at the Y or mtDNA.

If you're interested in a particular sample or set of samples being modeled, look it up on the annotation file of the Reich Dataset. They usually have some basic Y-haplogroup info on Male samples.

https://reichdata.hms.harvard.edu/pub/datasets/amh_repo/curated_releases/V50/V50.0/SHARE/public.dir/v50.0_1240K_public.anno

@Mnemonics

Do you have the HO-dataset? Please run the same model for Yemeni Jews if you can. Too lazy to do it myself lol...

Btw, if you're interested, I've got samples from the "Genomic History of the Middle-East" in eigenstrat format with around 1 million coverage for the 1240k dataset.

Take a look at the Saudi samples with less SSA-admixture (SaudiA).

https://drive.google.com/file/d/1_cwsHtWFkCr-XM8EfgiArphlMNgpdY5k/view?usp=sharing


There are a couple Yemenite_Jewish .DG samples.






left pops:
Jew_Yemenite.DG
Levant_BA
Kenya_PastoralN

right pops:
South_Africa_400BP.SG
Cameroon_SMA.SG
Mota
South_Africa_HG
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 Jew_Yemenite.DG 2
1 Levant_BA 3
2 Kenya_PastoralN 8
3 South_Africa_400BP.SG 4
4 Cameroon_SMA.SG 3
5 Mota 1
6 South_Africa_HG 4
7 ZlatyKun 1
8 BachoKiro 3
9 Eastern_Euro_UP.SG 3
10 DevilsGate_N.SG 4
11 Kazakhstan_Botai_Eneolithic.SG 3
12 Barcin_N.SG 3
13 Iran_N.SG 4
14 MA1 1
15 Morocco_EN.SG 3
16 CHG 2
17 Italy_Mesolithic.SG 3
18 Romania_IronGates_Mesolithic.SG 3
19 Russia_BA_Okunevo.SG 14
20 Germany_EarlyMedieval.SG 21
21 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150471 indivs: 106
number of blocks for block jackknife: 714
## ncols: 1150471
coverage: Jew_Yemenite.DG 1121336
coverage: Levant_BA 870086
coverage: Kenya_PastoralN 1061629
coverage: South_Africa_400BP.SG 1150180
coverage: Cameroon_SMA.SG 1149654
coverage: Mota 1134612
coverage: South_Africa_HG 1135129
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 623.434
numsnps used: 470822
codimension 1
f4info:
f4rank: 1 dof: 17 chisq: 23.380 tail: 0.137283226 dofdiff: 19 chisqdiff: -23.380 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG -0.039
Mota 0.037
South_Africa_HG -0.086
ZlatyKun 0.861
BachoKiro 0.832
Eastern_Euro_UP.SG 1.082
DevilsGate_N.SG 0.950
Kazakhstan_Botai_Eneolithic.SG 1.096
Barcin_N.SG 1.308
Iran_N.SG 1.166
MA1 1.085
Morocco_EN.SG 0.625
CHG 1.234
Italy_Mesolithic.SG 1.172
Romania_IronGates_Mesolithic.SG 1.183
Russia_BA_Okunevo.SG 1.089
Germany_EarlyMedieval.SG 1.277
Poland_Southeast_CordedWare.SG 1.237
A:
scale 132.102
Levant_BA 0.109
Kenya_PastoralN -1.410


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 17 chisqdiff: 23.380 taildiff: 0.137283226
B:
scale 893.300 93.543
Cameroon_SMA.SG 0.289 0.042
Mota 0.239 -0.033
South_Africa_HG 0.260 0.089
ZlatyKun 0.890 -0.860
BachoKiro 1.067 -0.825
Eastern_Euro_UP.SG 1.294 -1.077
DevilsGate_N.SG 0.976 -0.945
Kazakhstan_Botai_Eneolithic.SG 1.115 -1.099
Barcin_N.SG 1.030 -1.315
Iran_N.SG 1.227 -1.155
MA1 1.272 -1.086
Morocco_EN.SG 0.888 -0.618
CHG 1.186 -1.229
Italy_Mesolithic.SG 1.258 -1.170
Romania_IronGates_Mesolithic.SG 1.074 -1.188
Russia_BA_Okunevo.SG 0.946 -1.095
Germany_EarlyMedieval.SG 1.041 -1.280
Poland_Southeast_CordedWare.SG 0.946 -1.247
A:
scale 1.414 1.414
Levant_BA 1.414 0.000
Kenya_PastoralN 0.000 1.414


best coefficients: 0.928 0.072
Jackknife mean: 0.927858969 0.072141031
std. errors: 0.012 0.012

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


summ: Jew_Yemenite.DG 2 0.137283 0.928 0.072 139 -139 139

fixed pat wt dof chisq tail prob
00 0 17 23.380 0.137283 0.928 0.072
01 1 18 57.960 4.35069e-06 1.000 0.000
10 1 18 4183.569 0 0.000 1.000
best pat: 00 0.137283 - -
best pat: 01 4.35069e-06 chi(nested): 34.580 p-value for nested model: 4.09027e-09

coeffs: 0.928 0.072

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

details: Levant_BA Cameroon_SMA.SG 0.000324 1.916354
details: Kenya_PastoralN Cameroon_SMA.SG 0.000448 2.902760
dscore: Cameroon_SMA.SG f4: 0.000333 Z: 2.058401

details: Levant_BA Mota 0.000267 1.167541
details: Kenya_PastoralN Mota -0.000348 -1.647710
dscore: Mota f4: 0.000223 Z: 1.014637

details: Levant_BA South_Africa_HG 0.000291 1.737853
details: Kenya_PastoralN South_Africa_HG 0.000954 6.012059
dscore: South_Africa_HG f4: 0.000338 Z: 2.109003

details: Levant_BA ZlatyKun 0.000997 3.521312
details: Kenya_PastoralN ZlatyKun -0.009199 -31.261367
dscore: ZlatyKun f4: 0.000263 Z: 0.960649

details: Levant_BA BachoKiro 0.001195 4.660623
details: Kenya_PastoralN BachoKiro -0.008825 -34.593163
dscore: BachoKiro f4: 0.000474 Z: 1.917053

details: Levant_BA Eastern_Euro_UP.SG 0.001448 5.707467
details: Kenya_PastoralN Eastern_Euro_UP.SG -0.011512 -43.168621
dscore: Eastern_Euro_UP.SG f4: 0.000515 Z: 2.107414

details: Levant_BA DevilsGate_N.SG 0.001093 4.250497
details: Kenya_PastoralN DevilsGate_N.SG -0.010097 -40.095010
dscore: DevilsGate_N.SG f4: 0.000287 Z: 1.160676

details: Levant_BA Kazakhstan_Botai_Eneolithic.SG 0.001248 5.094405
details: Kenya_PastoralN Kazakhstan_Botai_Eneolithic.SG -0.011754 -45.398805
dscore: Kazakhstan_Botai_Eneolithic.SG f4: 0.000313 Z: 1.319379

details: Levant_BA Barcin_N.SG 0.001153 4.732148
details: Kenya_PastoralN Barcin_N.SG -0.014056 -53.377260
dscore: Barcin_N.SG f4: 0.000059 Z: 0.249408

details: Levant_BA Iran_N.SG 0.001374 5.972022
details: Kenya_PastoralN Iran_N.SG -0.012349 -49.422139
dscore: Iran_N.SG f4: 0.000386 Z: 1.733638

details: Levant_BA MA1 0.001424 4.520751
details: Kenya_PastoralN MA1 -0.011606 -38.629868
dscore: MA1 f4: 0.000486 Z: 1.606060

details: Levant_BA Morocco_EN.SG 0.000994 3.751526
details: Kenya_PastoralN Morocco_EN.SG -0.006605 -24.760311
dscore: Morocco_EN.SG f4: 0.000448 Z: 1.749453

details: Levant_BA CHG 0.001328 4.956868
details: Kenya_PastoralN CHG -0.013139 -47.842898
dscore: CHG f4: 0.000287 Z: 1.108553

details: Levant_BA Italy_Mesolithic.SG 0.001408 5.096107
details: Kenya_PastoralN Italy_Mesolithic.SG -0.012505 -43.704703
dscore: Italy_Mesolithic.SG f4: 0.000407 Z: 1.524340

details: Levant_BA Romania_IronGates_Mesolithic.SG 0.001202 4.539676
details: Kenya_PastoralN Romania_IronGates_Mesolithic.SG -0.012704 -46.279528
dscore: Romania_IronGates_Mesolithic.SG f4: 0.000201 Z: 0.786085

details: Levant_BA Russia_BA_Okunevo.SG 0.001058 4.690690
details: Kenya_PastoralN Russia_BA_Okunevo.SG -0.011704 -51.118939
dscore: Russia_BA_Okunevo.SG f4: 0.000140 Z: 0.643493

details: Levant_BA Germany_EarlyMedieval.SG 0.001165 5.288921
details: Kenya_PastoralN Germany_EarlyMedieval.SG -0.013683 -57.036142
dscore: Germany_EarlyMedieval.SG f4: 0.000097 Z: 0.454600

details: Levant_BA Poland_Southeast_CordedWare.SG 0.001059 4.499900
details: Kenya_PastoralN Poland_Southeast_CordedWare.SG -0.013331 -53.599263
dscore: Poland_Southeast_CordedWare.SG f4: 0.000023 Z: 0.101928

gendstat: South_Africa_400BP.SG Cameroon_SMA.SG 2.058
gendstat: South_Africa_400BP.SG Mota 1.015
gendstat: South_Africa_400BP.SG South_Africa_HG 2.109
gendstat: South_Africa_400BP.SG ZlatyKun 0.961
gendstat: South_Africa_400BP.SG BachoKiro 1.917
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG 2.107
gendstat: South_Africa_400BP.SG DevilsGate_N.SG 1.161
gendstat: South_Africa_400BP.SG Kazakhstan_Botai_Eneolithic.SG 1.319
gendstat: South_Africa_400BP.SG Barcin_N.SG 0.249
gendstat: South_Africa_400BP.SG Iran_N.SG 1.734
gendstat: South_Africa_400BP.SG MA1 1.606
gendstat: South_Africa_400BP.SG Morocco_EN.SG 1.749
gendstat: South_Africa_400BP.SG CHG 1.109
gendstat: South_Africa_400BP.SG Italy_Mesolithic.SG 1.524
gendstat: South_Africa_400BP.SG Romania_IronGates_Mesolithic.SG 0.786
gendstat: South_Africa_400BP.SG Russia_BA_Okunevo.SG 0.643
gendstat: South_Africa_400BP.SG Germany_EarlyMedieval.SG 0.455
gendstat: South_Africa_400BP.SG Poland_Southeast_CordedWare.SG 0.102
gendstat: Cameroon_SMA.SG Mota -0.477
gendstat: Cameroon_SMA.SG South_Africa_HG 0.035
gendstat: Cameroon_SMA.SG ZlatyKun -0.241
gendstat: Cameroon_SMA.SG BachoKiro 0.529
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG 0.699
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.173
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.080
gendstat: Cameroon_SMA.SG Barcin_N.SG -1.079
gendstat: Cameroon_SMA.SG Iran_N.SG 0.222
gendstat: Cameroon_SMA.SG MA1 0.495
gendstat: Cameroon_SMA.SG Morocco_EN.SG 0.424
gendstat: Cameroon_SMA.SG CHG -0.166
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG 0.267
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -0.491
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -0.834
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -1.044
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -1.251
gendstat: Mota South_Africa_HG 0.520
gendstat: Mota ZlatyKun 0.129
gendstat: Mota BachoKiro 0.863
gendstat: Mota Eastern_Euro_UP.SG 1.013
gendstat: Mota DevilsGate_N.SG 0.223
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG 0.328
gendstat: Mota Barcin_N.SG -0.576
gendstat: Mota Iran_N.SG 0.616
gendstat: Mota MA1 0.779
gendstat: Mota Morocco_EN.SG 0.754
gendstat: Mota CHG 0.218
gendstat: Mota Italy_Mesolithic.SG 0.603
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.076
gendstat: Mota Russia_BA_Okunevo.SG -0.328
gendstat: Mota Germany_EarlyMedieval.SG -0.498
gendstat: Mota Poland_Southeast_CordedWare.SG -0.739
gendstat: South_Africa_HG ZlatyKun -0.265
gendstat: South_Africa_HG BachoKiro 0.532
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.705
gendstat: South_Africa_HG DevilsGate_N.SG -0.201
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG -0.103
gendstat: South_Africa_HG Barcin_N.SG -1.123
gendstat: South_Africa_HG Iran_N.SG 0.199
gendstat: South_Africa_HG MA1 0.484
gendstat: South_Africa_HG Morocco_EN.SG 0.404
gendstat: South_Africa_HG CHG -0.189
gendstat: South_Africa_HG Italy_Mesolithic.SG 0.246
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -0.519
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -0.862
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -1.065
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG -1.303
gendstat: ZlatyKun BachoKiro 0.746
gendstat: ZlatyKun Eastern_Euro_UP.SG 0.893
gendstat: ZlatyKun DevilsGate_N.SG 0.088
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG 0.181
gendstat: ZlatyKun Barcin_N.SG -0.736
gendstat: ZlatyKun Iran_N.SG 0.455
gendstat: ZlatyKun MA1 0.652
gendstat: ZlatyKun Morocco_EN.SG 0.593
gendstat: ZlatyKun CHG 0.079
gendstat: ZlatyKun Italy_Mesolithic.SG 0.486
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG -0.210
gendstat: ZlatyKun Russia_BA_Okunevo.SG -0.486
gendstat: ZlatyKun Germany_EarlyMedieval.SG -0.658
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.892
gendstat: BachoKiro Eastern_Euro_UP.SG 0.168
gendstat: BachoKiro DevilsGate_N.SG -0.771
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -0.636
gendstat: BachoKiro Barcin_N.SG -1.628
gendstat: BachoKiro Iran_N.SG -0.355
gendstat: BachoKiro MA1 0.041
gendstat: BachoKiro Morocco_EN.SG -0.091
gendstat: BachoKiro CHG -0.659
gendstat: BachoKiro Italy_Mesolithic.SG -0.244
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -1.015
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.483
gendstat: BachoKiro Germany_EarlyMedieval.SG -1.648
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.928
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -0.931
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -0.903
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -2.005
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.553
gendstat: Eastern_Euro_UP.SG MA1 -0.102
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -0.247
gendstat: Eastern_Euro_UP.SG CHG -0.867
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG -0.456
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -1.302
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -1.843
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -2.053
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -2.291
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG 0.113
gendstat: DevilsGate_N.SG Barcin_N.SG -0.904
gendstat: DevilsGate_N.SG Iran_N.SG 0.424
gendstat: DevilsGate_N.SG MA1 0.663
gendstat: DevilsGate_N.SG Morocco_EN.SG 0.557
gendstat: DevilsGate_N.SG CHG -0.001
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.440
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG -0.339
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.747
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG -0.883
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG -1.147
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG -1.160
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 0.349
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 0.662
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG 0.502
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG -0.101
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.393
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG -0.493
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -1.054
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG -1.187
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG -1.515
gendstat: Barcin_N.SG Iran_N.SG 1.544
gendstat: Barcin_N.SG MA1 1.468
gendstat: Barcin_N.SG Morocco_EN.SG 1.426
gendstat: Barcin_N.SG CHG 0.937
gendstat: Barcin_N.SG Italy_Mesolithic.SG 1.435
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG 0.637
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG 0.397
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG 0.215
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG -0.193
gendstat: Iran_N.SG MA1 0.355
gendstat: Iran_N.SG Morocco_EN.SG 0.225
gendstat: Iran_N.SG CHG -0.416
gendstat: Iran_N.SG Italy_Mesolithic.SG 0.081
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -0.772
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -1.273
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -1.523
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -1.859
gendstat: MA1 Morocco_EN.SG -0.117
gendstat: MA1 CHG -0.637
gendstat: MA1 Italy_Mesolithic.SG -0.263
gendstat: MA1 Romania_IronGates_Mesolithic.SG -0.997
gendstat: MA1 Russia_BA_Okunevo.SG -1.373
gendstat: MA1 Germany_EarlyMedieval.SG -1.463
gendstat: MA1 Poland_Southeast_CordedWare.SG -1.722
gendstat: Morocco_EN.SG CHG -0.543
gendstat: Morocco_EN.SG Italy_Mesolithic.SG -0.134
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG -0.832
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG -1.221
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG -1.375
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG -1.601
gendstat: CHG Italy_Mesolithic.SG 0.428
gendstat: CHG Romania_IronGates_Mesolithic.SG -0.331
gendstat: CHG Russia_BA_Okunevo.SG -0.637
gendstat: CHG Germany_EarlyMedieval.SG -0.872
gendstat: CHG Poland_Southeast_CordedWare.SG -1.153
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -1.031
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -1.208
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -1.502
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -1.706
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.287
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG -0.536
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.858
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG -0.278
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG -0.679
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG -0.449

altvred
02-17-2022, 06:08 PM
There are a couple Yemenite_Jewish .DG samples.






left pops:
Jew_Yemenite.DG
Levant_BA
Kenya_PastoralN

right pops:
South_Africa_400BP.SG
Cameroon_SMA.SG
Mota
South_Africa_HG
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 Jew_Yemenite.DG 2
1 Levant_BA 3
2 Kenya_PastoralN 8
3 South_Africa_400BP.SG 4
4 Cameroon_SMA.SG 3
5 Mota 1
6 South_Africa_HG 4
7 ZlatyKun 1
8 BachoKiro 3
9 Eastern_Euro_UP.SG 3
10 DevilsGate_N.SG 4
11 Kazakhstan_Botai_Eneolithic.SG 3
12 Barcin_N.SG 3
13 Iran_N.SG 4
14 MA1 1
15 Morocco_EN.SG 3
16 CHG 2
17 Italy_Mesolithic.SG 3
18 Romania_IronGates_Mesolithic.SG 3
19 Russia_BA_Okunevo.SG 14
20 Germany_EarlyMedieval.SG 21
21 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150471 indivs: 106
number of blocks for block jackknife: 714
## ncols: 1150471
coverage: Jew_Yemenite.DG 1121336
coverage: Levant_BA 870086
coverage: Kenya_PastoralN 1061629
coverage: South_Africa_400BP.SG 1150180
coverage: Cameroon_SMA.SG 1149654
coverage: Mota 1134612
coverage: South_Africa_HG 1135129
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 623.434
numsnps used: 470822
codimension 1
f4info:
f4rank: 1 dof: 17 chisq: 23.380 tail: 0.137283226 dofdiff: 19 chisqdiff: -23.380 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG -0.039
Mota 0.037
South_Africa_HG -0.086
ZlatyKun 0.861
BachoKiro 0.832
Eastern_Euro_UP.SG 1.082
DevilsGate_N.SG 0.950
Kazakhstan_Botai_Eneolithic.SG 1.096
Barcin_N.SG 1.308
Iran_N.SG 1.166
MA1 1.085
Morocco_EN.SG 0.625
CHG 1.234
Italy_Mesolithic.SG 1.172
Romania_IronGates_Mesolithic.SG 1.183
Russia_BA_Okunevo.SG 1.089
Germany_EarlyMedieval.SG 1.277
Poland_Southeast_CordedWare.SG 1.237
A:
scale 132.102
Levant_BA 0.109
Kenya_PastoralN -1.410


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 17 chisqdiff: 23.380 taildiff: 0.137283226
B:
scale 893.300 93.543
Cameroon_SMA.SG 0.289 0.042
Mota 0.239 -0.033
South_Africa_HG 0.260 0.089
ZlatyKun 0.890 -0.860
BachoKiro 1.067 -0.825
Eastern_Euro_UP.SG 1.294 -1.077
DevilsGate_N.SG 0.976 -0.945
Kazakhstan_Botai_Eneolithic.SG 1.115 -1.099
Barcin_N.SG 1.030 -1.315
Iran_N.SG 1.227 -1.155
MA1 1.272 -1.086
Morocco_EN.SG 0.888 -0.618
CHG 1.186 -1.229
Italy_Mesolithic.SG 1.258 -1.170
Romania_IronGates_Mesolithic.SG 1.074 -1.188
Russia_BA_Okunevo.SG 0.946 -1.095
Germany_EarlyMedieval.SG 1.041 -1.280
Poland_Southeast_CordedWare.SG 0.946 -1.247
A:
scale 1.414 1.414
Levant_BA 1.414 0.000
Kenya_PastoralN 0.000 1.414


best coefficients: 0.928 0.072
Jackknife mean: 0.927858969 0.072141031
std. errors: 0.012 0.012

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


summ: Jew_Yemenite.DG 2 0.137283 0.928 0.072 139 -139 139

fixed pat wt dof chisq tail prob
00 0 17 23.380 0.137283 0.928 0.072
01 1 18 57.960 4.35069e-06 1.000 0.000
10 1 18 4183.569 0 0.000 1.000
best pat: 00 0.137283 - -
best pat: 01 4.35069e-06 chi(nested): 34.580 p-value for nested model: 4.09027e-09

coeffs: 0.928 0.072

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

details: Levant_BA Cameroon_SMA.SG 0.000324 1.916354
details: Kenya_PastoralN Cameroon_SMA.SG 0.000448 2.902760
dscore: Cameroon_SMA.SG f4: 0.000333 Z: 2.058401

details: Levant_BA Mota 0.000267 1.167541
details: Kenya_PastoralN Mota -0.000348 -1.647710
dscore: Mota f4: 0.000223 Z: 1.014637

details: Levant_BA South_Africa_HG 0.000291 1.737853
details: Kenya_PastoralN South_Africa_HG 0.000954 6.012059
dscore: South_Africa_HG f4: 0.000338 Z: 2.109003

details: Levant_BA ZlatyKun 0.000997 3.521312
details: Kenya_PastoralN ZlatyKun -0.009199 -31.261367
dscore: ZlatyKun f4: 0.000263 Z: 0.960649

details: Levant_BA BachoKiro 0.001195 4.660623
details: Kenya_PastoralN BachoKiro -0.008825 -34.593163
dscore: BachoKiro f4: 0.000474 Z: 1.917053

details: Levant_BA Eastern_Euro_UP.SG 0.001448 5.707467
details: Kenya_PastoralN Eastern_Euro_UP.SG -0.011512 -43.168621
dscore: Eastern_Euro_UP.SG f4: 0.000515 Z: 2.107414

details: Levant_BA DevilsGate_N.SG 0.001093 4.250497
details: Kenya_PastoralN DevilsGate_N.SG -0.010097 -40.095010
dscore: DevilsGate_N.SG f4: 0.000287 Z: 1.160676

details: Levant_BA Kazakhstan_Botai_Eneolithic.SG 0.001248 5.094405
details: Kenya_PastoralN Kazakhstan_Botai_Eneolithic.SG -0.011754 -45.398805
dscore: Kazakhstan_Botai_Eneolithic.SG f4: 0.000313 Z: 1.319379

details: Levant_BA Barcin_N.SG 0.001153 4.732148
details: Kenya_PastoralN Barcin_N.SG -0.014056 -53.377260
dscore: Barcin_N.SG f4: 0.000059 Z: 0.249408

details: Levant_BA Iran_N.SG 0.001374 5.972022
details: Kenya_PastoralN Iran_N.SG -0.012349 -49.422139
dscore: Iran_N.SG f4: 0.000386 Z: 1.733638

details: Levant_BA MA1 0.001424 4.520751
details: Kenya_PastoralN MA1 -0.011606 -38.629868
dscore: MA1 f4: 0.000486 Z: 1.606060

details: Levant_BA Morocco_EN.SG 0.000994 3.751526
details: Kenya_PastoralN Morocco_EN.SG -0.006605 -24.760311
dscore: Morocco_EN.SG f4: 0.000448 Z: 1.749453

details: Levant_BA CHG 0.001328 4.956868
details: Kenya_PastoralN CHG -0.013139 -47.842898
dscore: CHG f4: 0.000287 Z: 1.108553

details: Levant_BA Italy_Mesolithic.SG 0.001408 5.096107
details: Kenya_PastoralN Italy_Mesolithic.SG -0.012505 -43.704703
dscore: Italy_Mesolithic.SG f4: 0.000407 Z: 1.524340

details: Levant_BA Romania_IronGates_Mesolithic.SG 0.001202 4.539676
details: Kenya_PastoralN Romania_IronGates_Mesolithic.SG -0.012704 -46.279528
dscore: Romania_IronGates_Mesolithic.SG f4: 0.000201 Z: 0.786085

details: Levant_BA Russia_BA_Okunevo.SG 0.001058 4.690690
details: Kenya_PastoralN Russia_BA_Okunevo.SG -0.011704 -51.118939
dscore: Russia_BA_Okunevo.SG f4: 0.000140 Z: 0.643493

details: Levant_BA Germany_EarlyMedieval.SG 0.001165 5.288921
details: Kenya_PastoralN Germany_EarlyMedieval.SG -0.013683 -57.036142
dscore: Germany_EarlyMedieval.SG f4: 0.000097 Z: 0.454600

details: Levant_BA Poland_Southeast_CordedWare.SG 0.001059 4.499900
details: Kenya_PastoralN Poland_Southeast_CordedWare.SG -0.013331 -53.599263
dscore: Poland_Southeast_CordedWare.SG f4: 0.000023 Z: 0.101928

gendstat: South_Africa_400BP.SG Cameroon_SMA.SG 2.058
gendstat: South_Africa_400BP.SG Mota 1.015
gendstat: South_Africa_400BP.SG South_Africa_HG 2.109
gendstat: South_Africa_400BP.SG ZlatyKun 0.961
gendstat: South_Africa_400BP.SG BachoKiro 1.917
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG 2.107
gendstat: South_Africa_400BP.SG DevilsGate_N.SG 1.161
gendstat: South_Africa_400BP.SG Kazakhstan_Botai_Eneolithic.SG 1.319
gendstat: South_Africa_400BP.SG Barcin_N.SG 0.249
gendstat: South_Africa_400BP.SG Iran_N.SG 1.734
gendstat: South_Africa_400BP.SG MA1 1.606
gendstat: South_Africa_400BP.SG Morocco_EN.SG 1.749
gendstat: South_Africa_400BP.SG CHG 1.109
gendstat: South_Africa_400BP.SG Italy_Mesolithic.SG 1.524
gendstat: South_Africa_400BP.SG Romania_IronGates_Mesolithic.SG 0.786
gendstat: South_Africa_400BP.SG Russia_BA_Okunevo.SG 0.643
gendstat: South_Africa_400BP.SG Germany_EarlyMedieval.SG 0.455
gendstat: South_Africa_400BP.SG Poland_Southeast_CordedWare.SG 0.102
gendstat: Cameroon_SMA.SG Mota -0.477
gendstat: Cameroon_SMA.SG South_Africa_HG 0.035
gendstat: Cameroon_SMA.SG ZlatyKun -0.241
gendstat: Cameroon_SMA.SG BachoKiro 0.529
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG 0.699
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.173
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.080
gendstat: Cameroon_SMA.SG Barcin_N.SG -1.079
gendstat: Cameroon_SMA.SG Iran_N.SG 0.222
gendstat: Cameroon_SMA.SG MA1 0.495
gendstat: Cameroon_SMA.SG Morocco_EN.SG 0.424
gendstat: Cameroon_SMA.SG CHG -0.166
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG 0.267
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -0.491
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -0.834
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -1.044
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -1.251
gendstat: Mota South_Africa_HG 0.520
gendstat: Mota ZlatyKun 0.129
gendstat: Mota BachoKiro 0.863
gendstat: Mota Eastern_Euro_UP.SG 1.013
gendstat: Mota DevilsGate_N.SG 0.223
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG 0.328
gendstat: Mota Barcin_N.SG -0.576
gendstat: Mota Iran_N.SG 0.616
gendstat: Mota MA1 0.779
gendstat: Mota Morocco_EN.SG 0.754
gendstat: Mota CHG 0.218
gendstat: Mota Italy_Mesolithic.SG 0.603
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.076
gendstat: Mota Russia_BA_Okunevo.SG -0.328
gendstat: Mota Germany_EarlyMedieval.SG -0.498
gendstat: Mota Poland_Southeast_CordedWare.SG -0.739
gendstat: South_Africa_HG ZlatyKun -0.265
gendstat: South_Africa_HG BachoKiro 0.532
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.705
gendstat: South_Africa_HG DevilsGate_N.SG -0.201
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG -0.103
gendstat: South_Africa_HG Barcin_N.SG -1.123
gendstat: South_Africa_HG Iran_N.SG 0.199
gendstat: South_Africa_HG MA1 0.484
gendstat: South_Africa_HG Morocco_EN.SG 0.404
gendstat: South_Africa_HG CHG -0.189
gendstat: South_Africa_HG Italy_Mesolithic.SG 0.246
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -0.519
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -0.862
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -1.065
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG -1.303
gendstat: ZlatyKun BachoKiro 0.746
gendstat: ZlatyKun Eastern_Euro_UP.SG 0.893
gendstat: ZlatyKun DevilsGate_N.SG 0.088
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG 0.181
gendstat: ZlatyKun Barcin_N.SG -0.736
gendstat: ZlatyKun Iran_N.SG 0.455
gendstat: ZlatyKun MA1 0.652
gendstat: ZlatyKun Morocco_EN.SG 0.593
gendstat: ZlatyKun CHG 0.079
gendstat: ZlatyKun Italy_Mesolithic.SG 0.486
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG -0.210
gendstat: ZlatyKun Russia_BA_Okunevo.SG -0.486
gendstat: ZlatyKun Germany_EarlyMedieval.SG -0.658
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.892
gendstat: BachoKiro Eastern_Euro_UP.SG 0.168
gendstat: BachoKiro DevilsGate_N.SG -0.771
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -0.636
gendstat: BachoKiro Barcin_N.SG -1.628
gendstat: BachoKiro Iran_N.SG -0.355
gendstat: BachoKiro MA1 0.041
gendstat: BachoKiro Morocco_EN.SG -0.091
gendstat: BachoKiro CHG -0.659
gendstat: BachoKiro Italy_Mesolithic.SG -0.244
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -1.015
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.483
gendstat: BachoKiro Germany_EarlyMedieval.SG -1.648
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.928
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -0.931
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -0.903
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -2.005
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.553
gendstat: Eastern_Euro_UP.SG MA1 -0.102
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -0.247
gendstat: Eastern_Euro_UP.SG CHG -0.867
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG -0.456
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -1.302
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -1.843
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -2.053
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -2.291
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG 0.113
gendstat: DevilsGate_N.SG Barcin_N.SG -0.904
gendstat: DevilsGate_N.SG Iran_N.SG 0.424
gendstat: DevilsGate_N.SG MA1 0.663
gendstat: DevilsGate_N.SG Morocco_EN.SG 0.557
gendstat: DevilsGate_N.SG CHG -0.001
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.440
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG -0.339
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.747
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG -0.883
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG -1.147
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG -1.160
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 0.349
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 0.662
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG 0.502
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG -0.101
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.393
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG -0.493
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -1.054
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG -1.187
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG -1.515
gendstat: Barcin_N.SG Iran_N.SG 1.544
gendstat: Barcin_N.SG MA1 1.468
gendstat: Barcin_N.SG Morocco_EN.SG 1.426
gendstat: Barcin_N.SG CHG 0.937
gendstat: Barcin_N.SG Italy_Mesolithic.SG 1.435
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG 0.637
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG 0.397
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG 0.215
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG -0.193
gendstat: Iran_N.SG MA1 0.355
gendstat: Iran_N.SG Morocco_EN.SG 0.225
gendstat: Iran_N.SG CHG -0.416
gendstat: Iran_N.SG Italy_Mesolithic.SG 0.081
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -0.772
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -1.273
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -1.523
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -1.859
gendstat: MA1 Morocco_EN.SG -0.117
gendstat: MA1 CHG -0.637
gendstat: MA1 Italy_Mesolithic.SG -0.263
gendstat: MA1 Romania_IronGates_Mesolithic.SG -0.997
gendstat: MA1 Russia_BA_Okunevo.SG -1.373
gendstat: MA1 Germany_EarlyMedieval.SG -1.463
gendstat: MA1 Poland_Southeast_CordedWare.SG -1.722
gendstat: Morocco_EN.SG CHG -0.543
gendstat: Morocco_EN.SG Italy_Mesolithic.SG -0.134
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG -0.832
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG -1.221
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG -1.375
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG -1.601
gendstat: CHG Italy_Mesolithic.SG 0.428
gendstat: CHG Romania_IronGates_Mesolithic.SG -0.331
gendstat: CHG Russia_BA_Okunevo.SG -0.637
gendstat: CHG Germany_EarlyMedieval.SG -0.872
gendstat: CHG Poland_Southeast_CordedWare.SG -1.153
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -1.031
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -1.208
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -1.502
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -1.706
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.287
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG -0.536
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.858
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG -0.278
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG -0.679
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG -0.449




Thanks. Pretty similar overall to BedouinB. I guess not too surprising considering where they cluster on a PCA.

https://i.imgur.com/CQmY0XP.png

Mnemonics
02-17-2022, 07:40 PM
What? Haplotype as in Y-DNA/mtDNA?

It's qpAdm. It is a tool used to model autosomal ancestry; it has nothing to do with uniparental markers.

So nobody is omitting anything because qpAdm isn't even looking at the Y or mtDNA.

If you're interested in a particular sample or set of samples being modeled, look it up on the annotation file of the Reich Dataset. They usually have some basic Y-haplogroup info on Male samples.

https://reichdata.hms.harvard.edu/pub/datasets/amh_repo/curated_releases/V50/V50.0/SHARE/public.dir/v50.0_1240K_public.anno

@Mnemonics

Do you have the HO-dataset? Please run the same model for Yemeni Jews if you can. Too lazy to do it myself lol...

Btw, if you're interested, I've got samples from the "Genomic History of the Middle-East" in eigenstrat format with around 1 million coverage for the 1240k dataset.

Take a look at the Saudi samples with less SSA-admixture (SaudiA).

https://drive.google.com/file/d/1_cwsHtWFkCr-XM8EfgiArphlMNgpdY5k/view?usp=sharing


Those SaudiA samples seem to have a strong Mesopotamian shift, which probably is what is depleting their SSA. They also seem to need something more West African Shifted.

Without South_Africa_400BP in the model.






left pops:
SaudiA
BedouinB.DG
Jew_Iraqi.DG

right pops:
South_Africa_HG
Cameroon_SMA.SG
Mota
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 SaudiA 6
1 BedouinB.DG 2
2 Jew_Iraqi.DG 2
3 South_Africa_HG 4
4 Cameroon_SMA.SG 3
5 Mota 1
6 ZlatyKun 1
7 BachoKiro 3
8 Eastern_Euro_UP.SG 3
9 DevilsGate_N.SG 4
10 Kazakhstan_Botai_Eneolithic.SG 3
11 Barcin_N.SG 3
12 Iran_N.SG 4
13 MA1 1
14 Morocco_EN.SG 3
15 CHG 2
16 Italy_Mesolithic.SG 3
17 Romania_IronGates_Mesolithic.SG 3
18 Russia_BA_Okunevo.SG 14
19 Germany_EarlyMedieval.SG 21
20 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150445 indivs: 99
number of blocks for block jackknife: 714
## ncols: 1150445
coverage: SaudiA 955178
coverage: BedouinB.DG 1121254
coverage: Jew_Iraqi.DG 1121367
coverage: South_Africa_HG 1135129
coverage: Cameroon_SMA.SG 1149654
coverage: Mota 1134612
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 621.822
numsnps used: 447192
codimension 1
f4info:
f4rank: 1 dof: 16 chisq: 11.212 tail: 0.796189117 dofdiff: 18 chisqdiff: -11.212 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG 0.094
Mota 0.171
ZlatyKun 0.729
BachoKiro 0.783
Eastern_Euro_UP.SG 0.848
DevilsGate_N.SG 0.945
Kazakhstan_Botai_Eneolithic.SG 1.342
Barcin_N.SG 0.917
Iran_N.SG 1.613
MA1 1.037
Morocco_EN.SG -0.052
CHG 1.412
Italy_Mesolithic.SG 0.888
Romania_IronGates_Mesolithic.SG 0.915
Russia_BA_Okunevo.SG 1.110
Germany_EarlyMedieval.SG 1.092
Poland_Southeast_CordedWare.SG 1.358
A:
scale 1256.352
BedouinB.DG -0.881
Jew_Iraqi.DG 1.107


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 16 chisqdiff: 11.212 taildiff: 0.796189117
B:
scale 1398.583 1065.049
Cameroon_SMA.SG 0.202 0.257
Mota 0.037 0.280
ZlatyKun -0.407 0.926
BachoKiro -0.385 1.048
Eastern_Euro_UP.SG -0.727 0.911
DevilsGate_N.SG -0.953 0.920
Kazakhstan_Botai_Eneolithic.SG -1.324 1.297
Barcin_N.SG -1.052 0.799
Iran_N.SG -1.565 1.616
MA1 -1.048 1.002
Morocco_EN.SG -0.248 -0.200
CHG -1.554 1.271
Italy_Mesolithic.SG -0.772 0.998
Romania_IronGates_Mesolithic.SG -0.985 0.880
Russia_BA_Okunevo.SG -1.176 1.057
Germany_EarlyMedieval.SG -1.247 0.990
Poland_Southeast_CordedWare.SG -1.347 1.330
A:
scale 1.414 1.414
BedouinB.DG 1.414 0.000
Jew_Iraqi.DG 0.000 1.414


best coefficients: 0.557 0.443
Jackknife mean: 0.555320805 0.444679195
std. errors: 0.070 0.070

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


summ: SaudiA 2 0.796189 0.555 0.445 4858 -4858 4858

fixed pat wt dof chisq tail prob
00 0 16 11.212 0.796189 0.557 0.443
01 1 17 40.692 0.00103282 1.000 0.000
10 1 17 62.202 4.53022e-07 0.000 1.000
best pat: 00 0.796189 - -
best pat: 01 0.00103282 chi(nested): 29.479 p-value for nested model: 5.65176e-08

coeffs: 0.557 0.443

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

details: BedouinB.DG Cameroon_SMA.SG 0.000144 0.853408
details: Jew_Iraqi.DG Cameroon_SMA.SG 0.000242 1.504129
dscore: Cameroon_SMA.SG f4: 0.000187 Z: 1.438056

details: BedouinB.DG Mota 0.000026 0.120526
details: Jew_Iraqi.DG Mota 0.000263 1.219238
dscore: Mota f4: 0.000131 Z: 0.769001

details: BedouinB.DG ZlatyKun -0.000291 -0.983586
details: Jew_Iraqi.DG ZlatyKun 0.000869 3.072358
dscore: ZlatyKun f4: 0.000223 Z: 0.976828

details: BedouinB.DG BachoKiro -0.000275 -1.136714
details: Jew_Iraqi.DG BachoKiro 0.000984 3.927022
dscore: BachoKiro f4: 0.000283 Z: 1.462619

details: BedouinB.DG Eastern_Euro_UP.SG -0.000520 -1.912732
details: Jew_Iraqi.DG Eastern_Euro_UP.SG 0.000856 3.253421
dscore: Eastern_Euro_UP.SG f4: 0.000090 Z: 0.427584

details: BedouinB.DG DevilsGate_N.SG -0.000682 -2.561556
details: Jew_Iraqi.DG DevilsGate_N.SG 0.000864 3.287650
dscore: DevilsGate_N.SG f4: 0.000003 Z: 0.015493

details: BedouinB.DG Kazakhstan_Botai_Eneolithic.SG -0.000947 -3.644726
details: Jew_Iraqi.DG Kazakhstan_Botai_Eneolithic.SG 0.001217 4.829825
dscore: Kazakhstan_Botai_Eneolithic.SG f4: 0.000012 Z: 0.062439

details: BedouinB.DG Barcin_N.SG -0.000752 -2.737019
details: Jew_Iraqi.DG Barcin_N.SG 0.000750 3.041354
dscore: Barcin_N.SG f4: -0.000086 Z: -0.427726

details: BedouinB.DG Iran_N.SG -0.001119 -4.553794
details: Jew_Iraqi.DG Iran_N.SG 0.001517 6.462285
dscore: Iran_N.SG f4: 0.000049 Z: 0.262568

details: BedouinB.DG MA1 -0.000749 -2.308396
details: Jew_Iraqi.DG MA1 0.000940 2.998594
dscore: MA1 f4: -0.000001 Z: -0.002314

details: BedouinB.DG Morocco_EN.SG -0.000177 -0.649618
details: Jew_Iraqi.DG Morocco_EN.SG -0.000188 -0.686803
dscore: Morocco_EN.SG f4: -0.000182 Z: -0.867706

details: BedouinB.DG CHG -0.001111 -3.955875
details: Jew_Iraqi.DG CHG 0.001194 4.506456
dscore: CHG f4: -0.000090 Z: -0.435745

details: BedouinB.DG Italy_Mesolithic.SG -0.000552 -1.870786
details: Jew_Iraqi.DG Italy_Mesolithic.SG 0.000937 3.257735
dscore: Italy_Mesolithic.SG f4: 0.000108 Z: 0.476907

details: BedouinB.DG Romania_IronGates_Mesolithic.SG -0.000705 -2.554319
details: Jew_Iraqi.DG Romania_IronGates_Mesolithic.SG 0.000827 3.137806
dscore: Romania_IronGates_Mesolithic.SG f4: -0.000026 Z: -0.126128

details: BedouinB.DG Russia_BA_Okunevo.SG -0.000841 -3.641153
details: Jew_Iraqi.DG Russia_BA_Okunevo.SG 0.000992 4.303999
dscore: Russia_BA_Okunevo.SG f4: -0.000029 Z: -0.159261

details: BedouinB.DG Germany_EarlyMedieval.SG -0.000892 -3.638672
details: Jew_Iraqi.DG Germany_EarlyMedieval.SG 0.000929 4.189687
dscore: Germany_EarlyMedieval.SG f4: -0.000085 Z: -0.464289

details: BedouinB.DG Poland_Southeast_CordedWare.SG -0.000963 -3.965396
details: Jew_Iraqi.DG Poland_Southeast_CordedWare.SG 0.001248 5.426267
dscore: Poland_Southeast_CordedWare.SG f4: 0.000017 Z: 0.092370

gendstat: South_Africa_HG Cameroon_SMA.SG 1.438
gendstat: South_Africa_HG Mota 0.769
gendstat: South_Africa_HG ZlatyKun 0.977
gendstat: South_Africa_HG BachoKiro 1.463
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.428
gendstat: South_Africa_HG DevilsGate_N.SG 0.015
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG 0.062
gendstat: South_Africa_HG Barcin_N.SG -0.428
gendstat: South_Africa_HG Iran_N.SG 0.263
gendstat: South_Africa_HG MA1 -0.002
gendstat: South_Africa_HG Morocco_EN.SG -0.868
gendstat: South_Africa_HG CHG -0.436
gendstat: South_Africa_HG Italy_Mesolithic.SG 0.477
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -0.126
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -0.159
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -0.464
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG 0.092
gendstat: Cameroon_SMA.SG Mota -0.313
gendstat: Cameroon_SMA.SG ZlatyKun 0.153
gendstat: Cameroon_SMA.SG BachoKiro 0.459
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG -0.454
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.865
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.846
gendstat: Cameroon_SMA.SG Barcin_N.SG -1.286
gendstat: Cameroon_SMA.SG Iran_N.SG -0.709
gendstat: Cameroon_SMA.SG MA1 -0.753
gendstat: Cameroon_SMA.SG Morocco_EN.SG -1.686
gendstat: Cameroon_SMA.SG CHG -1.284
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG -0.342
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -1.002
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -1.123
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -1.408
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -0.875
gendstat: Mota ZlatyKun 0.379
gendstat: Mota BachoKiro 0.683
gendstat: Mota Eastern_Euro_UP.SG -0.174
gendstat: Mota DevilsGate_N.SG -0.554
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG -0.532
gendstat: Mota Barcin_N.SG -0.936
gendstat: Mota Iran_N.SG -0.394
gendstat: Mota MA1 -0.504
gendstat: Mota Morocco_EN.SG -1.357
gendstat: Mota CHG -0.988
gendstat: Mota Italy_Mesolithic.SG -0.093
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.669
gendstat: Mota Russia_BA_Okunevo.SG -0.757
gendstat: Mota Germany_EarlyMedieval.SG -1.020
gendstat: Mota Poland_Southeast_CordedWare.SG -0.539
gendstat: ZlatyKun BachoKiro 0.268
gendstat: ZlatyKun Eastern_Euro_UP.SG -0.565
gendstat: ZlatyKun DevilsGate_N.SG -0.966
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG -0.958
gendstat: ZlatyKun Barcin_N.SG -1.272
gendstat: ZlatyKun Iran_N.SG -0.799
gendstat: ZlatyKun MA1 -0.813
gendstat: ZlatyKun Morocco_EN.SG -1.587
gendstat: ZlatyKun CHG -1.323
gendstat: ZlatyKun Italy_Mesolithic.SG -0.465
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG -1.045
gendstat: ZlatyKun Russia_BA_Okunevo.SG -1.199
gendstat: ZlatyKun Germany_EarlyMedieval.SG -1.445
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.937
gendstat: BachoKiro Eastern_Euro_UP.SG -0.919
gendstat: BachoKiro DevilsGate_N.SG -1.447
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -1.388
gendstat: BachoKiro Barcin_N.SG -1.788
gendstat: BachoKiro Iran_N.SG -1.257
gendstat: BachoKiro MA1 -1.170
gendstat: BachoKiro Morocco_EN.SG -2.007
gendstat: BachoKiro CHG -1.755
gendstat: BachoKiro Italy_Mesolithic.SG -0.781
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -1.503
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.811
gendstat: BachoKiro Germany_EarlyMedieval.SG -2.024
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.438
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -0.428
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -0.410
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -0.918
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.210
gendstat: Eastern_Euro_UP.SG MA1 -0.382
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -1.202
gendstat: Eastern_Euro_UP.SG CHG -0.830
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG 0.090
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -0.600
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -0.702
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -1.044
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -0.416
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG 0.048
gendstat: DevilsGate_N.SG Barcin_N.SG -0.456
gendstat: DevilsGate_N.SG Iran_N.SG 0.245
gendstat: DevilsGate_N.SG MA1 -0.015
gendstat: DevilsGate_N.SG Morocco_EN.SG -0.805
gendstat: DevilsGate_N.SG CHG -0.436
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.479
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG -0.139
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.204
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG -0.487
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG 0.077
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG -0.531
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 0.214
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 -0.065
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG -0.901
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG -0.515
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.489
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG -0.226
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -0.332
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG -0.657
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG 0.032
gendstat: Barcin_N.SG Iran_N.SG 0.726
gendstat: Barcin_N.SG MA1 0.364
gendstat: Barcin_N.SG Morocco_EN.SG -0.427
gendstat: Barcin_N.SG CHG -0.018
gendstat: Barcin_N.SG Italy_Mesolithic.SG 0.987
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG 0.334
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG 0.357
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG 0.011
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG 0.683
gendstat: Iran_N.SG MA1 -0.224
gendstat: Iran_N.SG Morocco_EN.SG -1.014
gendstat: Iran_N.SG CHG -0.764
gendstat: Iran_N.SG Italy_Mesolithic.SG 0.278
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -0.389
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -0.501
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -0.864
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -0.206
gendstat: MA1 Morocco_EN.SG -0.690
gendstat: MA1 CHG -0.363
gendstat: MA1 Italy_Mesolithic.SG 0.455
gendstat: MA1 Romania_IronGates_Mesolithic.SG -0.118
gendstat: MA1 Russia_BA_Okunevo.SG -0.142
gendstat: MA1 Germany_EarlyMedieval.SG -0.411
gendstat: MA1 Poland_Southeast_CordedWare.SG 0.086
gendstat: Morocco_EN.SG CHG 0.384
gendstat: Morocco_EN.SG Italy_Mesolithic.SG 1.184
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG 0.702
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG 0.743
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG 0.484
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG 0.937
gendstat: CHG Italy_Mesolithic.SG 0.874
gendstat: CHG Romania_IronGates_Mesolithic.SG 0.304
gendstat: CHG Russia_BA_Okunevo.SG 0.334
gendstat: CHG Germany_EarlyMedieval.SG 0.030
gendstat: CHG Poland_Southeast_CordedWare.SG 0.615
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -0.794
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -0.751
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -1.158
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.505
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.016
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG -0.398
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG 0.275
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG -0.426
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG 0.354
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG 0.850







With South_Africa_400BP







left pops:
SaudiA
BedouinB.DG
Jew_Iraqi.DG

right pops:
South_Africa_HG
Cameroon_SMA.SG
South_Africa_400BP.SG
Mota
ZlatyKun
BachoKiro
Eastern_Euro_UP.SG
DevilsGate_N.SG
Kazakhstan_Botai_Eneolithic.SG
Barcin_N.SG
Iran_N.SG
MA1
Morocco_EN.SG
CHG
Italy_Mesolithic.SG
Romania_IronGates_Mesolithic.SG
Russia_BA_Okunevo.SG
Germany_EarlyMedieval.SG
Poland_Southeast_CordedWare.SG

0 SaudiA 6
1 BedouinB.DG 2
2 Jew_Iraqi.DG 2
3 South_Africa_HG 4
4 Cameroon_SMA.SG 3
5 South_Africa_400BP.SG 4
6 Mota 1
7 ZlatyKun 1
8 BachoKiro 3
9 Eastern_Euro_UP.SG 3
10 DevilsGate_N.SG 4
11 Kazakhstan_Botai_Eneolithic.SG 3
12 Barcin_N.SG 3
13 Iran_N.SG 4
14 MA1 1
15 Morocco_EN.SG 3
16 CHG 2
17 Italy_Mesolithic.SG 3
18 Romania_IronGates_Mesolithic.SG 3
19 Russia_BA_Okunevo.SG 14
20 Germany_EarlyMedieval.SG 21
21 Poland_Southeast_CordedWare.SG 13
jackknife block size: 0.050
snps: 1150450 indivs: 103
number of blocks for block jackknife: 714
## ncols: 1150450
coverage: SaudiA 955178
coverage: BedouinB.DG 1121254
coverage: Jew_Iraqi.DG 1121367
coverage: South_Africa_HG 1135129
coverage: Cameroon_SMA.SG 1149654
coverage: South_Africa_400BP.SG 1150180
coverage: Mota 1134612
coverage: ZlatyKun 1120502
coverage: BachoKiro 990561
coverage: Eastern_Euro_UP.SG 1150149
coverage: DevilsGate_N.SG 1134715
coverage: Kazakhstan_Botai_Eneolithic.SG 1149869
coverage: Barcin_N.SG 1150122
coverage: Iran_N.SG 1149571
coverage: MA1 796065
coverage: Morocco_EN.SG 653097
coverage: CHG 1134595
coverage: Italy_Mesolithic.SG 1141309
coverage: Romania_IronGates_Mesolithic.SG 1136384
coverage: Russia_BA_Okunevo.SG 1149711
coverage: Germany_EarlyMedieval.SG 1148477
coverage: Poland_Southeast_CordedWare.SG 1144291
Effective number of blocks: 621.897
numsnps used: 444038
codimension 1
f4info:
f4rank: 1 dof: 17 chisq: 25.326 tail: 0.087642507 dofdiff: 19 chisqdiff: -25.326 taildiff: 1
B:
scale 1.000
Cameroon_SMA.SG 0.102
South_Africa_400BP.SG 0.133
Mota 0.139
ZlatyKun 0.756
BachoKiro 0.795
Eastern_Euro_UP.SG 0.890
DevilsGate_N.SG 0.934
Kazakhstan_Botai_Eneolithic.SG 1.382
Barcin_N.SG 0.936
Iran_N.SG 1.703
MA1 1.020
Morocco_EN.SG -0.059
CHG 1.468
Italy_Mesolithic.SG 0.902
Romania_IronGates_Mesolithic.SG 0.917
Russia_BA_Okunevo.SG 1.169
Germany_EarlyMedieval.SG 1.118
Poland_Southeast_CordedWare.SG 1.388
A:
scale 1297.668
BedouinB.DG -0.851
Jew_Iraqi.DG 1.129


full rank
f4info:
f4rank: 2 dof: 0 chisq: 0.000 tail: 1 dofdiff: 17 chisqdiff: 25.326 taildiff: 0.087642507
B:
scale 1418.403 1087.581
Cameroon_SMA.SG 0.206 0.265
South_Africa_400BP.SG -0.515 -0.143
Mota 0.038 0.288
ZlatyKun -0.416 0.952
BachoKiro -0.393 1.078
Eastern_Euro_UP.SG -0.743 0.937
DevilsGate_N.SG -0.974 0.946
Kazakhstan_Botai_Eneolithic.SG -1.353 1.333
Barcin_N.SG -1.074 0.822
Iran_N.SG -1.599 1.661
MA1 -1.071 1.030
Morocco_EN.SG -0.253 -0.206
CHG -1.588 1.307
Italy_Mesolithic.SG -0.789 1.026
Romania_IronGates_Mesolithic.SG -1.007 0.905
Russia_BA_Okunevo.SG -1.201 1.087
Germany_EarlyMedieval.SG -1.274 1.018
Poland_Southeast_CordedWare.SG -1.376 1.367
A:
scale 1.414 1.414
BedouinB.DG 1.414 0.000
Jew_Iraqi.DG 0.000 1.414


best coefficients: 0.570 0.430
Jackknife mean: 0.566935754 0.433064246
std. errors: 0.085 0.085

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


summ: SaudiA 2 0.087643 0.567 0.433 7225 -7225 7225

fixed pat wt dof chisq tail prob
00 0 17 25.326 0.0876425 0.570 0.430
01 1 18 48.898 0.00011059 1.000 0.000
10 1 18 72.925 1.43674e-08 0.000 1.000
best pat: 00 0.0876425 - -
best pat: 01 0.00011059 chi(nested): 23.572 p-value for nested model: 1.20343e-06

coeffs: 0.570 0.430

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

details: BedouinB.DG Cameroon_SMA.SG 0.000145 0.852949
details: Jew_Iraqi.DG Cameroon_SMA.SG 0.000243 1.503825
dscore: Cameroon_SMA.SG f4: 0.000187 Z: 1.423038

details: BedouinB.DG South_Africa_400BP.SG -0.000363 -2.417435
details: Jew_Iraqi.DG South_Africa_400BP.SG -0.000131 -0.859768
dscore: South_Africa_400BP.SG f4: -0.000264 Z: -2.202786

details: BedouinB.DG Mota 0.000026 0.120514
details: Jew_Iraqi.DG Mota 0.000265 1.219135
dscore: Mota f4: 0.000129 Z: 0.748507

details: BedouinB.DG ZlatyKun -0.000293 -0.983551
details: Jew_Iraqi.DG ZlatyKun 0.000875 3.071859
dscore: ZlatyKun f4: 0.000209 Z: 0.906049

details: BedouinB.DG BachoKiro -0.000277 -1.136480
details: Jew_Iraqi.DG BachoKiro 0.000991 3.926246
dscore: BachoKiro f4: 0.000268 Z: 1.373213

details: BedouinB.DG Eastern_Euro_UP.SG -0.000524 -1.912454
details: Jew_Iraqi.DG Eastern_Euro_UP.SG 0.000862 3.253011
dscore: Eastern_Euro_UP.SG f4: 0.000072 Z: 0.338738

details: BedouinB.DG DevilsGate_N.SG -0.000686 -2.561453
details: Jew_Iraqi.DG DevilsGate_N.SG 0.000870 3.287506
dscore: DevilsGate_N.SG f4: -0.000018 Z: -0.086081

details: BedouinB.DG Kazakhstan_Botai_Eneolithic.SG -0.000954 -3.644336
details: Jew_Iraqi.DG Kazakhstan_Botai_Eneolithic.SG 0.001226 4.829242
dscore: Kazakhstan_Botai_Eneolithic.SG f4: -0.000017 Z: -0.085165

details: BedouinB.DG Barcin_N.SG -0.000757 -2.737116
details: Jew_Iraqi.DG Barcin_N.SG 0.000756 3.041351
dscore: Barcin_N.SG f4: -0.000107 Z: -0.524969

details: BedouinB.DG Iran_N.SG -0.001127 -4.553996
details: Jew_Iraqi.DG Iran_N.SG 0.001528 6.461914
dscore: Iran_N.SG f4: 0.000014 Z: 0.073355

details: BedouinB.DG MA1 -0.000755 -2.308265
details: Jew_Iraqi.DG MA1 0.000947 2.998225
dscore: MA1 f4: -0.000023 Z: -0.094952

details: BedouinB.DG Morocco_EN.SG -0.000179 -0.649544
details: Jew_Iraqi.DG Morocco_EN.SG -0.000189 -0.686734
dscore: Morocco_EN.SG f4: -0.000183 Z: -0.864891

details: BedouinB.DG CHG -0.001119 -3.955372
details: Jew_Iraqi.DG CHG 0.001202 4.506281
dscore: CHG f4: -0.000122 Z: -0.582840

details: BedouinB.DG Italy_Mesolithic.SG -0.000556 -1.870702
details: Jew_Iraqi.DG Italy_Mesolithic.SG 0.000943 3.257365
dscore: Italy_Mesolithic.SG f4: 0.000088 Z: 0.387395

details: BedouinB.DG Romania_IronGates_Mesolithic.SG -0.000710 -2.554088
details: Jew_Iraqi.DG Romania_IronGates_Mesolithic.SG 0.000833 3.137206
dscore: Romania_IronGates_Mesolithic.SG f4: -0.000047 Z: -0.224723

details: BedouinB.DG Russia_BA_Okunevo.SG -0.000847 -3.640994
details: Jew_Iraqi.DG Russia_BA_Okunevo.SG 0.000999 4.303239
dscore: Russia_BA_Okunevo.SG f4: -0.000054 Z: -0.295252

details: BedouinB.DG Germany_EarlyMedieval.SG -0.000898 -3.638559
details: Jew_Iraqi.DG Germany_EarlyMedieval.SG 0.000936 4.189115
dscore: Germany_EarlyMedieval.SG f4: -0.000110 Z: -0.595239

details: BedouinB.DG Poland_Southeast_CordedWare.SG -0.000970 -3.965327
details: Jew_Iraqi.DG Poland_Southeast_CordedWare.SG 0.001257 5.425342
dscore: Poland_Southeast_CordedWare.SG f4: -0.000013 Z: -0.068829

gendstat: South_Africa_HG Cameroon_SMA.SG 1.423
gendstat: South_Africa_HG South_Africa_400BP.SG -2.203
gendstat: South_Africa_HG Mota 0.749
gendstat: South_Africa_HG ZlatyKun 0.906
gendstat: South_Africa_HG BachoKiro 1.373
gendstat: South_Africa_HG Eastern_Euro_UP.SG 0.339
gendstat: South_Africa_HG DevilsGate_N.SG -0.086
gendstat: South_Africa_HG Kazakhstan_Botai_Eneolithic.SG -0.085
gendstat: South_Africa_HG Barcin_N.SG -0.525
gendstat: South_Africa_HG Iran_N.SG 0.073
gendstat: South_Africa_HG MA1 -0.095
gendstat: South_Africa_HG Morocco_EN.SG -0.865
gendstat: South_Africa_HG CHG -0.583
gendstat: South_Africa_HG Italy_Mesolithic.SG 0.387
gendstat: South_Africa_HG Romania_IronGates_Mesolithic.SG -0.225
gendstat: South_Africa_HG Russia_BA_Okunevo.SG -0.295
gendstat: South_Africa_HG Germany_EarlyMedieval.SG -0.595
gendstat: South_Africa_HG Poland_Southeast_CordedWare.SG -0.069
gendstat: Cameroon_SMA.SG South_Africa_400BP.SG -3.559
gendstat: Cameroon_SMA.SG Mota -0.323
gendstat: Cameroon_SMA.SG ZlatyKun 0.092
gendstat: Cameroon_SMA.SG BachoKiro 0.384
gendstat: Cameroon_SMA.SG Eastern_Euro_UP.SG -0.532
gendstat: Cameroon_SMA.SG DevilsGate_N.SG -0.953
gendstat: Cameroon_SMA.SG Kazakhstan_Botai_Eneolithic.SG -0.976
gendstat: Cameroon_SMA.SG Barcin_N.SG -1.368
gendstat: Cameroon_SMA.SG Iran_N.SG -0.880
gendstat: Cameroon_SMA.SG MA1 -0.835
gendstat: Cameroon_SMA.SG Morocco_EN.SG -1.674
gendstat: Cameroon_SMA.SG CHG -1.415
gendstat: Cameroon_SMA.SG Italy_Mesolithic.SG -0.421
gendstat: Cameroon_SMA.SG Romania_IronGates_Mesolithic.SG -1.088
gendstat: Cameroon_SMA.SG Russia_BA_Okunevo.SG -1.240
gendstat: Cameroon_SMA.SG Germany_EarlyMedieval.SG -1.521
gendstat: Cameroon_SMA.SG Poland_Southeast_CordedWare.SG -1.017
gendstat: South_Africa_400BP.SG Mota 2.390
gendstat: South_Africa_400BP.SG ZlatyKun 2.116
gendstat: South_Africa_400BP.SG BachoKiro 2.774
gendstat: South_Africa_400BP.SG Eastern_Euro_UP.SG 1.611
gendstat: South_Africa_400BP.SG DevilsGate_N.SG 1.221
gendstat: South_Africa_400BP.SG Kazakhstan_Botai_Eneolithic.SG 1.278
gendstat: South_Africa_400BP.SG Barcin_N.SG 0.762
gendstat: South_Africa_400BP.SG Iran_N.SG 1.493
gendstat: South_Africa_400BP.SG MA1 0.996
gendstat: South_Africa_400BP.SG Morocco_EN.SG 0.394
gendstat: South_Africa_400BP.SG CHG 0.694
gendstat: South_Africa_400BP.SG Italy_Mesolithic.SG 1.570
gendstat: South_Africa_400BP.SG Romania_IronGates_Mesolithic.SG 1.051
gendstat: South_Africa_400BP.SG Russia_BA_Okunevo.SG 1.176
gendstat: South_Africa_400BP.SG Germany_EarlyMedieval.SG 0.836
gendstat: South_Africa_400BP.SG Poland_Southeast_CordedWare.SG 1.366
gendstat: Mota ZlatyKun 0.327
gendstat: Mota BachoKiro 0.620
gendstat: Mota Eastern_Euro_UP.SG -0.237
gendstat: Mota DevilsGate_N.SG -0.628
gendstat: Mota Kazakhstan_Botai_Eneolithic.SG -0.646
gendstat: Mota Barcin_N.SG -1.004
gendstat: Mota Iran_N.SG -0.546
gendstat: Mota MA1 -0.577
gendstat: Mota Morocco_EN.SG -1.339
gendstat: Mota CHG -1.107
gendstat: Mota Italy_Mesolithic.SG -0.159
gendstat: Mota Romania_IronGates_Mesolithic.SG -0.740
gendstat: Mota Russia_BA_Okunevo.SG -0.855
gendstat: Mota Germany_EarlyMedieval.SG -1.116
gendstat: Mota Poland_Southeast_CordedWare.SG -0.661
gendstat: ZlatyKun BachoKiro 0.261
gendstat: ZlatyKun Eastern_Euro_UP.SG -0.575
gendstat: ZlatyKun DevilsGate_N.SG -0.986
gendstat: ZlatyKun Kazakhstan_Botai_Eneolithic.SG -1.016
gendstat: ZlatyKun Barcin_N.SG -1.287
gendstat: ZlatyKun Iran_N.SG -0.888
gendstat: ZlatyKun MA1 -0.836
gendstat: ZlatyKun Morocco_EN.SG -1.522
gendstat: ZlatyKun CHG -1.384
gendstat: ZlatyKun Italy_Mesolithic.SG -0.482
gendstat: ZlatyKun Romania_IronGates_Mesolithic.SG -1.063
gendstat: ZlatyKun Russia_BA_Okunevo.SG -1.239
gendstat: ZlatyKun Germany_EarlyMedieval.SG -1.482
gendstat: ZlatyKun Poland_Southeast_CordedWare.SG -0.998
gendstat: BachoKiro Eastern_Euro_UP.SG -0.924
gendstat: BachoKiro DevilsGate_N.SG -1.464
gendstat: BachoKiro Kazakhstan_Botai_Eneolithic.SG -1.448
gendstat: BachoKiro Barcin_N.SG -1.800
gendstat: BachoKiro Iran_N.SG -1.354
gendstat: BachoKiro MA1 -1.192
gendstat: BachoKiro Morocco_EN.SG -1.930
gendstat: BachoKiro CHG -1.818
gendstat: BachoKiro Italy_Mesolithic.SG -0.794
gendstat: BachoKiro Romania_IronGates_Mesolithic.SG -1.518
gendstat: BachoKiro Russia_BA_Okunevo.SG -1.853
gendstat: BachoKiro Germany_EarlyMedieval.SG -2.062
gendstat: BachoKiro Poland_Southeast_CordedWare.SG -1.506
gendstat: Eastern_Euro_UP.SG DevilsGate_N.SG -0.438
gendstat: Eastern_Euro_UP.SG Kazakhstan_Botai_Eneolithic.SG -0.465
gendstat: Eastern_Euro_UP.SG Barcin_N.SG -0.924
gendstat: Eastern_Euro_UP.SG Iran_N.SG -0.296
gendstat: Eastern_Euro_UP.SG MA1 -0.399
gendstat: Eastern_Euro_UP.SG Morocco_EN.SG -1.117
gendstat: Eastern_Euro_UP.SG CHG -0.886
gendstat: Eastern_Euro_UP.SG Italy_Mesolithic.SG 0.082
gendstat: Eastern_Euro_UP.SG Romania_IronGates_Mesolithic.SG -0.609
gendstat: Eastern_Euro_UP.SG Russia_BA_Okunevo.SG -0.736
gendstat: Eastern_Euro_UP.SG Germany_EarlyMedieval.SG -1.076
gendstat: Eastern_Euro_UP.SG Poland_Southeast_CordedWare.SG -0.479
gendstat: DevilsGate_N.SG Kazakhstan_Botai_Eneolithic.SG 0.004
gendstat: DevilsGate_N.SG Barcin_N.SG -0.452
gendstat: DevilsGate_N.SG Iran_N.SG 0.167
gendstat: DevilsGate_N.SG MA1 -0.023
gendstat: DevilsGate_N.SG Morocco_EN.SG -0.712
gendstat: DevilsGate_N.SG CHG -0.483
gendstat: DevilsGate_N.SG Italy_Mesolithic.SG 0.482
gendstat: DevilsGate_N.SG Romania_IronGates_Mesolithic.SG -0.138
gendstat: DevilsGate_N.SG Russia_BA_Okunevo.SG -0.229
gendstat: DevilsGate_N.SG Germany_EarlyMedieval.SG -0.506
gendstat: DevilsGate_N.SG Poland_Southeast_CordedWare.SG 0.027
gendstat: Kazakhstan_Botai_Eneolithic.SG Barcin_N.SG -0.482
gendstat: Kazakhstan_Botai_Eneolithic.SG Iran_N.SG 0.177
gendstat: Kazakhstan_Botai_Eneolithic.SG MA1 -0.032
gendstat: Kazakhstan_Botai_Eneolithic.SG Morocco_EN.SG -0.764
gendstat: Kazakhstan_Botai_Eneolithic.SG CHG -0.524
gendstat: Kazakhstan_Botai_Eneolithic.SG Italy_Mesolithic.SG 0.533
gendstat: Kazakhstan_Botai_Eneolithic.SG Romania_IronGates_Mesolithic.SG -0.175
gendstat: Kazakhstan_Botai_Eneolithic.SG Russia_BA_Okunevo.SG -0.295
gendstat: Kazakhstan_Botai_Eneolithic.SG Germany_EarlyMedieval.SG -0.624
gendstat: Kazakhstan_Botai_Eneolithic.SG Poland_Southeast_CordedWare.SG 0.028
gendstat: Barcin_N.SG Iran_N.SG 0.643
gendstat: Barcin_N.SG MA1 0.352
gendstat: Barcin_N.SG Morocco_EN.SG -0.336
gendstat: Barcin_N.SG CHG -0.071
gendstat: Barcin_N.SG Italy_Mesolithic.SG 0.987
gendstat: Barcin_N.SG Romania_IronGates_Mesolithic.SG 0.332
gendstat: Barcin_N.SG Russia_BA_Okunevo.SG 0.329
gendstat: Barcin_N.SG Germany_EarlyMedieval.SG -0.019
gendstat: Barcin_N.SG Poland_Southeast_CordedWare.SG 0.619
gendstat: Iran_N.SG MA1 -0.166
gendstat: Iran_N.SG Morocco_EN.SG -0.857
gendstat: Iran_N.SG CHG -0.738
gendstat: Iran_N.SG Italy_Mesolithic.SG 0.350
gendstat: Iran_N.SG Romania_IronGates_Mesolithic.SG -0.312
gendstat: Iran_N.SG Russia_BA_Okunevo.SG -0.432
gendstat: Iran_N.SG Germany_EarlyMedieval.SG -0.792
gendstat: Iran_N.SG Poland_Southeast_CordedWare.SG -0.169
gendstat: MA1 Morocco_EN.SG -0.602
gendstat: MA1 CHG -0.396
gendstat: MA1 Italy_Mesolithic.SG 0.465
gendstat: MA1 Romania_IronGates_Mesolithic.SG -0.108
gendstat: MA1 Russia_BA_Okunevo.SG -0.152
gendstat: MA1 Germany_EarlyMedieval.SG -0.419
gendstat: MA1 Poland_Southeast_CordedWare.SG 0.052
gendstat: Morocco_EN.SG CHG 0.254
gendstat: Morocco_EN.SG Italy_Mesolithic.SG 1.100
gendstat: Morocco_EN.SG Romania_IronGates_Mesolithic.SG 0.608
gendstat: Morocco_EN.SG Russia_BA_Okunevo.SG 0.622
gendstat: Morocco_EN.SG Germany_EarlyMedieval.SG 0.361
gendstat: Morocco_EN.SG Poland_Southeast_CordedWare.SG 0.795
gendstat: CHG Italy_Mesolithic.SG 0.921
gendstat: CHG Romania_IronGates_Mesolithic.SG 0.353
gendstat: CHG Russia_BA_Okunevo.SG 0.368
gendstat: CHG Germany_EarlyMedieval.SG 0.067
gendstat: CHG Poland_Southeast_CordedWare.SG 0.621
gendstat: Italy_Mesolithic.SG Romania_IronGates_Mesolithic.SG -0.796
gendstat: Italy_Mesolithic.SG Russia_BA_Okunevo.SG -0.775
gendstat: Italy_Mesolithic.SG Germany_EarlyMedieval.SG -1.183
gendstat: Italy_Mesolithic.SG Poland_Southeast_CordedWare.SG -0.559
gendstat: Romania_IronGates_Mesolithic.SG Russia_BA_Okunevo.SG -0.042
gendstat: Romania_IronGates_Mesolithic.SG Germany_EarlyMedieval.SG -0.423
gendstat: Romania_IronGates_Mesolithic.SG Poland_Southeast_CordedWare.SG 0.216
gendstat: Russia_BA_Okunevo.SG Germany_EarlyMedieval.SG -0.424
gendstat: Russia_BA_Okunevo.SG Poland_Southeast_CordedWare.SG 0.314
gendstat: Germany_EarlyMedieval.SG Poland_Southeast_CordedWare.SG 0.805

Gentica277282
02-17-2022, 08:15 PM
I tried to model Saudi with dzu and ana still a poor fit for them

Target: Saudi
Distance: 4.3824% / 0.04382357
67.0 Dzudzuana
26.4 Iran_Neolithic
5.4 Ancestral_North-African(DZU-55%-45%)
0.8 Early_European_Farmer
0.4 KEN_Kakapel_900BP

Mahra a better fit:

Target: Yemenite_Mahra
Distance: 3.6397% / 0.03639744
70.0 Dzudzuana
24.2 Iran_Neolithic
5.8 Ancestral_North-African(DZU-55%-45%)

Bedouin B same issue as Saudi:
Target: BedouinB
Distance: 4.2875% / 0.04287521
66.0 Dzudzuana
24.8 Iran_Neolithic
6.0 Ancestral_North-African(DZU-55%-45%)
2.2 Early_European_Farmer
1.0 KEN_Kakapel_900BP