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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-1808.4772	3649.6664	0.3154	0.3154
HKYIG	17	-1807.8610	3650.5242	0.2054	0.5209
HKYI	16	-1809.0572	3650.8265	0.1766	0.6975
GTRG	20	-1805.1111	3651.3275	0.1375	0.8350
GTRIG	21	-1804.5005	3652.2185	0.0881	0.9230
GTRI	20	-1805.7204	3652.5461	0.0747	0.9978
K2PG	13	-1817.6215	3661.7176	0.0008	0.9985
K2PIG	14	-1817.0096	3662.5676	0.0005	0.9990
K2PI	13	-1818.1594	3662.7933	0.0004	0.9995
SYMG	17	-1814.6067	3664.0154	0.0002	0.9997
SYMIG	18	-1814.0177	3664.9331	0.0002	0.9999
SYMI	17	-1815.1609	3665.1239	0.0001	1.0000
HKY	15	-1826.0261	3682.6797	0.0000	1.0000
GTR	19	-1822.5126	3684.0239	0.0000	1.0000
K2P	12	-1835.4038	3695.2139	0.0000	1.0000
SYM	16	-1832.1639	3697.0399	0.0000	1.0000
F81G	15	-1836.0508	3702.7291	0.0000	1.0000
F81IG	16	-1835.5406	3703.7933	0.0000	1.0000
F81I	15	-1836.6246	3703.8766	0.0000	1.0000
JC69G	12	-1844.8038	3714.0139	0.0000	1.0000
JC69IG	13	-1844.2944	3715.0635	0.0000	1.0000
JC69I	12	-1845.3715	3715.1493	0.0000	1.0000
F81	14	-1852.9576	3734.4636	0.0000	1.0000
JC69	11	-1861.7411	3745.8254	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1808.4772	3648.9543	0.2919	0.2919
HKYIG	17	-1807.8610	3649.7221	0.1988	0.4907
HKYI	16	-1809.0572	3650.1144	0.1634	0.6541
GTRG	20	-1805.1111	3650.2223	0.1548	0.8090
GTRIG	21	-1804.5005	3651.0011	0.1049	0.9139
GTRI	20	-1805.7204	3651.4409	0.0842	0.9981
K2PG	13	-1817.6215	3661.2430	0.0006	0.9987
K2PIG	14	-1817.0096	3662.0193	0.0004	0.9991
K2PI	13	-1818.1594	3662.3188	0.0004	0.9995
SYMG	17	-1814.6067	3663.2133	0.0002	0.9997
SYMIG	18	-1814.0177	3664.0355	0.0002	0.9999
SYMI	17	-1815.1609	3664.3218	0.0001	1.0000
HKY	15	-1826.0261	3682.0523	0.0000	1.0000
GTR	19	-1822.5126	3683.0252	0.0000	1.0000
K2P	12	-1835.4038	3694.8076	0.0000	1.0000
SYM	16	-1832.1639	3696.3279	0.0000	1.0000
F81G	15	-1836.0508	3702.1017	0.0000	1.0000
F81IG	16	-1835.5406	3703.0813	0.0000	1.0000
F81I	15	-1836.6246	3703.2492	0.0000	1.0000
JC69G	12	-1844.8038	3713.6077	0.0000	1.0000
JC69IG	13	-1844.2944	3714.5889	0.0000	1.0000
JC69I	12	-1845.3715	3714.7430	0.0000	1.0000
F81	14	-1852.9576	3733.9153	0.0000	1.0000
JC69	11	-1861.7411	3745.4821	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1817.6215	3721.8305	0.4260	0.4260
K2PI	13	-1818.1594	3722.9062	0.2488	0.6748
HKYG	16	-1808.4772	3723.5235	0.1827	0.8575
HKYI	16	-1809.0572	3724.6836	0.1023	0.9598
K2PIG	14	-1817.0096	3727.2673	0.0281	0.9879
HKYIG	17	-1807.8610	3728.9519	0.0121	1.0000
SYMG	17	-1814.6067	3742.4431	0.0000	1.0000
GTRG	20	-1805.1111	3743.4338	0.0000	1.0000
SYMI	17	-1815.1609	3743.5515	0.0000	1.0000
GTRI	20	-1805.7204	3744.6524	0.0000	1.0000
SYMIG	18	-1814.0177	3747.9259	0.0000	1.0000
GTRIG	21	-1804.5005	3748.8731	0.0000	1.0000
K2P	12	-1835.4038	3750.7345	0.0000	1.0000
HKY	15	-1826.0261	3751.9609	0.0000	1.0000
JC69G	12	-1844.8038	3769.5346	0.0000	1.0000
JC69I	12	-1845.3715	3770.6699	0.0000	1.0000
SYM	16	-1832.1639	3770.8971	0.0000	1.0000
GTR	19	-1822.5126	3771.5761	0.0000	1.0000
F81G	15	-1836.0508	3772.0103	0.0000	1.0000
F81I	15	-1836.6246	3773.1578	0.0000	1.0000
JC69IG	13	-1844.2944	3775.1764	0.0000	1.0000
F81IG	16	-1835.5406	3777.6505	0.0000	1.0000
JC69	11	-1861.7411	3796.7485	0.0000	1.0000
F81	14	-1852.9576	3799.1633	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/chr4_7559.phylip = 1 - 781;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7559.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/chr4_7559.phylip = 1 - 781;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7559.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_7559.phylip = 1 - 781;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7559.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

