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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-3454.2838	6948.5676	0.4722	0.4722
GTRG	20	-3454.4401	6948.8803	0.4038	0.8760
GTRIG	21	-3454.6244	6951.2489	0.1236	0.9995
HKYG	16	-3466.0060	6964.0121	0.0002	0.9997
HKYI	16	-3466.0846	6964.1691	0.0002	0.9999
HKYIG	17	-3466.1207	6966.2414	0.0001	1.0000
GTR	19	-3508.4621	7054.9243	0.0000	1.0000
HKY	15	-3518.4430	7066.8859	0.0000	1.0000
F81G	15	-3522.4890	7074.9780	0.0000	1.0000
F81I	15	-3522.5114	7075.0229	0.0000	1.0000
F81IG	16	-3522.6008	7077.2015	0.0000	1.0000
SYMG	17	-3521.9202	7077.8405	0.0000	1.0000
SYMI	17	-3522.3272	7078.6545	0.0000	1.0000
SYMIG	18	-3521.9651	7079.9302	0.0000	1.0000
K2PG	13	-3528.0843	7082.1685	0.0000	1.0000
K2PI	13	-3528.4072	7082.8144	0.0000	1.0000
K2PIG	14	-3528.1514	7084.3029	0.0000	1.0000
F81	14	-3570.4192	7168.8383	0.0000	1.0000
JC69G	12	-3573.9996	7171.9993	0.0000	1.0000
JC69I	12	-3574.1478	7172.2957	0.0000	1.0000
SYM	16	-3570.3091	7172.6181	0.0000	1.0000
JC69IG	13	-3573.9997	7173.9994	0.0000	1.0000
K2P	12	-3577.1042	7178.2083	0.0000	1.0000
JC69	11	-3620.3785	7262.7569	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3466.0060	7044.7306	0.4702	0.4702
HKYI	16	-3466.0846	7044.8876	0.4347	0.9049
GTRI	20	-3454.2838	7049.4657	0.0441	0.9490
GTRG	20	-3454.4401	7049.7784	0.0377	0.9867
HKYIG	17	-3466.1207	7052.0048	0.0124	0.9991
GTRIG	21	-3454.6244	7057.1919	0.0009	1.0000
HKY	15	-3518.4430	7142.5595	0.0000	1.0000
K2PG	13	-3528.0843	7147.7523	0.0000	1.0000
K2PI	13	-3528.4072	7148.3982	0.0000	1.0000
F81G	15	-3522.4890	7150.6516	0.0000	1.0000
F81I	15	-3522.5114	7150.6965	0.0000	1.0000
GTR	19	-3508.4621	7150.7775	0.0000	1.0000
K2PIG	14	-3528.1514	7154.9316	0.0000	1.0000
F81IG	16	-3522.6008	7157.9200	0.0000	1.0000
SYMG	17	-3521.9202	7163.6039	0.0000	1.0000
SYMI	17	-3522.3272	7164.4179	0.0000	1.0000
SYMIG	18	-3521.9651	7170.7385	0.0000	1.0000
JC69G	12	-3573.9996	7232.5381	0.0000	1.0000
JC69I	12	-3574.1478	7232.8346	0.0000	1.0000
K2P	12	-3577.1042	7238.7472	0.0000	1.0000
F81	14	-3570.4192	7239.4670	0.0000	1.0000
JC69IG	13	-3573.9997	7239.5831	0.0000	1.0000
SYM	16	-3570.3091	7253.3366	0.0000	1.0000
JC69	11	-3620.3785	7318.2509	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/chr1_5427.phylip = 1 - 1147;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5427.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/chr1_5427.phylip = 1 - 1147;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5427.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5427.phylip = 1 - 1147;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5427.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

