G25 model
qpAdm models:
qpAdm right pops
right = c('Mbuti.DG', 'Russia_Ust_Ishim.DG', 'Russia_DevilsCave_N.SG', 'Russia_MA1_HG.SG', 'Czech_Vestonice16', 'Morocco_Iberomaurusian', 'Serbia_IronGates_Mesolithic', 'Russia_Samara_EBA_Yamnaya', 'Anatolia_Neolithic_6000_BC', 'Jordan_PPNB', 'Iran_GanjDareh_N', 'Croatia_MBA', 'Levant_BA_1900_1500_BC', 'Mycenaean_Greek_1500_1050_BC', 'Turkey_Hellenistic_Halikarnassos', 'Latvia_BA', 'Russia_IA_Ingria.SG', 'Albania_1500_400_BC')
left = c('Aegean_Imperial_1_200ad', 'Levant_BA_1900_1500_BC_Sidon')
target = c('Greek_Cypriot')
results = qpadm(prefix, left, right, target, allsnps = TRUE)
results$weights
results$popdrop
right = c('Mbuti.DG', 'Russia_Ust_Ishim.DG', 'Russia_DevilsCave_N.SG', 'Russia_MA1_HG.SG', 'Czech_Vestonice16', 'Morocco_Iberomaurusian', 'Serbia_IronGates_Mesolithic', 'Russia_Samara_EBA_Yamnaya', 'Anatolia_Neolithic_6000_BC', 'Jordan_PPNB', 'Iran_GanjDareh_N', 'Croatia_MBA', 'Levant_BA_1900_1500_BC', 'Mycenaean_Greek_1500_1050_BC', 'Minoan_Crete_2600_1700_BC', 'Turkey_Hellenistic_Halikarnassos', 'Latvia_BA', 'Russia_IA_Ingria.SG', 'Albania_1500_400_BC', 'North_Macedonia_500_100_BC')
left = c('Aegean_Imperial_1_200ad', 'Levant_BA_1900_1500_BC_Sidon', 'Russia_Sunghir_Medieval.SG')
target = c('Greek_Cypriot')
--------------------------------------------------
P 0.157
target left weight se z
<chr> <chr> <dbl> <dbl> <dbl>
1 Greek_Cypriot Aegean_Imperial_1_200ad 0.680 0.0828 8.21
2 Greek_Cypriot Levant_BA_1900_1500_BC_Sidon 0.268 0.0654 4.10
3 Greek_Cypriot Russia_Sunghir_Medieval.SG 0.0523 0.0224 2.33
-------------------------------------------------------------------------------------------
right = c('Mbuti.DG', 'Russia_Ust_Ishim.DG', 'Russia_DevilsCave_N.SG', 'Russia_MA1_HG.SG', 'Czech_Vestonice16', 'Morocco_Iberomaurusian', 'Serbia_IronGates_Mesolithic', 'Russia_Samara_EBA_Yamnaya', 'Anatolia_Neolithic_6000_BC', 'Jordan_PPNB', 'Iran_GanjDareh_N', 'Croatia_MBA', 'Levant_BA_1900_1500_BC', 'Turkey_Alalakh_MLBA', 'Mycenaean_Greek_1500_1050_BC', 'Minoan_Crete_2600_1700_BC', 'Anatolia_Central_Gordion_780_544BC', 'Turkey_Hellenistic_Halikarnassos', 'Latvia_BA', 'Russia_IA_Ingria.SG', 'Albania_1500_400_BC', 'North_Macedonia_500_100_BC')
left = c('Aegean_Imperial_1_200ad', 'Levant_BA_1900_1500_BC_Sidon', 'Russia_Sunghir_Medieval.SG')
target = c('Greek_Cypriot')
-------------------------------------------------------------------------------------------
P 0.0386
target left weight se z
<chr> <chr> <dbl> <dbl> <dbl>
1 Greek_Cypriot Aegean_Imperial_1_200ad 0.725 0.0863 8.40
2 Greek_Cypriot Levant_BA_1900_1500_BC_Sidon 0.237 0.0676 3.50
3 Greek_Cypriot Russia_Sunghir_Medieval.SG 0.0384 0.0232 1.66
No comments:
Post a Comment