Post date: Jan 26, 2015 10:43:20 PM
#Concatenated all sims:
#To pull out abc sim file from each abc directory (one directory gets made per run) and concatenate them and add one header to the top of the file, no perl script is necessary. Instead:
tail -q -n +2 abcdir*/abc* > abccat.txt
#I overwrote the abccat.txt I made for the tests of correlations between sum stats and param, etc; oops.
wc -l abccat.txt
#total sims = 1164000
cat hdr.txt abccat.txt > abccathead.txt
mv abccathead.txt conc_spsims_27i15.abcout
#Uploaded Splbnds file to cluster.
#Ran R on the cluster:
use R
R
#see file on laptop: stygoparnus_abcanalyses_27i15.txt for R code
#copied two files back to laptop:
Sp.est.regression.csv
Sp.est.rejection.csv
#19Jan15: started copying conc_spsims_27i15.abcout and Spabc.Rwsp from USU cluster to laptop so I could make a combined migration rate figure with the other taxa. Also started copying hetabc.Rwsp from greenhouse to laptop. #I ended up having to get ZG to do this from campus, 1Feb15 and made the figure 2Feb15.
Also made a figure with abc-based thetas to compare to allele frequency thetas to see if the pattern is similar or different. See abctheta_fig_29jan15.R for code. For figure, see fig_abctheta_29jan15.pdf and compare to fig_theta_9dec14.pdf. Patterns are not identical. What to think?