
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:45 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/chr6_4126.phylip" (19 taxa, 578 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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
HKYI	40	-1274.1856	2634.4791	0.7018	0.7018
HKYIG	41	-1274.1651	2636.7555	0.2249	0.9267
K2PI	37	-1280.5024	2640.2123	0.0399	0.9666
GTRI	44	-1273.4223	2642.2742	0.0142	0.9808
K2PIG	38	-1280.5027	2642.5044	0.0127	0.9935
GTRIG	45	-1273.3962	2644.5744	0.0045	0.9980
SYMI	41	-1279.4180	2647.2614	0.0012	0.9992
SYMIG	42	-1278.7746	2648.3005	0.0007	0.9999
HKYG	40	-1283.1446	2652.3972	0.0001	1.0000
K2PG	37	-1289.8563	2658.9199	0.0000	1.0000
GTRG	44	-1282.2617	2659.9531	0.0000	1.0000
F81IG	40	-1287.6364	2661.3808	0.0000	1.0000
F81I	39	-1289.5242	2662.8476	0.0000	1.0000
F81G	39	-1290.2485	2664.2963	0.0000	1.0000
JC69I	36	-1294.1188	2665.1618	0.0000	1.0000
SYMG	41	-1289.1675	2666.7604	0.0000	1.0000
JC69IG	37	-1294.1191	2667.4456	0.0000	1.0000
JC69G	36	-1301.8531	2680.6304	0.0000	1.0000
HKY	39	-1316.6660	2717.1312	0.0000	1.0000
GTR	43	-1315.4495	2723.9852	0.0000	1.0000
K2P	36	-1324.1950	2725.3143	0.0000	1.0000
F81	38	-1321.9467	2725.3924	0.0000	1.0000
SYM	40	-1323.5800	2733.2679	0.0000	1.0000
JC69	35	-1335.5213	2745.6921	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1274.1856	2628.3711	0.6744	0.6744
HKYIG	41	-1274.1651	2630.3302	0.2532	0.9277
GTRI	44	-1273.4223	2634.8445	0.0265	0.9542
K2PI	37	-1280.5024	2635.0049	0.0245	0.9787
GTRIG	45	-1273.3962	2636.7925	0.0100	0.9887
K2PIG	38	-1280.5027	2637.0053	0.0090	0.9977
SYMI	41	-1279.4180	2640.8361	0.0013	0.9990
SYMIG	42	-1278.7746	2641.5491	0.0009	0.9999
HKYG	40	-1283.1446	2646.2892	0.0001	1.0000
GTRG	44	-1282.2617	2652.5235	0.0000	1.0000
K2PG	37	-1289.8563	2653.7125	0.0000	1.0000
F81IG	40	-1287.6364	2655.2728	0.0000	1.0000
F81I	39	-1289.5242	2657.0483	0.0000	1.0000
F81G	39	-1290.2485	2658.4970	0.0000	1.0000
JC69I	36	-1294.1188	2660.2376	0.0000	1.0000
SYMG	41	-1289.1675	2660.3351	0.0000	1.0000
JC69IG	37	-1294.1191	2662.2382	0.0000	1.0000
JC69G	36	-1301.8531	2675.7062	0.0000	1.0000
HKY	39	-1316.6660	2711.3320	0.0000	1.0000
GTR	43	-1315.4495	2716.8990	0.0000	1.0000
F81	38	-1321.9467	2719.8933	0.0000	1.0000
K2P	36	-1324.1950	2720.3900	0.0000	1.0000
SYM	40	-1323.5800	2727.1599	0.0000	1.0000
JC69	35	-1335.5213	2741.0426	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-1280.5024	2796.3091	0.9231	0.9231
K2PIG	38	-1280.5027	2802.6691	0.0384	0.9615
HKYI	40	-1274.1856	2802.7541	0.0368	0.9983
HKYIG	41	-1274.1651	2809.0727	0.0016	0.9999
K2PG	37	-1289.8563	2815.0167	0.0001	1.0000
JC69I	36	-1294.1188	2817.1822	0.0000	1.0000
SYMI	41	-1279.4180	2819.5786	0.0000	1.0000
HKYG	40	-1283.1446	2820.6721	0.0000	1.0000
JC69IG	37	-1294.1191	2823.5425	0.0000	1.0000
SYMIG	42	-1278.7746	2824.6512	0.0000	1.0000
GTRI	44	-1273.4223	2826.6658	0.0000	1.0000
F81I	39	-1289.5242	2827.0717	0.0000	1.0000
F81G	39	-1290.2485	2828.5204	0.0000	1.0000
F81IG	40	-1287.6364	2829.6557	0.0000	1.0000
JC69G	36	-1301.8531	2832.6508	0.0000	1.0000
GTRIG	45	-1273.3962	2832.9733	0.0000	1.0000
SYMG	41	-1289.1675	2839.0776	0.0000	1.0000
GTRG	44	-1282.2617	2844.3448	0.0000	1.0000
K2P	36	-1324.1950	2877.3347	0.0000	1.0000
HKY	39	-1316.6660	2881.3554	0.0000	1.0000
F81	38	-1321.9467	2885.5571	0.0000	1.0000
JC69	35	-1335.5213	2893.6277	0.0000	1.0000
SYM	40	-1323.5800	2901.5429	0.0000	1.0000
GTR	43	-1315.4495	2904.3607	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_4126.phylip = 1 - 578;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_4126.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_4126.phylip = 1 - 578;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_4126.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

