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

Input data from file "683-tomodel-phylip/chr5_3290.phylip" (7 taxa, 1224 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3792.4466	7624.8933	0.6990	0.6990
GTRIG	21	-3792.3706	7626.7412	0.2775	0.9765
GTRI	20	-3795.8406	7631.6811	0.0235	1.0000
SYMG	17	-3805.9597	7645.9194	0.0000	1.0000
SYMIG	18	-3805.9458	7647.8916	0.0000	1.0000
SYMI	17	-3808.1821	7650.3641	0.0000	1.0000
HKYG	16	-3809.8756	7651.7512	0.0000	1.0000
HKYIG	17	-3809.8675	7653.7351	0.0000	1.0000
HKYI	16	-3811.8212	7655.6424	0.0000	1.0000
K2PG	13	-3818.1457	7662.2915	0.0000	1.0000
K2PIG	14	-3818.1468	7664.2936	0.0000	1.0000
K2PI	13	-3820.0718	7666.1436	0.0000	1.0000
GTR	19	-3861.5951	7761.1902	0.0000	1.0000
HKY	15	-3874.3379	7778.6758	0.0000	1.0000
SYM	16	-3874.8105	7781.6209	0.0000	1.0000
K2P	12	-3881.7174	7787.4347	0.0000	1.0000
F81G	15	-3901.1880	7832.3760	0.0000	1.0000
F81I	15	-3901.7942	7833.5883	0.0000	1.0000
F81IG	16	-3901.1880	7834.3760	0.0000	1.0000
JC69G	12	-3909.3999	7842.7998	0.0000	1.0000
JC69I	12	-3909.9708	7843.9417	0.0000	1.0000
JC69IG	13	-3909.5282	7845.0563	0.0000	1.0000
F81	14	-3959.1839	7946.3678	0.0000	1.0000
JC69	11	-3967.4536	7956.9071	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3792.4466	7727.0909	0.5890	0.5890
K2PG	13	-3818.1457	7728.7199	0.2609	0.8499
K2PI	13	-3820.0718	7732.5721	0.0380	0.8879
SYMG	17	-3805.9597	7732.7873	0.0341	0.9220
HKYG	16	-3809.8756	7733.5093	0.0238	0.9458
GTRI	20	-3795.8406	7733.8787	0.0198	0.9656
GTRIG	21	-3792.3706	7734.0487	0.0182	0.9838
K2PIG	14	-3818.1468	7735.8319	0.0074	0.9912
SYMI	17	-3808.1821	7737.2321	0.0037	0.9949
HKYI	16	-3811.8212	7737.4005	0.0034	0.9983
SYMIG	18	-3805.9458	7739.8694	0.0010	0.9993
HKYIG	17	-3809.8675	7740.6030	0.0007	1.0000
K2P	12	-3881.7174	7848.7533	0.0000	1.0000
HKY	15	-3874.3379	7855.3240	0.0000	1.0000
GTR	19	-3861.5951	7858.2779	0.0000	1.0000
SYM	16	-3874.8105	7863.3790	0.0000	1.0000
JC69G	12	-3909.3999	7904.1183	0.0000	1.0000
JC69I	12	-3909.9708	7905.2602	0.0000	1.0000
F81G	15	-3901.1880	7909.0242	0.0000	1.0000
F81I	15	-3901.7942	7910.2365	0.0000	1.0000
JC69IG	13	-3909.5282	7911.4848	0.0000	1.0000
F81IG	16	-3901.1880	7916.1341	0.0000	1.0000
JC69	11	-3967.4536	8013.1158	0.0000	1.0000
F81	14	-3959.1839	8017.9062	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 683-tomodel-phylip/chr5_3290.phylip = 1 - 1224;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3290.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 683-tomodel-phylip/chr5_3290.phylip = 1 - 1224;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3290.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 (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_3290.phylip = 1 - 1224;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3290.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;


-------------------------------------------------------------
End of Output

