Rors and Jats

MODEL 1
target = Ror
left = c('Onge.SG', 'Turkmenistan_Gonur_BA_2', 'Russia_MLBA_Sintashta.SG', 'Kazakhstan_Kumsay_EBA')
right = c('Mbuti.DG', 'EHG', 'CHG', 'WSHG', 'Iran_GanjDareh_N', 'Turkey_N', 'Levant_PPN', 'Mongolia_North_N', 'Serbia_IronGates_Mesolithic', 'CHN_Tarim_EMBA1', 'China_Tianyuan', 'Dai.DG')
qpadm(prefix, left, right, target, allsnps = TRUE)

$weights
# A tibble: 4 × 5
  target left                     weight     se     z
  <chr>  <chr>                     <dbl>  <dbl> <dbl>
1 Ror    Onge.SG                  0.0871 0.0141  6.19
2 Ror    Turkmenistan_Gonur_BA_2  0.514  0.0247 20.8
3 Ror    Russia_MLBA_Sintashta.SG 0.346  0.0182 19.0
4 Ror    Kazakhstan_Kumsay_EBA    0.0528 0.0136  3.89

$rankdrop
# A tibble: 4 × 7
  f4rank   dof  chisq         p dofdiff chisqdiff   p_nested
   <int> <int>  <dbl>     <dbl>   <int>     <dbl>      <dbl>
1      3     8   11.0 2.00e-  1      10      779.  8.63e-161
2      2    18  790.  5.38e-156      12     2402.  0
3      1    30 3191.  0              14     4151.  0
4      0    44 7342.  0              NA       NA  NA

$popdrop
# A tibble: 15 × 15
   pat      wt   dof  chisq         p f4rank Onge.SG Turkmenistan_Gonur…¹ Russi…² Kazakh…³ feasi…⁴ best  dofdiff chisq…⁵ p_nes…⁶
   <chr> <dbl> <dbl>  <dbl>     <dbl>  <dbl>   <dbl>                <dbl>   <dbl>    <dbl> <lgl>   <lgl>   <dbl>   <dbl>   <dbl>
 1 0000      0     8   11.0 2.00e-  1      3  0.0871                0.514   0.346  0.0528  TRUE    NA         NA     NA       NA
 2 0001      1     9   25.3 2.64e-  3      2  0.0821                0.540   0.378 NA       TRUE    TRUE        0   -169.       1
 3 0010      1     9  195.  4.40e- 37      2 -0.0833                0.960  NA      0.123   FALSE   TRUE        0   -375.       1
 4 0100      1     9  569.  7.61e-117      2  0.283                NA       0.708  0.00834 TRUE    TRUE        0    528.       0
 5 1000      1     9   41.3 4.43e-  6      2 NA                     0.643   0.310  0.0465  TRUE    TRUE       NA     NA       NA
 6 0011      2    10  228.  2.80e- 43      1 -0.149                 1.15   NA     NA       FALSE   NA         NA     NA       NA
 7 0101      2    10  578.  1.02e-117      1  0.286                NA       0.714 NA       TRUE    NA         NA     NA       NA
 8 0110      2    10 2249.  0              1  0.268                NA      NA      0.732   TRUE    NA         NA     NA       NA
 9 1001      2    10   51.2 1.58e-  7      1 NA                     0.661   0.339 NA       TRUE    NA         NA     NA       NA
10 1010      2    10  218.  2.79e- 41      1 NA                     0.863  NA      0.137   TRUE    NA         NA     NA       NA
11 1100      2    10 1334.  1.94e-280      1 NA                    NA       1.28  -0.282   FALSE   NA         NA     NA       NA
12 0111      3    11 3539.  0              0  1                    NA      NA     NA       TRUE    NA         NA     NA       NA
13 1011      3    11  275.  1.60e- 52      0 NA                     1      NA     NA       TRUE    NA         NA     NA       NA
14 1101      3    11 1533.  2.6 e-322      0 NA                    NA       1     NA       TRUE    NA         NA     NA       NA
15 1110      3    11 2589.  0              0 NA                    NA      NA      1       TRUE    NA         NA     NA       NA
# … with abbreviated variable names ¹Turkmenistan_Gonur_BA_2, ²Russia_MLBA_Sintashta.SG, ³Kazakhstan_Kumsay_EBA, ⁴feasible,
#   ⁵chisqdiff, ⁶p_nested

MODEL 2
target = Jat_Haryana
left = c('Onge.SG', 'Turkmenistan_Gonur_BA_2', 'Russia_MLBA_Sintashta', 'Kazakhstan_Kumsay_EBA')
right = c('Mbuti.DG', 'EHG', 'CHG', 'WSHG', 'Iran_GanjDareh_N', 'Turkey_N', 'Levant_PPN', 'Mongolia_North_N', 'Serbia_IronGates_Mesolithic', 'CHN_Tarim_EMBA1', 'China_Tianyuan', 'Dai.DG')
qpadm(prefix, left, right, target, allsnps = TRUE)

$weights
# A tibble: 4 × 5
  target      left                    weight     se     z
  <chr>       <chr>                    <dbl>  <dbl> <dbl>
1 Jat_Haryana Onge.SG                 0.103  0.0249  4.14
2 Jat_Haryana Turkmenistan_Gonur_BA_2 0.512  0.0455 11.2
3 Jat_Haryana Russia_MLBA_Sintashta   0.318  0.0323  9.84
4 Jat_Haryana Kazakhstan_Kumsay_EBA   0.0668 0.0238  2.81

$rankdrop
# A tibble: 4 × 7
  f4rank   dof  chisq         p dofdiff chisqdiff   p_nested
   <int> <int>  <dbl>     <dbl>   <int>     <dbl>      <dbl>
1      3     8   12.4 1.34e-  1      10      575.  4.31e-117
2      2    18  587.  4.28e-113      12     2565.  0
3      1    30 3152.  0              14     4862.  0
4      0    44 8014.  0              NA       NA  NA

$popdrop
# A tibble: 15 × 15
   pat      wt   dof  chisq         p f4rank Onge.SG Turkmenistan_Gonur_…¹ Russi…² Kazak…³ feasi…⁴ best  dofdiff chisq…⁵ p_nes…⁶
   <chr> <dbl> <dbl>  <dbl>     <dbl>  <dbl>   <dbl>                 <dbl>   <dbl>   <dbl> <lgl>   <lgl>   <dbl>   <dbl>   <dbl>
 1 0000      0     8   12.4 1.34e-  1      3  0.103                  0.512   0.318  0.0668 TRUE    NA         NA    NA        NA
 2 0001      1     9   33.5 1.09e-  4      2  0.0975                 0.548   0.355 NA      TRUE    TRUE        0   -91.6       1
 3 0010      1     9  125.  1.21e- 22      2 -0.0511                 0.936  NA      0.115  FALSE   TRUE        0  -165.        1
 4 0100      1     9  290.  2.72e- 57      2  0.314                 NA       0.555  0.131  TRUE    TRUE        0   250.        0
 5 1000      1     9   40.3 6.78e-  6      2 NA                      0.668   0.270  0.0622 TRUE    TRUE       NA    NA        NA
 6 0011      2    10  156.  2.60e- 28      1 -0.109                  1.11   NA     NA      FALSE   NA         NA    NA        NA
 7 0101      2    10  537.  4.48e-109      1  0.332                 NA       0.668 NA      TRUE    NA         NA    NA        NA
 8 0110      2    10 1782.  0              1  0.360                 NA      NA      0.640  TRUE    NA         NA    NA        NA
 9 1001      2    10   61.3 2.05e-  9      1 NA                      0.705   0.295 NA      TRUE    NA         NA    NA        NA
10 1010      2    10  140.  3.81e- 25      1 NA                      0.875  NA      0.125  TRUE    NA         NA    NA        NA
11 1100      2    10 1367.  1.50e-287      1 NA                     NA       0.750  0.250  TRUE    NA         NA    NA        NA
12 0111      3    11 3076.  0              0  1                     NA      NA     NA      TRUE    NA         NA    NA        NA
13 1011      3    11  183.  2.85e- 33      0 NA                      1      NA     NA      TRUE    NA         NA    NA        NA
14 1101      3    11 2136.  0              0 NA                     NA       1     NA      TRUE    NA         NA    NA        NA
15 1110      3    11 2389.  0              0 NA                     NA      NA      1      TRUE    NA         NA    NA        NA
# … with abbreviated variable names ¹Turkmenistan_Gonur_BA_2, ²Russia_MLBA_Sintashta, ³Kazakhstan_Kumsay_EBA, ⁴feasible,
#   ⁵chisqdiff, ⁶p_nested