
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:48 02/17/2011

Input data from file "183-phylip/chr4_17854.phylip" (29 taxa, 217 characters)

Minimum AIC  model: GTRG
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	57	-729.5974	1614.7796	0.4522	0.4522
K2PG	57	-729.7415	1615.0679	0.3915	0.8437
K2PIG	58	-729.7411	1618.7988	0.0606	0.9043
HKYG	60	-726.2532	1619.4294	0.0442	0.9486
HKYI	60	-726.4033	1619.7296	0.0381	0.9866
HKYIG	61	-726.0383	1622.8766	0.0079	0.9945
GTRG	64	-721.6278	1625.9924	0.0017	0.9962
SYMI	61	-727.8819	1626.5639	0.0012	0.9974
SYMG	61	-728.0228	1626.8455	0.0011	0.9985
GTRI	64	-722.2726	1627.2820	0.0009	0.9994
GTRIG	65	-721.0153	1628.8518	0.0004	0.9998
SYMIG	62	-727.9179	1630.5631	0.0002	0.9999
JC69I	56	-741.1229	1634.1459	0.0000	1.0000
JC69G	56	-741.2833	1634.4666	0.0000	1.0000
JC69IG	57	-741.2829	1638.1507	0.0000	1.0000
F81I	59	-737.6451	1638.3858	0.0000	1.0000
F81G	59	-737.7832	1638.6619	0.0000	1.0000
F81IG	60	-737.7829	1642.4889	0.0000	1.0000
K2P	56	-751.0792	1654.0585	0.0000	1.0000
HKY	59	-748.2385	1659.5725	0.0000	1.0000
GTR	63	-742.6472	1664.0002	0.0000	1.0000
SYM	60	-749.3320	1665.5870	0.0000	1.0000
JC69	55	-762.3988	1673.0584	0.0000	1.0000
F81	58	-758.3893	1676.0951	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-721.6278	1571.2555	0.2295	0.2295
GTRIG	65	-721.0153	1572.0306	0.1558	0.3853
HKYG	60	-726.2532	1572.5063	0.1228	0.5081
GTRI	64	-722.2726	1572.5452	0.1204	0.6285
HKYI	60	-726.4033	1572.8065	0.1057	0.7342
K2PI	57	-729.5974	1573.1947	0.0870	0.8212
K2PG	57	-729.7415	1573.4830	0.0754	0.8966
HKYIG	61	-726.0383	1574.0766	0.0560	0.9526
K2PIG	58	-729.7411	1575.4823	0.0277	0.9803
SYMI	61	-727.8819	1577.7639	0.0089	0.9891
SYMG	61	-728.0228	1578.0455	0.0077	0.9968
SYMIG	62	-727.9179	1579.8358	0.0031	1.0000
F81I	59	-737.6451	1593.2902	0.0000	1.0000
F81G	59	-737.7832	1593.5663	0.0000	1.0000
JC69I	56	-741.1229	1594.2459	0.0000	1.0000
JC69G	56	-741.2833	1594.5666	0.0000	1.0000
F81IG	60	-737.7829	1595.5658	0.0000	1.0000
JC69IG	57	-741.2829	1596.5658	0.0000	1.0000
GTR	63	-742.6472	1611.2943	0.0000	1.0000
K2P	56	-751.0792	1614.1585	0.0000	1.0000
HKY	59	-748.2385	1614.4769	0.0000	1.0000
SYM	60	-749.3320	1618.6640	0.0000	1.0000
F81	58	-758.3893	1632.7786	0.0000	1.0000
JC69	55	-762.3988	1634.7976	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-729.5974	1765.8489	0.5148	0.5148
K2PG	57	-729.7415	1766.1372	0.4457	0.9606
K2PIG	58	-729.7411	1771.5163	0.0303	0.9908
HKYG	60	-726.2532	1775.3002	0.0046	0.9954
HKYI	60	-726.4033	1775.6004	0.0039	0.9993
HKYIG	61	-726.0383	1780.2503	0.0004	0.9997
JC69I	56	-741.1229	1783.5201	0.0001	0.9998
JC69G	56	-741.2833	1783.8409	0.0001	0.9999
SYMI	61	-727.8819	1783.9376	0.0001	0.9999
SYMG	61	-728.0228	1784.2193	0.0001	1.0000
GTRG	64	-721.6278	1787.5690	0.0000	1.0000
GTRI	64	-722.2726	1788.8586	0.0000	1.0000
JC69IG	57	-741.2829	1789.2200	0.0000	1.0000
SYMIG	62	-727.9179	1789.3894	0.0000	1.0000
GTRIG	65	-721.0153	1791.7240	0.0000	1.0000
F81I	59	-737.6451	1792.7042	0.0000	1.0000
F81G	59	-737.7832	1792.9803	0.0000	1.0000
F81IG	60	-737.7829	1798.3596	0.0000	1.0000
K2P	56	-751.0792	1803.4327	0.0000	1.0000
HKY	59	-748.2385	1813.8909	0.0000	1.0000
JC69	55	-762.3988	1820.6919	0.0000	1.0000
SYM	60	-749.3320	1821.4578	0.0000	1.0000
GTR	63	-742.6472	1824.2279	0.0000	1.0000
F81	58	-758.3893	1828.8127	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_17854.phylip = 1 - 217;
 Partition Dummy = 1:183-phylip/chr4_17854.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_17854.phylip = 1 - 217;
 Partition Dummy = 1:183-phylip/chr4_17854.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_17854.phylip = 1 - 217;
 Partition Dummy = 1:183-phylip/chr4_17854.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

