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

Input data from file "591-tomodel-phylip/chr5_3405.phylip" (7 taxa, 1151 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2284.0288	4608.0576	0.3416	0.3416
GTRI	20	-2284.1400	4608.2800	0.3057	0.6473
GTRIG	21	-2283.9314	4609.8627	0.1385	0.7858
HKYG	16	-2289.3268	4610.6535	0.0933	0.8791
HKYI	16	-2289.4681	4610.9363	0.0810	0.9601
HKYIG	17	-2289.1795	4612.3589	0.0398	0.9999
GTR	19	-2293.3336	4624.6672	0.0001	1.0000
HKY	15	-2298.6794	4627.3589	0.0000	1.0000
SYMG	17	-2317.0569	4668.1138	0.0000	1.0000
SYMI	17	-2317.1225	4668.2451	0.0000	1.0000
SYMIG	18	-2317.0111	4670.0223	0.0000	1.0000
K2PG	13	-2323.6509	4673.3019	0.0000	1.0000
K2PI	13	-2323.7436	4673.4872	0.0000	1.0000
K2PIG	14	-2323.5586	4675.1173	0.0000	1.0000
F81G	15	-2324.4691	4678.9382	0.0000	1.0000
F81I	15	-2324.5455	4679.0910	0.0000	1.0000
F81IG	16	-2324.4187	4680.8373	0.0000	1.0000
SYM	16	-2326.4003	4684.8005	0.0000	1.0000
K2P	12	-2333.0962	4690.1924	0.0000	1.0000
F81	14	-2333.1526	4694.3052	0.0000	1.0000
JC69G	12	-2357.0040	4738.0079	0.0000	1.0000
JC69I	12	-2357.0641	4738.1282	0.0000	1.0000
JC69IG	13	-2356.9700	4739.9399	0.0000	1.0000
JC69	11	-2365.7811	4753.5622	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2289.3268	4691.4277	0.5248	0.5248
HKYI	16	-2289.4681	4691.7105	0.4556	0.9804
HKYIG	17	-2289.1795	4698.1815	0.0179	0.9983
HKY	15	-2298.6794	4703.0847	0.0015	0.9998
GTRG	20	-2284.0288	4709.0254	0.0001	0.9999
GTRI	20	-2284.1400	4709.2477	0.0001	1.0000
GTRIG	21	-2283.9314	4715.8789	0.0000	1.0000
GTR	19	-2293.3336	4720.5865	0.0000	1.0000
K2PG	13	-2323.6509	4738.9309	0.0000	1.0000
K2PI	13	-2323.7436	4739.1163	0.0000	1.0000
K2PIG	14	-2323.5586	4745.7947	0.0000	1.0000
K2P	12	-2333.0962	4750.7731	0.0000	1.0000
SYMG	17	-2317.0569	4753.9363	0.0000	1.0000
SYMI	17	-2317.1225	4754.0676	0.0000	1.0000
F81G	15	-2324.4691	4754.6640	0.0000	1.0000
F81I	15	-2324.5455	4754.8168	0.0000	1.0000
SYMIG	18	-2317.0111	4760.8932	0.0000	1.0000
F81IG	16	-2324.4187	4761.6115	0.0000	1.0000
F81	14	-2333.1526	4764.9826	0.0000	1.0000
SYM	16	-2326.4003	4765.5747	0.0000	1.0000
JC69G	12	-2357.0040	4798.5886	0.0000	1.0000
JC69I	12	-2357.0641	4798.7088	0.0000	1.0000
JC69IG	13	-2356.9700	4805.5690	0.0000	1.0000
JC69	11	-2365.7811	4809.0944	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3405.phylip = 1 - 1151;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3405.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3405.phylip = 1 - 1151;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3405.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3405.phylip = 1 - 1151;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3405.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

