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

Input data from file "591-tomodel-phylip/chr7_10603.phylip" (7 taxa, 1111 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	-3125.0801	6282.1601	0.5550	0.5550
HKYIG	17	-3124.9413	6283.8826	0.2346	0.7896
HKYI	16	-3126.5541	6285.1082	0.1271	0.9167
GTRG	20	-3123.4750	6286.9499	0.0506	0.9673
GTRIG	21	-3123.3326	6288.6653	0.0215	0.9888
GTRI	20	-3124.9816	6289.9633	0.0112	1.0000
K2PG	13	-3151.8944	6329.7888	0.0000	1.0000
K2PIG	14	-3151.7718	6331.5436	0.0000	1.0000
K2PI	13	-3153.2854	6332.5708	0.0000	1.0000
SYMG	17	-3149.9568	6333.9136	0.0000	1.0000
SYMIG	18	-3149.8386	6335.6771	0.0000	1.0000
SYMI	17	-3151.3305	6336.6611	0.0000	1.0000
HKY	15	-3153.5141	6337.0282	0.0000	1.0000
GTR	19	-3151.9890	6341.9780	0.0000	1.0000
F81G	15	-3168.0695	6366.1390	0.0000	1.0000
F81IG	16	-3167.9794	6367.9587	0.0000	1.0000
F81I	15	-3169.2856	6368.5712	0.0000	1.0000
K2P	12	-3179.3957	6382.7914	0.0000	1.0000
SYM	16	-3177.3687	6386.7374	0.0000	1.0000
JC69G	12	-3190.9874	6405.9749	0.0000	1.0000
JC69IG	13	-3190.8993	6407.7987	0.0000	1.0000
JC69I	12	-3192.1871	6408.3742	0.0000	1.0000
F81	14	-3194.6429	6417.2859	0.0000	1.0000
JC69	11	-3217.2690	6456.5380	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3125.0801	6362.3684	0.7915	0.7915
HKYI	16	-3126.5541	6365.3165	0.1812	0.9727
HKYIG	17	-3124.9413	6369.1038	0.0273	1.0000
GTRG	20	-3123.4750	6387.2102	0.0000	1.0000
GTRI	20	-3124.9816	6390.2236	0.0000	1.0000
GTRIG	21	-3123.3326	6393.9386	0.0000	1.0000
K2PG	13	-3151.8944	6394.9580	0.0000	1.0000
K2PI	13	-3153.2854	6397.7400	0.0000	1.0000
K2PIG	14	-3151.7718	6401.7259	0.0000	1.0000
HKY	15	-3153.5141	6412.2235	0.0000	1.0000
SYMG	17	-3149.9568	6419.1349	0.0000	1.0000
SYMI	17	-3151.3305	6421.8823	0.0000	1.0000
SYMIG	18	-3149.8386	6425.9114	0.0000	1.0000
GTR	19	-3151.9890	6437.2253	0.0000	1.0000
F81G	15	-3168.0695	6441.3343	0.0000	1.0000
K2P	12	-3179.3957	6442.9476	0.0000	1.0000
F81I	15	-3169.2856	6443.7664	0.0000	1.0000
F81IG	16	-3167.9794	6448.1670	0.0000	1.0000
JC69G	12	-3190.9874	6466.1310	0.0000	1.0000
SYM	16	-3177.3687	6466.9457	0.0000	1.0000
JC69I	12	-3192.1871	6468.5304	0.0000	1.0000
JC69IG	13	-3190.8993	6472.9679	0.0000	1.0000
F81	14	-3194.6429	6487.4681	0.0000	1.0000
JC69	11	-3217.2690	6511.6812	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/chr7_10603.phylip = 1 - 1111;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10603.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/chr7_10603.phylip = 1 - 1111;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10603.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/chr7_10603.phylip = 1 - 1111;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10603.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

