
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:17 02/24/2011

Input data from file "/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip" (19 taxa, 422 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	44	-1507.1371	3112.7782	0.3016	0.3016
SYMG	41	-1511.3692	3113.8016	0.1808	0.4823
K2PG	37	-1516.4425	3114.2079	0.1475	0.6299
GTRIG	45	-1507.1013	3115.2132	0.0893	0.7192
HKYG	40	-1513.7546	3116.1181	0.0568	0.7759
SYMIG	42	-1511.3269	3116.1841	0.0549	0.8309
GTRI	44	-1508.9694	3116.4427	0.0483	0.8791
K2PI	37	-1517.5654	3116.4536	0.0480	0.9271
SYMI	41	-1512.9440	3116.9512	0.0374	0.9646
HKYIG	41	-1513.6899	3118.4430	0.0178	0.9823
HKYI	40	-1515.0455	3118.6999	0.0156	0.9979
K2PIG	38	-1519.4952	3122.7294	0.0021	1.0000
GTR	43	-1522.7176	3141.4458	0.0000	1.0000
SYM	40	-1527.1019	3142.8128	0.0000	1.0000
K2P	36	-1533.8458	3146.6110	0.0000	1.0000
HKY	39	-1530.8545	3147.8766	0.0000	1.0000
JC69G	36	-1553.9219	3186.7634	0.0000	1.0000
JC69I	36	-1554.9385	3188.7965	0.0000	1.0000
JC69IG	37	-1553.8746	3189.0721	0.0000	1.0000
F81G	39	-1551.8906	3189.9487	0.0000	1.0000
F81I	39	-1552.9623	3192.0922	0.0000	1.0000
F81IG	40	-1551.8454	3192.2996	0.0000	1.0000
JC69	35	-1568.5857	3213.6999	0.0000	1.0000
F81	38	-1566.5945	3216.9279	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1507.1371	3102.2742	0.4423	0.4423
GTRIG	45	-1507.1013	3104.2026	0.1686	0.6109
SYMG	41	-1511.3692	3104.7384	0.1290	0.7399
GTRI	44	-1508.9694	3105.9387	0.0708	0.8107
SYMIG	42	-1511.3269	3106.6538	0.0495	0.8602
K2PG	37	-1516.4425	3106.8849	0.0441	0.9043
HKYG	40	-1513.7546	3107.5092	0.0323	0.9366
SYMI	41	-1512.9440	3107.8880	0.0267	0.9633
K2PI	37	-1517.5654	3109.1307	0.0143	0.9777
HKYIG	41	-1513.6899	3109.3799	0.0127	0.9904
HKYI	40	-1515.0455	3110.0910	0.0089	0.9992
K2PIG	38	-1519.4952	3114.9905	0.0008	1.0000
GTR	43	-1522.7176	3131.4353	0.0000	1.0000
SYM	40	-1527.1019	3134.2039	0.0000	1.0000
K2P	36	-1533.8458	3139.6915	0.0000	1.0000
HKY	39	-1530.8545	3139.7090	0.0000	1.0000
JC69G	36	-1553.9219	3179.8439	0.0000	1.0000
JC69IG	37	-1553.8746	3181.7492	0.0000	1.0000
F81G	39	-1551.8906	3181.7811	0.0000	1.0000
JC69I	36	-1554.9385	3181.8770	0.0000	1.0000
F81IG	40	-1551.8454	3183.6907	0.0000	1.0000
F81I	39	-1552.9623	3183.9246	0.0000	1.0000
JC69	35	-1568.5857	3207.1714	0.0000	1.0000
F81	38	-1566.5945	3209.1890	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1516.4425	3256.5501	0.7513	0.7513
K2PI	37	-1517.5654	3258.7959	0.2444	0.9957
K2PIG	38	-1519.4952	3268.7007	0.0017	0.9975
HKYG	40	-1513.7546	3269.3094	0.0013	0.9987
SYMG	41	-1511.3692	3270.5836	0.0007	0.9994
HKYI	40	-1515.0455	3271.8912	0.0004	0.9998
SYMI	41	-1512.9440	3273.7332	0.0001	0.9999
HKYIG	41	-1513.6899	3275.2251	0.0001	1.0000
SYMIG	42	-1511.3269	3276.5440	0.0000	1.0000
GTRG	44	-1507.1371	3280.2545	0.0000	1.0000
GTRI	44	-1508.9694	3283.9190	0.0000	1.0000
K2P	36	-1533.8458	3285.3117	0.0000	1.0000
GTRIG	45	-1507.1013	3286.2278	0.0000	1.0000
SYM	40	-1527.1019	3296.0041	0.0000	1.0000
HKY	39	-1530.8545	3297.4642	0.0000	1.0000
GTR	43	-1522.7176	3305.3705	0.0000	1.0000
JC69G	36	-1553.9219	3325.4641	0.0000	1.0000
JC69I	36	-1554.9385	3327.4972	0.0000	1.0000
JC69IG	37	-1553.8746	3331.4144	0.0000	1.0000
F81G	39	-1551.8906	3339.5363	0.0000	1.0000
F81I	39	-1552.9623	3341.6798	0.0000	1.0000
F81IG	40	-1551.8454	3345.4909	0.0000	1.0000
JC69	35	-1568.5857	3348.7466	0.0000	1.0000
F81	38	-1566.5945	3362.8992	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip = 1 - 422;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip = 1 - 422;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip = 1 - 422;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_72.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

