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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2989.4865	6018.9730	0.6903	0.6903
GTRIG	21	-2989.4714	6020.9429	0.2578	0.9481
GTRI	20	-2992.0987	6024.1974	0.0506	0.9987
HKYG	16	-3000.1922	6032.3845	0.0008	0.9996
HKYIG	17	-3000.1831	6034.3661	0.0003	0.9999
HKYI	16	-3002.1316	6036.2632	0.0001	1.0000
GTR	19	-3009.5555	6057.1110	0.0000	1.0000
HKY	15	-3020.3524	6070.7049	0.0000	1.0000
K2PG	13	-3069.3911	6164.7823	0.0000	1.0000
K2PI	13	-3070.2590	6166.5180	0.0000	1.0000
K2PIG	14	-3069.3793	6166.7585	0.0000	1.0000
SYMG	17	-3066.6180	6167.2359	0.0000	1.0000
SYMI	17	-3067.3949	6168.7897	0.0000	1.0000
SYMIG	18	-3066.5964	6169.1928	0.0000	1.0000
K2P	12	-3091.2191	6206.4381	0.0000	1.0000
SYM	16	-3087.6772	6207.3543	0.0000	1.0000
F81G	15	-3088.8591	6207.7183	0.0000	1.0000
F81IG	16	-3088.8620	6209.7240	0.0000	1.0000
F81I	15	-3090.0263	6210.0526	0.0000	1.0000
F81	14	-3109.1669	6246.3338	0.0000	1.0000
JC69G	12	-3145.4781	6314.9562	0.0000	1.0000
JC69I	12	-3145.8180	6315.6359	0.0000	1.0000
JC69IG	13	-3145.4853	6316.9707	0.0000	1.0000
JC69	11	-3165.7029	6353.4058	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3000.1922	6110.8605	0.8161	0.8161
HKYI	16	-3002.1316	6114.7393	0.1173	0.9334
GTRG	20	-2989.4865	6117.0680	0.0366	0.9701
HKYIG	17	-3000.1831	6117.7469	0.0261	0.9961
GTRI	20	-2992.0987	6122.2924	0.0027	0.9988
GTRIG	21	-2989.4714	6123.9426	0.0012	1.0000
HKY	15	-3020.3524	6144.2761	0.0000	1.0000
GTR	19	-3009.5555	6150.3013	0.0000	1.0000
K2PG	13	-3069.3911	6228.5440	0.0000	1.0000
K2PI	13	-3070.2590	6230.2797	0.0000	1.0000
K2PIG	14	-3069.3793	6235.4250	0.0000	1.0000
SYMG	17	-3066.6180	6250.6167	0.0000	1.0000
SYMI	17	-3067.3949	6252.1705	0.0000	1.0000
SYMIG	18	-3066.5964	6257.4784	0.0000	1.0000
K2P	12	-3091.2191	6265.2951	0.0000	1.0000
F81G	15	-3088.8591	6281.2895	0.0000	1.0000
F81I	15	-3090.0263	6283.6238	0.0000	1.0000
SYM	16	-3087.6772	6285.8304	0.0000	1.0000
F81IG	16	-3088.8620	6288.2001	0.0000	1.0000
F81	14	-3109.1669	6315.0003	0.0000	1.0000
JC69G	12	-3145.4781	6373.8132	0.0000	1.0000
JC69I	12	-3145.8180	6374.4929	0.0000	1.0000
JC69IG	13	-3145.4853	6380.7324	0.0000	1.0000
JC69	11	-3165.7029	6407.3581	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr12_5739.phylip = 1 - 997;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5739.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr12_5739.phylip = 1 - 997;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5739.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

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

