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

Input data from file "591-tomodel-phylip/chr3_5470.phylip" (7 taxa, 657 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	-1722.9526	3487.2260	0.4584	0.4584
GTRG	20	-1723.0795	3487.4798	0.4038	0.8622
GTRIG	21	-1723.1136	3489.6824	0.1342	0.9965
HKYI	16	-1732.8816	3498.6132	0.0015	0.9980
HKYG	16	-1732.9289	3498.7077	0.0015	0.9995
HKYIG	17	-1732.9550	3500.8677	0.0005	1.0000
GTR	19	-1734.1855	3507.5640	0.0000	1.0000
HKY	15	-1744.2336	3519.2160	0.0000	1.0000
F81G	15	-1760.1448	3551.0383	0.0000	1.0000
F81I	15	-1760.1470	3551.0429	0.0000	1.0000
F81IG	16	-1760.1623	3553.1746	0.0000	1.0000
SYMI	17	-1765.7833	3566.5243	0.0000	1.0000
SYMG	17	-1765.8135	3566.5848	0.0000	1.0000
SYMIG	18	-1765.8390	3568.7502	0.0000	1.0000
F81	14	-1770.4341	3569.5224	0.0000	1.0000
K2PG	13	-1771.8790	3570.3241	0.0000	1.0000
K2PI	13	-1771.9174	3570.4008	0.0000	1.0000
K2PIG	14	-1771.8974	3572.4489	0.0000	1.0000
SYM	16	-1776.9575	3586.7651	0.0000	1.0000
K2P	12	-1783.0779	3590.6404	0.0000	1.0000
JC69G	12	-1796.4356	3617.3557	0.0000	1.0000
JC69I	12	-1796.4899	3617.4643	0.0000	1.0000
JC69IG	13	-1796.4504	3619.4668	0.0000	1.0000
JC69	11	-1807.0866	3636.5826	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1722.9526	3485.9052	0.4545	0.4545
GTRG	20	-1723.0795	3486.1590	0.4003	0.8549
GTRIG	21	-1723.1136	3488.2273	0.1423	0.9972
HKYI	16	-1732.8816	3497.7632	0.0012	0.9984
HKYG	16	-1732.9289	3497.8577	0.0012	0.9996
HKYIG	17	-1732.9550	3499.9099	0.0004	1.0000
GTR	19	-1734.1855	3506.3709	0.0000	1.0000
HKY	15	-1744.2336	3518.4671	0.0000	1.0000
F81G	15	-1760.1448	3550.2895	0.0000	1.0000
F81I	15	-1760.1470	3550.2941	0.0000	1.0000
F81IG	16	-1760.1623	3552.3246	0.0000	1.0000
SYMI	17	-1765.7833	3565.5666	0.0000	1.0000
SYMG	17	-1765.8135	3565.6270	0.0000	1.0000
SYMIG	18	-1765.8390	3567.6781	0.0000	1.0000
F81	14	-1770.4341	3568.8682	0.0000	1.0000
K2PG	13	-1771.8790	3569.7580	0.0000	1.0000
K2PI	13	-1771.9174	3569.8347	0.0000	1.0000
K2PIG	14	-1771.8974	3571.7947	0.0000	1.0000
SYM	16	-1776.9575	3585.9151	0.0000	1.0000
K2P	12	-1783.0779	3590.1559	0.0000	1.0000
JC69G	12	-1796.4356	3616.8712	0.0000	1.0000
JC69I	12	-1796.4899	3616.9798	0.0000	1.0000
JC69IG	13	-1796.4504	3618.9007	0.0000	1.0000
JC69	11	-1807.0866	3636.1733	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1732.8816	3569.5662	0.4804	0.4804
HKYG	16	-1732.9289	3569.6607	0.4583	0.9387
GTRI	20	-1722.9526	3575.6589	0.0228	0.9616
GTRG	20	-1723.0795	3575.9127	0.0201	0.9817
HKYIG	17	-1732.9550	3576.2005	0.0174	0.9991
GTRIG	21	-1723.1136	3582.4686	0.0008	0.9998
HKY	15	-1744.2336	3585.7824	0.0001	1.0000
GTR	19	-1734.1855	3591.6369	0.0000	1.0000
F81G	15	-1760.1448	3617.6048	0.0000	1.0000
F81I	15	-1760.1470	3617.6093	0.0000	1.0000
F81IG	16	-1760.1623	3624.1275	0.0000	1.0000
K2PG	13	-1771.8790	3628.0979	0.0000	1.0000
K2PI	13	-1771.9174	3628.1746	0.0000	1.0000
F81	14	-1770.4341	3631.6958	0.0000	1.0000
K2PIG	14	-1771.8974	3634.6223	0.0000	1.0000
SYMI	17	-1765.7833	3641.8572	0.0000	1.0000
SYMG	17	-1765.8135	3641.9177	0.0000	1.0000
K2P	12	-1783.0779	3644.0081	0.0000	1.0000
SYMIG	18	-1765.8390	3648.4564	0.0000	1.0000
SYM	16	-1776.9575	3657.7180	0.0000	1.0000
JC69G	12	-1796.4356	3670.7234	0.0000	1.0000
JC69I	12	-1796.4899	3670.8320	0.0000	1.0000
JC69IG	13	-1796.4504	3677.2406	0.0000	1.0000
JC69	11	-1807.0866	3685.5378	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/chr3_5470.phylip = 1 - 657;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5470.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/chr3_5470.phylip = 1 - 657;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5470.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 591-tomodel-phylip/chr3_5470.phylip = 1 - 657;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5470.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

