
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:09 04/04/2011

Input data from file "683-tomodel-phylip/chr12_5903.phylip" (7 taxa, 1014 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	-3466.5525	6973.1049	0.4510	0.4510
GTRG	20	-3466.9898	6973.9796	0.2912	0.7422
GTRIG	21	-3466.4033	6974.8067	0.1926	0.9348
HKYI	16	-3473.1709	6978.3418	0.0329	0.9677
HKYG	16	-3473.7124	6979.4248	0.0191	0.9868
HKYIG	17	-3473.0847	6980.1693	0.0132	1.0000
SYMI	17	-3480.7374	6995.4747	0.0000	1.0000
SYMG	17	-3481.1630	6996.3259	0.0000	1.0000
SYMIG	18	-3480.6009	6997.2019	0.0000	1.0000
K2PI	13	-3485.7474	6997.4949	0.0000	1.0000
K2PG	13	-3486.2923	6998.5845	0.0000	1.0000
K2PIG	14	-3485.6734	6999.3469	0.0000	1.0000
GTR	19	-3493.1083	7024.2167	0.0000	1.0000
HKY	15	-3499.9608	7029.9215	0.0000	1.0000
SYM	16	-3507.6844	7047.3688	0.0000	1.0000
K2P	12	-3512.7852	7049.5705	0.0000	1.0000
F81I	15	-3578.9284	7187.8568	0.0000	1.0000
F81IG	16	-3578.9333	7189.8665	0.0000	1.0000
F81G	15	-3580.0231	7190.0462	0.0000	1.0000
JC69I	12	-3592.8929	7209.7859	0.0000	1.0000
JC69IG	13	-3592.8964	7211.7929	0.0000	1.0000
JC69G	12	-3593.9300	7211.8599	0.0000	1.0000
F81	14	-3601.0972	7230.1943	0.0000	1.0000
JC69	11	-3615.2326	7252.4652	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3473.1709	7057.0884	0.5564	0.5564
HKYG	16	-3473.7124	7058.1713	0.3238	0.8802
K2PI	13	-3485.7474	7061.4764	0.0620	0.9422
K2PG	13	-3486.2923	7062.5661	0.0360	0.9782
HKYIG	17	-3473.0847	7063.8375	0.0190	0.9972
K2PIG	14	-3485.6734	7068.2501	0.0021	0.9993
GTRI	20	-3466.5525	7071.5381	0.0004	0.9997
GTRG	20	-3466.9898	7072.4128	0.0003	1.0000
GTRIG	21	-3466.4033	7078.1615	0.0000	1.0000
SYMI	17	-3480.7374	7079.1429	0.0000	1.0000
SYMG	17	-3481.1630	7079.9941	0.0000	1.0000
SYMIG	18	-3480.6009	7085.7917	0.0000	1.0000
HKY	15	-3499.9608	7103.7464	0.0000	1.0000
K2P	12	-3512.7852	7108.6304	0.0000	1.0000
GTR	19	-3493.1083	7117.7282	0.0000	1.0000
SYM	16	-3507.6844	7126.1153	0.0000	1.0000
F81I	15	-3578.9284	7261.6817	0.0000	1.0000
F81G	15	-3580.0231	7263.8711	0.0000	1.0000
F81IG	16	-3578.9333	7268.6130	0.0000	1.0000
JC69I	12	-3592.8929	7268.8458	0.0000	1.0000
JC69G	12	-3593.9300	7270.9198	0.0000	1.0000
JC69IG	13	-3592.8964	7275.7744	0.0000	1.0000
F81	14	-3601.0972	7299.0975	0.0000	1.0000
JC69	11	-3615.2326	7306.6034	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 683-tomodel-phylip/chr12_5903.phylip = 1 - 1014;
 Partition Dummy = 1:683-tomodel-phylip/chr12_5903.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 683-tomodel-phylip/chr12_5903.phylip = 1 - 1014;
 Partition Dummy = 1:683-tomodel-phylip/chr12_5903.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 683-tomodel-phylip/chr12_5903.phylip = 1 - 1014;
 Partition Dummy = 1:683-tomodel-phylip/chr12_5903.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

