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

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


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	-684.0872	1463.5016	0.4042	0.4042
HKYI	40	-684.6241	1464.5753	0.2363	0.6405
HKY	39	-686.3356	1465.1829	0.1744	0.8149
HKYIG	41	-684.0867	1466.3425	0.0977	0.9125
K2PG	37	-690.9849	1468.9283	0.0268	0.9393
K2PI	37	-691.4021	1469.7628	0.0177	0.9570
GTRG	44	-682.1102	1471.0776	0.0092	0.9662
K2P	36	-693.5655	1471.3512	0.0080	0.9741
K2PIG	38	-690.9882	1471.6986	0.0067	0.9808
GTR	43	-683.9950	1471.9236	0.0060	0.9868
GTRI	44	-682.6215	1472.1001	0.0055	0.9923
SYMG	41	-687.7991	1473.7673	0.0024	0.9947
GTRIG	45	-682.1099	1474.0285	0.0021	0.9968
SYMI	41	-688.2199	1474.6088	0.0016	0.9984
SYM	40	-690.0266	1475.3803	0.0011	0.9994
SYMIG	42	-687.7987	1476.6352	0.0006	1.0000
F81G	39	-701.8514	1496.2144	0.0000	1.0000
F81I	39	-702.2914	1497.0944	0.0000	1.0000
F81	38	-704.0685	1497.8592	0.0000	1.0000
F81IG	40	-701.8510	1499.0291	0.0000	1.0000
JC69G	36	-708.4118	1501.0437	0.0000	1.0000
JC69I	36	-708.8118	1501.8437	0.0000	1.0000
JC69	35	-710.6045	1502.7159	0.0000	1.0000
JC69IG	37	-708.4113	1503.7812	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-684.0872	1448.1745	0.3816	0.3816
HKYI	40	-684.6241	1449.2482	0.2231	0.6047
HKYIG	41	-684.0867	1450.1735	0.1405	0.7451
HKY	39	-686.3356	1450.6713	0.1095	0.8546
GTRG	44	-682.1102	1452.2204	0.0505	0.9051
GTRI	44	-682.6215	1453.2430	0.0303	0.9354
GTR	43	-683.9950	1453.9899	0.0208	0.9562
GTRIG	45	-682.1099	1454.2199	0.0186	0.9748
K2PG	37	-690.9849	1455.9698	0.0077	0.9825
K2PI	37	-691.4021	1456.8042	0.0051	0.9876
SYMG	41	-687.7991	1457.5982	0.0034	0.9910
K2PIG	38	-690.9882	1457.9764	0.0028	0.9939
SYMI	41	-688.2199	1458.4398	0.0023	0.9961
K2P	36	-693.5655	1459.1310	0.0016	0.9977
SYMIG	42	-687.7987	1459.5975	0.0013	0.9990
SYM	40	-690.0266	1460.0532	0.0010	1.0000
F81G	39	-701.8514	1481.7028	0.0000	1.0000
F81I	39	-702.2914	1482.5828	0.0000	1.0000
F81IG	40	-701.8510	1483.7020	0.0000	1.0000
F81	38	-704.0685	1484.1370	0.0000	1.0000
JC69G	36	-708.4118	1488.8236	0.0000	1.0000
JC69I	36	-708.8118	1489.6235	0.0000	1.0000
JC69IG	37	-708.4113	1490.8227	0.0000	1.0000
JC69	35	-710.6045	1491.2091	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	36	-693.5655	1586.6165	0.3231	0.3231
K2PG	37	-690.9849	1586.9965	0.2672	0.5904
K2PI	37	-691.4021	1587.8310	0.1761	0.7664
HKY	39	-686.3356	1588.7805	0.1095	0.8759
HKYG	40	-684.0872	1589.8250	0.0650	0.9409
HKYI	40	-684.6241	1590.8987	0.0380	0.9789
K2PIG	38	-690.9882	1592.5444	0.0167	0.9956
HKYIG	41	-684.0867	1595.3653	0.0041	0.9996
SYM	40	-690.0266	1601.7038	0.0002	0.9998
SYMG	41	-687.7991	1602.7900	0.0001	0.9999
SYMI	41	-688.2199	1603.6316	0.0001	1.0000
GTR	43	-683.9950	1606.2643	0.0000	1.0000
GTRG	44	-682.1102	1608.0360	0.0000	1.0000
SYMIG	42	-687.7987	1608.3305	0.0000	1.0000
GTRI	44	-682.6215	1609.0586	0.0000	1.0000
GTRIG	45	-682.1099	1613.5767	0.0000	1.0000
JC69	35	-710.6045	1615.1533	0.0000	1.0000
JC69G	36	-708.4118	1616.3090	0.0000	1.0000
JC69I	36	-708.8118	1617.1090	0.0000	1.0000
F81	38	-704.0685	1618.7050	0.0000	1.0000
F81G	39	-701.8514	1619.8121	0.0000	1.0000
F81I	39	-702.2914	1620.6921	0.0000	1.0000
JC69IG	37	-708.4113	1621.8495	0.0000	1.0000
F81IG	40	-701.8510	1625.3525	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/chr13_1159.phylip = 1 - 255;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_1159.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/chr13_1159.phylip = 1 - 255;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_1159.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 (K2P)]

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


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

