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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2155.4024	4453.3503	0.6206	0.6206
GTRIG	65	-2155.0636	4455.1534	0.2519	0.8726
HKYG	60	-2162.2953	4457.2990	0.0862	0.9588
HKYIG	61	-2161.9241	4459.0030	0.0368	0.9955
GTRI	64	-2160.5121	4463.5697	0.0037	0.9993
HKYI	60	-2167.8232	4468.3548	0.0003	0.9996
K2PG	57	-2171.9116	4469.2429	0.0002	0.9998
K2PIG	58	-2171.6053	4471.0515	0.0001	0.9999
SYMG	61	-2168.5114	4472.1775	0.0001	1.0000
SYMIG	62	-2168.2033	4474.0163	0.0000	1.0000
K2PI	57	-2176.7492	4478.9181	0.0000	1.0000
SYMI	61	-2173.6637	4482.4822	0.0000	1.0000
GTR	63	-2186.2585	4512.5903	0.0000	1.0000
HKY	59	-2196.0707	4522.4117	0.0000	1.0000
F81G	59	-2198.6439	4527.5582	0.0000	1.0000
F81IG	60	-2198.3755	4529.4592	0.0000	1.0000
K2P	56	-2203.3924	4529.7916	0.0000	1.0000
SYM	60	-2199.2218	4531.1519	0.0000	1.0000
F81I	59	-2203.7513	4537.7730	0.0000	1.0000
JC69G	56	-2211.9552	4546.9173	0.0000	1.0000
JC69IG	57	-2211.6888	4548.7973	0.0000	1.0000
JC69I	56	-2217.0989	4557.2048	0.0000	1.0000
F81	58	-2229.5315	4586.9038	0.0000	1.0000
JC69	55	-2242.4602	4605.5227	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2155.4024	4438.8048	0.6220	0.6220
GTRIG	65	-2155.0636	4440.1272	0.3211	0.9431
HKYG	60	-2162.2953	4444.5906	0.0345	0.9776
HKYIG	61	-2161.9241	4445.8483	0.0184	0.9960
GTRI	64	-2160.5121	4449.0242	0.0038	0.9998
HKYI	60	-2167.8232	4455.6465	0.0001	0.9999
K2PG	57	-2171.9116	4457.8232	0.0000	0.9999
SYMG	61	-2168.5114	4459.0228	0.0000	1.0000
K2PIG	58	-2171.6053	4459.2107	0.0000	1.0000
SYMIG	62	-2168.2033	4460.4066	0.0000	1.0000
K2PI	57	-2176.7492	4467.4984	0.0000	1.0000
SYMI	61	-2173.6637	4469.3274	0.0000	1.0000
GTR	63	-2186.2585	4498.5170	0.0000	1.0000
HKY	59	-2196.0707	4510.1413	0.0000	1.0000
F81G	59	-2198.6439	4515.2879	0.0000	1.0000
F81IG	60	-2198.3755	4516.7509	0.0000	1.0000
SYM	60	-2199.2218	4518.4436	0.0000	1.0000
K2P	56	-2203.3924	4518.7847	0.0000	1.0000
F81I	59	-2203.7513	4525.5026	0.0000	1.0000
JC69G	56	-2211.9552	4535.9104	0.0000	1.0000
JC69IG	57	-2211.6888	4537.3776	0.0000	1.0000
JC69I	56	-2217.0989	4546.1979	0.0000	1.0000
F81	58	-2229.5315	4575.0629	0.0000	1.0000
JC69	55	-2242.4602	4594.9203	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2171.9116	4711.8591	0.4866	0.4866
HKYG	60	-2162.2953	4711.9968	0.4542	0.9409
K2PIG	58	-2171.6053	4717.7033	0.0262	0.9670
HKYIG	61	-2161.9241	4717.7112	0.0261	0.9931
K2PI	57	-2176.7492	4721.5343	0.0039	0.9970
HKYI	60	-2167.8232	4723.0527	0.0018	0.9988
GTRG	64	-2155.4024	4724.0381	0.0011	0.9999
GTRIG	65	-2155.0636	4729.8172	0.0001	1.0000
SYMG	61	-2168.5114	4730.8857	0.0000	1.0000
GTRI	64	-2160.5121	4734.2575	0.0000	1.0000
SYMIG	62	-2168.2033	4736.7263	0.0000	1.0000
SYMI	61	-2173.6637	4741.1904	0.0000	1.0000
K2P	56	-2203.3924	4768.3638	0.0000	1.0000
HKY	59	-2196.0707	4773.0907	0.0000	1.0000
F81G	59	-2198.6439	4778.2373	0.0000	1.0000
GTR	63	-2186.2585	4779.2935	0.0000	1.0000
F81IG	60	-2198.3755	4784.1571	0.0000	1.0000
JC69G	56	-2211.9552	4785.4895	0.0000	1.0000
SYM	60	-2199.2218	4785.8498	0.0000	1.0000
F81I	59	-2203.7513	4788.4520	0.0000	1.0000
JC69IG	57	-2211.6888	4791.4135	0.0000	1.0000
JC69I	56	-2217.0989	4795.7770	0.0000	1.0000
F81	58	-2229.5315	4833.5556	0.0000	1.0000
JC69	55	-2242.4602	4840.0427	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/chr5_10260.phylip = 1 - 637;
 Partition Dummy = 1:183-phylip/chr5_10260.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10260.phylip = 1 - 637;
 Partition Dummy = 1:183-phylip/chr5_10260.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 BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10260.phylip = 1 - 637;
 Partition Dummy = 1:183-phylip/chr5_10260.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

