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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-1714.0917	3458.1835	0.2312	0.2312
HKYI	16	-1713.3320	3458.6640	0.1818	0.4131
HKYG	16	-1713.3819	3458.7638	0.1730	0.5861
GTR	19	-1710.7198	3459.4396	0.1234	0.7095
GTRI	20	-1709.9527	3459.9053	0.0978	0.8072
GTRG	20	-1710.0053	3460.0107	0.0927	0.9000
HKYIG	17	-1713.3798	3460.7597	0.0638	0.9637
GTRIG	21	-1709.9522	3461.9044	0.0360	0.9997
F81	14	-1722.8490	3473.6981	0.0001	0.9998
F81I	15	-1722.1415	3474.2829	0.0001	0.9999
F81G	15	-1722.1787	3474.3573	0.0001	1.0000
F81IG	16	-1722.1764	3476.3528	0.0000	1.0000
SYMI	17	-1764.6630	3563.3260	0.0000	1.0000
SYM	16	-1765.6935	3563.3870	0.0000	1.0000
SYMG	17	-1764.7028	3563.4057	0.0000	1.0000
SYMIG	18	-1764.7058	3565.4116	0.0000	1.0000
K2P	12	-1771.5096	3567.0192	0.0000	1.0000
K2PI	13	-1770.6186	3567.2371	0.0000	1.0000
K2PG	13	-1770.6561	3567.3123	0.0000	1.0000
K2PIG	14	-1770.6570	3569.3141	0.0000	1.0000
JC69	11	-1779.6701	3581.3401	0.0000	1.0000
JC69I	12	-1778.9025	3581.8050	0.0000	1.0000
JC69G	12	-1778.9338	3581.8677	0.0000	1.0000
JC69IG	13	-1778.9330	3583.8659	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1714.0917	3530.8397	0.8742	0.8742
HKYI	16	-1713.3320	3536.1640	0.0610	0.9352
HKYG	16	-1713.3819	3536.2638	0.0581	0.9933
F81	14	-1722.8490	3541.5106	0.0042	0.9975
HKYIG	17	-1713.3798	3543.1034	0.0019	0.9994
F81I	15	-1722.1415	3546.9391	0.0003	0.9997
F81G	15	-1722.1787	3547.0135	0.0003	1.0000
GTR	19	-1710.7198	3551.4708	0.0000	1.0000
F81IG	16	-1722.1764	3553.8528	0.0000	1.0000
GTRI	20	-1709.9527	3556.7803	0.0000	1.0000
GTRG	20	-1710.0053	3556.8857	0.0000	1.0000
GTRIG	21	-1709.9522	3563.6231	0.0000	1.0000
K2P	12	-1771.5096	3625.1442	0.0000	1.0000
K2PI	13	-1770.6186	3630.2059	0.0000	1.0000
K2PG	13	-1770.6561	3630.2810	0.0000	1.0000
JC69	11	-1779.6701	3634.6214	0.0000	1.0000
K2PIG	14	-1770.6570	3637.1266	0.0000	1.0000
JC69I	12	-1778.9025	3639.9300	0.0000	1.0000
JC69G	12	-1778.9338	3639.9927	0.0000	1.0000
SYM	16	-1765.6935	3640.8870	0.0000	1.0000
SYMI	17	-1764.6630	3645.6697	0.0000	1.0000
SYMG	17	-1764.7028	3645.7494	0.0000	1.0000
JC69IG	13	-1778.9330	3646.8347	0.0000	1.0000
SYMIG	18	-1764.7058	3652.5991	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10633.phylip = 1 - 938;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10633.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10633.phylip = 1 - 938;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10633.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10633.phylip = 1 - 938;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10633.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

