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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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	-1726.3292	3493.9567	0.6285	0.6285
GTRG	20	-1727.2873	3495.8728	0.2411	0.8697
GTRIG	21	-1727.2874	3498.0051	0.0830	0.9527
HKYI	16	-1733.7545	3500.3447	0.0258	0.9785
HKYG	16	-1734.2340	3501.3036	0.0160	0.9944
HKYIG	17	-1734.2340	3503.4096	0.0056	1.0000
GTR	19	-1746.8720	3532.9169	0.0000	1.0000
HKY	15	-1752.4606	3535.6574	0.0000	1.0000
SYMI	17	-1761.7034	3558.3484	0.0000	1.0000
K2PI	13	-1765.9744	3558.5055	0.0000	1.0000
K2PG	13	-1766.1834	3558.9235	0.0000	1.0000
SYMG	17	-1762.2527	3559.4469	0.0000	1.0000
F81I	15	-1764.4321	3559.6003	0.0000	1.0000
F81G	15	-1764.9075	3560.5512	0.0000	1.0000
K2PIG	14	-1766.1834	3561.0100	0.0000	1.0000
SYMIG	18	-1762.2527	3561.5593	0.0000	1.0000
F81IG	16	-1764.9074	3562.6505	0.0000	1.0000
K2P	12	-1783.1592	3590.7947	0.0000	1.0000
F81	14	-1782.2135	3593.0703	0.0000	1.0000
SYM	16	-1780.4361	3593.7079	0.0000	1.0000
JC69I	12	-1791.4600	3607.3963	0.0000	1.0000
JC69G	12	-1791.8490	3608.1744	0.0000	1.0000
JC69IG	13	-1791.8490	3610.2545	0.0000	1.0000
JC69	11	-1808.5134	3639.4292	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1726.3292	3492.6584	0.6310	0.6310
GTRG	20	-1727.2873	3494.5745	0.2421	0.8730
GTRIG	21	-1727.2874	3496.5748	0.0890	0.9621
HKYI	16	-1733.7545	3499.5090	0.0205	0.9826
HKYG	16	-1734.2340	3500.4680	0.0127	0.9953
HKYIG	17	-1734.2340	3502.4680	0.0047	1.0000
GTR	19	-1746.8720	3531.7441	0.0000	1.0000
HKY	15	-1752.4606	3534.9212	0.0000	1.0000
SYMI	17	-1761.7034	3557.4068	0.0000	1.0000
K2PI	13	-1765.9744	3557.9489	0.0000	1.0000
K2PG	13	-1766.1834	3558.3669	0.0000	1.0000
SYMG	17	-1762.2527	3558.5053	0.0000	1.0000
F81I	15	-1764.4321	3558.8641	0.0000	1.0000
F81G	15	-1764.9075	3559.8150	0.0000	1.0000
K2PIG	14	-1766.1834	3560.3668	0.0000	1.0000
SYMIG	18	-1762.2527	3560.5053	0.0000	1.0000
F81IG	16	-1764.9074	3561.8148	0.0000	1.0000
K2P	12	-1783.1592	3590.3184	0.0000	1.0000
F81	14	-1782.2135	3592.4271	0.0000	1.0000
SYM	16	-1780.4361	3592.8723	0.0000	1.0000
JC69I	12	-1791.4600	3606.9199	0.0000	1.0000
JC69G	12	-1791.8490	3607.6980	0.0000	1.0000
JC69IG	13	-1791.8490	3609.6980	0.0000	1.0000
JC69	11	-1808.5134	3639.0268	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1733.7545	3571.5777	0.6067	0.6067
HKYG	16	-1734.2340	3572.5366	0.3756	0.9823
HKYIG	17	-1734.2340	3579.0409	0.0145	0.9968
GTRI	20	-1726.3292	3582.7442	0.0023	0.9991
GTRG	20	-1727.2873	3584.6603	0.0009	1.0000
GTRIG	21	-1727.2874	3591.1649	0.0000	1.0000
HKY	15	-1752.4606	3602.4855	0.0000	1.0000
K2PI	13	-1765.9744	3616.5046	0.0000	1.0000
K2PG	13	-1766.1834	3616.9226	0.0000	1.0000
GTR	19	-1746.8720	3617.3256	0.0000	1.0000
K2PIG	14	-1766.1834	3623.4268	0.0000	1.0000
F81I	15	-1764.4321	3626.4285	0.0000	1.0000
F81G	15	-1764.9075	3627.3793	0.0000	1.0000
F81IG	16	-1764.9074	3633.8835	0.0000	1.0000
SYMI	17	-1761.7034	3633.9797	0.0000	1.0000
SYMG	17	-1762.2527	3635.0782	0.0000	1.0000
SYMIG	18	-1762.2527	3641.5825	0.0000	1.0000
K2P	12	-1783.1592	3644.3699	0.0000	1.0000
F81	14	-1782.2135	3655.4871	0.0000	1.0000
JC69I	12	-1791.4600	3660.9714	0.0000	1.0000
JC69G	12	-1791.8490	3661.7495	0.0000	1.0000
SYM	16	-1780.4361	3664.9409	0.0000	1.0000
JC69IG	13	-1791.8490	3668.2537	0.0000	1.0000
JC69	11	-1808.5134	3688.5739	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 683-tomodel-phylip/chr5_1757.phylip = 1 - 668;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1757.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 683-tomodel-phylip/chr5_1757.phylip = 1 - 668;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1757.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 (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_1757.phylip = 1 - 668;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1757.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

