#take off header from prior sample file. Do for each of the 8 prior sample files corresponding to M1 ... M8
 grep -v "PRI.numTauClass" prior_M7 > prior_M7_headless

# run accept/reject script for each set of prior replicates (ie from each prior model M1 .... M8) to filter each set of 500,000 simulations down to 10,000 each

#obs_ssKansas7.txt is the observed summary statistic  vector calculated  from Oaks et al. (Dryad doi:10.5061/dryad.5s07m)
# using the msbayes function 
#$obsSumStats.pl -T obs_ssKansas7.table batch_table1_prior.txt > obs_ssKansas7.txt
#obs_ssKansas7HEADLESS.txt is the same observed summary statistic  vector without the header and with a dummy first column

rm ~/ABC_Kansas/Posterior_Kansas7_072913_M7headless_10K
cd $HOME/msbayes/trunk/src
./msReject ~/ABC_Kansas/obs_ssKansas7HEADLESS.txt   ~/ABC_Kansas/priorKansas7_072913_M7headless 0.002 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 > ~/ABC_Kansas/Posterior_Kansas7_072913_M7headless_10K
cd ~/ABC_Kansas

# add extra column to prior sample file (prior_M7_headless) with a value corresponding to model indicator value (eg 7 for model 7)

perl -pi -e 's/0\t1\t/7\t1\t/g' prior_M7_headless
perl -pi -e 's/0\t2\t/7\t2\t/g' prior_M7_headless
perl -pi -e 's/0\t3\t/7\t3\t/g' prior_M7_headless
perl -pi -e 's/0\t4\t/7\t4\t/g' prior_M7_headless
perl -pi -e 's/0\t5\t/7\t5\t/g' prior_M7_headless
perl -pi -e 's/0\t6\t/7\t6\t/g' prior_M7_headless
perl -pi -e 's/0\t7\t/7\t7\t/g' prior_M7_headless
perl -pi -e 's/0\t8\t/7\t8\t/g' prior_M7_headless
perl -pi -e 's/0\t9\t/7\t9\t/g' prior_M7_headless
perl -pi -e 's/0\t10\t/7\t10\t/g' prior_M7_headless
perl -pi -e 's/0\t11\t/7\t11\t/g' prior_M7_headless
perl -pi -e 's/0\t12\t/7\t12\t/g' prior_M7_headless
perl -pi -e 's/0\t13\t/7\t13\t/g' prior_M7_headless
perl -pi -e 's/0\t14\t/7\t14\t/g' prior_M7_headless
perl -pi -e 's/0\t15\t/7\t15\t/g' prior_M7_headless
perl -pi -e 's/0\t16\t/7\t16\t/g' prior_M7_headless
perl -pi -e 's/0\t17\t/7\t17\t/g' prior_M7_headless
perl -pi -e 's/0\t18\t/7\t18\t/g' prior_M7_headless
perl -pi -e 's/0\t19\t/7\t19\t/g' prior_M7_headless
perl -pi -e 's/0\t20\t/7\t20\t/g' prior_M7_headless
perl -pi -e 's/0\t21\t/7\t21\t/g' prior_M7_headless
perl -pi -e 's/0\t22\t/7\t22\t/g' prior_M7_headless



#concatenate the the eight sets of prior replicates (ie from each prior model M1 .... M8)

cat Posterior_M1_headless Posterior_M2_headless Posterior_M3_headless Posterior_M4_headless Posterior_M5_headless Posterior_M6_headless Posterior_M7_headless Posterior_M8_headless > Posterior_M12345678

# set the proportion Z.ZZ to run accept/reject script on Posterior_M12345678 to distill posterior sample down to final desired number of accepted prior simulations to estimate the posterior strength of each model prior

cd $HOME/msbayes/trunk/src
rm Posterior_M12345678_Z.ZZ
./msReject obs_ssKansas7HEADLESS.txt Posterior_M12345678 Z.ZZ 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 > Posterior_M12345678_Z.ZZ