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

Input data from file "591-tomodel-phylip/chr3_5877.phylip" (7 taxa, 910 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	-2700.0180	5432.0359	0.4388	0.4388
GTRG	20	-2696.7644	5433.5287	0.2080	0.6468
HKYIG	17	-2699.9460	5433.8921	0.1735	0.8203
GTRIG	21	-2696.7000	5435.4000	0.0816	0.9019
HKYI	16	-2701.9228	5435.8455	0.0653	0.9672
GTRI	20	-2698.6112	5437.2224	0.0328	1.0000
K2PG	13	-2721.8149	5469.6298	0.0000	1.0000
K2PIG	14	-2721.7280	5471.4561	0.0000	1.0000
K2PI	13	-2723.6520	5473.3040	0.0000	1.0000
SYMG	17	-2720.3786	5474.7573	0.0000	1.0000
SYMIG	18	-2720.3037	5476.6073	0.0000	1.0000
SYMI	17	-2722.2035	5478.4071	0.0000	1.0000
HKY	15	-2725.4795	5480.9590	0.0000	1.0000
GTR	19	-2721.8804	5481.7609	0.0000	1.0000
K2P	12	-2748.2150	5520.4300	0.0000	1.0000
SYM	16	-2745.9389	5523.8778	0.0000	1.0000
F81G	15	-2755.5720	5541.1440	0.0000	1.0000
F81IG	16	-2755.5392	5543.0783	0.0000	1.0000
F81I	15	-2757.1268	5544.2536	0.0000	1.0000
JC69G	12	-2774.7754	5573.5508	0.0000	1.0000
JC69IG	13	-2774.7446	5575.4893	0.0000	1.0000
JC69I	12	-2776.2947	5576.5895	0.0000	1.0000
F81	14	-2778.5523	5585.1046	0.0000	1.0000
JC69	11	-2797.8720	5617.7440	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2700.0180	5509.0511	0.8442	0.8442
HKYI	16	-2701.9228	5512.8607	0.1257	0.9699
HKYIG	17	-2699.9460	5515.7206	0.0301	1.0000
GTRG	20	-2696.7644	5529.7976	0.0000	1.0000
K2PG	13	-2721.8149	5532.2045	0.0000	1.0000
GTRI	20	-2698.6112	5533.4913	0.0000	1.0000
K2PI	13	-2723.6520	5535.8788	0.0000	1.0000
GTRIG	21	-2696.7000	5536.4824	0.0000	1.0000
K2PIG	14	-2721.7280	5538.8443	0.0000	1.0000
HKY	15	-2725.4795	5553.1606	0.0000	1.0000
SYMG	17	-2720.3786	5556.5858	0.0000	1.0000
SYMI	17	-2722.2035	5560.2356	0.0000	1.0000
SYMIG	18	-2720.3037	5563.2493	0.0000	1.0000
GTR	19	-2721.8804	5573.2163	0.0000	1.0000
K2P	12	-2748.2150	5578.1914	0.0000	1.0000
SYM	16	-2745.9389	5600.8930	0.0000	1.0000
F81G	15	-2755.5720	5613.3456	0.0000	1.0000
F81I	15	-2757.1268	5616.4553	0.0000	1.0000
F81IG	16	-2755.5392	5620.0935	0.0000	1.0000
JC69G	12	-2774.7754	5631.3121	0.0000	1.0000
JC69I	12	-2776.2947	5634.3508	0.0000	1.0000
JC69IG	13	-2774.7446	5638.0640	0.0000	1.0000
F81	14	-2778.5523	5652.4928	0.0000	1.0000
JC69	11	-2797.8720	5670.6919	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_5877.phylip = 1 - 910;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5877.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_5877.phylip = 1 - 910;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5877.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_5877.phylip = 1 - 910;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5877.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

