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

Minimum AIC  model: F81G
Minimum AICc model: F81G
Minimum BIC  model: F81G


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
F81G	39	-729.5323	1548.0894	0.3523	0.3523
F81I	39	-729.8300	1548.6848	0.2616	0.6140
HKYG	40	-729.3722	1550.3756	0.1123	0.7263
F81IG	40	-729.4739	1550.5789	0.1015	0.8278
HKYI	40	-729.6723	1550.9759	0.0832	0.9110
F81	38	-732.9786	1552.3939	0.0410	0.9520
HKYIG	41	-729.3112	1552.8787	0.0321	0.9841
HKY	39	-732.8427	1554.7102	0.0129	0.9970
GTRG	44	-728.4171	1559.0788	0.0014	0.9984
GTRI	44	-728.7625	1559.7696	0.0010	0.9994
GTRIG	45	-728.3448	1561.6355	0.0004	0.9998
GTR	43	-731.9044	1563.3715	0.0002	1.0000
JC69G	36	-753.5160	1588.3466	0.0000	1.0000
JC69I	36	-753.8156	1588.9458	0.0000	1.0000
JC69IG	37	-753.4666	1590.7998	0.0000	1.0000
K2PG	37	-753.5033	1590.8732	0.0000	1.0000
K2PI	37	-753.8036	1591.4739	0.0000	1.0000
JC69	35	-756.8365	1592.4535	0.0000	1.0000
K2PIG	38	-753.4532	1593.3430	0.0000	1.0000
K2P	36	-756.8274	1594.9694	0.0000	1.0000
SYMG	41	-750.8696	1595.9953	0.0000	1.0000
SYMI	41	-751.1942	1596.6446	0.0000	1.0000
SYMIG	42	-750.8086	1598.5173	0.0000	1.0000
SYM	40	-754.3058	1600.2428	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
F81G	39	-729.5323	1537.0647	0.3105	0.3105
F81I	39	-729.8300	1537.6601	0.2305	0.5410
HKYG	40	-729.3722	1538.7444	0.1340	0.6750
F81IG	40	-729.4739	1538.9477	0.1211	0.7961
HKYI	40	-729.6723	1539.3447	0.0993	0.8954
HKYIG	41	-729.3112	1540.6225	0.0524	0.9478
F81	38	-732.9786	1541.9572	0.0269	0.9747
HKY	39	-732.8427	1543.6854	0.0113	0.9860
GTRG	44	-728.4171	1544.8342	0.0064	0.9924
GTRI	44	-728.7625	1545.5250	0.0045	0.9969
GTRIG	45	-728.3448	1546.6896	0.0025	0.9995
GTR	43	-731.9044	1549.8088	0.0005	1.0000
JC69G	36	-753.5160	1579.0320	0.0000	1.0000
JC69I	36	-753.8156	1579.6311	0.0000	1.0000
JC69IG	37	-753.4666	1580.9331	0.0000	1.0000
K2PG	37	-753.5033	1581.0065	0.0000	1.0000
K2PI	37	-753.8036	1581.6073	0.0000	1.0000
K2PIG	38	-753.4532	1582.9064	0.0000	1.0000
JC69	35	-756.8365	1583.6731	0.0000	1.0000
SYMG	41	-750.8696	1583.7391	0.0000	1.0000
SYMI	41	-751.1942	1584.3883	0.0000	1.0000
SYMIG	42	-750.8086	1585.6173	0.0000	1.0000
K2P	36	-756.8274	1585.6547	0.0000	1.0000
SYM	40	-754.3058	1588.6116	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81G	39	-729.5323	1684.3931	0.3955	0.3955
F81I	39	-729.8300	1684.9885	0.2937	0.6892
F81	38	-732.9786	1685.5080	0.2265	0.9157
HKYG	40	-729.3722	1689.8505	0.0258	0.9416
F81IG	40	-729.4739	1690.0538	0.0233	0.9649
HKYI	40	-729.6723	1690.4508	0.0191	0.9840
HKY	39	-732.8427	1691.0139	0.0144	0.9985
HKYIG	41	-729.3112	1695.5062	0.0015	1.0000
GTRG	44	-728.4171	1711.0509	0.0000	1.0000
GTRI	44	-728.7625	1711.7417	0.0000	1.0000
GTR	43	-731.9044	1712.2478	0.0000	1.0000
JC69G	36	-753.5160	1715.0274	0.0000	1.0000
JC69I	36	-753.8156	1715.6266	0.0000	1.0000
JC69	35	-756.8365	1715.8909	0.0000	1.0000
GTRIG	45	-728.3448	1716.6840	0.0000	1.0000
JC69IG	37	-753.4666	1720.7063	0.0000	1.0000
K2PG	37	-753.5033	1720.7797	0.0000	1.0000
K2PI	37	-753.8036	1721.3804	0.0000	1.0000
K2P	36	-756.8274	1721.6502	0.0000	1.0000
K2PIG	38	-753.4532	1726.4572	0.0000	1.0000
SYMG	41	-750.8696	1738.6228	0.0000	1.0000
SYMI	41	-751.1942	1739.2721	0.0000	1.0000
SYM	40	-754.3058	1739.7177	0.0000	1.0000
SYMIG	42	-750.8086	1744.2787	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (F81G)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2007.phylip = 1 - 323;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2007.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (F81G)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2007.phylip = 1 - 323;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2007.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (F81G)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2007.phylip = 1 - 323;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2007.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

