
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:25 04/04/2011

Input data from file "591-tomodel-phylip/chr10_4599.phylip" (7 taxa, 518 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1754.3748	3550.4397	0.6179	0.6179
GTRIG	21	-1754.2918	3552.4465	0.2266	0.8445
GTRI	20	-1756.5611	3554.8123	0.0694	0.9139
HKYG	16	-1761.1675	3555.4209	0.0512	0.9651
HKYI	16	-1762.1928	3557.4714	0.0184	0.9835
HKYIG	17	-1761.2301	3557.6843	0.0165	1.0000
K2PG	13	-1777.7791	3582.2805	0.0000	1.0000
K2PI	13	-1778.7013	3584.1249	0.0000	1.0000
K2PIG	14	-1777.7854	3584.4058	0.0000	1.0000
SYMG	17	-1776.0447	3587.3134	0.0000	1.0000
SYMIG	18	-1776.0529	3589.4765	0.0000	1.0000
SYMI	17	-1777.1707	3589.5654	0.0000	1.0000
GTR	19	-1786.7421	3613.0103	0.0000	1.0000
HKY	15	-1796.8189	3624.5940	0.0000	1.0000
F81G	15	-1798.2926	3627.5414	0.0000	1.0000
F81I	15	-1798.9522	3628.8606	0.0000	1.0000
F81IG	16	-1798.3339	3629.7536	0.0000	1.0000
JC69G	12	-1812.8150	3650.2478	0.0000	1.0000
JC69I	12	-1813.2620	3651.1418	0.0000	1.0000
JC69IG	13	-1812.8918	3652.5059	0.0000	1.0000
K2P	12	-1815.4361	3655.4901	0.0000	1.0000
SYM	16	-1811.8749	3656.8357	0.0000	1.0000
F81	14	-1831.1579	3691.1508	0.0000	1.0000
JC69	11	-1847.0226	3716.5669	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1754.3748	3548.7496	0.6186	0.6186
GTRIG	21	-1754.2918	3550.5836	0.2473	0.8659
GTRI	20	-1756.5611	3553.1222	0.0695	0.9354
HKYG	16	-1761.1675	3554.3351	0.0379	0.9733
HKYI	16	-1762.1928	3556.3855	0.0136	0.9869
HKYIG	17	-1761.2301	3556.4603	0.0131	1.0000
K2PG	13	-1777.7791	3581.5582	0.0000	1.0000
K2PI	13	-1778.7013	3583.4027	0.0000	1.0000
K2PIG	14	-1777.7854	3583.5708	0.0000	1.0000
SYMG	17	-1776.0447	3586.0894	0.0000	1.0000
SYMIG	18	-1776.0529	3588.1058	0.0000	1.0000
SYMI	17	-1777.1707	3588.3414	0.0000	1.0000
GTR	19	-1786.7421	3611.4842	0.0000	1.0000
HKY	15	-1796.8189	3623.6378	0.0000	1.0000
F81G	15	-1798.2926	3626.5852	0.0000	1.0000
F81I	15	-1798.9522	3627.9045	0.0000	1.0000
F81IG	16	-1798.3339	3628.6678	0.0000	1.0000
JC69G	12	-1812.8150	3649.6300	0.0000	1.0000
JC69I	12	-1813.2620	3650.5239	0.0000	1.0000
JC69IG	13	-1812.8918	3651.7837	0.0000	1.0000
K2P	12	-1815.4361	3654.8722	0.0000	1.0000
SYM	16	-1811.8749	3655.7499	0.0000	1.0000
F81	14	-1831.1579	3690.3158	0.0000	1.0000
JC69	11	-1847.0226	3716.0451	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1761.1675	3622.3347	0.7118	0.7118
HKYI	16	-1762.1928	3624.3851	0.2553	0.9671
HKYIG	17	-1761.2301	3628.7099	0.0294	0.9965
GTRG	20	-1754.3748	3633.7491	0.0024	0.9989
K2PG	13	-1777.7791	3636.8079	0.0005	0.9994
GTRI	20	-1756.5611	3638.1217	0.0003	0.9997
K2PI	13	-1778.7013	3638.6524	0.0002	0.9999
GTRIG	21	-1754.2918	3639.8331	0.0001	1.0000
K2PIG	14	-1777.7854	3643.0704	0.0000	1.0000
SYMG	17	-1776.0447	3658.3390	0.0000	1.0000
SYMI	17	-1777.1707	3660.5910	0.0000	1.0000
SYMIG	18	-1776.0529	3664.6053	0.0000	1.0000
HKY	15	-1796.8189	3687.3874	0.0000	1.0000
F81G	15	-1798.2926	3690.3348	0.0000	1.0000
F81I	15	-1798.9522	3691.6541	0.0000	1.0000
GTR	19	-1786.7421	3692.2337	0.0000	1.0000
F81IG	16	-1798.3339	3696.6674	0.0000	1.0000
JC69G	12	-1812.8150	3700.6297	0.0000	1.0000
JC69I	12	-1813.2620	3701.5236	0.0000	1.0000
K2P	12	-1815.4361	3705.8719	0.0000	1.0000
JC69IG	13	-1812.8918	3707.0333	0.0000	1.0000
SYM	16	-1811.8749	3723.7495	0.0000	1.0000
F81	14	-1831.1579	3749.8155	0.0000	1.0000
JC69	11	-1847.0226	3762.7948	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 591-tomodel-phylip/chr10_4599.phylip = 1 - 518;
 Partition Dummy = 1:591-tomodel-phylip/chr10_4599.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 591-tomodel-phylip/chr10_4599.phylip = 1 - 518;
 Partition Dummy = 1:591-tomodel-phylip/chr10_4599.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr10_4599.phylip = 1 - 518;
 Partition Dummy = 1:591-tomodel-phylip/chr10_4599.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

