
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:43 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/chr7_10536.phylip" (19 taxa, 311 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	-779.3627	1650.8736	0.4473	0.4473
HKYI	40	-779.6280	1651.4041	0.3431	0.7904
HKYIG	41	-779.1084	1653.0198	0.1530	0.9434
GTRG	44	-776.7985	1656.4843	0.0271	0.9704
GTRI	44	-777.2053	1657.2978	0.0180	0.9884
GTRIG	45	-776.3799	1658.3825	0.0105	0.9989
HKY	39	-787.6955	1664.9040	0.0004	0.9993
K2PG	37	-790.9020	1666.1044	0.0002	0.9995
K2PI	37	-791.1764	1666.6532	0.0002	0.9997
F81G	39	-789.1937	1667.9004	0.0001	0.9998
K2PIG	38	-790.7599	1668.4169	0.0001	0.9999
F81I	39	-789.5175	1668.5478	0.0001	0.9999
GTR	43	-784.7781	1669.7286	0.0000	1.0000
F81IG	40	-789.0907	1670.3296	0.0000	1.0000
SYMG	41	-789.1355	1673.0740	0.0000	1.0000
SYMI	41	-789.4981	1673.7991	0.0000	1.0000
SYMIG	42	-788.9911	1675.4599	0.0000	1.0000
K2P	36	-798.4282	1678.5790	0.0000	1.0000
F81	38	-795.9469	1678.7909	0.0000	1.0000
JC69G	36	-800.0617	1681.8461	0.0000	1.0000
JC69I	36	-800.3669	1682.4564	0.0000	1.0000
SYM	40	-795.8995	1683.9472	0.0000	1.0000
JC69IG	37	-799.9889	1684.2782	0.0000	1.0000
JC69	35	-806.5960	1692.3556	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-779.3627	1638.7254	0.3607	0.3607
HKYI	40	-779.6280	1639.2560	0.2767	0.6374
HKYIG	41	-779.1084	1640.2169	0.1711	0.8085
GTRG	44	-776.7985	1641.5971	0.0858	0.8943
GTRI	44	-777.2053	1642.4105	0.0571	0.9514
GTRIG	45	-776.3799	1642.7598	0.0480	0.9994
HKY	39	-787.6955	1653.3911	0.0002	0.9996
GTR	43	-784.7781	1655.5563	0.0001	0.9997
K2PG	37	-790.9020	1655.8041	0.0001	0.9998
K2PI	37	-791.1764	1656.3528	0.0001	0.9998
F81G	39	-789.1937	1656.3875	0.0001	0.9999
F81I	39	-789.5175	1657.0349	0.0000	0.9999
K2PIG	38	-790.7599	1657.5199	0.0000	1.0000
F81IG	40	-789.0907	1658.1814	0.0000	1.0000
SYMG	41	-789.1355	1660.2710	0.0000	1.0000
SYMI	41	-789.4981	1660.9961	0.0000	1.0000
SYMIG	42	-788.9911	1661.9823	0.0000	1.0000
F81	38	-795.9469	1667.8939	0.0000	1.0000
K2P	36	-798.4282	1668.8564	0.0000	1.0000
SYM	40	-795.8995	1671.7991	0.0000	1.0000
JC69G	36	-800.0617	1672.1235	0.0000	1.0000
JC69I	36	-800.3669	1672.7338	0.0000	1.0000
JC69IG	37	-799.9889	1673.9779	0.0000	1.0000
JC69	35	-806.5960	1683.1920	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-779.3627	1788.3172	0.5143	0.5143
HKYI	40	-779.6280	1788.8477	0.3944	0.9087
HKYIG	41	-779.1084	1793.5484	0.0376	0.9463
K2PG	37	-790.9020	1794.1764	0.0275	0.9738
K2PI	37	-791.1764	1794.7252	0.0209	0.9947
HKY	39	-787.6955	1799.2430	0.0022	0.9968
K2PIG	38	-790.7599	1799.6320	0.0018	0.9986
F81G	39	-789.1937	1802.2394	0.0005	0.9991
F81I	39	-789.5175	1802.8868	0.0004	0.9995
K2P	36	-798.4282	1803.4889	0.0003	0.9997
GTRG	44	-776.7985	1806.1479	0.0001	0.9998
JC69G	36	-800.0617	1806.7560	0.0001	0.9999
GTRI	44	-777.2053	1806.9614	0.0000	0.9999
JC69I	36	-800.3669	1807.3664	0.0000	0.9999
F81IG	40	-789.0907	1807.7732	0.0000	1.0000
F81	38	-795.9469	1810.0060	0.0000	1.0000
GTRIG	45	-776.3799	1811.0505	0.0000	1.0000
JC69IG	37	-799.9889	1812.3502	0.0000	1.0000
SYMG	41	-789.1355	1813.6025	0.0000	1.0000
JC69	35	-806.5960	1814.0848	0.0000	1.0000
SYMI	41	-789.4981	1814.3276	0.0000	1.0000
GTR	43	-784.7781	1816.3674	0.0000	1.0000
SYMIG	42	-788.9911	1819.0536	0.0000	1.0000
SYM	40	-795.8995	1821.3908	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/chr7_10536.phylip = 1 - 311;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10536.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/chr7_10536.phylip = 1 - 311;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10536.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/chr7_10536.phylip = 1 - 311;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10536.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

