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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-1262.6695	2566.6895	0.3613	0.3613
GTRG	20	-1262.7211	2566.7926	0.3431	0.7044
GTRIG	21	-1262.7454	2568.9788	0.1150	0.8194
HKYI	16	-1268.5945	2570.0580	0.0670	0.8865
HKYG	16	-1268.6514	2570.1718	0.0633	0.9498
GTR	19	-1266.4960	2572.2119	0.0228	0.9726
HKYIG	17	-1268.6606	2572.3004	0.0219	0.9945
HKY	15	-1272.1442	2575.0539	0.0055	1.0000
F81I	15	-1287.8415	2606.4486	0.0000	1.0000
F81G	15	-1287.9344	2606.6343	0.0000	1.0000
F81IG	16	-1287.9343	2608.7376	0.0000	1.0000
F81	14	-1291.5852	2611.8391	0.0000	1.0000
K2PI	13	-1332.9531	2692.4850	0.0000	1.0000
K2PG	13	-1333.0333	2692.6453	0.0000	1.0000
SYMI	17	-1329.0869	2693.1530	0.0000	1.0000
SYMG	17	-1329.1579	2693.2950	0.0000	1.0000
K2PIG	14	-1333.0534	2694.7755	0.0000	1.0000
SYMIG	18	-1329.1745	2695.4453	0.0000	1.0000
K2P	12	-1335.8553	2696.2059	0.0000	1.0000
SYM	16	-1331.9243	2696.7176	0.0000	1.0000
JC69I	12	-1347.5048	2719.5048	0.0000	1.0000
JC69G	12	-1347.6030	2719.7012	0.0000	1.0000
JC69IG	13	-1347.6263	2721.8313	0.0000	1.0000
JC69	11	-1350.3745	2723.1674	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1262.6695	2565.3391	0.3710	0.3710
GTRG	20	-1262.7211	2565.4421	0.3524	0.7234
GTRIG	21	-1262.7454	2567.4909	0.1265	0.8499
HKYI	16	-1268.5945	2569.1890	0.0541	0.9040
HKYG	16	-1268.6514	2569.3028	0.0511	0.9552
GTR	19	-1266.4960	2570.9920	0.0220	0.9771
HKYIG	17	-1268.6606	2571.3212	0.0186	0.9958
HKY	15	-1272.1442	2574.2883	0.0042	1.0000
F81I	15	-1287.8415	2605.6830	0.0000	1.0000
F81G	15	-1287.9344	2605.8687	0.0000	1.0000
F81IG	16	-1287.9343	2607.8686	0.0000	1.0000
F81	14	-1291.5852	2611.1703	0.0000	1.0000
K2PI	13	-1332.9531	2691.9063	0.0000	1.0000
K2PG	13	-1333.0333	2692.0666	0.0000	1.0000
SYMI	17	-1329.0869	2692.1738	0.0000	1.0000
SYMG	17	-1329.1579	2692.3158	0.0000	1.0000
K2PIG	14	-1333.0534	2694.1068	0.0000	1.0000
SYMIG	18	-1329.1745	2694.3491	0.0000	1.0000
K2P	12	-1335.8553	2695.7107	0.0000	1.0000
SYM	16	-1331.9243	2695.8486	0.0000	1.0000
JC69I	12	-1347.5048	2719.0095	0.0000	1.0000
JC69G	12	-1347.6030	2719.2059	0.0000	1.0000
JC69IG	13	-1347.6263	2721.2526	0.0000	1.0000
JC69	11	-1350.3745	2722.7490	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1268.5945	2640.6473	0.3687	0.3687
HKYG	16	-1268.6514	2640.7611	0.3483	0.7169
HKY	15	-1272.1442	2641.2805	0.2686	0.9855
HKYIG	17	-1268.6606	2647.2456	0.0136	0.9992
GTRI	20	-1262.6695	2654.6620	0.0003	0.9995
GTRG	20	-1262.7211	2654.7650	0.0003	0.9998
GTR	19	-1266.4960	2655.8487	0.0002	1.0000
GTRIG	21	-1262.7454	2661.2799	0.0000	1.0000
F81I	15	-1287.8415	2672.6752	0.0000	1.0000
F81G	15	-1287.9344	2672.8609	0.0000	1.0000
F81	14	-1291.5852	2673.6963	0.0000	1.0000
F81IG	16	-1287.9343	2679.3270	0.0000	1.0000
K2P	12	-1335.8553	2749.3044	0.0000	1.0000
K2PI	13	-1332.9531	2749.9662	0.0000	1.0000
K2PG	13	-1333.0333	2750.1265	0.0000	1.0000
K2PIG	14	-1333.0534	2756.6328	0.0000	1.0000
SYM	16	-1331.9243	2767.3069	0.0000	1.0000
SYMI	17	-1329.0869	2768.0982	0.0000	1.0000
SYMG	17	-1329.1579	2768.2403	0.0000	1.0000
JC69	11	-1350.3745	2771.8766	0.0000	1.0000
JC69I	12	-1347.5048	2772.6033	0.0000	1.0000
JC69G	12	-1347.6030	2772.7997	0.0000	1.0000
SYMIG	18	-1329.1745	2774.7397	0.0000	1.0000
JC69IG	13	-1347.6263	2779.3125	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10859.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10859.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10859.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10859.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10859.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10859.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

