Page 3 of 3 FirstFirst 123
Results 21 to 28 of 28

Thread: How did Villabruna get its ANE related ancestry?

  1. #21
    Registered Users
    Posts
    1,267
    Sex
    Omitted

    What uniparental markers would you suggest accompany the West-Asian/Basal ancestry in Sunghir/Kostenki/Gravettians/WHG and why?
    Collection of 14,000 d-stats: Hidden Content Part 2: Hidden Content Part 3: Hidden Content PM me for d-stats, qpadm, qpgraph, or f3-outgroup nmonte models.

  2. #22
    Registered Users
    Posts
    1,106
    Sex
    Omitted

    Quote Originally Posted by Chad Rohlfsen View Post
    https://populationgenomics.blog/2019...cient-eurasia/

    Look at the above. Best fit is Ust-Ishim as early West Eurasian, with later West Eurasians having flow from East Eurasians. This groups splits, with the group leading to MA1 admixing into Native Americans and the European branch forming the bulk of "WHG".

    I do have Yana now, so I can go back to this and also add Iron Gates, etc..
    OK, that would be the deepest split. But after the LGM we have three clearly distinctive branches. How did the ANE branch reach Villabruna? There are two ways to explain how WHG came into existence. A merger locally in a south European refuge or as some sort of repopulation from further east. I think that ANE in Villabruna needs an explanation.

  3. #23
    Registered Users
    Posts
    1,106
    Sex
    Omitted

    Quote Originally Posted by Chad Rohlfsen View Post
    https://populationgenomics.blog/2019...cient-eurasia/

    Look at the above. Best fit is Ust-Ishim as early West Eurasian, with later West Eurasians having flow from East Eurasians. This groups splits, with the group leading to MA1 admixing into Native Americans and the European branch forming the bulk of "WHG".

    I do have Yana now, so I can go back to this and also add Iron Gates, etc..
    Another thing. What do you think of the claim in this paper?

    https://genome.cshlp.org/content/27/1/1.full

    We found that Siberian and East Asian populations shared 38% of their ancestry with a 45,000-yr-old Ustí-Ishim individual who was previously believed to have no modern-day descendants.

  4. #24
    Registered Users
    Posts
    155
    Sex
    Location
    Haifa
    Ethnicity
    Jewish & Slavic
    Nationality
    Israel & Russia
    Y-DNA
    I2a1b - S17250
    mtDNA
    R0a4

    Israel Israel Jerusalem Russian Federation Serbia Montenegro
    Quote Originally Posted by epoch View Post
    OK, that would be the deepest split. But after the LGM we have three clearly distinctive branches. How did the ANE branch reach Villabruna? There are two ways to explain how WHG came into existence. A merger locally in a south European refuge or as some sort of repopulation from further east. I think that ANE in Villabruna needs an explanation.
    I guess R1b in Villabruna is a clear sign of ANE ancestry. In mesolithic Balkans there is more R1b than I2. So, I guess ANE got to Villabruna through Balkans, and we will see an ANE ancestry cline from Pontic-Caspian Steppe to Balkans and further to Appenines.
    We need some samples from what is now Pontic-Caspian Steppe and from Balkans, older than 12 000 BCE to say something more definite.

  5. The Following User Says Thank You to artemv For This Useful Post:

     parasar (06-12-2019)

  6. #25
    Registered Users
    Posts
    1,267
    Sex
    Omitted

    To follow up from earlier. A hybrid Eastern population (Ust_Ishim + East-Eurasian-proper) mixing into paleolithic Europe turns out some pretty weird results, probably not workable. Think the main issue is that Ust_Ishim just doesn't have enough drift, which ends up requiring paleolithic Europeans to have exhuberant amounts of Ust_Ishim ancestry (on the order of 70-75%).

    EDIT: Trying a variation of Chad's idea now. Ust_Ishim as early West Eurasian and Goyet being the lone paleo-European without basal.
    Last edited by Kale; 06-12-2019 at 04:39 AM.
    Collection of 14,000 d-stats: Hidden Content Part 2: Hidden Content Part 3: Hidden Content PM me for d-stats, qpadm, qpgraph, or f3-outgroup nmonte models.

  7. The Following 3 Users Say Thank You to Kale For This Useful Post:

     epoch (06-12-2019),  palamede (06-12-2019),  parasar (06-12-2019)

  8. #26
    Registered Users
    Posts
    1,106
    Sex
    Omitted

    Quote Originally Posted by Kale View Post
    To follow up from earlier. A hybrid Eastern population (Ust_Ishim + East-Eurasian-proper) mixing into paleolithic Europe turns out some pretty weird results, probably not workable. Think the main issue is that Ust_Ishim just doesn't have enough drift, which ends up requiring paleolithic Europeans to have exhuberant amounts of Ust_Ishim ancestry (on the order of 70-75%).

    EDIT: Trying a variation of Chad's idea now. Ust_Ishim as early West Eurasian and Goyet being the lone paleo-European without basal.
    Remember the difference between ElMiron and Chan_Meso you thought you saw?

    Quote Originally Posted by Kale View Post
    It's interesting to compare ElMiron and Chan_Meso

    They don't differ significantly in relation to... Barcin_N, Levant_N, Iran_N, Taforalt, MA1, Native Americans, East Eurasians, Ust_Ishim, GoyetQ116-1, Gravettians, Kostenki14, Sunghir
    Maybe significant differences to (I'm looking at outgroup F3 rather than D-stats so forgive the uncertainty)... CHG and EHG, with both closer to Chan
    And WHG are significantly more related to Chan, even ones with minimal Goyet ancestry.

    Basically they look identical in every way except Chan's WHG is more 'WHG' than ElMiron's WHG, if that makes any sense?
    Maybe to build on what Epoch said, ElMiron got a shot of proto-WHG Solutrean, whereas Chan got WHG-proper?
    It is quite interesting to see if ANE is the difference. If El_miron is half WHG it should also have substantial ANE admixture. If it's not there we have found a clue for dating the admixture.

  9. #27
    Registered Users
    Posts
    1,267
    Sex
    Omitted

    ANE amongst WHG increases over time and proximity to Eastern Europe. It was gradual, not punctuated. Compared to Villabruna, ElMiron is both earlier, farther from Eastern Europe, and mixed with Goyet. It would be difficult to determine if it has ANE ancestry because if it does, it's <5%, but I would venture to say probably.
    Collection of 14,000 d-stats: Hidden Content Part 2: Hidden Content Part 3: Hidden Content PM me for d-stats, qpadm, qpgraph, or f3-outgroup nmonte models.

  10. The Following User Says Thank You to Kale For This Useful Post:

     Pribislav (06-12-2019)

  11. #28
    Registered Users
    Posts
    1,267
    Sex
    Omitted

    Well, Chad's basic premise of Ust_Ishim as early West Eurasian and Goyet being the lone paleo-European without basal spliced well enough into the other working phylogeny I had.
    To see the graph, paste the contents of the dotfile spoiler at Viz.js

    dotfile
     

    digraph G {
    size = "7.5,10" ;
    labelloc = "t" ;
    label = "61119 :: Sun Gra Taf Bar 0.004486 -0.000909 -0.005395 0.001266 -4.263
    " ;

    Mbuti_DG [ label = "Mbuti_DG" ] ;
    Mota [ label = "Mota" ] ;
    Taforalt [ label = "Taforalt" ] ;
    Barcin_N [ label = "Barcin_N" ] ;
    Onge_DG [ label = "Onge_DG" ] ;
    Ust_Ishim [ label = "Ust_Ishim" ] ;
    Dai_DG [ label = "Dai_DG" ] ;
    GoyetQ116_1 [ label = "GoyetQ116_1" ] ;
    Kostenki14 [ label = "Kostenki14" ] ;
    Sunghir [ label = "Sunghir" ] ;
    Gravettian [ label = "Gravettian" ] ;
    Brazil_Sumidouro_10100BP [ label = "Brazil_Sumidouro_10100BP" ] ;
    Baikal_EN [ label = "Baikal_EN" ] ;
    EHG [ label = "EHG" ] ;
    Loschbour [ label = "Loschbour" ] ;
    Iran_N [ label = "Iran_N" ] ;
    CHG [ label = "CHG" ] ;
    Root -> Mbuti_DG [ label = "36" ] ;
    Root -> A [ label = "36" ] ;
    A -> Mota [ label = "494" ] ;
    A -> B [ label = "48" ] ;
    B -> ANA [ label = "120" ] ;
    B -> C [ label = "30" ] ;
    C -> Basal1 [ label = "30" ] ;
    C -> Eurasia [ label = "28" ] ;
    Eurasia -> E1 [ label = "24" ] ;
    Eurasia -> W1 [ label = "3" ] ;
    pre_IBM2 -> Taforalt [ label = "146" ] ;
    ANF2 -> Barcin_N [ label = "39" ] ;
    E1 -> Onge_DG [ label = "91" ] ;
    E1 -> E2 [ label = "30" ] ;
    W1 -> Ust_Ishim [ label = "9" ] ;
    W1 -> W2 [ label = "38" ] ;
    W3 -> W4 [ label = "7" ] ;
    NE1 -> NE2 [ label = "47" ] ;
    E2 -> Dai_DG [ label = "15" ] ;
    E2 -> E3 [ label = "12" ] ;
    W2 -> W2a [ label = "216" ] ;
    W4 -> W5 [ label = "12" ] ;
    W4 -> WHG0 [ label = "26" ] ;
    NE2 -> NE3 [ label = "19" ] ;
    NE2 -> NE4 [ label = "28" ] ;
    W2a -> GoyetQ116_1 [ label = "220" ] ;
    W5 -> Kostenki14 [ label = "425" ] ;
    W5 -> W6 [ label = "30" ] ;
    pre_ANE -> ANE [ label = "114" ] ;
    W6 -> Sunghir [ label = "66" ] ;
    pre_GRV -> Gravettian [ label = "282" ] ;
    America -> Brazil_Sumidouro_10100BP [ label = "185" ] ;
    pre_BKN -> Baikal_EN [ label = "47" ] ;
    WHG1 -> WHG2 [ label = "10" ] ;
    WHG3 -> WHG4 [ label = "27" ] ;
    pre_EHG2 -> EHG [ label = "117" ] ;
    WHG2 -> WHG5 [ label = "43" ] ;
    WHG4 -> Loschbour [ label = "57" ] ;
    INF2 -> Iran_N [ label = "56" ] ;
    pre_CHG2 -> CHG [ label = "256" ] ;
    ANA -> pre_IBM2 [ style=dotted, label = "70%" ] ;
    pre_IBM1 -> pre_IBM2 [ style=dotted, label = "30%" ] ;
    ANA -> ANF2 [ style=dotted, label = "10%" ] ;
    ANF1 -> ANF2 [ style=dotted, label = "90%" ] ;
    Basal1 -> W3 [ style=dotted, label = "15%" ] ;
    W2 -> W3 [ style=dotted, label = "85%" ] ;
    Basal1 -> NE1 [ style=dotted, label = "69%" ] ;
    E1 -> NE1 [ style=dotted, label = "31%" ] ;
    E2 -> pre_ANE [ style=dotted, label = "30%" ] ;
    W4 -> pre_ANE [ style=dotted, label = "70%" ] ;
    W6 -> pre_GRV [ style=dotted, label = "41%" ] ;
    WHG0 -> pre_GRV [ style=dotted, label = "59%" ] ;
    ANE -> America [ style=dotted, label = "37%" ] ;
    E3 -> America [ style=dotted, label = "63%" ] ;
    ANE -> pre_BKN [ style=dotted, label = "12%" ] ;
    E3 -> pre_BKN [ style=dotted, label = "88%" ] ;
    ANE -> WHG1 [ style=dotted, label = "19%" ] ;
    WHG0 -> WHG1 [ style=dotted, label = "81%" ] ;
    W2a -> WHG3 [ style=dotted, label = "4%" ] ;
    WHG2 -> WHG3 [ style=dotted, label = "96%" ] ;
    W2a -> pre_IBM1 [ style=dotted, label = "5%" ] ;
    WHG2 -> pre_IBM1 [ style=dotted, label = "95%" ] ;
    NE4 -> ANF1 [ style=dotted, label = "28%" ] ;
    WHG5 -> ANF1 [ style=dotted, label = "72%" ] ;
    NE4 -> pre_EHG2 [ style=dotted, label = "17%" ] ;
    pre_EHG1 -> pre_EHG2 [ style=dotted, label = "83%" ] ;
    ANE -> INF1 [ style=dotted, label = "13%" ] ;
    NE3 -> INF1 [ style=dotted, label = "87%" ] ;
    ANE -> pre_CHG1 [ style=dotted, label = "16%" ] ;
    NE3 -> pre_CHG1 [ style=dotted, label = "84%" ] ;
    ANE -> pre_EHG1 [ style=dotted, label = "56%" ] ;
    WHG4 -> pre_EHG1 [ style=dotted, label = "44%" ] ;
    INF1 -> INF2 [ style=dotted, label = "77%" ] ;
    WHG5 -> INF2 [ style=dotted, label = "23%" ] ;
    pre_CHG1 -> pre_CHG2 [ style=dotted, label = "62%" ] ;
    WHG5 -> pre_CHG2 [ style=dotted, label = "38%" ] ;
    }


    terminal output - further information
     

    parameter file: parqpgraph
    ### THE INPUT PARAMETERS
    ##PARAMETER NAME: VALUE
    DIR: /home/user/AdmixTools-master/bin
    S1: dataset19d
    indivname: /home/user/AdmixTools-master/bin/dataset19d.ind
    snpname: /home/user/AdmixTools-master/bin/dataset19d.snp
    genotypename: /home/user/AdmixTools-master/bin/dataset19d.geno
    outpop: Mbuti.DG
    blgsize: 0.05
    lsqmode: YES
    diag: .0001
    hires: YES
    initmix: 1000
    precision: .0001
    zthresh: 3.0
    terse: NO
    useallsnps: NO
    ## qpGraph version: 6450
    outpop: Mbuti.DG
    seed: 1616053385
    simple lsqmode
    Mbuti.DG
    Mota
    Onge.DG
    Dai.DG
    Ust_Ishim
    GoyetQ116-1
    Kostenki14
    Sunghir
    Gravettian
    Brazil_Sumidouro_10100BP
    Baikal_EN
    Loschbour
    Taforalt
    Barcin_N
    Iran_N
    EHG
    CHG

    graph: 61119


    numedge: 40 numancestor: 1
    0 Mbuti.DG
    1 Mota
    2 Onge.DG
    3 Dai.DG
    4 Ust_Ishim
    5 GoyetQ116-1
    6 Kostenki14
    7 Sunghir
    8 Gravettian
    9 Brazil_Sumidouro_10100BP
    10 Baikal_EN
    11 Loschbour
    12 Taforalt
    13 Barcin_N
    14 Iran_N
    15 EHG
    16 CHG
    root label: :Mbuti.DG
    nodata: rs61780219
    nodata: rs71629813
    nodata: rs12077803
    nodata: Affx-50004616
    nodata: Affx-50006276
    nodata: rs2799108
    nodata: rs17159675
    nodata: rs61816207
    nodata: Affx-50009874
    nodata: Affx-50009934
    nodata: rs4500359
    nodata: Affx-50009936
    nodata: Affx-50009937
    nodata: Affx-50066782
    nodata: rs142953672
    nodata: Affx-50066783
    nodata: Affx-50069565
    nodata: Affx-50066086
    nodata: Affx-50066092
    nodata: rs451002
    nodata: rs62162900
    nodata: rs71224406
    nodata: Affx-50066302
    nodata: Affx-50070296
    nodata: Affx-50070298
    nodata: Affx-50070299
    nodata: rs1347870
    nodata: Affx-50070300
    nodata: rs62164875
    nodata: Affx-50070303
    nodata: Affx-50086450
    nodata: rs138840766
    nodata: Affx-50086444
    nodata: rs141141725
    nodata: Affx-50099670
    nodata: Affx-50099671
    nodata: rs150428
    nodata: Affx-50099676
    nodata: rs1061581
    nodata: Affx-50112420
    nodata: rs148349653
    nodata: rs17201511
    nodata: rs139740532
    nodata: rs77335078
    nodata: Affx-50119664
    nodata: rs2311899
    nodata: rs2109936
    nodata: rs801050
    nodata: Affx-50120102
    nodata: Affx-50120105
    nodata: rs34472859
    nodata: Affx-50120106
    nodata: rs1637142
    nodata: rs186407955
    nodata: rs75647144
    nodata: Affx-50116060
    nodata: Affx-50116077
    nodata: Affx-50116081
    nodata: Affx-50124647
    nodata: Affx-50122316
    nodata: Affx-50122318
    nodata: rs11545441
    nodata: Affx-50122419
    nodata: Affx-50127044
    nodata: Affx-50127888
    nodata: rs61843996
    nodata: rs10508545
    nodata: Affx-50001783
    nodata: rs71497203
    nodata: rs145676940
    nodata: rs71497209
    nodata: rs3858193
    nodata: Affx-50001788
    nodata: rs35701218
    nodata: Affx-50001790
    nodata: Affx-50001799
    nodata: rs45470101
    nodata: Affx-50001803
    nodata: rs61844042
    nodata: Affx-50002744
    nodata: rs2942947
    nodata: Affx-50002746
    nodata: rs146410326
    nodata: Affx-50002771
    nodata: rs147158483
    nodata: rs150050324
    nodata: Affx-50011536
    nodata: rs79351613
    nodata: rs71473119
    nodata: Affx-50036749
    nodata: Affx-50036772
    nodata: Affx-50036779
    nodata: Affx-50037259
    nodata: rs145128593
    nodata: Affx-50037264
    nodata: Affx-50037270
    nodata: rs145771210
    nodata: Affx-50038810
    nodata: Affx-50040402
    nodata: Affx-50040413
    nodata: rs76539946
    nodata: Affx-50040139
    nodata: Affx-50047810
    nodata: Affx-50047802
    nodata: Affx-50047800
    nodata: Affx-50051857
    nodata: Affx-50064549
    nodata: Affx-50064552
    nodata: snp_24_3221179
    nodata: snp_24_3709585
    nodata: snp_24_3882879
    nodata: L381
    nodata: L1231
    nodata: rs186809987
    nodata: rs182663965
    nodata: 92R7_1
    nodata: snp_24_6424628
    nodata: Z3177
    nodata: Z3223
    nodata: 92R7_2
    nodata: rs201194017
    nodata: P68
    nodata: rs7892898
    nodata: P9.2
    nodata: V58
    nodata: rs199952688
    nodata: L353_1
    nodata: P16_1
    nodata: P16_2
    nodata: L353_2
    nodata: P112
    nodata: rs76320386
    nodata: Z3467
    nodata: Z3468
    nodata: snp_24_22272938
    nodata: Z3474
    nodata: P19_1
    nodata: P34_1
    nodata: P19_2
    nodata: P34_2
    nodata: rs55880099
    nodata: P8_2
    nodata: P19_3
    nodata: P34_3
    nodata: P8_3
    nodata: P19_4
    nodata: P18_1
    nodata: P7_1
    nodata: rs150173
    nodata: P12_1
    nodata: P65_1
    nodata: P66_1
    nodata: P67_1
    nodata: P35_1
    nodata: P33_1
    nodata: P21_1
    nodata: P18_2
    nodata: P7_2
    nodata: P12_2
    nodata: P65_2
    nodata: P66_2
    nodata: P67_2
    nodata: P35_2
    nodata: P33_2
    nodata: P21_2
    nodata: P21_3
    nodata: P33_3
    nodata: P35_3
    nodata: P67_3
    nodata: P66_3
    nodata: P65_3
    nodata: P7_3
    nodata: P18_3
    outpop: Mbuti.DG
    population: 0 Mbuti.DG 5
    population: 1 Mota 1
    population: 2 Onge.DG 2
    population: 3 Dai.DG 5
    population: 4 Ust_Ishim 1
    population: 5 GoyetQ116-1 1
    population: 6 Kostenki14 1
    population: 7 Sunghir 4
    population: 8 Gravettian 3
    population: 9 Brazil_Sumidouro_10100BP 4
    population: 10 Baikal_EN 13
    population: 11 Loschbour 1
    population: 12 Taforalt 4
    population: 13 Barcin_N 10
    population: 14 Iran_N 7
    population: 15 EHG 4
    population: 16 CHG 2
    before setwt numsnps: 1233013 outpop: Mbuti.DG
    setwt numsnps: 358121
    number of blocks for moving block jackknife: 710
    snps: 358121 indivs: 68
    lambdascale: 0.450
    fst:
    M M O D U G K S G B B L T B I E C
    M 0 504 234 193 146 530 528 252 409 331 222 252 271 195 206 273 371
    M 504 0 675 628 577 1000 1000 694 868 780 659 690 714 629 642 715 825
    O 234 675 0 152 134 628 628 277 477 350 197 273 353 215 226 299 432
    D 193 628 152 0 86 573 575 225 424 262 86 220 303 161 172 241 378
    U 146 577 134 86 0 512 514 164 360 262 125 164 249 108 127 200 329
    G 530 1000 628 573 512 0 1000 626 837 745 607 616 691 566 594 660 798
    K 528 1000 628 575 514 1000 0 613 836 746 609 623 692 566 590 661 797
    S 252 694 277 225 164 626 613 0 439 390 256 248 362 207 239 296 440
    G 409 868 477 424 360 837 836 439 0 591 456 446 545 403 438 498 642
    B 331 780 350 262 262 745 746 390 591 0 286 378 459 321 336 384 542
    B 222 659 197 86 125 607 609 256 456 286 0 249 334 193 207 267 410
    L 252 690 273 220 164 616 623 248 446 378 249 0 351 177 227 251 426
    T 271 714 353 303 249 691 692 362 545 459 334 351 0 279 314 383 505
    B 195 629 215 161 108 566 566 207 403 321 193 177 279 0 144 220 347
    I 206 642 226 172 127 594 590 239 438 336 207 227 314 144 0 248 355
    E 273 715 299 241 200 660 661 296 498 384 267 251 383 220 248 0 452
    C 371 825 432 378 329 798 797 440 642 542 410 426 505 347 355 452 0

    f2:
    M M O D U G K S G B B L T B I E C
    M 0 567 292 246 189 654 645 315 505 413 280 317 335 248 260 345 460
    M 567 0 714 670 612 1077 1071 740 927 831 703 736 753 672 683 768 884
    O 292 714 0 135 126 593 592 262 452 317 177 257 357 206 215 284 412
    D 246 670 135 0 81 542 545 214 403 228 73 208 308 155 164 227 361
    U 189 612 126 81 0 481 481 156 339 247 118 156 253 106 123 193 318
    G 654 1077 593 542 481 0 887 556 739 696 573 541 690 527 567 603 751
    K 645 1071 592 545 481 887 0 524 726 698 575 554 691 524 556 605 744
    S 315 740 262 214 156 556 524 0 373 365 242 219 362 191 226 270 410
    G 505 927 452 403 339 739 726 373 0 551 431 386 539 370 414 451 597
    B 413 831 317 228 247 696 698 365 551 0 243 350 463 305 318 343 508
    B 280 703 177 73 118 573 575 242 431 243 0 233 338 186 197 248 390
    L 317 736 257 208 156 541 554 219 386 350 233 0 349 159 212 216 390
    T 335 753 357 308 253 690 691 362 539 463 338 349 0 275 317 384 507
    B 248 672 206 155 106 527 524 191 370 305 186 159 275 0 133 202 315
    I 260 683 215 164 123 567 556 226 414 318 197 212 317 133 0 231 318
    E 345 768 284 227 193 603 605 270 451 343 248 216 384 202 231 0 414
    C 460 884 412 361 318 751 744 410 597 508 390 390 507 315 318 414 0

    ff3:
    M M O D U G K S G B B L T B I E C
    M 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    M 0 567 72 72 72 72 70 71 73 74 72 74 75 72 72 72 72
    O 0 72 292 201 177 176 173 173 173 194 198 176 135 167 169 177 170
    D 0 72 201 246 177 179 173 174 174 215 227 178 137 169 171 182 173
    U 0 72 177 177 189 181 176 174 177 177 176 175 136 166 163 171 165
    G 0 72 176 179 181 654 206 207 210 185 180 215 150 188 174 198 181
    K 0 70 173 173 176 206 645 218 212 180 175 204 145 184 175 193 181
    S 0 71 173 174 174 207 218 315 224 181 177 207 144 186 175 195 183
    G 0 73 173 174 177 210 212 224 505 183 177 218 151 192 176 200 184
    B 0 74 194 215 177 185 180 181 183 413 225 190 142 178 178 207 182
    B 0 72 198 227 176 180 175 177 177 225 280 182 139 171 172 189 175
    L 0 74 176 178 175 215 204 207 218 190 182 317 152 203 183 223 194
    T 0 75 135 137 136 150 145 144 151 142 139 152 335 154 139 148 144
    B 0 72 167 169 166 188 184 186 192 178 171 203 154 248 188 196 197
    I 0 72 169 171 163 174 175 175 176 178 172 183 139 188 260 187 201
    E 0 72 177 182 171 198 193 195 200 207 189 223 148 196 187 345 196
    C 0 72 170 173 165 181 181 183 184 182 175 194 144 197 201 196 460

    ff3sig*10:
    M M O D U G K S G B B L T B I E C
    M 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    M 0 24 18 17 19 21 20 17 19 19 17 19 19 17 17 18 19
    O 0 18 32 23 25 26 26 22 23 23 22 25 21 21 21 23 23
    D 0 17 23 24 25 24 25 21 22 23 22 24 20 20 20 22 22
    U 0 19 25 25 52 29 29 24 26 25 24 28 22 22 22 24 25
    G 0 21 26 24 29 27 33 26 27 26 24 29 23 23 22 26 26
    K 0 20 26 25 29 33 28 27 28 27 24 28 24 23 23 25 26
    S 0 17 22 21 24 26 27 26 24 23 21 25 20 21 20 22 23
    G 0 19 23 22 26 27 28 24 25 23 22 26 22 22 21 23 24
    B 0 19 23 23 25 26 27 23 23 26 23 25 21 20 21 23 23
    B 0 17 22 22 24 24 24 21 22 23 22 23 19 20 20 22 22
    L 0 19 25 24 28 29 28 25 26 25 23 45 22 22 22 25 26
    T 0 19 21 20 22 23 24 20 22 21 19 22 24 20 19 21 21
    B 0 17 21 20 22 23 23 21 22 20 20 22 20 21 19 21 21
    I 0 17 21 20 22 22 23 20 21 21 20 22 19 19 20 21 22
    E 0 18 23 22 24 26 25 22 23 23 22 25 21 21 21 25 23
    C 0 19 23 22 25 26 26 23 24 23 22 26 21 21 22 23 26

    starting analysis
    number of initial random trials: 1850
    initmixwts: :::Kostenki14 W3 Basal1 W2 0.116 0.884
    initmixwts: :Gravettian pre_GRV W6 WHG0 0.478 0.522
    initmixwts: : pre_ANE E2 W4 0.240 0.760
    initmixwts: :Brazil_Sumidouro_10100BP America ANE E3 0.351 0.649
    initmixwts: :Baikal_EN pre_BKN ANE E3 0.110 0.890
    initmixwts: :: WHG1 ANE WHG0 0.399 0.601
    initmixwts: ::Loschbour WHG3 W2a WHG2 0.111 0.889
    initmixwts: pre_IBM1 W2a WHG2 0.125 0.875
    initmixwts: :Taforalt pre_IBM2 ANA pre_IBM1 0.743 0.257
    initmixwts: :: NE1 Basal1 E1 0.729 0.271
    initmixwts: ANF1 NE4 WHG5 0.025 0.975
    initmixwts: :Barcin_N ANF2 ANA ANF1 0.216 0.784
    initmixwts: pre_EHG1 ANE WHG4 0.825 0.175
    initmixwts: :EHG pre_EHG2 NE4 pre_EHG1 0.102 0.898
    initmixwts: INF1 ANE NE3 0.308 0.692
    initmixwts: :Iran_N INF2 INF1 WHG5 0.830 0.170
    initmixwts: pre_CHG1 ANE NE3 0.430 0.570
    initmixwts: :CHG pre_CHG2 pre_CHG1 WHG5 0.776 0.224


    initial score: 32378.903
    init vg:
    0.116 0.478 0.240 0.351 0.110 ...
    0.399 0.111 0.125 0.743 0.729 ...
    0.025 0.216 0.825 0.102 0.308 ...
    0.830 0.430 0.776
    gslsetup called
    gslans: 3897 17309.563879
    0.150 0.408 0.298 0.366 0.123 ...
    0.187 0.039 0.053 0.700 0.691 ...
    0.284 0.097 0.556 0.170 0.129 ...
    0.774 0.163 0.616
    final vg:
    0.150 0.408 0.298 0.366 0.123 ...
    0.187 0.039 0.053 0.700 0.691 ...
    0.284 0.097 0.556 0.170 0.129 ...
    0.774 0.163 0.616

    final score: 17309.564 dof: 78 nominal p-val 0.000000
    ff3fit:
    M M O D U G K S G B B L T B I E C
    M 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    M 0 566 72 72 72 72 72 72 72 72 72 72 72 72 72 72 72
    O 0 72 292 201 177 177 173 173 173 194 198 175 137 167 169 176 171
    D 0 72 201 246 177 177 173 173 173 216 226 176 137 168 171 181 172
    U 0 72 177 177 189 180 175 175 175 177 177 176 137 166 164 173 167
    G 0 72 177 177 180 654 208 208 208 185 180 215 150 186 173 198 181
    K 0 72 173 173 175 208 644 219 212 182 176 205 146 185 174 194 181
    S 0 72 173 173 175 208 219 315 224 182 176 205 146 185 174 194 181
    G 0 72 173 173 175 208 212 224 505 182 176 217 149 193 177 198 186
    B 0 72 194 216 177 185 182 182 182 412 225 191 142 178 177 206 182
    B 0 72 198 226 177 180 176 176 176 225 280 181 139 171 173 190 175
    L 0 72 175 176 176 215 205 205 217 191 181 317 154 204 183 223 194
    T 0 72 137 137 137 150 146 146 149 142 139 154 335 153 139 148 142
    B 0 72 167 168 166 186 185 185 193 178 171 204 153 248 187 195 197
    I 0 72 169 171 164 173 174 174 177 177 173 183 139 187 260 188 202
    E 0 72 176 181 173 198 194 194 198 206 190 223 148 195 188 345 195
    C 0 72 171 172 167 181 181 181 186 182 175 194 142 197 202 195 460

    ff3diff:
    M M O D U G K S G B B L T B I E C
    M 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    M 0 -1 0 1 0 0 2 2 0 -2 0 -2 -3 1 0 1 1
    O 0 0 0 -1 0 1 1 1 1 0 1 -2 1 0 0 0 0
    D 0 1 -1 0 1 -2 0 -1 -1 1 -1 -1 0 -1 -1 -1 0
    U 0 0 0 1 0 -1 -1 1 -2 0 2 1 1 0 1 2 1
    G 0 0 1 -2 -1 0 2 1 -2 0 0 0 0 -2 0 0 -1
    K 0 2 1 0 -1 2 -1 1 0 2 1 1 1 1 -1 1 0
    S 0 2 1 -1 1 1 1 0 0 1 0 -1 1 -1 -1 -1 -1
    G 0 0 1 -1 -2 -2 0 0 0 -2 -1 -1 -2 2 1 -1 2
    B 0 -2 0 1 0 0 2 1 -2 0 0 1 -1 0 0 -1 -1
    B 0 0 1 -1 2 0 1 0 -1 0 0 -1 0 0 1 1 0
    L 0 -2 -2 -1 1 0 1 -1 -1 1 -1 0 2 0 0 0 0
    T 0 -3 1 0 1 0 1 1 -2 -1 0 2 0 -1 0 0 -2
    B 0 1 0 -1 0 -2 1 -1 2 0 0 0 -1 0 -1 -1 0
    I 0 0 0 -1 1 0 -1 -1 1 0 1 0 0 -1 0 1 0
    E 0 1 0 -1 2 0 1 -1 -1 -1 1 0 0 -1 1 0 -1
    C 0 1 0 0 1 -1 0 -1 2 -1 0 0 -2 0 0 -1 0

    mix: :::Kostenki14 W3 Basal1 W2 0.150 0.850
    mix: :Gravettian pre_GRV W6 WHG0 0.408 0.592
    mix: : pre_ANE E2 W4 0.298 0.702
    mix: :Brazil_Sumidouro_10100BP America ANE E3 0.366 0.634
    mix: :Baikal_EN pre_BKN ANE E3 0.123 0.877
    mix: :: WHG1 ANE WHG0 0.187 0.813
    mix: ::Loschbour WHG3 W2a WHG2 0.039 0.961
    mix: pre_IBM1 W2a WHG2 0.053 0.947
    mix: :Taforalt pre_IBM2 ANA pre_IBM1 0.700 0.300
    mix: :: NE1 Basal1 E1 0.691 0.309
    mix: ANF1 NE4 WHG5 0.284 0.716
    mix: :Barcin_N ANF2 ANA ANF1 0.097 0.903
    mix: pre_EHG1 ANE WHG4 0.556 0.444
    mix: :EHG pre_EHG2 NE4 pre_EHG1 0.170 0.830
    mix: INF1 ANE NE3 0.129 0.871
    mix: :Iran_N INF2 INF1 WHG5 0.774 0.226
    mix: pre_CHG1 ANE NE3 0.163 0.837
    mix: :CHG pre_CHG2 pre_CHG1 WHG5 0.616 0.384

    Mbuti.DG A Mota B ANA C Basal1 Eurasia E1 W1 Onge.DG E2 Dai.DG E3 Ust_Ishim W2 W2a GoyetQ116-1 W4 W5 WHG0 Kostenki14 W6 Sunghir Gravettian ANE Brazil_Sumidouro_10100BP Baikal_EN WHG2 WHG5 WHG4 Loschbour Taforalt NE2 NE3 NE4 Barcin_N EHG Iran_N CHG
    0.0362 0.0362 0.4939 0.0478 0.1197 0.0296 0.0302 0.0276 0.0235 0.0025 0.0909 0.0302 0.0149 0.0118 0.0091 0.0385 0.2156 0.2196 0.0071 0.0116 0.0257 0.4253 0.0304 0.0659 0.2816 0.1138 0.1849 0.0475 0.0099 0.0435 0.0274 0.0572 0.1463 0.0470 0.0187 0.0277 0.0393 0.1173 0.0557 0.2564


    fst: fitted estim diff std. err Z score
    Mbu Mot f2: 0.566302 0.567094 0.000792 0.002366 0.335
    Mbu Ong f2: 0.291829 0.291975 0.000146 0.003204 0.046
    Mbu Dai f2: 0.246024 0.246048 0.000024 0.002383 0.010
    Mbu Ust f2: 0.188930 0.188944 0.000015 0.005230 0.003
    Mbu Goy f2: 0.653520 0.653872 0.000352 0.002726 0.129
    Mbu Kos f2: 0.643957 0.644639 0.000682 0.002774 0.246
    Mbu Sun f2: 0.314998 0.315104 0.000106 0.002561 0.041
    Mbu Gra f2: 0.504711 0.505163 0.000452 0.002476 0.182
    Mbu Bra f2: 0.412332 0.412629 0.000297 0.002645 0.112
    Mbu Bai f2: 0.280071 0.280148 0.000076 0.002228 0.034
    Mbu Los f2: 0.317216 0.317308 0.000092 0.004469 0.021
    Mbu Taf f2: 0.335137 0.335371 0.000235 0.002429 0.097
    Mbu Bar f2: 0.248028 0.248091 0.000063 0.002061 0.031
    Mbu Ira f2: 0.260232 0.260321 0.000089 0.002026 0.044
    Mbu EHG f2: 0.344828 0.345016 0.000188 0.002460 0.076
    Mbu CHG f2: 0.459993 0.460404 0.000411 0.002569 0.160
    Mot Ong f2: 0.713334 0.714390 0.001056 0.003810 0.277
    Mot Dai f2: 0.667529 0.670036 0.002507 0.002854 0.879
    Mot Ust f2: 0.610435 0.611658 0.001223 0.005385 0.227
    Mot Goy f2: 1.075025 1.076656 0.001631 0.003545 0.460
    Mot Kos f2: 1.065462 1.071092 0.005630 0.003521 1.599
    Mot Sun f2: 0.736503 0.740458 0.003955 0.003221 1.228
    Mot Gra f2: 0.926216 0.926500 0.000284 0.003247 0.087
    Mot Bra f2: 0.833837 0.831375 -0.002462 0.003201 -0.769
    Mot Bai f2: 0.701576 0.703060 0.001484 0.002709 0.548
    Mot Los f2: 0.738721 0.736060 -0.002661 0.004842 -0.550
    Mot Taf f2: 0.756642 0.752630 -0.004011 0.002971 -1.350
    Mot Bar f2: 0.669533 0.671628 0.002094 0.002570 0.815
    Mot Ira f2: 0.681737 0.682667 0.000930 0.002588 0.360
    Mot EHG f2: 0.766333 0.768474 0.002141 0.002880 0.743
    Mot CHG f2: 0.881498 0.884430 0.002932 0.003343 0.877
    Ong Dai f2: 0.136091 0.135060 -0.001031 0.002566 -0.402
    Ong Ust f2: 0.126067 0.126135 0.000067 0.005332 0.013
    Ong Goy f2: 0.590657 0.592911 0.002254 0.003946 0.571
    Ong Kos f2: 0.589357 0.591599 0.002243 0.003959 0.566
    Ong Sun f2: 0.260398 0.261756 0.001358 0.003142 0.432
    Ong Gra f2: 0.450111 0.451820 0.001709 0.003309 0.516
    Ong Bra f2: 0.316622 0.317346 0.000724 0.003355 0.216
    Ong Bai f2: 0.174907 0.176952 0.002045 0.002630 0.778
    Ong Los f2: 0.259327 0.256508 -0.002819 0.004955 -0.569
    Ong Taf f2: 0.353763 0.356747 0.002985 0.003391 0.880
    Ong Bar f2: 0.205636 0.205855 0.000219 0.002765 0.079
    Ong Ira f2: 0.213593 0.214549 0.000957 0.002754 0.347
    Ong EHG f2: 0.283999 0.283628 -0.000371 0.003098 -0.120
    Ong CHG f2: 0.410415 0.411894 0.001479 0.003435 0.431
    Dai Ust f2: 0.080262 0.081312 0.001051 0.004830 0.218
    Dai Goy f2: 0.544852 0.542162 -0.002690 0.003104 -0.866
    Dai Kos f2: 0.543551 0.544958 0.001406 0.002977 0.472
    Dai Sun f2: 0.214592 0.213616 -0.000976 0.002136 -0.457
    Dai Gra f2: 0.404305 0.402796 -0.001509 0.002387 -0.632
    Dai Bra f2: 0.225888 0.228266 0.002378 0.001996 1.192
    Dai Bai f2: 0.073839 0.072707 -0.001133 0.000990 -1.144
    Dai Los f2: 0.210280 0.208146 -0.002134 0.004324 -0.493
    Dai Taf f2: 0.306998 0.307501 0.000504 0.002441 0.206
    Dai Bar f2: 0.157650 0.155292 -0.002358 0.001604 -1.471
    Dai Ira f2: 0.165230 0.164241 -0.000990 0.001580 -0.626
    Dai EHG f2: 0.228683 0.227458 -0.001225 0.002042 -0.600
    Dai CHG f2: 0.361500 0.361372 -0.000128 0.002413 -0.053
    Ust Goy f2: 0.482745 0.480657 -0.002088 0.006092 -0.343
    Ust Kos f2: 0.482195 0.481134 -0.001062 0.005639 -0.188
    Ust Sun f2: 0.153236 0.155689 0.002453 0.005083 0.483
    Ust Gra f2: 0.342949 0.339395 -0.003554 0.005147 -0.691
    Ust Bra f2: 0.247598 0.247153 -0.000445 0.005140 -0.087
    Ust Bai f2: 0.114654 0.117993 0.003339 0.004703 0.710
    Ust Los f2: 0.154888 0.156417 0.001529 0.006303 0.243
    Ust Taf f2: 0.250385 0.252903 0.002518 0.005277 0.477
    Ust Bar f2: 0.105568 0.105588 0.000020 0.004708 0.004
    Ust Ira f2: 0.121297 0.123081 0.001784 0.004704 0.379
    Ust EHG f2: 0.188093 0.192627 0.004534 0.005012 0.905
    Ust CHG f2: 0.315589 0.318419 0.002830 0.005224 0.542
    Goy Kos f2: 0.881359 0.886971 0.005612 0.005349 1.049
    Goy Sun f2: 0.552400 0.555707 0.003307 0.003622 0.913
    Goy Gra f2: 0.742113 0.738900 -0.003213 0.004280 -0.751
    Goy Bra f2: 0.695370 0.696486 0.001116 0.003948 0.283
    Goy Bai f2: 0.573605 0.573099 -0.000507 0.003154 -0.161
    Goy Los f2: 0.540262 0.540862 0.000600 0.005538 0.108
    Goy Taf f2: 0.689300 0.689710 0.000410 0.003470 0.118
    Goy Bar f2: 0.530215 0.526779 -0.003436 0.002798 -1.228
    Goy Ira f2: 0.567347 0.566910 -0.000437 0.002892 -0.151
    Goy EHG f2: 0.602269 0.603311 0.001042 0.003406 0.306
    Goy CHG f2: 0.751813 0.751297 -0.000515 0.003928 -0.131
    Kos Sun f2: 0.521549 0.523923 0.002374 0.003315 0.716
    Kos Gra f2: 0.725064 0.725801 0.000737 0.004031 0.183
    Kos Bra f2: 0.692461 0.697804 0.005343 0.003913 1.366
    Kos Bai f2: 0.571766 0.575150 0.003385 0.003068 1.103
    Kos Los f2: 0.550610 0.553828 0.003218 0.005311 0.606
    Kos Taf f2: 0.687610 0.690565 0.002955 0.003291 0.898
    Kos Bar f2: 0.521229 0.524001 0.002772 0.002586 1.072
    Kos Ira f2: 0.556216 0.555627 -0.000590 0.002586 -0.228
    Kos EHG f2: 0.601208 0.604627 0.003419 0.003369 1.015
    Kos CHG f2: 0.741646 0.743545 0.001899 0.003584 0.530
    Sun Gra f2: 0.371347 0.372553 0.001206 0.002624 0.459
    Sun Bra f2: 0.363502 0.365464 0.001962 0.002732 0.718
    Sun Bai f2: 0.242806 0.242172 -0.000635 0.002101 -0.302
    Sun Los f2: 0.221651 0.219327 -0.002324 0.004328 -0.537
    Sun Taf f2: 0.358651 0.361681 0.003030 0.002678 1.132
    Sun Bar f2: 0.192270 0.191172 -0.001098 0.001846 -0.595
    Sun Ira f2: 0.227257 0.226062 -0.001196 0.001993 -0.600
    Sun EHG f2: 0.272249 0.270264 -0.001985 0.002150 -0.923
    Sun CHG f2: 0.412687 0.410497 -0.002189 0.002667 -0.821
    Gra Bra f2: 0.553215 0.550953 -0.002262 0.003225 -0.702
    Gra Bai f2: 0.432520 0.430855 -0.001665 0.002353 -0.707
    Gra Los f2: 0.387557 0.386485 -0.001073 0.004787 -0.224
    Gra Taf f2: 0.541317 0.538572 -0.002745 0.002883 -0.952
    Gra Bar f2: 0.365964 0.369881 0.003917 0.001995 1.964
    Gra Ira f2: 0.411371 0.413860 0.002489 0.002154 1.155
    Gra EHG f2: 0.453184 0.450993 -0.002191 0.002598 -0.843
    Gra CHG f2: 0.592878 0.597444 0.004566 0.003001 1.521
    Bra Bai f2: 0.242509 0.243453 0.000944 0.001902 0.497
    Bra Los f2: 0.346785 0.350118 0.003333 0.004663 0.715
    Bra Taf f2: 0.464414 0.463398 -0.001016 0.002897 -0.351
    Bra Bar f2: 0.305037 0.305215 0.000178 0.002019 0.088
    Bra Ira f2: 0.317923 0.317694 -0.000229 0.002100 -0.109
    Bra EHG f2: 0.344468 0.343171 -0.001297 0.002429 -0.534
    Bra CHG f2: 0.508852 0.508352 -0.000501 0.002949 -0.170
    Bai Los f2: 0.234335 0.233176 -0.001159 0.004260 -0.272
    Bai Taf f2: 0.338064 0.338434 0.000370 0.002353 0.157
    Bai Bar f2: 0.185353 0.185596 0.000242 0.001460 0.166
    Bai Ira f2: 0.194713 0.196653 0.001940 0.001515 1.281
    Bai EHG f2: 0.245792 0.248139 0.002348 0.001844 1.273
    Bai CHG f2: 0.389192 0.389846 0.000654 0.002305 0.284
    Los Taf f2: 0.344595 0.348598 0.004003 0.004663 0.859
    Los Bar f2: 0.157905 0.158616 0.000711 0.003904 0.182
    Los Ira f2: 0.212403 0.212387 -0.000017 0.004127 -0.004
    Los EHG f2: 0.215365 0.216043 0.000678 0.004261 0.159
    Los CHG f2: 0.388932 0.390031 0.001098 0.004696 0.234
    Taf Bar f2: 0.276602 0.274644 -0.001958 0.002107 -0.930
    Taf Ira f2: 0.317664 0.317354 -0.000311 0.002274 -0.137
    Taf EHG f2: 0.383851 0.384280 0.000430 0.002607 0.165
    Taf CHG f2: 0.510481 0.506779 -0.003703 0.002868 -1.291
    Bar Ira f2: 0.134112 0.133011 -0.001101 0.001042 -1.057
    Bar EHG f2: 0.203142 0.201788 -0.001353 0.001405 -0.964
    Bar CHG f2: 0.314285 0.315432 0.001147 0.001871 0.613
    Ira EHG f2: 0.228650 0.230951 0.002301 0.001689 1.362
    Ira CHG f2: 0.316972 0.318162 0.001190 0.001826 0.652
    EHG CHG f2: 0.414507 0.413899 -0.000608 0.002314 -0.263


    outliers:
    Fit Obs Diff Std. error Z
    Mot Sun Taf Bar 0.039636 0.044753 0.005117 0.001660 3.082
    Mot Los Dai Bra 0.014901 0.009684 -0.005218 0.001728 -3.020
    Mot Los Sun Taf -0.051403 -0.058550 -0.007147 0.002206 -3.240
    Mot Taf Los Bar -0.000597 0.004761 0.005358 0.001666 3.217
    Mot Taf Los CHG -0.011555 -0.004905 0.006649 0.002103 3.162
    Mot CHG Sun Gra 0.004761 -0.000452 -0.005213 0.001733 -3.008
    Dai Gra Bra Bar 0.059507 0.054049 -0.005458 0.001721 -3.171
    Dai Gra Taf Bar 0.013002 0.008241 -0.004762 0.001437 -3.313
    Ust Gra Bar EHG -0.002347 0.002964 0.005311 0.001507 3.524
    Ust Bar Sun Gra 0.008010 0.002498 -0.005511 0.001610 -3.423
    Ust Bar Gra Taf -0.001601 0.004373 0.005974 0.001938 3.082
    Ust Bar Gra EHG 0.003983 0.010662 0.006679 0.001867 3.578
    Ust CHG Sun Gra 0.004761 -0.001620 -0.006381 0.001981 -3.221
    Ust CHG Gra Taf -0.005084 0.002086 0.007170 0.002185 3.282
    Ust CHG Gra EHG 0.011757 0.018389 0.006631 0.002172 3.053
    Goy Gra Bar EHG -0.007583 -0.002290 0.005293 0.001585 3.339
    Goy Bar Sun Gra 0.008010 0.002242 -0.005768 0.001906 -3.026
    Goy CHG Sun Gra 0.004761 -0.001877 -0.006638 0.002094 -3.170
    Sun Gra Bra Bar 0.008010 0.003390 -0.004619 0.001328 -3.479
    Sun Gra Bra Ira 0.002800 -0.001154 -0.003954 0.001313 -3.011
    Sun Gra Bra CHG 0.004761 -0.000729 -0.005490 0.001630 -3.368
    Sun Gra Taf Bar 0.004486 -0.000909 -0.005395 0.001266 -4.263
    Sun Gra Taf Ira -0.000724 -0.005454 -0.004730 0.001326 -3.567
    Sun Gra Taf CHG 0.001237 -0.005028 -0.006265 0.001653 -3.791
    Sun CHG Taf Bar 0.014908 0.010419 -0.004489 0.001416 -3.170
    Gra Bai Bar EHG 0.013391 0.009284 -0.004107 0.001275 -3.221
    Gra Los Taf Bar 0.005669 0.010646 0.004977 0.001588 3.133
    Gra Taf Bai CHG -0.006030 -0.000878 0.005152 0.001706 3.021
    Gra Taf Los Bar 0.023200 0.028675 0.005475 0.001811 3.024
    Gra Taf Los CHG 0.019717 0.026389 0.006672 0.002045 3.262
    Gra Taf Bar EHG -0.010015 -0.014263 -0.004248 0.001371 -3.099
    Gra Taf EHG CHG 0.006532 0.011976 0.005445 0.001681 3.239
    Gra Ira Bai Bar -0.002977 0.001334 0.004312 0.001293 3.334
    Gra Ira Bar EHG -0.003659 -0.008414 -0.004755 0.001282 -3.710
    Gra EHG Bai Bar -0.011953 -0.007311 0.004642 0.001333 3.481
    Gra EHG Taf Bar 0.002678 0.006901 0.004222 0.001367 3.088


    worst f-stat: Sun Gra Taf Bar 0.004486 -0.000909 -0.005395 0.001266 -4.263

    vertex Root 0
    vertex Mbuti.DG 0
    vertex A 0
    vertex Mota 0
    vertex B 0
    vertex ANA 0
    vertex C 0
    vertex Basal1 0
    vertex Eurasia 0
    vertex pre_IBM2 0
    vertex ANF2 0
    vertex Taforalt 0
    vertex Barcin_N 0
    vertex E1 0
    vertex W1 0
    vertex W3 0
    vertex NE1 0
    vertex Onge.DG 0
    vertex Ust_Ishim 0
    vertex E2 0
    vertex W2 0
    vertex W4 0
    vertex NE2 0
    vertex Dai.DG 0
    vertex E3 0
    vertex W2a 0
    vertex W5 0
    vertex WHG0 0
    vertex pre_ANE 0
    vertex NE3 0
    vertex NE4 0
    vertex GoyetQ116-1 0
    vertex Kostenki14 0
    vertex W6 0
    vertex pre_GRV 0
    vertex ANE 0
    vertex America 0
    vertex pre_BKN 0
    vertex WHG1 0
    vertex WHG3 0
    vertex pre_IBM1 0
    vertex ANF1 0
    vertex pre_EHG2 0
    vertex INF1 0
    vertex pre_CHG1 0
    vertex Sunghir 0
    vertex Gravettian 0
    vertex Brazil_Sumidouro_10100BP 0
    vertex Baikal_EN 0
    vertex EHG 0
    vertex WHG2 0
    vertex WHG4 0
    vertex pre_EHG1 0
    vertex INF2 0
    vertex pre_CHG2 0
    vertex Loschbour 0
    vertex Iran_N 0
    vertex CHG 0
    vertex WHG5 0
    label Mbuti.DG Mbuti.DG 5
    label Mota Mota 1
    label Taforalt Taforalt 4
    label Barcin_N Barcin_N 10
    label Onge.DG Onge.DG 2
    label Ust_Ishim Ust_Ishim 1
    label Dai.DG Dai.DG 5
    label GoyetQ116-1 GoyetQ116-1 1
    label Kostenki14 Kostenki14 1
    label Sunghir Sunghir 4
    label Gravettian Gravettian 3
    label Brazil_Sumidouro_10100BP Brazil_Sumidouro_10100BP 4
    label Baikal_EN Baikal_EN 13
    label EHG EHG 4
    label Loschbour Loschbour 1
    label Iran_N Iran_N 7
    label CHG CHG 2
    ledge Mbuti.DG Root Mbuti.DG 0.036199
    redge A Root A 0.036199
    ledge Mota A Mota 0.493903
    redge B A B 0.047760
    ledge ANA B ANA 0.119658
    redge C B C 0.029599
    ledge Basal1 C Basal1 0.030186
    redge Eurasia C Eurasia 0.027588
    ledge E1 Eurasia E1 0.023535
    redge W1 Eurasia W1 0.002506
    ledge Taforalt pre_IBM2 Taforalt 0.146302
    ledge Barcin_N ANF2 Barcin_N 0.039275
    ledge Onge.DG E1 Onge.DG 0.090948
    redge E2 E1 E2 0.030237
    ledge Ust_Ishim W1 Ust_Ishim 0.009078
    redge W2 W1 W2 0.038474
    ledge W4 W3 W4 0.007051
    ledge NE2 NE1 NE2 0.047042
    ledge Dai.DG E2 Dai.DG 0.014906
    redge E3 E2 E3 0.011791
    ledge W2a W2 W2a 0.215642
    ledge W5 W4 W5 0.011647
    redge WHG0 W4 WHG0 0.025725
    ledge NE3 NE2 NE3 0.018702
    redge NE4 NE2 NE4 0.027726
    ledge GoyetQ116-1 W2a GoyetQ116-1 0.219552
    ledge Kostenki14 W5 Kostenki14 0.425254
    redge W6 W5 W6 0.030377
    ledge ANE pre_ANE ANE 0.113802
    ledge Sunghir W6 Sunghir 0.065918
    ledge Gravettian pre_GRV Gravettian 0.281646
    ledge Brazil_Sumidouro_10100BP America Brazil_Sumidouro_10100BP 0.184929
    ledge Baikal_EN pre_BKN Baikal_EN 0.047464
    ledge WHG2 WHG1 WHG2 0.009920
    ledge WHG4 WHG3 WHG4 0.027364
    ledge EHG pre_EHG2 EHG 0.117315
    ledge WHG5 WHG2 WHG5 0.043499
    ledge Loschbour WHG4 Loschbour 0.057190
    ledge Iran_N INF2 Iran_N 0.055738
    ledge CHG pre_CHG2 CHG 0.256383
    admix pre_IBM2 ANA pre_IBM1 0.699596 0.300404
    admix ANF2 ANA ANF1 0.096937 0.903063
    admix W3 Basal1 W2 0.149738 0.850262
    admix NE1 Basal1 E1 0.690571 0.309429
    admix pre_ANE E2 W4 0.297854 0.702146
    admix pre_GRV W6 WHG0 0.407501 0.592499
    admix America ANE E3 0.366094 0.633906
    admix pre_BKN ANE E3 0.122742 0.877258
    admix WHG1 ANE WHG0 0.187285 0.812715
    admix WHG3 W2a WHG2 0.039061 0.960939
    admix pre_IBM1 W2a WHG2 0.053070 0.946930
    admix ANF1 NE4 WHG5 0.284015 0.715985
    admix pre_EHG2 NE4 pre_EHG1 0.169667 0.830333
    admix INF1 ANE NE3 0.128734 0.871266
    admix pre_CHG1 ANE NE3 0.163470 0.836530
    admix pre_EHG1 ANE WHG4 0.555926 0.444074
    admix INF2 INF1 WHG5 0.773973 0.226027
    admix pre_CHG2 pre_CHG1 WHG5 0.615660 0.384340
    ## end of qpGraph time: 4299.447 seconds memory: 1396.388 Mbytes


    The worst outlier is 4.263, not markedly different from the ~3.7 with the phylogeny being a Goyet-related population feeding into the root of all East Eurasians (including Ust_Ishim)
    In both cases though, the only real cause of either model not reaching below the generally accepted passing level of worst Z-score <3, is Gravettians. Taforalt not close enough, Near Eastern populations too close.
    So either way works, the question is what can be done to rule 1 or the other out?
    Until that can be answered I'm going to be biased and say I like the Goyet into East-Eurasian version better, it avoids the need for Basal Eurasian.
    Until we find a pure one, I'll be skeptical of Basal Eurasian.

    Also, here ANE is shown as deriving it's Eastern ancestry from a Dai-like population.
    It works equally well to bump that to a more basal Onge-like population, then have Dai receive 5% ANE.
    Last edited by Kale; 06-13-2019 at 03:01 AM.
    Collection of 14,000 d-stats: Hidden Content Part 2: Hidden Content Part 3: Hidden Content PM me for d-stats, qpadm, qpgraph, or f3-outgroup nmonte models.

  12. The Following User Says Thank You to Kale For This Useful Post:

     Ryukendo (06-13-2019)

Page 3 of 3 FirstFirst 123

Similar Threads

  1. Certainly Related
    By curiousII in forum Autosomal (auDNA)
    Replies: 0
    Last Post: 06-08-2019, 07:43 AM
  2. Replies: 7
    Last Post: 06-14-2018, 05:54 AM
  3. Replies: 4
    Last Post: 02-22-2018, 05:25 PM
  4. We're Related app
    By LauraHolland in forum AncestryDNA
    Replies: 5
    Last Post: 06-30-2017, 11:37 AM
  5. Replies: 6
    Last Post: 07-10-2016, 08:32 AM

Posting Permissions

  • You may not post new threads
  • You may not post replies
  • You may not post attachments
  • You may not edit your posts
  •