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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2805.3058	5650.6115	0.5772	0.5772
GTRIG	21	-2805.2851	5652.5702	0.2168	0.7940
GTRG	20	-2806.3363	5652.6727	0.2059	0.9999
HKYI	16	-2818.9614	5669.9229	0.0000	1.0000
HKYG	16	-2819.3874	5670.7748	0.0000	1.0000
HKYIG	17	-2818.9175	5671.8350	0.0000	1.0000
GTR	19	-2840.1364	5718.2728	0.0000	1.0000
HKY	15	-2852.4599	5734.9197	0.0000	1.0000
SYMI	17	-2860.5711	5755.1422	0.0000	1.0000
SYMG	17	-2860.5862	5755.1725	0.0000	1.0000
K2PG	13	-2865.0280	5756.0561	0.0000	1.0000
SYMIG	18	-2860.3646	5756.7291	0.0000	1.0000
K2PI	13	-2865.3991	5756.7982	0.0000	1.0000
K2PIG	14	-2865.0515	5758.1029	0.0000	1.0000
F81G	15	-2869.0294	5768.0588	0.0000	1.0000
F81I	15	-2869.1190	5768.2380	0.0000	1.0000
F81IG	16	-2869.0624	5770.1248	0.0000	1.0000
SYM	16	-2893.7432	5819.4864	0.0000	1.0000
K2P	12	-2898.3233	5820.6466	0.0000	1.0000
F81	14	-2899.8766	5827.7532	0.0000	1.0000
JC69G	12	-2907.0453	5838.0906	0.0000	1.0000
JC69I	12	-2907.4824	5838.9647	0.0000	1.0000
JC69IG	13	-2907.0537	5840.1074	0.0000	1.0000
JC69	11	-2938.0944	5898.1888	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2818.9614	5747.3030	0.3274	0.3274
GTRI	20	-2805.3058	5747.3367	0.3219	0.6493
HKYG	16	-2819.3874	5748.1549	0.2138	0.8632
GTRG	20	-2806.3363	5749.3979	0.1149	0.9780
HKYIG	17	-2818.9175	5754.0514	0.0112	0.9892
GTRIG	21	-2805.2851	5754.1317	0.0108	1.0000
HKY	15	-2852.4599	5807.4636	0.0000	1.0000
GTR	19	-2840.1364	5810.1617	0.0000	1.0000
K2PG	13	-2865.0280	5818.9274	0.0000	1.0000
K2PI	13	-2865.3991	5819.6696	0.0000	1.0000
K2PIG	14	-2865.0515	5825.8106	0.0000	1.0000
SYMI	17	-2860.5711	5837.3586	0.0000	1.0000
SYMG	17	-2860.5862	5837.3889	0.0000	1.0000
F81G	15	-2869.0294	5840.6027	0.0000	1.0000
F81I	15	-2869.1190	5840.7819	0.0000	1.0000
SYMIG	18	-2860.3646	5843.7818	0.0000	1.0000
F81IG	16	-2869.0624	5847.5050	0.0000	1.0000
K2P	12	-2898.3233	5878.6817	0.0000	1.0000
F81	14	-2899.8766	5895.4608	0.0000	1.0000
JC69G	12	-2907.0453	5896.1257	0.0000	1.0000
SYM	16	-2893.7432	5896.8666	0.0000	1.0000
JC69I	12	-2907.4824	5896.9998	0.0000	1.0000
JC69IG	13	-2907.0537	5902.9788	0.0000	1.0000
JC69	11	-2938.0944	5951.3877	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/chr2_27234.phylip = 1 - 931;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27234.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/chr2_27234.phylip = 1 - 931;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27234.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/chr2_27234.phylip = 1 - 931;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27234.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

