Post date: Feb 16, 2014 12:15:43 AM
I ran treemix with the -se option to get p-values for migration events on dorc, but these times out. I am trying again on my computer treating each SNPs as independent or each set of 100 SNPs as independent (-k 100). Here is the code (I am using migration events = 4, 6, and 7, all with rep 0, which is the higher of the two likelihoods for #m = 6).
treemix -i trmx_in.txt.gz -g trmxout_mle_rootanna__m4rep0.vertices.gz trmxout_mle_rootanna__m4rep0.edges.gz -se -o trmxout_mle_rootanna_m4rep0_se &
treemix -i trmx_in.txt.gz -g trmxout_mle_rootanna__m6rep0.vertices.gz trmxout_mle_rootanna__m6rep0.edges.gz -se -o trmxout_mle_rootanna_m6rep0_se &
treemix -i trmx_in.txt.gz -g trmxout_mle_rootanna__m7rep0.vertices.gz trmxout_mle_rootanna__m7rep0.edges.gz -se -o trmxout_mle_rootanna_m7rep0_se &
treemix -i trmx_in.txt.gz -k 100 -g trmxout_mle_rootanna__m4rep0.vertices.gz trmxout_mle_rootanna__m4rep0.edges.gz -se -o trmxout_mle_rootanna_m4rep0_sek100 &
treemix -i trmx_in.txt.gz -k 100 -g trmxout_mle_rootanna__m6rep0.vertices.gz trmxout_mle_rootanna__m6rep0.edges.gz -se -o trmxout_mle_rootanna_m6rep0_sek100 &
treemix -i trmx_in.txt.gz -k 100 -g trmxout_mle_rootanna__m7rep0.vertices.gz trmxout_mle_rootanna__m7rep0.edges.gz -se -o trmxout_mle_rootanna_m7rep0_sek100 &