multiple_regression_v2

[14 files] Data and R code used for multiple regression analysis of neutral diversity. The raw unprocessed data are given in the files "set80.Zupdated.realign80.GQ30.4D.lowCUB.autoScafs.PiDxyGC.w100m150s100.csv" (pi, dxy and gc content for 100 kb windows); "set80.Zupdated.realign80.GQ30.4D.autoScafs.PiDxyGC.w100m250s100.csv" (the same but only for low-CUB genes); "mel1_hec1_wal1_hecu1_era1.realign80.GQ30.ALLSITES.cons.codeml_nW.w100m500.csv" (Paml analysis results, where branch 8..1 refers to the branch leading to H. melpomene); and "MK_gene_a.csv" (gene by gene a estimates). The "model data" files give the processed data that is imported by the R scripts to run the model, as described in the text.