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

Input data from file "591-tomodel-phylip/chr17_2981.phylip" (7 taxa, 1003 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	-3055.4050	6142.8100	0.4718	0.4718
HKYI	16	-3055.9875	6143.9750	0.2635	0.7353
HKYIG	17	-3055.5424	6145.0847	0.1513	0.8866
GTRG	20	-3053.4296	6146.8592	0.0623	0.9489
GTRI	20	-3054.2226	6148.4452	0.0282	0.9771
GTRIG	21	-3053.4296	6148.8592	0.0229	1.0000
HKY	15	-3103.1875	6236.3750	0.0000	1.0000
GTR	19	-3101.0689	6240.1377	0.0000	1.0000
F81I	15	-3119.2630	6268.5261	0.0000	1.0000
F81G	15	-3119.3713	6268.7425	0.0000	1.0000
F81IG	16	-3119.5470	6271.0940	0.0000	1.0000
K2PG	13	-3137.2814	6300.5628	0.0000	1.0000
K2PI	13	-3137.7188	6301.4376	0.0000	1.0000
K2PIG	14	-3137.3850	6302.7699	0.0000	1.0000
SYMI	17	-3134.7271	6303.4542	0.0000	1.0000
SYMG	17	-3135.8850	6305.7700	0.0000	1.0000
SYMIG	18	-3135.8925	6307.7850	0.0000	1.0000
F81	14	-3164.9872	6357.9743	0.0000	1.0000
JC69G	12	-3188.1926	6400.3851	0.0000	1.0000
JC69I	12	-3188.4344	6400.8688	0.0000	1.0000
JC69IG	13	-3188.3188	6402.6375	0.0000	1.0000
K2P	12	-3189.5663	6403.1327	0.0000	1.0000
SYM	16	-3187.8261	6407.6522	0.0000	1.0000
JC69	11	-3237.0959	6496.1919	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3055.4050	6221.3820	0.6305	0.6305
HKYI	16	-3055.9875	6222.5470	0.3521	0.9826
HKYIG	17	-3055.5424	6228.5675	0.0174	1.0000
GTRG	20	-3053.4296	6245.0742	0.0000	1.0000
GTRI	20	-3054.2226	6246.6602	0.0000	1.0000
GTRIG	21	-3053.4296	6251.9850	0.0000	1.0000
HKY	15	-3103.1875	6310.0363	0.0000	1.0000
GTR	19	-3101.0689	6333.4420	0.0000	1.0000
F81I	15	-3119.2630	6342.1873	0.0000	1.0000
F81G	15	-3119.3713	6342.4038	0.0000	1.0000
F81IG	16	-3119.5470	6349.6660	0.0000	1.0000
K2PG	13	-3137.2814	6364.4026	0.0000	1.0000
K2PI	13	-3137.7188	6365.2774	0.0000	1.0000
K2PIG	14	-3137.3850	6371.5204	0.0000	1.0000
SYMI	17	-3134.7271	6386.9370	0.0000	1.0000
SYMG	17	-3135.8850	6389.2527	0.0000	1.0000
SYMIG	18	-3135.8925	6396.1785	0.0000	1.0000
F81	14	-3164.9872	6426.7248	0.0000	1.0000
JC69G	12	-3188.1926	6459.3141	0.0000	1.0000
JC69I	12	-3188.4344	6459.7978	0.0000	1.0000
K2P	12	-3189.5663	6462.0617	0.0000	1.0000
JC69IG	13	-3188.3188	6466.4773	0.0000	1.0000
SYM	16	-3187.8261	6486.2242	0.0000	1.0000
JC69	11	-3237.0959	6550.2101	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/chr17_2981.phylip = 1 - 1003;
 Partition Dummy = 1:591-tomodel-phylip/chr17_2981.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/chr17_2981.phylip = 1 - 1003;
 Partition Dummy = 1:591-tomodel-phylip/chr17_2981.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/chr17_2981.phylip = 1 - 1003;
 Partition Dummy = 1:591-tomodel-phylip/chr17_2981.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

