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

Input data from file "591-tomodel-phylip/chr5_11206.phylip" (7 taxa, 829 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	-2626.6265	5294.2926	0.6126	0.6126
GTRIG	21	-2626.6265	5296.3980	0.2138	0.8264
GTRI	20	-2628.2484	5297.5364	0.1210	0.9473
HKYG	16	-2633.7035	5300.0769	0.0340	0.9813
HKYIG	17	-2633.7035	5302.1616	0.0120	0.9933
HKYI	16	-2635.3273	5303.3246	0.0067	1.0000
GTR	19	-2640.1818	5319.3031	0.0000	1.0000
HKY	15	-2647.3154	5325.2213	0.0000	1.0000
SYMG	17	-2672.1658	5379.0863	0.0000	1.0000
K2PG	13	-2677.1838	5380.8143	0.0000	1.0000
SYMIG	18	-2672.1658	5381.1760	0.0000	1.0000
SYMI	17	-2673.5445	5381.8435	0.0000	1.0000
K2PIG	14	-2677.1838	5382.8836	0.0000	1.0000
K2PI	13	-2678.4633	5383.3732	0.0000	1.0000
SYM	16	-2685.7442	5404.1584	0.0000	1.0000
K2P	12	-2690.4711	5405.3246	0.0000	1.0000
F81G	15	-2707.7004	5445.9913	0.0000	1.0000
F81I	15	-2708.5380	5447.6664	0.0000	1.0000
F81IG	16	-2707.7107	5448.0914	0.0000	1.0000
F81	14	-2718.3728	5465.2615	0.0000	1.0000
JC69G	12	-2743.3320	5511.0464	0.0000	1.0000
JC69I	12	-2744.0500	5512.4824	0.0000	1.0000
JC69IG	13	-2743.3456	5513.1379	0.0000	1.0000
JC69	11	-2754.0655	5530.4541	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2626.6265	5293.2530	0.6107	0.6107
GTRIG	21	-2626.6265	5295.2530	0.2247	0.8353
GTRI	20	-2628.2484	5296.4968	0.1206	0.9559
HKYG	16	-2633.7035	5299.4069	0.0282	0.9841
HKYIG	17	-2633.7035	5301.4069	0.0104	0.9944
HKYI	16	-2635.3273	5302.6547	0.0055	1.0000
GTR	19	-2640.1818	5318.3637	0.0000	1.0000
HKY	15	-2647.3154	5324.6309	0.0000	1.0000
SYMG	17	-2672.1658	5378.3317	0.0000	1.0000
SYMIG	18	-2672.1658	5380.3315	0.0000	1.0000
K2PG	13	-2677.1838	5380.3677	0.0000	1.0000
SYMI	17	-2673.5445	5381.0889	0.0000	1.0000
K2PIG	14	-2677.1838	5382.3676	0.0000	1.0000
K2PI	13	-2678.4633	5382.9266	0.0000	1.0000
SYM	16	-2685.7442	5403.4885	0.0000	1.0000
K2P	12	-2690.4711	5404.9423	0.0000	1.0000
F81G	15	-2707.7004	5445.4009	0.0000	1.0000
F81I	15	-2708.5380	5447.0760	0.0000	1.0000
F81IG	16	-2707.7107	5447.4214	0.0000	1.0000
F81	14	-2718.3728	5464.7456	0.0000	1.0000
JC69G	12	-2743.3320	5510.6641	0.0000	1.0000
JC69I	12	-2744.0500	5512.1000	0.0000	1.0000
JC69IG	13	-2743.3456	5512.6912	0.0000	1.0000
JC69	11	-2754.0655	5530.1309	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2633.7035	5374.9305	0.8104	0.8104
HKYI	16	-2635.3273	5378.1782	0.1598	0.9701
HKYIG	17	-2633.7035	5381.6507	0.0281	0.9983
GTRG	20	-2626.6265	5387.6574	0.0014	0.9996
GTRI	20	-2628.2484	5390.9012	0.0003	0.9999
GTRIG	21	-2626.6265	5394.3776	0.0000	1.0000
HKY	15	-2647.3154	5395.4342	0.0000	1.0000
GTR	19	-2640.1818	5408.0478	0.0000	1.0000
K2PG	13	-2677.1838	5441.7305	0.0000	1.0000
K2PI	13	-2678.4633	5444.2895	0.0000	1.0000
K2PIG	14	-2677.1838	5448.4507	0.0000	1.0000
SYMG	17	-2672.1658	5458.5754	0.0000	1.0000
SYMI	17	-2673.5445	5461.3327	0.0000	1.0000
K2P	12	-2690.4711	5461.5849	0.0000	1.0000
SYMIG	18	-2672.1658	5465.2955	0.0000	1.0000
SYM	16	-2685.7442	5479.0120	0.0000	1.0000
F81G	15	-2707.7004	5516.2042	0.0000	1.0000
F81I	15	-2708.5380	5517.8793	0.0000	1.0000
F81IG	16	-2707.7107	5522.9450	0.0000	1.0000
F81	14	-2718.3728	5530.8286	0.0000	1.0000
JC69G	12	-2743.3320	5567.3067	0.0000	1.0000
JC69I	12	-2744.0500	5568.7426	0.0000	1.0000
JC69IG	13	-2743.3456	5574.0541	0.0000	1.0000
JC69	11	-2754.0655	5582.0533	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/chr5_11206.phylip = 1 - 829;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11206.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/chr5_11206.phylip = 1 - 829;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11206.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/chr5_11206.phylip = 1 - 829;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11206.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

