
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:38 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/chr5_12037.phylip" (19 taxa, 234 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	-712.9252	1514.1974	0.4330	0.4330
K2PI	37	-712.9502	1514.2473	0.4223	0.8553
K2PIG	38	-712.9229	1517.0459	0.1042	0.9595
HKYI	40	-711.9131	1520.8211	0.0158	0.9753
HKYG	40	-711.9222	1520.8393	0.0156	0.9910
HKYIG	41	-711.8832	1523.7038	0.0037	0.9947
SYMI	41	-712.3600	1524.6575	0.0023	0.9970
SYMG	41	-712.3602	1524.6579	0.0023	0.9993
SYMIG	42	-712.3301	1527.5712	0.0005	0.9999
GTRI	44	-711.5382	1532.0288	0.0001	0.9999
GTRG	44	-711.5423	1532.0369	0.0001	1.0000
GTRIG	45	-711.5059	1535.0331	0.0000	1.0000
JC69G	36	-725.6452	1536.8132	0.0000	1.0000
JC69I	36	-726.2389	1538.0007	0.0000	1.0000
JC69IG	37	-725.6450	1539.6370	0.0000	1.0000
F81I	39	-724.6884	1543.4593	0.0000	1.0000
F81G	39	-724.6890	1543.4604	0.0000	1.0000
F81IG	40	-724.6884	1546.3717	0.0000	1.0000
K2P	36	-741.0999	1567.7226	0.0000	1.0000
SYM	40	-738.7257	1574.4462	0.0000	1.0000
HKY	39	-740.4301	1574.9427	0.0000	1.0000
GTR	43	-738.4260	1582.7679	0.0000	1.0000
JC69	35	-752.1609	1587.0490	0.0000	1.0000
F81	38	-751.5451	1594.2902	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-712.9252	1499.8504	0.3631	0.3631
K2PI	37	-712.9502	1499.9004	0.3541	0.7172
K2PIG	38	-712.9229	1501.8459	0.1339	0.8511
HKYI	40	-711.9131	1503.8263	0.0497	0.9009
HKYG	40	-711.9222	1503.8445	0.0493	0.9501
HKYIG	41	-711.8832	1505.7663	0.0189	0.9690
SYMI	41	-712.3600	1506.7200	0.0117	0.9807
SYMG	41	-712.3602	1506.7204	0.0117	0.9924
SYMIG	42	-712.3301	1508.6602	0.0044	0.9968
GTRI	44	-711.5382	1511.0764	0.0013	0.9982
GTRG	44	-711.5423	1511.0846	0.0013	0.9995
GTRIG	45	-711.5059	1513.0118	0.0005	1.0000
JC69G	36	-725.6452	1523.2903	0.0000	1.0000
JC69I	36	-726.2389	1524.4779	0.0000	1.0000
JC69IG	37	-725.6450	1525.2900	0.0000	1.0000
F81I	39	-724.6884	1527.3769	0.0000	1.0000
F81G	39	-724.6890	1527.3779	0.0000	1.0000
F81IG	40	-724.6884	1529.3769	0.0000	1.0000
K2P	36	-741.0999	1554.1998	0.0000	1.0000
SYM	40	-738.7257	1557.4513	0.0000	1.0000
HKY	39	-740.4301	1558.8602	0.0000	1.0000
GTR	43	-738.4260	1562.8521	0.0000	1.0000
JC69	35	-752.1609	1574.3217	0.0000	1.0000
F81	38	-751.5451	1579.0902	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-712.9252	1627.6973	0.4896	0.4896
K2PI	37	-712.9502	1627.7473	0.4775	0.9671
K2PIG	38	-712.9229	1633.1481	0.0321	0.9992
HKYI	40	-711.9131	1642.0391	0.0004	0.9995
HKYG	40	-711.9222	1642.0573	0.0004	0.9999
HKYIG	41	-711.8832	1647.4345	0.0000	0.9999
JC69G	36	-725.6452	1647.6819	0.0000	1.0000
SYMI	41	-712.3600	1648.3882	0.0000	1.0000
SYMG	41	-712.3602	1648.3886	0.0000	1.0000
JC69I	36	-726.2389	1648.8694	0.0000	1.0000
JC69IG	37	-725.6450	1653.1369	0.0000	1.0000
SYMIG	42	-712.3301	1653.7837	0.0000	1.0000
F81I	39	-724.6884	1662.1344	0.0000	1.0000
F81G	39	-724.6890	1662.1355	0.0000	1.0000
GTRI	44	-711.5382	1663.1106	0.0000	1.0000
GTRG	44	-711.5423	1663.1187	0.0000	1.0000
F81IG	40	-724.6884	1667.5897	0.0000	1.0000
GTRIG	45	-711.5059	1668.5013	0.0000	1.0000
K2P	36	-741.0999	1678.5913	0.0000	1.0000
HKY	39	-740.4301	1693.6177	0.0000	1.0000
JC69	35	-752.1609	1695.2580	0.0000	1.0000
SYM	40	-738.7257	1695.6642	0.0000	1.0000
F81	38	-751.5451	1710.3924	0.0000	1.0000
GTR	43	-738.4260	1711.4309	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/chr5_12037.phylip = 1 - 234;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12037.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/chr5_12037.phylip = 1 - 234;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12037.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/chr5_12037.phylip = 1 - 234;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12037.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

