library(grur);library(radiator)


setwd("~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/06.Population_genetics/combined")

het.prob <- radiator::detect_het_outliers(
  data = "combined.recode.vcf", strata = "~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/strata_lom_ddr.tsv", nreps = 10000)

het.prob <- radiator::detect_het_outliers(
  data = "combined.recode.vcf", strata = "~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/strata_lom_epi.tsv", nreps = 10000)

het.prob <- radiator::detect_het_outliers(
  data = "combined.recode.vcf", strata = "~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/strata_eck_ddr.tsv", nreps = 10000)

het.prob <- radiator::detect_het_outliers(
  data = "combined.recode.vcf", strata = "~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/strata_eck_epi.tsv", nreps = 10000)

het.prob <- radiator::detect_het_outliers(
  data = "combined.recode.vcf", strata = "~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/strata_sloy_ddr.tsv", nreps = 10000)

het.prob <- radiator::detect_het_outliers(
  data = "combined.recode.vcf", strata = "~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/strata_sloy_epi.tsv", nreps = 10000)