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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2290.0734	4621.5352	0.5183	0.5183
GTRI	20	-2290.6890	4622.7665	0.2800	0.7983
GTRIG	21	-2290.0783	4623.6863	0.1768	0.9751
GTR	19	-2294.1956	4627.6453	0.0244	0.9996
HKYG	16	-2302.0407	4636.9746	0.0002	0.9998
HKYI	16	-2302.6543	4638.2019	0.0001	0.9999
HKYIG	17	-2302.0407	4639.0880	0.0001	1.0000
HKY	15	-2306.0904	4642.9678	0.0000	1.0000
SYMG	17	-2332.1684	4699.3433	0.0000	1.0000
K2PG	13	-2336.4106	4699.4159	0.0000	1.0000
SYMI	17	-2333.1351	4701.2769	0.0000	1.0000
K2PI	13	-2337.3550	4701.3048	0.0000	1.0000
SYMIG	18	-2332.1691	4701.4650	0.0000	1.0000
K2PIG	14	-2336.4106	4701.5086	0.0000	1.0000
SYM	16	-2336.7839	4706.4610	0.0000	1.0000
K2P	12	-2341.3458	4707.2006	0.0000	1.0000
F81G	15	-2346.7832	4724.3532	0.0000	1.0000
F81I	15	-2347.4618	4725.7104	0.0000	1.0000
F81IG	16	-2346.7832	4726.4597	0.0000	1.0000
F81	14	-2350.7456	4730.1786	0.0000	1.0000
JC69G	12	-2373.9619	4772.4327	0.0000	1.0000
JC69I	12	-2374.8470	4774.2030	0.0000	1.0000
JC69IG	13	-2373.9619	4774.5187	0.0000	1.0000
JC69	11	-2378.5165	4779.4629	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2290.0734	4620.1468	0.5125	0.5125
GTRI	20	-2290.6890	4621.3781	0.2769	0.7894
GTRIG	21	-2290.0783	4622.1565	0.1876	0.9771
GTR	19	-2294.1956	4626.3912	0.0226	0.9997
HKYG	16	-2302.0407	4636.0814	0.0002	0.9998
HKYI	16	-2302.6543	4637.3087	0.0001	0.9999
HKYIG	17	-2302.0407	4638.0814	0.0001	1.0000
HKY	15	-2306.0904	4642.1809	0.0000	1.0000
SYMG	17	-2332.1684	4698.3367	0.0000	1.0000
K2PG	13	-2336.4106	4698.8211	0.0000	1.0000
SYMI	17	-2333.1351	4700.2703	0.0000	1.0000
SYMIG	18	-2332.1691	4700.3381	0.0000	1.0000
K2PI	13	-2337.3550	4700.7101	0.0000	1.0000
K2PIG	14	-2336.4106	4700.8212	0.0000	1.0000
SYM	16	-2336.7839	4705.5678	0.0000	1.0000
K2P	12	-2341.3458	4706.6916	0.0000	1.0000
F81G	15	-2346.7832	4723.5663	0.0000	1.0000
F81I	15	-2347.4618	4724.9235	0.0000	1.0000
F81IG	16	-2346.7832	4725.5665	0.0000	1.0000
F81	14	-2350.7456	4729.4912	0.0000	1.0000
JC69G	12	-2373.9619	4771.9237	0.0000	1.0000
JC69I	12	-2374.8470	4773.6940	0.0000	1.0000
JC69IG	13	-2373.9619	4773.9239	0.0000	1.0000
JC69	11	-2378.5165	4779.0329	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2302.0407	4707.1110	0.3552	0.3552
HKYI	16	-2302.6543	4708.3383	0.1923	0.5475
HKY	15	-2306.0904	4708.7711	0.1549	0.7023
GTRG	20	-2290.0734	4708.9338	0.1428	0.8451
GTRI	20	-2290.6890	4710.1651	0.0771	0.9222
GTR	19	-2294.1956	4710.7388	0.0579	0.9801
HKYIG	17	-2302.0407	4713.5504	0.0142	0.9943
GTRIG	21	-2290.0783	4715.3829	0.0057	1.0000
K2PG	13	-2336.4106	4756.5327	0.0000	1.0000
K2PI	13	-2337.3550	4758.4216	0.0000	1.0000
K2P	12	-2341.3458	4759.9638	0.0000	1.0000
K2PIG	14	-2336.4106	4762.9721	0.0000	1.0000
SYMG	17	-2332.1684	4773.8057	0.0000	1.0000
SYMI	17	-2333.1351	4775.7392	0.0000	1.0000
SYM	16	-2336.7839	4776.5974	0.0000	1.0000
SYMIG	18	-2332.1691	4780.2464	0.0000	1.0000
F81G	15	-2346.7832	4790.1566	0.0000	1.0000
F81I	15	-2347.4618	4791.5138	0.0000	1.0000
F81	14	-2350.7456	4791.6421	0.0000	1.0000
F81IG	16	-2346.7832	4796.5961	0.0000	1.0000
JC69G	12	-2373.9619	4825.1959	0.0000	1.0000
JC69I	12	-2374.8470	4826.9662	0.0000	1.0000
JC69	11	-2378.5165	4827.8658	0.0000	1.0000
JC69IG	13	-2373.9619	4831.6354	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/chr4_6508.phylip = 1 - 626;
 Partition Dummy = 1:591-tomodel-phylip/chr4_6508.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/chr4_6508.phylip = 1 - 626;
 Partition Dummy = 1:591-tomodel-phylip/chr4_6508.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/chr4_6508.phylip = 1 - 626;
 Partition Dummy = 1:591-tomodel-phylip/chr4_6508.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

