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

Input data from file "683-tomodel-phylip/chr5_1737.phylip" (7 taxa, 1001 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2852.3603	5736.7206	0.3792	0.3792
GTRG	20	-2848.6899	5737.3799	0.2727	0.6520
HKYIG	17	-2852.1467	5738.2933	0.1727	0.8247
GTRIG	21	-2848.4434	5738.8867	0.1284	0.9531
HKYI	16	-2855.0051	5742.0101	0.0269	0.9800
GTRI	20	-2851.3047	5742.6095	0.0200	1.0000
HKY	15	-2878.4745	5786.9491	0.0000	1.0000
GTR	19	-2876.0676	5790.1353	0.0000	1.0000
K2PG	13	-2901.9862	5829.9723	0.0000	1.0000
K2PIG	14	-2901.7608	5831.5215	0.0000	1.0000
K2PI	13	-2905.1177	5836.2354	0.0000	1.0000
SYMG	17	-2901.3740	5836.7479	0.0000	1.0000
SYMIG	18	-2901.1373	5838.2747	0.0000	1.0000
SYMI	17	-2904.6446	5843.2893	0.0000	1.0000
F81G	15	-2909.1294	5848.2589	0.0000	1.0000
F81IG	16	-2909.0189	5850.0378	0.0000	1.0000
F81I	15	-2911.5753	5853.1507	0.0000	1.0000
K2P	12	-2926.2436	5876.4873	0.0000	1.0000
SYM	16	-2925.5205	5883.0410	0.0000	1.0000
F81	14	-2931.7222	5891.4444	0.0000	1.0000
JC69G	12	-2950.0530	5924.1059	0.0000	1.0000
JC69IG	13	-2949.9371	5925.8742	0.0000	1.0000
JC69I	12	-2952.7118	5929.4235	0.0000	1.0000
JC69	11	-2972.1639	5966.3278	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2852.3603	5815.2607	0.9007	0.9007
HKYI	16	-2855.0051	5820.5502	0.0640	0.9647
HKYIG	17	-2852.1467	5821.7422	0.0353	1.0000
GTRG	20	-2848.6899	5835.5550	0.0000	1.0000
GTRI	20	-2851.3047	5840.7846	0.0000	1.0000
GTRIG	21	-2848.4434	5841.9706	0.0000	1.0000
HKY	15	-2878.4745	5860.5804	0.0000	1.0000
GTR	19	-2876.0676	5883.4016	0.0000	1.0000
K2PG	13	-2901.9862	5893.7861	0.0000	1.0000
K2PI	13	-2905.1177	5900.0492	0.0000	1.0000
K2PIG	14	-2901.7608	5900.2441	0.0000	1.0000
SYMG	17	-2901.3740	5920.1968	0.0000	1.0000
F81G	15	-2909.1294	5921.8902	0.0000	1.0000
SYMIG	18	-2901.1373	5926.6323	0.0000	1.0000
SYMI	17	-2904.6446	5926.7381	0.0000	1.0000
F81I	15	-2911.5753	5926.7820	0.0000	1.0000
F81IG	16	-2909.0189	5928.5779	0.0000	1.0000
K2P	12	-2926.2436	5935.3923	0.0000	1.0000
F81	14	-2931.7222	5960.1669	0.0000	1.0000
SYM	16	-2925.5205	5961.5810	0.0000	1.0000
JC69G	12	-2950.0530	5983.0110	0.0000	1.0000
JC69I	12	-2952.7118	5988.3286	0.0000	1.0000
JC69IG	13	-2949.9371	5989.6880	0.0000	1.0000
JC69	11	-2972.1639	6020.3241	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_1737.phylip = 1 - 1001;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1737.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_1737.phylip = 1 - 1001;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1737.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_1737.phylip = 1 - 1001;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1737.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

