
Originally Posted by
afbarwaaqo
could there be additional iran_n ancestry hidden when natufian is used? natufian gets 15% iran_n on this model
That's quite unlikely considering Iran N and Natufian are very differentiated ancestries.

Originally Posted by
afbarwaaqo
can you post how natufian is modeled on qpadm?
As Gentica said, there are some models in Lazaridis et al. official paper. Since you are interested, here are some models I did by myself:
right = c('Mbuti.HO', 'Ami.HO', 'Mixe.HO', 'Russia_Kostenki14', 'Russia_Samara_HG', 'Italy_North_Villabruna_HG')
Code:
> results$weights
# A tibble: 2 x 5
target left weight se z
<chr> <chr> <dbl> <dbl> <dbl>
1 Levant_Natufian Turkey_Epipaleolithic 0.678 0.0445 15.2
2 Levant_Natufian Morocco_Iberomaurusian 0.322 0.0445 7.24
>
> results$popdrop
# A tibble: 3 x 13
pat wt dof chisq p f4rank Turkey_Epipaleolithic Morocco_Iberomaurusian feasible best dofdiff chisqdiff p_nested
<chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <lgl> <lgl> <dbl> <dbl> <dbl>
1 00 0 4 1.40 8.44e- 1 1 0.678 0.322 TRUE NA NA NA NA
2 01 1 5 53.5 2.66e-10 0 1 NA TRUE TRUE 0 -261. 1
3 10 1 5 315. 6.17e-66 0 NA 1 TRUE TRUE NA NA NA
>
Keep in mind Pinarbasi HG is WHG-admixed so it's not the best proxy indeed. Here's another model using a variety of Neolithic Anatolians (excepting Boncuklu N):
Code:
> results$weights
# A tibble: 2 x 5
target left weight se z
<chr> <chr> <dbl> <dbl> <dbl>
1 Levant_Natufian Turkey_N 0.709 0.0362 19.6
2 Levant_Natufian Morocco_Iberomaurusian 0.291 0.0362 8.03
>
> results$popdrop
# A tibble: 3 x 13
pat wt dof chisq p f4rank Turkey_N Morocco_Iberomaurusian feasible best dofdiff chisqdiff p_nested
<chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <lgl> <lgl> <dbl> <dbl> <dbl>
1 00 0 4 2.31 6.79e- 1 1 0.709 0.291 TRUE NA NA NA NA
2 01 1 5 122. 1.27e-24 0 1 NA TRUE TRUE 0 -300. 1
3 10 1 5 422. 5.78e-89 0 NA 1 TRUE TRUE NA NA NA