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

Input data from file "591-tomodel-phylip/chr12_5852.phylip" (7 taxa, 984 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3410.1731	6860.3461	0.3659	0.3659
GTRI	20	-3410.4437	6860.8874	0.2792	0.6451
GTRIG	21	-3409.7082	6861.4165	0.2143	0.8594
HKYG	16	-3415.9191	6863.8382	0.0638	0.9232
HKYI	16	-3416.2996	6864.5993	0.0436	0.9668
HKYIG	17	-3415.5741	6865.1481	0.0332	1.0000
GTR	19	-3459.3750	6956.7501	0.0000	1.0000
HKY	15	-3465.8643	6961.7285	0.0000	1.0000
F81I	15	-3471.9654	6973.9308	0.0000	1.0000
F81G	15	-3472.0854	6974.1707	0.0000	1.0000
F81IG	16	-3471.6139	6975.2279	0.0000	1.0000
K2PG	13	-3476.8624	6979.7249	0.0000	1.0000
K2PIG	14	-3476.8251	6981.6501	0.0000	1.0000
K2PI	13	-3478.5158	6983.0316	0.0000	1.0000
SYMG	17	-3475.2297	6984.4594	0.0000	1.0000
SYMIG	18	-3475.1687	6986.3374	0.0000	1.0000
SYMI	17	-3477.1404	6988.2807	0.0000	1.0000
F81	14	-3518.3012	7064.6025	0.0000	1.0000
JC69G	12	-3521.9646	7067.9293	0.0000	1.0000
JC69I	12	-3522.7980	7069.5961	0.0000	1.0000
JC69IG	13	-3521.8515	7069.7030	0.0000	1.0000
K2P	12	-3529.4565	7082.9130	0.0000	1.0000
SYM	16	-3528.0081	7088.0162	0.0000	1.0000
JC69	11	-3570.7888	7163.5777	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3415.9191	6942.1042	0.5783	0.5783
HKYI	16	-3416.2996	6942.8653	0.3953	0.9736
HKYIG	17	-3415.5741	6948.3057	0.0260	0.9997
GTRG	20	-3410.1731	6958.1786	0.0002	0.9998
GTRI	20	-3410.4437	6958.7199	0.0001	1.0000
GTRIG	21	-3409.7082	6964.1406	0.0000	1.0000
HKY	15	-3465.8643	7035.1029	0.0000	1.0000
K2PG	13	-3476.8624	7043.3160	0.0000	1.0000
K2PI	13	-3478.5158	7046.6228	0.0000	1.0000
F81I	15	-3471.9654	7047.3052	0.0000	1.0000
F81G	15	-3472.0854	7047.5451	0.0000	1.0000
GTR	19	-3459.3750	7049.6910	0.0000	1.0000
K2PIG	14	-3476.8251	7050.1329	0.0000	1.0000
F81IG	16	-3471.6139	7053.4939	0.0000	1.0000
SYMG	17	-3475.2297	7067.6170	0.0000	1.0000
SYMI	17	-3477.1404	7071.4384	0.0000	1.0000
SYMIG	18	-3475.1687	7074.3866	0.0000	1.0000
JC69G	12	-3521.9646	7126.6288	0.0000	1.0000
JC69I	12	-3522.7980	7128.2956	0.0000	1.0000
F81	14	-3518.3012	7133.0852	0.0000	1.0000
JC69IG	13	-3521.8515	7133.2941	0.0000	1.0000
K2P	12	-3529.4565	7141.6125	0.0000	1.0000
SYM	16	-3528.0081	7166.2822	0.0000	1.0000
JC69	11	-3570.7888	7217.3856	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 591-tomodel-phylip/chr12_5852.phylip = 1 - 984;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5852.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 591-tomodel-phylip/chr12_5852.phylip = 1 - 984;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5852.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr12_5852.phylip = 1 - 984;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5852.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

