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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1755.8445	3544.5706	0.4953	0.4953
HKYI	16	-1756.3553	3545.5922	0.2972	0.7925
HKYIG	17	-1755.8100	3546.6136	0.1783	0.9709
GTRG	20	-1755.0119	3551.3941	0.0163	0.9872
GTRI	20	-1755.8991	3553.1685	0.0067	0.9939
GTRIG	21	-1754.9325	3553.3748	0.0061	1.0000
F81G	15	-1770.6394	3572.0555	0.0000	1.0000
HKY	15	-1770.8422	3572.4611	0.0000	1.0000
F81I	15	-1771.3089	3573.3944	0.0000	1.0000
F81IG	16	-1770.5992	3574.0800	0.0000	1.0000
K2PG	13	-1775.4932	3577.5735	0.0000	1.0000
GTR	19	-1769.4516	3578.1410	0.0000	1.0000
K2PI	13	-1776.2010	3578.9892	0.0000	1.0000
K2PIG	14	-1775.4461	3579.5708	0.0000	1.0000
SYMG	17	-1774.9553	3584.9040	0.0000	1.0000
SYMI	17	-1775.6261	3586.2457	0.0000	1.0000
SYMIG	18	-1774.9065	3586.9252	0.0000	1.0000
F81	14	-1784.2267	3597.1319	0.0000	1.0000
JC69G	12	-1788.0695	3600.6415	0.0000	1.0000
JC69I	12	-1788.6769	3601.8561	0.0000	1.0000
JC69IG	13	-1788.0387	3602.6645	0.0000	1.0000
K2P	12	-1789.1320	3602.7664	0.0000	1.0000
SYM	16	-1788.6759	3610.2335	0.0000	1.0000
JC69	11	-1801.3391	3625.1027	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1755.8445	3543.6889	0.4861	0.4861
HKYI	16	-1756.3553	3544.7106	0.2917	0.7778
HKYIG	17	-1755.8100	3545.6201	0.1851	0.9629
GTRG	20	-1755.0119	3550.0238	0.0205	0.9834
GTRI	20	-1755.8991	3551.7982	0.0084	0.9918
GTRIG	21	-1754.9325	3551.8650	0.0082	1.0000
F81G	15	-1770.6394	3571.2788	0.0000	1.0000
HKY	15	-1770.8422	3571.6844	0.0000	1.0000
F81I	15	-1771.3089	3572.6177	0.0000	1.0000
F81IG	16	-1770.5992	3573.1983	0.0000	1.0000
GTR	19	-1769.4516	3576.9032	0.0000	1.0000
K2PG	13	-1775.4932	3576.9864	0.0000	1.0000
K2PI	13	-1776.2010	3578.4021	0.0000	1.0000
K2PIG	14	-1775.4461	3578.8923	0.0000	1.0000
SYMG	17	-1774.9553	3583.9105	0.0000	1.0000
SYMI	17	-1775.6261	3585.2522	0.0000	1.0000
SYMIG	18	-1774.9065	3585.8130	0.0000	1.0000
F81	14	-1784.2267	3596.4534	0.0000	1.0000
JC69G	12	-1788.0695	3600.1391	0.0000	1.0000
JC69I	12	-1788.6769	3601.3537	0.0000	1.0000
JC69IG	13	-1788.0387	3602.0774	0.0000	1.0000
K2P	12	-1789.1320	3602.2640	0.0000	1.0000
SYM	16	-1788.6759	3609.3519	0.0000	1.0000
JC69	11	-1801.3391	3624.6782	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1755.8445	3614.9217	0.6093	0.6093
HKYI	16	-1756.3553	3615.9433	0.3656	0.9749
HKYIG	17	-1755.8100	3621.3049	0.0250	0.9999
K2PG	13	-1775.4932	3634.8630	0.0000	1.0000
K2PI	13	-1776.2010	3636.2787	0.0000	1.0000
F81G	15	-1770.6394	3638.0596	0.0000	1.0000
HKY	15	-1770.8422	3638.4652	0.0000	1.0000
GTRG	20	-1755.0119	3639.0647	0.0000	1.0000
F81I	15	-1771.3089	3639.3985	0.0000	1.0000
GTRI	20	-1755.8991	3640.8391	0.0000	1.0000
K2PIG	14	-1775.4461	3641.2210	0.0000	1.0000
F81IG	16	-1770.5992	3644.4311	0.0000	1.0000
GTRIG	21	-1754.9325	3645.3580	0.0000	1.0000
JC69G	12	-1788.0695	3653.5637	0.0000	1.0000
JC69I	12	-1788.6769	3654.7783	0.0000	1.0000
K2P	12	-1789.1320	3655.6886	0.0000	1.0000
F81	14	-1784.2267	3658.7821	0.0000	1.0000
SYMG	17	-1774.9553	3659.5954	0.0000	1.0000
JC69IG	13	-1788.0387	3659.9541	0.0000	1.0000
SYMI	17	-1775.6261	3660.9370	0.0000	1.0000
GTR	19	-1769.4516	3661.4922	0.0000	1.0000
SYMIG	18	-1774.9065	3665.9499	0.0000	1.0000
JC69	11	-1801.3391	3673.6508	0.0000	1.0000
SYM	16	-1788.6759	3680.5846	0.0000	1.0000



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

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

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


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

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


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

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

