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

Minimum AIC  model: HKYG
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	-739.5564	1565.4462	0.3313	0.3313
K2PI	37	-739.7476	1565.8286	0.2736	0.6049
HKYG	40	-736.5045	1567.5867	0.1136	0.7185
K2PIG	38	-739.3478	1567.7528	0.1045	0.8230
HKYI	40	-736.6525	1567.8828	0.0980	0.9210
HKYIG	41	-736.3706	1570.1163	0.0321	0.9531
SYMG	41	-737.0964	1571.5677	0.0155	0.9686
SYMI	41	-737.3056	1571.9861	0.0126	0.9812
GTRG	44	-733.7013	1573.3212	0.0065	0.9876
GTRI	44	-733.8826	1573.6838	0.0054	0.9930
SYMIG	42	-736.9074	1574.0120	0.0046	0.9976
GTRIG	45	-733.5698	1575.9579	0.0017	0.9993
JC69G	36	-748.0758	1579.7847	0.0003	0.9996
JC69I	36	-748.3053	1580.2438	0.0002	0.9998
JC69IG	37	-747.8227	1581.9787	0.0001	0.9999
F81G	39	-745.3328	1582.4710	0.0001	0.9999
F81I	39	-745.5364	1582.8781	0.0001	1.0000
F81IG	40	-745.1259	1584.8296	0.0000	1.0000
K2P	36	-755.4902	1594.6137	0.0000	1.0000
HKY	39	-752.3500	1596.5054	0.0000	1.0000
SYM	40	-752.5001	1599.5779	0.0000	1.0000
GTR	43	-748.9269	1600.8988	0.0000	1.0000
JC69	35	-763.4515	1607.8596	0.0000	1.0000
F81	38	-760.6092	1610.2757	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-736.5045	1553.0089	0.1789	0.1789
K2PG	37	-739.5564	1553.1129	0.1699	0.3488
HKYI	40	-736.6525	1553.3050	0.1543	0.5031
K2PI	37	-739.7476	1553.4953	0.1403	0.6434
K2PIG	38	-739.3478	1554.6956	0.0770	0.7204
HKYIG	41	-736.3706	1554.7413	0.0752	0.7956
GTRG	44	-733.7013	1555.4027	0.0541	0.8497
GTRI	44	-733.8826	1555.7652	0.0451	0.8948
SYMG	41	-737.0964	1556.1927	0.0364	0.9312
SYMI	41	-737.3056	1556.6111	0.0295	0.9608
GTRIG	45	-733.5698	1557.1397	0.0227	0.9834
SYMIG	42	-736.9074	1557.8147	0.0162	0.9996
JC69G	36	-748.0758	1568.1516	0.0001	0.9997
JC69I	36	-748.3053	1568.6107	0.0001	0.9998
F81G	39	-745.3328	1568.6657	0.0001	0.9999
F81I	39	-745.5364	1569.0728	0.0001	0.9999
JC69IG	37	-747.8227	1569.6454	0.0000	1.0000
F81IG	40	-745.1259	1570.2518	0.0000	1.0000
HKY	39	-752.3500	1582.7001	0.0000	1.0000
K2P	36	-755.4902	1582.9805	0.0000	1.0000
GTR	43	-748.9269	1583.8537	0.0000	1.0000
SYM	40	-752.5001	1585.0001	0.0000	1.0000
JC69	35	-763.4515	1596.9031	0.0000	1.0000
F81	38	-760.6092	1597.2184	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-739.5564	1685.7022	0.5216	0.5216
K2PI	37	-739.7476	1686.0846	0.4308	0.9523
K2PIG	38	-739.3478	1690.8684	0.0394	0.9917
HKYG	40	-736.5045	1696.3488	0.0025	0.9943
HKYI	40	-736.6525	1696.6448	0.0022	0.9965
JC69G	36	-748.0758	1697.1574	0.0017	0.9982
JC69I	36	-748.3053	1697.6165	0.0013	0.9995
HKYIG	41	-736.3706	1701.6646	0.0002	0.9997
JC69IG	37	-747.8227	1702.2347	0.0001	0.9998
SYMG	41	-737.0964	1703.1160	0.0001	0.9999
SYMI	41	-737.3056	1703.5345	0.0001	1.0000
SYMIG	42	-736.9074	1708.3216	0.0000	1.0000
F81G	39	-745.3328	1708.4220	0.0000	1.0000
F81I	39	-745.5364	1708.8292	0.0000	1.0000
K2P	36	-755.4902	1711.9864	0.0000	1.0000
GTRG	44	-733.7013	1713.0765	0.0000	1.0000
GTRI	44	-733.8826	1713.4391	0.0000	1.0000
F81IG	40	-745.1259	1713.5916	0.0000	1.0000
GTRIG	45	-733.5698	1718.3970	0.0000	1.0000
JC69	35	-763.4515	1722.3255	0.0000	1.0000
HKY	39	-752.3500	1722.4564	0.0000	1.0000
SYM	40	-752.5001	1728.3400	0.0000	1.0000
F81	38	-760.6092	1733.3913	0.0000	1.0000
GTR	43	-748.9269	1737.9441	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/chr2_29129.phylip = 1 - 266;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_29129.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 (K2PG)]

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

