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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1657.8990	3357.1819	0.4291	0.4291
GTRG	20	-1657.9158	3357.2154	0.4220	0.8511
GTRIG	21	-1657.9071	3359.3389	0.1459	0.9970
HKYG	16	-1667.8621	3368.6145	0.0014	0.9984
HKYI	16	-1668.3156	3369.5215	0.0009	0.9993
HKYIG	17	-1667.8310	3370.6654	0.0005	0.9998
GTR	19	-1666.7115	3372.6730	0.0002	1.0000
HKY	15	-1678.9473	3388.6789	0.0000	1.0000
SYMG	17	-1677.5979	3390.1990	0.0000	1.0000
SYMI	17	-1677.6029	3390.2091	0.0000	1.0000
SYMIG	18	-1677.5921	3392.3073	0.0000	1.0000
K2PG	13	-1687.0931	3400.7791	0.0000	1.0000
K2PI	13	-1687.4561	3401.5049	0.0000	1.0000
K2PIG	14	-1687.0818	3402.8488	0.0000	1.0000
SYM	16	-1685.1678	3403.2259	0.0000	1.0000
K2P	12	-1696.5237	3417.5547	0.0000	1.0000
F81G	15	-1698.9504	3428.6851	0.0000	1.0000
F81I	15	-1699.1612	3429.1067	0.0000	1.0000
F81IG	16	-1698.9440	3430.7783	0.0000	1.0000
F81	14	-1709.1891	3447.0633	0.0000	1.0000
JC69G	12	-1716.2519	3457.0112	0.0000	1.0000
JC69I	12	-1716.4312	3457.3696	0.0000	1.0000
JC69IG	13	-1716.2495	3459.0918	0.0000	1.0000
JC69	11	-1726.1013	3474.6313	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1657.8990	3355.7980	0.4248	0.4248
GTRG	20	-1657.9158	3355.8315	0.4178	0.8426
GTRIG	21	-1657.9071	3357.8142	0.1550	0.9976
HKYG	16	-1667.8621	3367.7242	0.0011	0.9987
HKYI	16	-1668.3156	3368.6311	0.0007	0.9994
HKYIG	17	-1667.8310	3369.6621	0.0004	0.9998
GTR	19	-1666.7115	3371.4230	0.0002	1.0000
HKY	15	-1678.9473	3387.8946	0.0000	1.0000
SYMG	17	-1677.5979	3389.1957	0.0000	1.0000
SYMI	17	-1677.6029	3389.2059	0.0000	1.0000
SYMIG	18	-1677.5921	3391.1841	0.0000	1.0000
K2PG	13	-1687.0931	3400.1862	0.0000	1.0000
K2PI	13	-1687.4561	3400.9121	0.0000	1.0000
K2PIG	14	-1687.0818	3402.1636	0.0000	1.0000
SYM	16	-1685.1678	3402.3356	0.0000	1.0000
K2P	12	-1696.5237	3417.0474	0.0000	1.0000
F81G	15	-1698.9504	3427.9008	0.0000	1.0000
F81I	15	-1699.1612	3428.3223	0.0000	1.0000
F81IG	16	-1698.9440	3429.8880	0.0000	1.0000
F81	14	-1709.1891	3446.3782	0.0000	1.0000
JC69G	12	-1716.2519	3456.5039	0.0000	1.0000
JC69I	12	-1716.4312	3456.8623	0.0000	1.0000
JC69IG	13	-1716.2495	3458.4990	0.0000	1.0000
JC69	11	-1726.1013	3474.2027	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1667.8621	3438.8048	0.5599	0.5599
HKYI	16	-1668.3156	3439.7118	0.3557	0.9156
GTRI	20	-1657.8990	3444.6488	0.0301	0.9457
GTRG	20	-1657.9158	3444.6823	0.0296	0.9754
HKYIG	17	-1667.8310	3445.1853	0.0230	0.9984
GTRIG	21	-1657.9071	3451.1075	0.0012	0.9996
HKY	15	-1678.9473	3454.5327	0.0002	0.9998
GTR	19	-1666.7115	3455.8313	0.0001	0.9999
K2PG	13	-1687.0931	3457.9393	0.0000	1.0000
K2PI	13	-1687.4561	3458.6651	0.0000	1.0000
K2PIG	14	-1687.0818	3464.3592	0.0000	1.0000
SYMG	17	-1677.5979	3464.7189	0.0000	1.0000
SYMI	17	-1677.6029	3464.7290	0.0000	1.0000
K2P	12	-1696.5237	3470.3579	0.0000	1.0000
SYMIG	18	-1677.5921	3471.1498	0.0000	1.0000
SYM	16	-1685.1678	3473.4162	0.0000	1.0000
F81G	15	-1698.9504	3494.5389	0.0000	1.0000
F81I	15	-1699.1612	3494.9604	0.0000	1.0000
F81IG	16	-1698.9440	3500.9686	0.0000	1.0000
F81	14	-1709.1891	3508.5737	0.0000	1.0000
JC69G	12	-1716.2519	3509.8143	0.0000	1.0000
JC69I	12	-1716.4312	3510.1728	0.0000	1.0000
JC69IG	13	-1716.2495	3516.2520	0.0000	1.0000
JC69	11	-1726.1013	3523.0706	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_7434.phylip = 1 - 628;
 Partition Dummy = 1:591-tomodel-phylip/chr9_7434.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_7434.phylip = 1 - 628;
 Partition Dummy = 1:591-tomodel-phylip/chr9_7434.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

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

