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

Input data from file "683-tomodel-phylip/chr3_20422.phylip" (7 taxa, 998 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	-2887.8297	5807.6595	0.4055	0.4055
HKYI	16	-2888.0126	5808.0252	0.3377	0.7432
HKYIG	17	-2887.8335	5809.6669	0.1486	0.8918
GTRG	20	-2886.0761	5812.1521	0.0429	0.9347
GTRI	20	-2886.2743	5812.5485	0.0352	0.9699
GTRIG	21	-2885.7584	5813.5167	0.0217	0.9915
HKY	15	-2893.0304	5816.0609	0.0061	0.9976
GTR	19	-2889.9668	5817.9336	0.0024	1.0000
K2PI	13	-2925.5864	5877.1728	0.0000	1.0000
K2PG	13	-2925.6726	5877.3451	0.0000	1.0000
K2PIG	14	-2925.5638	5879.1276	0.0000	1.0000
SYMG	17	-2924.0196	5882.0392	0.0000	1.0000
SYMI	17	-2924.2710	5882.5421	0.0000	1.0000
SYMIG	18	-2923.8955	5883.7909	0.0000	1.0000
K2P	12	-2931.3313	5886.6627	0.0000	1.0000
SYM	16	-2929.9716	5891.9431	0.0000	1.0000
F81I	15	-2994.9870	6019.9740	0.0000	1.0000
F81G	15	-2995.4352	6020.8703	0.0000	1.0000
F81IG	16	-2994.9981	6021.9962	0.0000	1.0000
F81	14	-2998.7396	6025.4793	0.0000	1.0000
JC69I	12	-3024.9967	6073.9933	0.0000	1.0000
JC69G	12	-3025.4995	6074.9989	0.0000	1.0000
JC69IG	13	-3025.0081	6076.0163	0.0000	1.0000
JC69	11	-3029.5284	6081.0568	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2887.8297	5886.1515	0.4905	0.4905
HKYI	16	-2888.0126	5886.5172	0.4086	0.8991
HKY	15	-2893.0304	5889.6472	0.0854	0.9845
HKYIG	17	-2887.8335	5893.0647	0.0155	1.0000
GTRG	20	-2886.0761	5910.2672	0.0000	1.0000
GTRI	20	-2886.2743	5910.6636	0.0000	1.0000
GTR	19	-2889.9668	5911.1429	0.0000	1.0000
GTRIG	21	-2885.7584	5916.5375	0.0000	1.0000
K2PI	13	-2925.5864	5940.9476	0.0000	1.0000
K2PG	13	-2925.6726	5941.1199	0.0000	1.0000
K2P	12	-2931.3313	5945.5317	0.0000	1.0000
K2PIG	14	-2925.5638	5947.8081	0.0000	1.0000
SYMG	17	-2924.0196	5965.4370	0.0000	1.0000
SYMI	17	-2924.2710	5965.9399	0.0000	1.0000
SYM	16	-2929.9716	5970.4352	0.0000	1.0000
SYMIG	18	-2923.8955	5972.0945	0.0000	1.0000
F81I	15	-2994.9870	6093.5603	0.0000	1.0000
F81	14	-2998.7396	6094.1598	0.0000	1.0000
F81G	15	-2995.4352	6094.4566	0.0000	1.0000
F81IG	16	-2994.9981	6100.4883	0.0000	1.0000
JC69I	12	-3024.9967	6132.8624	0.0000	1.0000
JC69G	12	-3025.4995	6133.8680	0.0000	1.0000
JC69	11	-3029.5284	6135.0201	0.0000	1.0000
JC69IG	13	-3025.0081	6139.7911	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 683-tomodel-phylip/chr3_20422.phylip = 1 - 998;
 Partition Dummy = 1:683-tomodel-phylip/chr3_20422.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 683-tomodel-phylip/chr3_20422.phylip = 1 - 998;
 Partition Dummy = 1:683-tomodel-phylip/chr3_20422.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 683-tomodel-phylip/chr3_20422.phylip = 1 - 998;
 Partition Dummy = 1:683-tomodel-phylip/chr3_20422.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

