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

Input data from file "591-tomodel-phylip/chr2_23651.phylip" (7 taxa, 407 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	-1368.4132	2779.0026	0.5239	0.5239
GTRIG	21	-1368.4132	2781.2263	0.1723	0.6962
GTRI	20	-1369.5732	2781.3226	0.1642	0.8604
HKYG	16	-1374.7298	2782.8544	0.0763	0.9367
HKYI	16	-1375.4422	2784.2793	0.0374	0.9742
HKYIG	17	-1374.7297	2785.0326	0.0257	0.9999
GTR	19	-1378.0046	2795.9730	0.0001	1.0000
HKY	15	-1384.7496	2800.7269	0.0000	1.0000
F81G	15	-1400.3514	2831.9305	0.0000	1.0000
F81I	15	-1400.4509	2832.1294	0.0000	1.0000
F81IG	16	-1400.3686	2834.1321	0.0000	1.0000
K2PG	13	-1404.4119	2835.7501	0.0000	1.0000
K2PI	13	-1404.9831	2836.8924	0.0000	1.0000
K2PIG	14	-1404.4200	2837.9115	0.0000	1.0000
SYMG	17	-1402.8031	2841.1794	0.0000	1.0000
SYMI	17	-1403.3973	2842.3678	0.0000	1.0000
SYMIG	18	-1402.8111	2843.3852	0.0000	1.0000
F81	14	-1409.1269	2847.3251	0.0000	1.0000
K2P	12	-1415.7024	2856.1967	0.0000	1.0000
SYM	16	-1413.7724	2860.9396	0.0000	1.0000
JC69G	12	-1425.4310	2875.6539	0.0000	1.0000
JC69I	12	-1425.5334	2875.8587	0.0000	1.0000
JC69IG	13	-1425.4465	2877.8192	0.0000	1.0000
JC69	11	-1435.0919	2892.8522	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1368.4132	2776.8264	0.5362	0.5362
GTRIG	21	-1368.4132	2778.8263	0.1973	0.7335
GTRI	20	-1369.5732	2779.1465	0.1681	0.9016
HKYG	16	-1374.7298	2781.4596	0.0529	0.9545
HKYI	16	-1375.4422	2782.8844	0.0259	0.9804
HKYIG	17	-1374.7297	2783.4594	0.0195	0.9999
GTR	19	-1378.0046	2794.0092	0.0001	1.0000
HKY	15	-1384.7496	2799.4993	0.0000	1.0000
F81G	15	-1400.3514	2830.7028	0.0000	1.0000
F81I	15	-1400.4509	2830.9017	0.0000	1.0000
F81IG	16	-1400.3686	2832.7372	0.0000	1.0000
K2PG	13	-1404.4119	2834.8239	0.0000	1.0000
K2PI	13	-1404.9831	2835.9662	0.0000	1.0000
K2PIG	14	-1404.4200	2836.8400	0.0000	1.0000
SYMG	17	-1402.8031	2839.6061	0.0000	1.0000
SYMI	17	-1403.3973	2840.7945	0.0000	1.0000
SYMIG	18	-1402.8111	2841.6223	0.0000	1.0000
F81	14	-1409.1269	2846.2537	0.0000	1.0000
K2P	12	-1415.7024	2855.4048	0.0000	1.0000
SYM	16	-1413.7724	2859.5448	0.0000	1.0000
JC69G	12	-1425.4310	2874.8620	0.0000	1.0000
JC69I	12	-1425.5334	2875.0669	0.0000	1.0000
JC69IG	13	-1425.4465	2876.8930	0.0000	1.0000
JC69	11	-1435.0919	2892.1838	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1374.7298	2845.6006	0.6470	0.6470
HKYI	16	-1375.4422	2847.0254	0.3173	0.9644
HKYIG	17	-1374.7297	2851.6092	0.0321	0.9965
GTRG	20	-1368.4132	2857.0027	0.0022	0.9986
GTRI	20	-1369.5732	2859.3227	0.0007	0.9993
HKY	15	-1384.7496	2859.6315	0.0006	0.9999
GTRIG	21	-1368.4132	2863.0114	0.0001	1.0000
GTR	19	-1378.0046	2870.1767	0.0000	1.0000
K2PG	13	-1404.4119	2886.9385	0.0000	1.0000
K2PI	13	-1404.9831	2888.0808	0.0000	1.0000
F81G	15	-1400.3514	2890.8350	0.0000	1.0000
F81I	15	-1400.4509	2891.0339	0.0000	1.0000
K2PIG	14	-1404.4200	2892.9634	0.0000	1.0000
F81IG	16	-1400.3686	2896.8782	0.0000	1.0000
F81	14	-1409.1269	2902.3771	0.0000	1.0000
K2P	12	-1415.7024	2903.5105	0.0000	1.0000
SYMG	17	-1402.8031	2907.7560	0.0000	1.0000
SYMI	17	-1403.3973	2908.9444	0.0000	1.0000
SYMIG	18	-1402.8111	2913.7809	0.0000	1.0000
JC69G	12	-1425.4310	2922.9678	0.0000	1.0000
JC69I	12	-1425.5334	2923.1726	0.0000	1.0000
SYM	16	-1413.7724	2923.6858	0.0000	1.0000
JC69IG	13	-1425.4465	2929.0076	0.0000	1.0000
JC69	11	-1435.0919	2936.2807	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/chr2_23651.phylip = 1 - 407;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23651.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/chr2_23651.phylip = 1 - 407;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23651.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/chr2_23651.phylip = 1 - 407;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23651.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

