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

Input data from file "591-tomodel-phylip/chr3_17832.phylip" (7 taxa, 1147 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	-2672.7432	5377.4864	0.2704	0.2704
HKYI	16	-2672.8649	5377.7299	0.2394	0.5098
HKYIG	17	-2672.3553	5378.7106	0.1466	0.6563
GTRG	20	-2669.3974	5378.7947	0.1406	0.7969
GTRI	20	-2669.4878	5378.9756	0.1284	0.9253
GTRIG	21	-2669.0296	5380.0593	0.0747	1.0000
HKY	15	-2706.1050	5442.2100	0.0000	1.0000
F81G	15	-2707.1072	5444.2144	0.0000	1.0000
F81I	15	-2707.1641	5444.3282	0.0000	1.0000
GTR	19	-2703.2411	5444.4823	0.0000	1.0000
F81IG	16	-2707.0263	5446.0526	0.0000	1.0000
K2PG	13	-2717.5596	5461.1192	0.0000	1.0000
K2PI	13	-2717.6003	5461.2007	0.0000	1.0000
SYMG	17	-2713.8770	5461.7540	0.0000	1.0000
SYMI	17	-2713.8865	5461.7731	0.0000	1.0000
K2PIG	14	-2717.3784	5462.7569	0.0000	1.0000
SYMIG	18	-2713.8841	5463.7681	0.0000	1.0000
F81	14	-2739.6197	5507.2394	0.0000	1.0000
JC69G	12	-2748.2564	5520.5127	0.0000	1.0000
JC69I	12	-2748.2865	5520.5731	0.0000	1.0000
JC69IG	13	-2748.2342	5522.4685	0.0000	1.0000
K2P	12	-2750.5530	5525.1060	0.0000	1.0000
SYM	16	-2746.9373	5525.8746	0.0000	1.0000
JC69	11	-2780.5202	5583.0405	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2672.7432	5458.2048	0.5184	0.5184
HKYI	16	-2672.8649	5458.4483	0.4590	0.9774
HKYIG	17	-2672.3553	5464.4740	0.0226	1.0000
GTRG	20	-2669.3974	5479.6928	0.0000	1.0000
GTRI	20	-2669.4878	5479.8737	0.0000	1.0000
GTRIG	21	-2669.0296	5486.0023	0.0000	1.0000
HKY	15	-2706.1050	5517.8836	0.0000	1.0000
F81G	15	-2707.1072	5519.8880	0.0000	1.0000
F81I	15	-2707.1641	5520.0017	0.0000	1.0000
K2PG	13	-2717.5596	5526.7030	0.0000	1.0000
F81IG	16	-2707.0263	5526.7711	0.0000	1.0000
K2PI	13	-2717.6003	5526.7844	0.0000	1.0000
K2PIG	14	-2717.3784	5533.3856	0.0000	1.0000
GTR	19	-2703.2411	5540.3355	0.0000	1.0000
SYMG	17	-2713.8770	5547.5174	0.0000	1.0000
SYMI	17	-2713.8865	5547.5364	0.0000	1.0000
SYMIG	18	-2713.8841	5554.5764	0.0000	1.0000
F81	14	-2739.6197	5577.8681	0.0000	1.0000
JC69G	12	-2748.2564	5581.0516	0.0000	1.0000
JC69I	12	-2748.2865	5581.1119	0.0000	1.0000
K2P	12	-2750.5530	5585.6449	0.0000	1.0000
JC69IG	13	-2748.2342	5588.0523	0.0000	1.0000
SYM	16	-2746.9373	5606.5931	0.0000	1.0000
JC69	11	-2780.5202	5638.5344	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/chr3_17832.phylip = 1 - 1147;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17832.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/chr3_17832.phylip = 1 - 1147;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17832.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/chr3_17832.phylip = 1 - 1147;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17832.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

