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

Input data from file "591-tomodel-phylip/chr12_5908.phylip" (7 taxa, 943 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	-3372.3136	6776.6272	0.4714	0.4714
HKYI	16	-3373.2530	6778.5060	0.1842	0.6556
HKYIG	17	-3372.3142	6778.6284	0.1733	0.8289
GTRG	20	-3369.9452	6779.8904	0.0922	0.9211
GTRI	20	-3370.6838	6781.3677	0.0441	0.9652
GTRIG	21	-3369.9193	6781.8387	0.0348	1.0000
HKY	15	-3394.4266	6818.8532	0.0000	1.0000
GTR	19	-3392.5179	6823.0357	0.0000	1.0000
K2PG	13	-3423.8317	6873.6633	0.0000	1.0000
K2PI	13	-3423.9343	6873.8685	0.0000	1.0000
SYMG	17	-3420.4879	6874.9758	0.0000	1.0000
K2PIG	14	-3423.8347	6875.6694	0.0000	1.0000
SYMI	17	-3421.3628	6876.7256	0.0000	1.0000
SYMIG	18	-3420.5028	6877.0055	0.0000	1.0000
F81I	15	-3428.7461	6887.4921	0.0000	1.0000
F81IG	16	-3427.9467	6887.8934	0.0000	1.0000
F81G	15	-3428.9788	6887.9575	0.0000	1.0000
K2P	12	-3444.8277	6913.6553	0.0000	1.0000
SYM	16	-3441.5194	6915.0388	0.0000	1.0000
F81	14	-3448.8761	6925.7522	0.0000	1.0000
JC69G	12	-3469.3321	6962.6642	0.0000	1.0000
JC69I	12	-3470.1292	6964.2584	0.0000	1.0000
JC69IG	13	-3469.3482	6964.6965	0.0000	1.0000
JC69	11	-3488.8867	6999.7735	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3372.3136	6854.2123	0.7025	0.7025
HKYI	16	-3373.2530	6856.0911	0.2746	0.9771
HKYIG	17	-3372.3142	6861.0625	0.0229	1.0000
GTRG	20	-3369.9452	6876.8718	0.0000	1.0000
GTRI	20	-3370.6838	6878.3490	0.0000	1.0000
GTRIG	21	-3369.9193	6883.6691	0.0000	1.0000
HKY	15	-3394.4266	6891.5892	0.0000	1.0000
GTR	19	-3392.5179	6915.1680	0.0000	1.0000
K2PG	13	-3423.8317	6936.7012	0.0000	1.0000
K2PI	13	-3423.9343	6936.9064	0.0000	1.0000
K2PIG	14	-3423.8347	6943.5563	0.0000	1.0000
SYMG	17	-3420.4879	6957.4099	0.0000	1.0000
SYMI	17	-3421.3628	6959.1597	0.0000	1.0000
F81I	15	-3428.7461	6960.2281	0.0000	1.0000
F81G	15	-3428.9788	6960.6935	0.0000	1.0000
SYMIG	18	-3420.5028	6964.2887	0.0000	1.0000
F81IG	16	-3427.9467	6965.4785	0.0000	1.0000
K2P	12	-3444.8277	6971.8441	0.0000	1.0000
SYM	16	-3441.5194	6992.6239	0.0000	1.0000
F81	14	-3448.8761	6993.6392	0.0000	1.0000
JC69G	12	-3469.3321	7020.8530	0.0000	1.0000
JC69I	12	-3470.1292	7022.4472	0.0000	1.0000
JC69IG	13	-3469.3482	7027.7343	0.0000	1.0000
JC69	11	-3488.8867	7053.1132	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/chr12_5908.phylip = 1 - 943;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5908.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/chr12_5908.phylip = 1 - 943;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5908.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/chr12_5908.phylip = 1 - 943;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5908.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

