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

Input data from file "591-tomodel-phylip/chr8_4340.phylip" (7 taxa, 871 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	-2541.1957	5114.3913	0.4551	0.4551
HKYIG	17	-2540.7299	5115.4597	0.2667	0.7218
GTRG	20	-2538.3525	5116.7049	0.1431	0.8649
GTRIG	21	-2537.8956	5117.7911	0.0831	0.9481
HKYI	16	-2543.6288	5119.2575	0.0399	0.9880
GTRI	20	-2540.8327	5121.6654	0.0120	1.0000
K2PG	13	-2558.0363	5142.0725	0.0000	1.0000
K2PIG	14	-2557.6336	5143.2672	0.0000	1.0000
SYMG	17	-2555.0203	5144.0406	0.0000	1.0000
SYMIG	18	-2554.6905	5145.3810	0.0000	1.0000
K2PI	13	-2560.0322	5146.0643	0.0000	1.0000
SYMI	17	-2556.8396	5147.6792	0.0000	1.0000
HKY	15	-2578.3637	5186.7274	0.0000	1.0000
GTR	19	-2575.5858	5189.1717	0.0000	1.0000
F81G	15	-2587.1692	5204.3385	0.0000	1.0000
F81IG	16	-2586.8337	5205.6675	0.0000	1.0000
F81I	15	-2589.2135	5208.4269	0.0000	1.0000
K2P	12	-2594.6242	5213.2483	0.0000	1.0000
SYM	16	-2591.0167	5214.0335	0.0000	1.0000
JC69G	12	-2602.1748	5228.3496	0.0000	1.0000
JC69IG	13	-2601.9040	5229.8081	0.0000	1.0000
JC69I	12	-2603.6437	5231.2874	0.0000	1.0000
F81	14	-2622.8070	5273.6141	0.0000	1.0000
JC69	11	-2636.9448	5295.8897	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2541.1957	5190.7056	0.8747	0.8747
HKYI	16	-2543.6288	5195.5718	0.0768	0.9515
HKYIG	17	-2540.7299	5196.5436	0.0472	0.9987
K2PG	13	-2558.0363	5204.0779	0.0011	0.9998
K2PI	13	-2560.0322	5208.0697	0.0001	0.9999
K2PIG	14	-2557.6336	5210.0422	0.0001	1.0000
GTRG	20	-2538.3525	5212.0978	0.0000	1.0000
GTRI	20	-2540.8327	5217.0582	0.0000	1.0000
GTRIG	21	-2537.8956	5217.9536	0.0000	1.0000
SYMG	17	-2555.0203	5225.1245	0.0000	1.0000
SYMI	17	-2556.8396	5228.7631	0.0000	1.0000
SYMIG	18	-2554.6905	5231.2346	0.0000	1.0000
HKY	15	-2578.3637	5258.2720	0.0000	1.0000
K2P	12	-2594.6242	5270.4840	0.0000	1.0000
F81G	15	-2587.1692	5275.8831	0.0000	1.0000
GTR	19	-2575.5858	5279.7949	0.0000	1.0000
F81I	15	-2589.2135	5279.9716	0.0000	1.0000
F81IG	16	-2586.8337	5281.9817	0.0000	1.0000
JC69G	12	-2602.1748	5285.5853	0.0000	1.0000
JC69I	12	-2603.6437	5288.5231	0.0000	1.0000
SYM	16	-2591.0167	5290.3478	0.0000	1.0000
JC69IG	13	-2601.9040	5291.8134	0.0000	1.0000
F81	14	-2622.8070	5340.3890	0.0000	1.0000
JC69	11	-2636.9448	5348.3557	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 591-tomodel-phylip/chr8_4340.phylip = 1 - 871;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4340.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 591-tomodel-phylip/chr8_4340.phylip = 1 - 871;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4340.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 591-tomodel-phylip/chr8_4340.phylip = 1 - 871;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4340.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

