
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:28 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/chr1_14993.phylip" (19 taxa, 878 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-3253.3108	6590.5404	0.7065	0.7065
HKYIG	41	-3253.3069	6592.7335	0.2360	0.9425
HKYI	40	-3256.1795	6596.2777	0.0401	0.9826
GTRG	44	-3252.8820	6598.5178	0.0131	0.9957
GTRIG	45	-3252.8778	6600.7315	0.0043	1.0000
K2PG	37	-3267.7179	6612.7834	0.0000	1.0000
K2PIG	38	-3267.7029	6614.9385	0.0000	1.0000
GTRI	44	-3262.6050	6617.9640	0.0000	1.0000
SYMIG	42	-3266.5412	6621.4082	0.0000	1.0000
SYMG	41	-3273.1548	6632.4293	0.0000	1.0000
K2PI	37	-3277.9021	6633.1517	0.0000	1.0000
SYMI	41	-3276.8311	6639.7818	0.0000	1.0000
HKY	39	-3310.1686	6702.0604	0.0000	1.0000
GTR	43	-3309.7766	6710.0904	0.0000	1.0000
K2P	36	-3324.1004	6723.3684	0.0000	1.0000
SYM	40	-3322.9248	6729.7683	0.0000	1.0000
F81G	39	-3338.9191	6759.5613	0.0000	1.0000
F81IG	40	-3338.9262	6761.7711	0.0000	1.0000
F81I	39	-3341.1509	6764.0250	0.0000	1.0000
JC69G	36	-3350.9369	6777.0415	0.0000	1.0000
JC69IG	37	-3350.9455	6779.2387	0.0000	1.0000
JC69I	36	-3353.2985	6781.7647	0.0000	1.0000
F81	38	-3383.3150	6846.1627	0.0000	1.0000
JC69	35	-3395.2412	6863.4753	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3253.3108	6586.6217	0.6828	0.6828
HKYIG	41	-3253.3069	6588.6138	0.2522	0.9349
HKYI	40	-3256.1795	6592.3589	0.0388	0.9737
GTRG	44	-3252.8820	6593.7639	0.0192	0.9929
GTRIG	45	-3252.8778	6595.7556	0.0071	1.0000
K2PG	37	-3267.7179	6609.4358	0.0000	1.0000
K2PIG	38	-3267.7029	6611.4057	0.0000	1.0000
GTRI	44	-3262.6050	6613.2101	0.0000	1.0000
SYMIG	42	-3266.5412	6617.0825	0.0000	1.0000
SYMG	41	-3273.1548	6628.3097	0.0000	1.0000
K2PI	37	-3277.9021	6629.8041	0.0000	1.0000
SYMI	41	-3276.8311	6635.6622	0.0000	1.0000
HKY	39	-3310.1686	6698.3373	0.0000	1.0000
GTR	43	-3309.7766	6705.5532	0.0000	1.0000
K2P	36	-3324.1004	6720.2007	0.0000	1.0000
SYM	40	-3322.9248	6725.8496	0.0000	1.0000
F81G	39	-3338.9191	6755.8382	0.0000	1.0000
F81IG	40	-3338.9262	6757.8523	0.0000	1.0000
F81I	39	-3341.1509	6760.3018	0.0000	1.0000
JC69G	36	-3350.9369	6773.8738	0.0000	1.0000
JC69IG	37	-3350.9455	6775.8911	0.0000	1.0000
JC69I	36	-3353.2985	6778.5971	0.0000	1.0000
F81	38	-3383.3150	6842.6299	0.0000	1.0000
JC69	35	-3395.2412	6860.4824	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3253.3108	6777.7275	0.9045	0.9045
HKYI	40	-3256.1795	6783.4648	0.0514	0.9559
HKYIG	41	-3253.3069	6784.4974	0.0306	0.9865
K2PG	37	-3267.7179	6786.2087	0.0130	0.9996
K2PIG	38	-3267.7029	6792.9563	0.0004	1.0000
GTRG	44	-3252.8820	6803.9804	0.0000	1.0000
K2PI	37	-3277.9021	6806.5770	0.0000	1.0000
GTRIG	45	-3252.8778	6810.7497	0.0000	1.0000
SYMIG	42	-3266.5412	6817.7436	0.0000	1.0000
GTRI	44	-3262.6050	6823.4265	0.0000	1.0000
SYMG	41	-3273.1548	6824.1932	0.0000	1.0000
SYMI	41	-3276.8311	6831.5457	0.0000	1.0000
HKY	39	-3310.1686	6884.6655	0.0000	1.0000
K2P	36	-3324.1004	6892.1960	0.0000	1.0000
GTR	43	-3309.7766	6910.9920	0.0000	1.0000
SYM	40	-3322.9248	6916.9554	0.0000	1.0000
F81G	39	-3338.9191	6942.1664	0.0000	1.0000
JC69G	36	-3350.9369	6945.8691	0.0000	1.0000
F81I	39	-3341.1509	6946.6301	0.0000	1.0000
F81IG	40	-3338.9262	6948.9582	0.0000	1.0000
JC69I	36	-3353.2985	6950.5923	0.0000	1.0000
JC69IG	37	-3350.9455	6952.6640	0.0000	1.0000
F81	38	-3383.3150	7024.1805	0.0000	1.0000
JC69	35	-3395.2412	7027.7001	0.0000	1.0000



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

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

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


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

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


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

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


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

