
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:15 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_24880.phylip" (19 taxa, 272 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	-738.4194	1571.0379	0.4657	0.4657
HKYI	40	-738.5676	1571.3344	0.4015	0.8672
HKYIG	41	-738.3632	1573.7003	0.1230	0.9902
GTRG	44	-737.9453	1581.3356	0.0027	0.9929
GTRI	44	-738.0132	1581.4714	0.0025	0.9954
K2PG	37	-748.0627	1582.1425	0.0018	0.9972
K2PI	37	-748.2090	1582.4352	0.0016	0.9988
GTRIG	45	-737.9184	1584.1554	0.0007	0.9994
K2PIG	38	-747.9921	1584.7053	0.0005	0.9999
HKY	39	-749.1409	1589.7300	0.0000	1.0000
SYMG	41	-747.2420	1591.4579	0.0000	1.0000
SYMI	41	-747.3474	1591.6686	0.0000	1.0000
SYMIG	42	-747.1753	1594.1235	0.0000	1.0000
GTR	43	-748.8433	1600.2831	0.0000	1.0000
K2P	36	-758.7712	1600.8785	0.0000	1.0000
F81G	39	-754.7470	1600.9423	0.0000	1.0000
F81I	39	-755.8637	1603.1756	0.0000	1.0000
F81IG	40	-754.7282	1603.6555	0.0000	1.0000
JC69G	36	-763.4432	1610.2227	0.0000	1.0000
SYM	40	-758.1113	1610.4218	0.0000	1.0000
JC69I	36	-764.5378	1612.4117	0.0000	1.0000
JC69IG	37	-763.3876	1612.7923	0.0000	1.0000
F81	38	-765.5474	1619.8158	0.0000	1.0000
JC69	35	-774.1221	1628.9221	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-738.4194	1556.8388	0.4306	0.4306
HKYI	40	-738.5676	1557.1353	0.3712	0.8018
HKYIG	41	-738.3632	1558.7264	0.1676	0.9694
GTRG	44	-737.9453	1563.8906	0.0127	0.9820
GTRI	44	-738.0132	1564.0265	0.0118	0.9939
GTRIG	45	-737.9184	1565.8369	0.0048	0.9987
K2PG	37	-748.0627	1570.1254	0.0006	0.9992
K2PI	37	-748.2090	1570.4181	0.0005	0.9997
K2PIG	38	-747.9921	1571.9842	0.0002	0.9999
HKY	39	-749.1409	1576.2817	0.0000	0.9999
SYMG	41	-747.2420	1576.4840	0.0000	1.0000
SYMI	41	-747.3474	1576.6947	0.0000	1.0000
SYMIG	42	-747.1753	1578.3506	0.0000	1.0000
GTR	43	-748.8433	1583.6866	0.0000	1.0000
F81G	39	-754.7470	1587.4940	0.0000	1.0000
F81IG	40	-754.7282	1589.4564	0.0000	1.0000
K2P	36	-758.7712	1589.5423	0.0000	1.0000
F81I	39	-755.8637	1589.7274	0.0000	1.0000
SYM	40	-758.1113	1596.2227	0.0000	1.0000
JC69G	36	-763.4432	1598.8865	0.0000	1.0000
JC69IG	37	-763.3876	1600.7752	0.0000	1.0000
JC69I	36	-764.5378	1601.0755	0.0000	1.0000
F81	38	-765.5474	1607.0947	0.0000	1.0000
JC69	35	-774.1221	1618.2442	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-738.4194	1701.0708	0.4019	0.4019
HKYI	40	-738.5676	1701.3673	0.3465	0.7484
K2PG	37	-748.0627	1703.5401	0.1169	0.8654
K2PI	37	-748.2090	1703.8327	0.1010	0.9664
HKYIG	41	-738.3632	1706.5642	0.0258	0.9922
K2PIG	38	-747.9921	1709.0047	0.0076	0.9998
HKY	39	-749.1409	1716.9080	0.0001	0.9999
K2P	36	-758.7712	1719.3512	0.0000	1.0000
GTRG	44	-737.9453	1722.5459	0.0000	1.0000
GTRI	44	-738.0132	1722.6818	0.0000	1.0000
SYMG	41	-747.2420	1724.3219	0.0000	1.0000
SYMI	41	-747.3474	1724.5326	0.0000	1.0000
GTRIG	45	-737.9184	1728.0980	0.0000	1.0000
F81G	39	-754.7470	1728.1203	0.0000	1.0000
JC69G	36	-763.4432	1728.6954	0.0000	1.0000
SYMIG	42	-747.1753	1729.7942	0.0000	1.0000
F81I	39	-755.8637	1730.3536	0.0000	1.0000
JC69I	36	-764.5378	1730.8844	0.0000	1.0000
F81IG	40	-754.7282	1733.6885	0.0000	1.0000
JC69IG	37	-763.3876	1734.1899	0.0000	1.0000
GTR	43	-748.8433	1738.7361	0.0000	1.0000
SYM	40	-758.1113	1740.4548	0.0000	1.0000
F81	38	-765.5474	1744.1152	0.0000	1.0000
JC69	35	-774.1221	1744.4472	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_24880.phylip = 1 - 272;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24880.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/chr2_24880.phylip = 1 - 272;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24880.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/chr2_24880.phylip = 1 - 272;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24880.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

