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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	37	-1270.0331	2621.3511	0.3597	0.3597
K2PI	37	-1270.3372	2621.9594	0.2654	0.6252
K2PIG	38	-1269.9591	2623.6169	0.1159	0.7410
HKYG	40	-1268.0446	2624.6532	0.0690	0.8100
SYMG	41	-1267.1242	2625.2642	0.0508	0.8609
HKYI	40	-1268.3632	2625.2903	0.0502	0.9111
SYMI	41	-1267.5926	2626.2008	0.0318	0.9429
HKYIG	41	-1267.9677	2626.9512	0.0219	0.9648
SYMIG	42	-1267.0507	2627.5817	0.0160	0.9808
GTRG	44	-1265.0247	2628.4979	0.0101	0.9908
GTRI	44	-1265.5502	2629.5489	0.0060	0.9968
GTRIG	45	-1264.9407	2630.8337	0.0031	1.0000
K2P	36	-1280.6068	2640.0973	0.0000	1.0000
SYM	40	-1277.3791	2643.3221	0.0000	1.0000
HKY	39	-1278.7789	2643.6828	0.0000	1.0000
GTR	43	-1275.7352	2647.4284	0.0000	1.0000
JC69G	36	-1292.9330	2664.7497	0.0000	1.0000
JC69I	36	-1293.1978	2665.2793	0.0000	1.0000
JC69IG	37	-1292.8625	2667.0100	0.0000	1.0000
F81G	39	-1291.0215	2668.1681	0.0000	1.0000
F81I	39	-1291.3048	2668.7346	0.0000	1.0000
F81IG	40	-1290.9506	2670.4652	0.0000	1.0000
JC69	35	-1303.5723	2683.6395	0.0000	1.0000
F81	38	-1301.8568	2687.4123	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-1270.0331	2614.0661	0.2592	0.2592
K2PI	37	-1270.3372	2614.6744	0.1913	0.4505
K2PIG	38	-1269.9591	2615.9182	0.1027	0.5532
HKYG	40	-1268.0446	2616.0892	0.0943	0.6475
SYMG	41	-1267.1242	2616.2485	0.0871	0.7345
HKYI	40	-1268.3632	2616.7264	0.0686	0.8031
SYMI	41	-1267.5926	2617.1851	0.0545	0.8576
HKYIG	41	-1267.9677	2617.9355	0.0375	0.8950
GTRG	44	-1265.0247	2618.0494	0.0354	0.9304
SYMIG	42	-1267.0507	2618.1014	0.0345	0.9649
GTRI	44	-1265.5502	2619.1004	0.0209	0.9858
GTRIG	45	-1264.9407	2619.8813	0.0142	1.0000
K2P	36	-1280.6068	2633.2136	0.0000	1.0000
SYM	40	-1277.3791	2634.7582	0.0000	1.0000
HKY	39	-1278.7789	2635.5578	0.0000	1.0000
GTR	43	-1275.7352	2637.4705	0.0000	1.0000
JC69G	36	-1292.9330	2657.8660	0.0000	1.0000
JC69I	36	-1293.1978	2658.3956	0.0000	1.0000
JC69IG	37	-1292.8625	2659.7251	0.0000	1.0000
F81G	39	-1291.0215	2660.0431	0.0000	1.0000
F81I	39	-1291.3048	2660.6096	0.0000	1.0000
F81IG	40	-1290.9506	2661.9012	0.0000	1.0000
JC69	35	-1303.5723	2677.1447	0.0000	1.0000
F81	38	-1301.8568	2679.7136	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1270.0331	2763.9063	0.5580	0.5580
K2PI	37	-1270.3372	2764.5145	0.4116	0.9696
K2PIG	38	-1269.9591	2769.8081	0.0292	0.9988
HKYG	40	-1268.0446	2778.0786	0.0005	0.9992
HKYI	40	-1268.3632	2778.7157	0.0003	0.9996
K2P	36	-1280.6068	2779.0040	0.0003	0.9999
SYMG	41	-1267.1242	2782.2876	0.0001	0.9999
SYMI	41	-1267.5926	2783.2242	0.0000	1.0000
HKYIG	41	-1267.9677	2783.9745	0.0000	1.0000
SYMIG	42	-1267.0507	2788.1902	0.0000	1.0000
HKY	39	-1278.7789	2793.4974	0.0000	1.0000
GTRG	44	-1265.0247	2796.2377	0.0000	1.0000
SYM	40	-1277.3791	2796.7475	0.0000	1.0000
GTRI	44	-1265.5502	2797.2886	0.0000	1.0000
GTRIG	45	-1264.9407	2802.1193	0.0000	1.0000
JC69G	36	-1292.9330	2803.6564	0.0000	1.0000
JC69I	36	-1293.1978	2804.1860	0.0000	1.0000
JC69IG	37	-1292.8625	2809.5652	0.0000	1.0000
GTR	43	-1275.7352	2811.6090	0.0000	1.0000
F81G	39	-1291.0215	2817.9827	0.0000	1.0000
F81I	39	-1291.3048	2818.5492	0.0000	1.0000
JC69	35	-1303.5723	2818.8854	0.0000	1.0000
F81IG	40	-1290.9506	2823.8905	0.0000	1.0000
F81	38	-1301.8568	2833.6035	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6301.phylip = 1 - 424;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6301.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;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6301.phylip = 1 - 424;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6301.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;


[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/chr7_6301.phylip = 1 - 424;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6301.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

