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

Input data from file "591-tomodel-phylip/chr2_22071.phylip" (7 taxa, 956 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	-2985.8544	6003.7087	0.5838	0.5838
HKYIG	17	-2985.7514	6005.5028	0.2381	0.8219
GTRG	20	-2983.5693	6007.1386	0.1051	0.9270
GTRIG	21	-2983.4883	6008.9766	0.0419	0.9689
HKYI	16	-2988.9793	6009.9587	0.0257	0.9945
GTRI	20	-2986.5256	6013.0513	0.0055	1.0000
HKY	15	-3019.4258	6068.8516	0.0000	1.0000
K2PG	13	-3022.5737	6071.1474	0.0000	1.0000
GTR	19	-3017.1232	6072.2464	0.0000	1.0000
SYMG	17	-3019.3372	6072.6744	0.0000	1.0000
K2PIG	14	-3022.4867	6072.9733	0.0000	1.0000
SYMIG	18	-3019.2596	6074.5192	0.0000	1.0000
K2PI	13	-3025.5354	6077.0708	0.0000	1.0000
SYMI	17	-3022.3620	6078.7241	0.0000	1.0000
F81G	15	-3044.0514	6118.1027	0.0000	1.0000
F81IG	16	-3044.0080	6120.0160	0.0000	1.0000
F81I	15	-3046.2194	6122.4388	0.0000	1.0000
K2P	12	-3056.6266	6137.2532	0.0000	1.0000
SYM	16	-3053.9960	6139.9920	0.0000	1.0000
JC69G	12	-3075.1001	6174.2002	0.0000	1.0000
JC69IG	13	-3075.0574	6176.1148	0.0000	1.0000
F81	14	-3074.8249	6177.6498	0.0000	1.0000
JC69I	12	-3077.5539	6179.1078	0.0000	1.0000
JC69	11	-3106.4905	6234.9809	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2985.8544	6081.5128	0.9261	0.9261
HKYI	16	-2988.9793	6087.7628	0.0407	0.9668
HKYIG	17	-2985.7514	6088.1697	0.0332	1.0000
GTRG	20	-2983.5693	6104.3937	0.0000	1.0000
GTRI	20	-2986.5256	6110.3064	0.0000	1.0000
GTRIG	21	-2983.4883	6111.0945	0.0000	1.0000
K2PG	13	-3022.5737	6134.3632	0.0000	1.0000
K2PI	13	-3025.5354	6140.2867	0.0000	1.0000
K2PIG	14	-3022.4867	6141.0519	0.0000	1.0000
HKY	15	-3019.4258	6141.7930	0.0000	1.0000
SYMG	17	-3019.3372	6155.3412	0.0000	1.0000
SYMI	17	-3022.3620	6161.3909	0.0000	1.0000
SYMIG	18	-3019.2596	6162.0489	0.0000	1.0000
GTR	19	-3017.1232	6164.6388	0.0000	1.0000
F81G	15	-3044.0514	6191.0441	0.0000	1.0000
F81I	15	-3046.2194	6195.3801	0.0000	1.0000
K2P	12	-3056.6266	6195.6063	0.0000	1.0000
F81IG	16	-3044.0080	6197.8202	0.0000	1.0000
SYM	16	-3053.9960	6217.7961	0.0000	1.0000
JC69G	12	-3075.1001	6232.5533	0.0000	1.0000
JC69I	12	-3077.5539	6237.4609	0.0000	1.0000
JC69IG	13	-3075.0574	6239.3307	0.0000	1.0000
F81	14	-3074.8249	6245.7285	0.0000	1.0000
JC69	11	-3106.4905	6288.4713	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/chr2_22071.phylip = 1 - 956;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22071.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/chr2_22071.phylip = 1 - 956;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22071.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/chr2_22071.phylip = 1 - 956;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22071.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

