### written by Kevin Schneider ### ### last edited: 29.07.2020 ### ### Kevin.Schneider@glasgow.ac.uk ### ### KevinSchneider@gmx.at ### require(reshape2) require(ggplot2) require(rcartocolor) require(plyr) require(stringr) require(ggpubr) require(cowplot) setwd("/Volumes/LaCie/demographic_simulation_paper/fvec_ms_and_vcf_simulation_output_20200709/") path <- "/Volumes/LaCie/demographic_simulation_paper/fvec_ms_and_vcf_simulation_output_20200709/" path2 <- "/Volumes/LaCie/demographic_simulation_paper/ms_and_vcf_simulation_output_20200709/" path3 <- "/Volumes/LaCie/demographic_simulation_paper/ms_and_vcf_simulation_output_20200709/rds_dir/" # specifying metric input files scenario_vec <- c("scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.80K.fvec") scenario_vecb <- c("scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p3.80K.fvec") beta_vec <- c("scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.rds", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.rds") fst_vec <- c("scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst", "scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample.80001.noNA.fst.windowed.weir.fst") hscan_vec <- c("hscan_scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario1_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario2_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario3_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_1_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_01_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_hard_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_constant_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_unequal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms", "hscan_scenario4_sim_sel0_1_increasing_equal_popsize_migrate0_001_soft_sweep.*.sample_p2.ms") title_vec <- c("scenario 1 hard sweep (constant, unequal popsize; migrate 0.1)", "scenario 1 hard sweep (constant, equal popsize; migrate 0.1)", "scenario 1 hard sweep (growing, unequal popsize; migrate 0.1)", "scenario 1 hard sweep (growing, equal popsize; migrate 0.1)", "scenario 1 soft sweep (constant, unequal popsize; migrate 0.1)", "scenario 1 soft sweep (constant, equal popsize; migrate 0.1)", "scenario 1 soft sweep (growing, unequal popsize; migrate 0.1)", "scenario 1 soft sweep (growing, equal popsize; migrate 0.1)", "scenario 1 hard sweep (constant, unequal popsize; migrate 0.01)", "scenario 1 hard sweep (constant, equal popsize; migrate 0.01)", "scenario 1 hard sweep (growing, unequal popsize; migrate 0.01)", "scenario 1 hard sweep (growing, equal popsize; migrate 0.01)", "scenario 1 soft sweep (constant, unequal popsize; migrate 0.01)", "scenario 1 soft sweep (constant, equal popsize; migrate 0.01)", "scenario 1 soft sweep (growing, unequal popsize; migrate 0.01)", "scenario 1 soft sweep (growing, equal popsize; migrate 0.01)", "scenario 1 hard sweep (constant, unequal popsize; migrate 0.001)", "scenario 1 hard sweep (constant, equal popsize; migrate 0.001)", "scenario 1 hard sweep (growing, unequal popsize; migrate 0.001)", "scenario 1 hard sweep (growing, equal popsize; migrate 0.001)", "scenario 1 soft sweep (constant, unequal popsize; migrate 0.001)", "scenario 1 soft sweep (constant, equal popsize; migrate 0.001)", "scenario 1 soft sweep (growing, unequal popsize; migrate 0.001)", "scenario 1 soft sweep (growing, equal popsize; migrate 0.001)", "scenario 2 hard sweep (constant, unequal popsize; migrate 0.1)", "scenario 2 hard sweep (constant, equal popsize; migrate 0.1)", "scenario 2 hard sweep (growing, unequal popsize; migrate 0.1)", "scenario 2 hard sweep (growing, equal popsize; migrate 0.1)", "scenario 2 soft sweep (constant, unequal popsize; migrate 0.1)", "scenario 2 soft sweep (constant, equal popsize; migrate 0.1)", "scenario 2 soft sweep (growing, unequal popsize; migrate 0.1)", "scenario 2 soft sweep (growing, equal popsize; migrate 0.1)", "scenario 2 hard sweep (constant, unequal popsize; migrate 0.01)", "scenario 2 hard sweep (constant, equal popsize; migrate 0.01)", "scenario 2 hard sweep (growing, unequal popsize; migrate 0.01)", "scenario 2 hard sweep (growing, equal popsize; migrate 0.01)", "scenario 2 soft sweep (constant, unequal popsize; migrate 0.01)", "scenario 2 soft sweep (constant, equal popsize; migrate 0.01)", "scenario 2 soft sweep (growing, unequal popsize; migrate 0.01)", "scenario 2 soft sweep (growing, equal popsize; migrate 0.01)", "scenario 2 hard sweep (constant, unequal popsize; migrate 0.001)", "scenario 2 hard sweep (constant, equal popsize; migrate 0.001)", "scenario 2 hard sweep (growing, unequal popsize; migrate 0.001)", "scenario 2 hard sweep (growing, equal popsize; migrate 0.001)", "scenario 2 soft sweep (constant, unequal popsize; migrate 0.001)", "scenario 2 soft sweep (constant, equal popsize; migrate 0.001)", "scenario 2 soft sweep (growing, unequal popsize; migrate 0.001)", "scenario 2 soft sweep (growing, equal popsize; migrate 0.001)", "scenario 3 hard sweep (constant, unequal popsize; migrate 0.1)", "scenario 3 hard sweep (constant, equal popsize; migrate 0.1)", "scenario 3 hard sweep (growing, unequal popsize; migrate 0.1)", "scenario 3 hard sweep (growing, equal popsize; migrate 0.1)", "scenario 3 soft sweep (constant, unequal popsize; migrate 0.1)", "scenario 3 soft sweep (constant, equal popsize; migrate 0.1)", "scenario 3 soft sweep (growing, unequal popsize; migrate 0.1)", "scenario 3 soft sweep (growing, equal popsize; migrate 0.1)", "scenario 3 hard sweep (constant, unequal popsize; migrate 0.01)", "scenario 3 hard sweep (constant, equal popsize; migrate 0.01)", "scenario 3 hard sweep (growing, unequal popsize; migrate 0.01)", "scenario 3 hard sweep (growing, equal popsize; migrate 0.01)", "scenario 3 soft sweep (constant, unequal popsize; migrate 0.01)", "scenario 3 soft sweep (constant, equal popsize; migrate 0.01)", "scenario 3 soft sweep (growing, unequal popsize; migrate 0.01)", "scenario 3 soft sweep (growing, equal popsize; migrate 0.01)", "scenario 3 hard sweep (constant, unequal popsize; migrate 0.001)", "scenario 3 hard sweep (constant, equal popsize; migrate 0.001)", "scenario 3 hard sweep (growing, unequal popsize; migrate 0.001)", "scenario 3 hard sweep (growing, equal popsize; migrate 0.001)", "scenario 3 soft sweep (constant, unequal popsize; migrate 0.001)", "scenario 3 soft sweep (constant, equal popsize; migrate 0.001)", "scenario 3 soft sweep (growing, unequal popsize; migrate 0.001)", "scenario 3 soft sweep (growing, equal popsize; migrate 0.001)", "scenario 4 hard sweep (constant, unequal popsize; migrate 0.1)", "scenario 4 hard sweep (constant, equal popsize; migrate 0.1)", "scenario 4 hard sweep (growing, unequal popsize; migrate 0.1)", "scenario 4 hard sweep (growing, equal popsize; migrate 0.1)", "scenario 4 soft sweep (constant, unequal popsize; migrate 0.1)", "scenario 4 soft sweep (constant, equal popsize; migrate 0.1)", "scenario 4 soft sweep (growing, unequal popsize; migrate 0.1)", "scenario 4 soft sweep (growing, equal popsize; migrate 0.1)", "scenario 4 hard sweep (constant, unequal popsize; migrate 0.01)", "scenario 4 hard sweep (constant, equal popsize; migrate 0.01)", "scenario 4 hard sweep (growing, unequal popsize; migrate 0.01)", "scenario 4 hard sweep (growing, equal popsize; migrate 0.01)", "scenario 4 soft sweep (constant, unequal popsize; migrate 0.01)", "scenario 4 soft sweep (constant, equal popsize; migrate 0.01)", "scenario 4 soft sweep (growing, unequal popsize; migrate 0.01)", "scenario 4 soft sweep (growing, equal popsize; migrate 0.01)", "scenario 4 hard sweep (constant, unequal popsize; migrate 0.001)", "scenario 4 hard sweep (constant, equal popsize; migrate 0.001)", "scenario 4 hard sweep (growing, unequal popsize; migrate 0.001)", "scenario 4 hard sweep (growing, equal popsize; migrate 0.001)", "scenario 4 soft sweep (constant, unequal popsize; migrate 0.001)", "scenario 4 soft sweep (constant, equal popsize; migrate 0.001)", "scenario 4 soft sweep (growing, unequal popsize; migrate 0.001)", "scenario 4 soft sweep (growing, equal popsize; migrate 0.001)") # parsing and changing format of scenario names filename_vec <- sapply(title_vec, function(x) gsub(" ", "_", x)) filename_vec <- sapply(filename_vec, function(x) str_remove(x, "\\(")) filename_vec <- sapply(filename_vec, function(x) str_remove(x, "\\)")) filename_vec <- sapply(filename_vec, function(x) str_remove(x, ",")) filename_vec <- sapply(filename_vec, function(x) str_remove(x, ";")) filename_vec <- sapply(filename_vec, function(x) gsub("0\\.", "0_", x)) plot_list <- list() for (i in 1:length(title_vec)) { glob_vec <- Sys.glob(file.path(paste(path, scenario_vec[i], sep = ""))) glob_vecb <- Sys.glob(file.path(paste(path, scenario_vecb[i], sep = ""))) glob_vec2 <- Sys.glob(file.path(paste(path2, fst_vec[i], sep = ""))) glob_vec3 <- Sys.glob(file.path(paste(path3, beta_vec[i], sep = ""))) glob_vec4 <- Sys.glob(file.path(paste(path2, hscan_vec[i], ".windows", sep = ""))) df_list <- list() for (j in 1:100) { # reading in metric input files fv <- read.table(glob_vec[j], header = TRUE) fv2 <- read.table(glob_vecb[j], header = TRUE) statvars <- colsplit(colnames(fv), "_win", c("metric", "window")) fv <- unlist(sapply(fv, unlist)) fv <- as.data.frame(matrix(fv, nrow = 50)) rownames(fv) <- 1:50 colnames(fv) <- unique(statvars[, 1]) # fv2 <- unlist(sapply(fv2, unlist)) fv2 <- as.data.frame(matrix(fv2, nrow = 50)) rownames(fv2) <- 1:50 colnames(fv2) <- unique(statvars[, 1]) # fst_frame <- read.table(glob_vec2[j], header = TRUE) beta_frame <- readRDS(glob_vec3[j]) hscan_frame <- read.table(glob_vec4[j], header = TRUE) fv <- data.frame(fv, "H-scan" = hscan_frame[1:50, 2] / 1000000) fv <- data.frame(fv, "delta_pi" = abs(fv[1:50, 1] - fv2[1:50, 1]) * 10) fv <- data.frame(fv, "FST" = fst_frame[1:50, 6]) fv <- data.frame(fv, "BbarGD" = beta_frame[1:50, 2:4]) fv <- data.frame(fv, "window_num" = ((1:50) * 4000050 / 50) - 40000) colnames(fv)[1:12] <- c("pi", "theta_hat_W", "TajD", "var_gkl", "skew_gkl", "kurt_gkl", "nDiplos", "J1", "J12", "J2/J1", "diplo_ZNS", "diplo_omega") df_list[[j]] <- fv } # calculating means of metrics for each scenario df_mean <- as.data.frame(aaply(laply(df_list, as.matrix), c(2, 3), mean, na.rm = TRUE)) df2plot <- melt(df_mean, id = c("window_num")) df2plot1 <- df2plot[which(df2plot$variable %in% levels(df2plot$variable)[1:6]), ] df2plot2 <- df2plot[which(df2plot$variable %in% levels(df2plot$variable)[7:13]), ] df2plot3 <- df2plot[which(df2plot$variable %in% levels(df2plot$variable)[14:15]), ] df2plot4 <- df2plot[which(df2plot$variable %in% levels(df2plot$variable)[16:18]), ] my_plot1 <- ggplot(df2plot1) + geom_line(aes(x = window_num, y = value, colour = variable), size = 1.5) + scale_colour_manual(values = carto_pal(12, "Bold")[1:12]) + xlab(label = "Genomic Position [bp]") + ylab(label = "Summary Statistic") + coord_cartesian(ylim = c(0, 0.06)) + theme_bw() + theme(legend.title = element_blank(), legend.text = element_text(size = 18), axis.text = element_text(size = 18), axis.title = element_text(size = 20), plot.title = element_text(size = 24)) + geom_rect(data = data.frame(xmin = 2000026, xmax = 4000050, ymin = -Inf, ymax = Inf), aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax), fill = "grey", alpha = 0.2) + geom_vline(xintercept = 1000013, linetype = "dotted", color = "black", size = 1.5) + geom_vline(xintercept = 3000038, linetype = "dotted", color = "black", size = 1.5) plot_list[[i]] <- my_plot1 my_plot2 <- ggplot(df2plot2) + geom_line(aes(x = window_num, y = value, colour = variable), size = 1.5) + scale_colour_manual(values = carto_pal(12, "Bold")[1:12]) + xlab(label = "Genomic Position [bp]") + ylab(label = "Summary Statistic") + coord_cartesian(ylim = c(0, 0.06)) + theme_bw() + theme(legend.title = element_blank(), legend.text = element_text(size = 18), axis.text = element_text(size = 18), axis.title = element_text(size = 20), plot.title = element_text(size = 24)) + geom_rect(data = data.frame(xmin = 2000026, xmax = 4000050, ymin = -Inf, ymax = Inf), aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax), fill = "grey", alpha = 0.2) + geom_vline(xintercept = 1000013, linetype = "dotted", color = "black", size = 1.5) + geom_vline(xintercept = 3000038, linetype = "dotted", color = "black", size = 1.5) plot_list[[i]] <- my_plot2 my_plot3 <- ggplot(df2plot3) + geom_line(aes(x = window_num, y = value, colour = variable), size = 1.5) + scale_colour_manual(values = carto_pal(12, "Bold")[1:12]) + xlab(label = "Genomic Position [bp]") + ylab(label = "Summary Statistic") + coord_cartesian(ylim = c(0, 0.5)) + theme_bw() + theme(legend.title = element_blank(), legend.text = element_text(size = 18), axis.text = element_text(size = 18), axis.title = element_text(size = 20), plot.title = element_text(size = 24)) + geom_rect(data = data.frame(xmin = 2000026, xmax = 4000050, ymin = -Inf, ymax = Inf), aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax), fill = "grey", alpha = 0.2) + geom_vline(xintercept = 1000013, linetype = "dotted", color = "black", size = 1.5) + geom_vline(xintercept = 3000038, linetype = "dotted", color = "black", size = 1.5) plot_list[[i]] <- my_plot3 my_plot4 <- ggplot(df2plot4) + geom_line(aes(x = window_num, y = value, colour = variable), size = 1.5) + scale_colour_manual(values = carto_pal(12, "Bold")[1:12]) + xlab(label = "Genomic Position [bp]") + ylab(label = "Summary Statistic") + coord_cartesian(ylim = c(1, 1.01)) + theme_bw() + theme(legend.title = element_blank(), legend.text = element_text(size = 18), axis.text = element_text(size = 18), axis.title = element_text(size = 20), plot.title = element_text(size = 24)) + geom_rect(data = data.frame(xmin = 2000026, xmax = 4000050, ymin = -Inf, ymax = Inf), aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax), fill = "grey", alpha = 0.2) + geom_vline(xintercept = 1000013, linetype = "dotted", color = "black", size = 1.5) + geom_vline(xintercept = 3000038, linetype = "dotted", color = "black", size = 1.5) plot_list[[i]] <- my_plot4 comb_plot <- cowplot::plot_grid(my_plot1, my_plot2, my_plot3, my_plot4, ncol = 2, align = "v") comb_plot_ann <- annotate_figure(comb_plot, top = text_grob(title_vec[i], size = 30), fig.lab.pos = "top") cowplot::save_plot(paste("/Volumes/LaCie/demographic_simulation_paper/fvec_ms_and_vcf_simulation_output_20200709/", filename_vec[[i]], ".png", sep = ""), comb_plot_ann, base_width = 20, base_height = 15, limitsize = FALSE) }