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

Input data from file "683-tomodel-phylip/chr5_10135.phylip" (7 taxa, 990 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	-2918.5762	5877.1524	0.4376	0.4376
GTRG	20	-2918.9920	5877.9839	0.2888	0.7264
GTRIG	21	-2919.0390	5880.0779	0.1014	0.8277
HKYI	16	-2924.3240	5880.6480	0.0762	0.9039
HKYG	16	-2924.3894	5880.7789	0.0714	0.9753
HKYIG	17	-2924.4519	5882.9037	0.0247	1.0000
GTR	19	-2946.4709	5930.9418	0.0000	1.0000
HKY	15	-2950.6227	5931.2453	0.0000	1.0000
K2PI	13	-2979.2327	5984.4653	0.0000	1.0000
K2PG	13	-2979.4344	5984.8688	0.0000	1.0000
SYMI	17	-2975.5541	5985.1082	0.0000	1.0000
SYMG	17	-2976.0809	5986.1617	0.0000	1.0000
K2PIG	14	-2979.5252	5987.0504	0.0000	1.0000
SYMIG	18	-2976.1739	5988.3478	0.0000	1.0000
F81I	15	-2999.4746	6028.9491	0.0000	1.0000
F81G	15	-3000.3111	6030.6223	0.0000	1.0000
F81IG	16	-2999.4903	6030.9806	0.0000	1.0000
SYM	16	-3002.5617	6037.1233	0.0000	1.0000
K2P	12	-3007.0246	6038.0493	0.0000	1.0000
F81	14	-3024.8751	6077.7502	0.0000	1.0000
JC69I	12	-3045.2786	6114.5573	0.0000	1.0000
JC69G	12	-3046.1184	6116.2368	0.0000	1.0000
JC69IG	13	-3045.2959	6116.5918	0.0000	1.0000
JC69	11	-3070.7326	6163.4652	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2924.3240	5959.0113	0.5089	0.5089
HKYG	16	-2924.3894	5959.1422	0.4766	0.9855
HKYIG	17	-2924.4519	5966.1647	0.0142	0.9997
GTRI	20	-2918.5762	5975.1065	0.0002	0.9999
GTRG	20	-2918.9920	5975.9380	0.0001	1.0000
GTRIG	21	-2919.0390	5982.9297	0.0000	1.0000
HKY	15	-2950.6227	6004.7109	0.0000	1.0000
GTR	19	-2946.4709	6023.9982	0.0000	1.0000
K2PI	13	-2979.2327	6048.1355	0.0000	1.0000
K2PG	13	-2979.4344	6048.5390	0.0000	1.0000
K2PIG	14	-2979.5252	6055.6183	0.0000	1.0000
SYMI	17	-2975.5541	6068.3692	0.0000	1.0000
SYMG	17	-2976.0809	6069.4227	0.0000	1.0000
SYMIG	18	-2976.1739	6076.5065	0.0000	1.0000
K2P	12	-3007.0246	6096.8217	0.0000	1.0000
F81I	15	-2999.4746	6102.4147	0.0000	1.0000
F81G	15	-3000.3111	6104.0878	0.0000	1.0000
F81IG	16	-2999.4903	6109.3439	0.0000	1.0000
SYM	16	-3002.5617	6115.4866	0.0000	1.0000
F81	14	-3024.8751	6146.3181	0.0000	1.0000
JC69I	12	-3045.2786	6173.3297	0.0000	1.0000
JC69G	12	-3046.1184	6175.0093	0.0000	1.0000
JC69IG	13	-3045.2959	6180.2620	0.0000	1.0000
JC69	11	-3070.7326	6217.3400	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/chr5_10135.phylip = 1 - 990;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10135.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/chr5_10135.phylip = 1 - 990;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10135.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/chr5_10135.phylip = 1 - 990;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10135.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

