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

Input data from file "683-tomodel-phylip/chr9_6388.phylip" (7 taxa, 1095 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	-3519.9048	7071.8096	0.3378	0.3378
GTRG	20	-3516.3636	7072.7273	0.2135	0.5513
HKYI	16	-3520.7831	7073.5662	0.1404	0.6917
HKYIG	17	-3519.8417	7073.6834	0.1324	0.8241
GTRI	20	-3517.2339	7074.4678	0.0894	0.9135
GTRIG	21	-3516.2669	7074.5338	0.0865	1.0000
HKY	15	-3553.7738	7137.5476	0.0000	1.0000
GTR	19	-3550.7311	7139.4622	0.0000	1.0000
K2PG	13	-3558.7395	7143.4789	0.0000	1.0000
K2PIG	14	-3558.6838	7145.3676	0.0000	1.0000
K2PI	13	-3560.2833	7146.5666	0.0000	1.0000
SYMG	17	-3556.5804	7147.1607	0.0000	1.0000
SYMI	17	-3557.4151	7148.8303	0.0000	1.0000
SYMIG	18	-3556.5412	7149.0824	0.0000	1.0000
K2P	12	-3591.9526	7207.9052	0.0000	1.0000
SYM	16	-3589.7314	7211.4628	0.0000	1.0000
F81G	15	-3609.3909	7248.7819	0.0000	1.0000
F81I	15	-3609.7952	7249.5904	0.0000	1.0000
F81IG	16	-3609.3072	7250.6145	0.0000	1.0000
JC69G	12	-3639.8372	7303.6743	0.0000	1.0000
JC69I	12	-3640.3371	7304.6743	0.0000	1.0000
JC69IG	13	-3639.7950	7305.5901	0.0000	1.0000
F81	14	-3639.2950	7306.5900	0.0000	1.0000
JC69	11	-3669.7133	7361.4266	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3519.9048	7151.7858	0.6907	0.6907
HKYI	16	-3520.7831	7153.5423	0.2870	0.9777
HKYIG	17	-3519.8417	7158.6580	0.0222	1.0000
GTRG	20	-3516.3636	7172.6975	0.0000	1.0000
GTRI	20	-3517.2339	7174.4380	0.0000	1.0000
GTRIG	21	-3516.2669	7179.5025	0.0000	1.0000
K2PG	13	-3558.7395	7208.4596	0.0000	1.0000
K2PI	13	-3560.2833	7211.5472	0.0000	1.0000
HKY	15	-3553.7738	7212.5252	0.0000	1.0000
K2PIG	14	-3558.6838	7215.3467	0.0000	1.0000
SYMG	17	-3556.5804	7232.1354	0.0000	1.0000
SYMI	17	-3557.4151	7233.8049	0.0000	1.0000
GTR	19	-3550.7311	7234.4339	0.0000	1.0000
SYMIG	18	-3556.5412	7239.0555	0.0000	1.0000
K2P	12	-3591.9526	7267.8873	0.0000	1.0000
SYM	16	-3589.7314	7291.4390	0.0000	1.0000
F81G	15	-3609.3909	7323.7595	0.0000	1.0000
F81I	15	-3609.7952	7324.5680	0.0000	1.0000
F81IG	16	-3609.3072	7330.5906	0.0000	1.0000
JC69G	12	-3639.8372	7363.6564	0.0000	1.0000
JC69I	12	-3640.3371	7364.6564	0.0000	1.0000
JC69IG	13	-3639.7950	7370.5707	0.0000	1.0000
F81	14	-3639.2950	7376.5691	0.0000	1.0000
JC69	11	-3669.7133	7416.4102	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/chr9_6388.phylip = 1 - 1095;
 Partition Dummy = 1:683-tomodel-phylip/chr9_6388.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/chr9_6388.phylip = 1 - 1095;
 Partition Dummy = 1:683-tomodel-phylip/chr9_6388.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/chr9_6388.phylip = 1 - 1095;
 Partition Dummy = 1:683-tomodel-phylip/chr9_6388.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

