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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2P


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	-1349.5819	2740.9978	0.3807	0.3807
GTRG	20	-1349.5951	2741.0243	0.3757	0.7564
GTRIG	21	-1349.5853	2743.1925	0.1271	0.8835
GTR	19	-1352.0558	2743.7675	0.0953	0.9788
HKYI	16	-1357.8189	2748.8154	0.0076	0.9865
HKYG	16	-1357.8436	2748.8646	0.0075	0.9939
HKYIG	17	-1357.8315	2750.9905	0.0026	0.9965
HKY	15	-1360.5480	2752.1326	0.0015	0.9979
SYMG	17	-1359.0715	2753.4706	0.0007	0.9987
SYMI	17	-1359.1621	2753.6517	0.0007	0.9994
SYMIG	18	-1359.0795	2755.6460	0.0003	0.9996
SYM	16	-1361.2847	2755.7469	0.0002	0.9999
K2PG	13	-1365.9518	2758.6864	0.0001	0.9999
K2PI	13	-1366.0078	2758.7984	0.0001	1.0000
K2PIG	14	-1365.9526	2760.8103	0.0000	1.0000
K2P	12	-1368.5717	2761.8129	0.0000	1.0000
F81I	15	-1386.1363	2803.3093	0.0000	1.0000
F81G	15	-1386.2104	2803.4574	0.0000	1.0000
F81IG	16	-1386.2023	2805.5821	0.0000	1.0000
F81	14	-1388.6351	2806.1753	0.0000	1.0000
JC69I	12	-1393.0551	2810.7797	0.0000	1.0000
JC69G	12	-1393.1587	2810.9870	0.0000	1.0000
JC69IG	13	-1393.1565	2813.0959	0.0000	1.0000
JC69	11	-1395.4381	2813.4415	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1349.5819	2739.1638	0.3812	0.3812
GTRG	20	-1349.5951	2739.1903	0.3762	0.7575
GTRIG	21	-1349.5853	2741.1706	0.1398	0.8973
GTR	19	-1352.0558	2742.1117	0.0873	0.9846
HKYI	16	-1357.8189	2747.6379	0.0055	0.9901
HKYG	16	-1357.8436	2747.6871	0.0054	0.9954
HKYIG	17	-1357.8315	2749.6630	0.0020	0.9974
HKY	15	-1360.5480	2751.0959	0.0010	0.9984
SYMG	17	-1359.0715	2752.1431	0.0006	0.9990
SYMI	17	-1359.1621	2752.3242	0.0005	0.9995
SYMIG	18	-1359.0795	2754.1591	0.0002	0.9997
SYM	16	-1361.2847	2754.5694	0.0002	0.9999
K2PG	13	-1365.9518	2757.9036	0.0000	1.0000
K2PI	13	-1366.0078	2758.0156	0.0000	1.0000
K2PIG	14	-1365.9526	2759.9051	0.0000	1.0000
K2P	12	-1368.5717	2761.1433	0.0000	1.0000
F81I	15	-1386.1363	2802.2726	0.0000	1.0000
F81G	15	-1386.2104	2802.4207	0.0000	1.0000
F81IG	16	-1386.2023	2804.4046	0.0000	1.0000
F81	14	-1388.6351	2805.2701	0.0000	1.0000
JC69I	12	-1393.0551	2810.1102	0.0000	1.0000
JC69G	12	-1393.1587	2810.3174	0.0000	1.0000
JC69IG	13	-1393.1565	2812.3131	0.0000	1.0000
JC69	11	-1395.4381	2812.8762	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-1368.5717	2811.2037	0.3359	0.3359
K2PG	13	-1365.9518	2812.1357	0.2108	0.5467
K2PI	13	-1366.0078	2812.2477	0.1993	0.7460
HKY	15	-1360.5480	2813.6714	0.0978	0.8438
HKYI	16	-1357.8189	2814.3851	0.0685	0.9122
HKYG	16	-1357.8436	2814.4343	0.0668	0.9790
K2PIG	14	-1365.9526	2818.3089	0.0096	0.9887
HKYIG	17	-1357.8315	2820.5819	0.0031	0.9917
SYM	16	-1361.2847	2821.3166	0.0021	0.9939
GTR	19	-1352.0558	2821.3740	0.0021	0.9960
GTRI	20	-1349.5819	2822.5978	0.0011	0.9971
GTRG	20	-1349.5951	2822.6243	0.0011	0.9982
SYMG	17	-1359.0715	2823.0620	0.0009	0.9991
SYMI	17	-1359.1621	2823.2431	0.0008	0.9999
GTRIG	21	-1349.5853	2828.7763	0.0001	1.0000
SYMIG	18	-1359.0795	2829.2497	0.0000	1.0000
JC69	11	-1395.4381	2858.7649	0.0000	1.0000
JC69I	12	-1393.0551	2860.1706	0.0000	1.0000
JC69G	12	-1393.1587	2860.3778	0.0000	1.0000
F81	14	-1388.6351	2863.6739	0.0000	1.0000
F81I	15	-1386.1363	2864.8481	0.0000	1.0000
F81G	15	-1386.2104	2864.9962	0.0000	1.0000
JC69IG	13	-1393.1565	2866.5452	0.0000	1.0000
F81IG	16	-1386.2023	2871.1518	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/chr3_17607.phylip = 1 - 479;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17607.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/chr3_17607.phylip = 1 - 479;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17607.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 (K2P)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_17607.phylip = 1 - 479;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17607.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

