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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3743.1973	7526.3946	0.4608	0.4608
GTRI	20	-3743.3841	7526.7682	0.3823	0.8430
GTRIG	21	-3743.2786	7528.5572	0.1563	0.9993
HKYG	16	-3754.6050	7541.2101	0.0003	0.9996
HKYI	16	-3754.7764	7541.5529	0.0002	0.9998
HKYIG	17	-3754.4099	7542.8198	0.0001	0.9999
SYMG	17	-3756.0030	7546.0061	0.0000	1.0000
SYMI	17	-3756.0838	7546.1676	0.0000	1.0000
SYMIG	18	-3755.7742	7547.5483	0.0000	1.0000
K2PG	13	-3763.7631	7553.5263	0.0000	1.0000
K2PI	13	-3763.9707	7553.9414	0.0000	1.0000
K2PIG	14	-3763.5678	7555.1356	0.0000	1.0000
GTR	19	-3790.3857	7618.7714	0.0000	1.0000
HKY	15	-3800.3804	7630.7607	0.0000	1.0000
SYM	16	-3801.2535	7634.5071	0.0000	1.0000
K2P	12	-3808.8094	7641.6188	0.0000	1.0000
F81I	15	-3835.0687	7700.1374	0.0000	1.0000
F81G	15	-3835.4313	7700.8625	0.0000	1.0000
F81IG	16	-3834.9560	7701.9120	0.0000	1.0000
JC69I	12	-3842.1749	7708.3498	0.0000	1.0000
JC69G	12	-3842.5279	7709.0558	0.0000	1.0000
JC69IG	13	-3842.0549	7710.1099	0.0000	1.0000
F81	14	-3877.6029	7783.2057	0.0000	1.0000
JC69	11	-3884.5717	7791.1434	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3763.7631	7619.2341	0.4284	0.4284
K2PI	13	-3763.9707	7619.6492	0.3481	0.7766
HKYG	16	-3754.6050	7622.0813	0.1032	0.8797
HKYI	16	-3754.7764	7622.4241	0.0869	0.9667
K2PIG	14	-3763.5678	7625.8979	0.0153	0.9820
GTRG	20	-3743.1973	7627.4836	0.0069	0.9889
GTRI	20	-3743.3841	7627.8572	0.0057	0.9947
HKYIG	17	-3754.4099	7628.7454	0.0037	0.9983
SYMG	17	-3756.0030	7631.9317	0.0007	0.9991
SYMI	17	-3756.0838	7632.0932	0.0007	0.9998
GTRIG	21	-3743.2786	7634.7006	0.0002	1.0000
SYMIG	18	-3755.7742	7638.5284	0.0000	1.0000
K2P	12	-3808.8094	7702.2722	0.0000	1.0000
HKY	15	-3800.3804	7706.5775	0.0000	1.0000
GTR	19	-3790.3857	7714.8060	0.0000	1.0000
SYM	16	-3801.2535	7715.3783	0.0000	1.0000
JC69I	12	-3842.1749	7769.0032	0.0000	1.0000
JC69G	12	-3842.5279	7769.7092	0.0000	1.0000
JC69IG	13	-3842.0549	7775.8177	0.0000	1.0000
F81I	15	-3835.0687	7775.9541	0.0000	1.0000
F81G	15	-3835.4313	7776.6793	0.0000	1.0000
F81IG	16	-3834.9560	7782.7832	0.0000	1.0000
JC69	11	-3884.5717	7846.7424	0.0000	1.0000
F81	14	-3877.6029	7853.9680	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_11359.phylip = 1 - 1158;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11359.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_11359.phylip = 1 - 1158;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11359.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_11359.phylip = 1 - 1158;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11359.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

