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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1562.4647	3166.6127	0.5692	0.5692
GTRIG	21	-1562.4673	3168.7900	0.1916	0.7608
GTRI	20	-1564.1644	3170.0121	0.1040	0.8648
HKYG	16	-1568.6803	3170.4421	0.0839	0.9487
HKYIG	17	-1568.6778	3172.5746	0.0289	0.9776
HKYI	16	-1570.0015	3173.0846	0.0224	0.9999
GTR	19	-1572.7350	3184.9900	0.0001	1.0000
HKY	15	-1580.4194	3191.7912	0.0000	1.0000
F81G	15	-1590.7337	3212.4197	0.0000	1.0000
F81IG	16	-1590.7336	3214.5487	0.0000	1.0000
F81I	15	-1591.8348	3214.6221	0.0000	1.0000
K2PG	13	-1594.3145	3215.3483	0.0000	1.0000
K2PIG	14	-1594.2921	3217.4158	0.0000	1.0000
K2PI	13	-1595.7887	3218.2969	0.0000	1.0000
SYMG	17	-1592.1621	3219.5433	0.0000	1.0000
SYMIG	18	-1592.1462	3221.6576	0.0000	1.0000
SYMI	17	-1593.7821	3222.7833	0.0000	1.0000
F81	14	-1602.2311	3233.2938	0.0000	1.0000
K2P	12	-1606.8623	3238.3401	0.0000	1.0000
SYM	16	-1603.9951	3241.0717	0.0000	1.0000
JC69G	12	-1613.4438	3251.5029	0.0000	1.0000
JC69IG	13	-1613.4330	3253.5854	0.0000	1.0000
JC69I	12	-1614.6859	3253.9873	0.0000	1.0000
JC69	11	-1625.3956	3273.3110	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1562.4647	3164.9293	0.5786	0.5786
GTRIG	21	-1562.4673	3166.9346	0.2123	0.7910
GTRI	20	-1564.1644	3168.3288	0.1057	0.8967
HKYG	16	-1568.6803	3169.3606	0.0631	0.9598
HKYIG	17	-1568.6778	3171.3555	0.0233	0.9831
HKYI	16	-1570.0015	3172.0030	0.0168	0.9999
GTR	19	-1572.7350	3183.4700	0.0001	1.0000
HKY	15	-1580.4194	3190.8388	0.0000	1.0000
F81G	15	-1590.7337	3211.4673	0.0000	1.0000
F81IG	16	-1590.7336	3213.4671	0.0000	1.0000
F81I	15	-1591.8348	3213.6697	0.0000	1.0000
K2PG	13	-1594.3145	3214.6289	0.0000	1.0000
K2PIG	14	-1594.2921	3216.5842	0.0000	1.0000
K2PI	13	-1595.7887	3217.5775	0.0000	1.0000
SYMG	17	-1592.1621	3218.3242	0.0000	1.0000
SYMIG	18	-1592.1462	3220.2923	0.0000	1.0000
SYMI	17	-1593.7821	3221.5641	0.0000	1.0000
F81	14	-1602.2311	3232.4621	0.0000	1.0000
K2P	12	-1606.8623	3237.7247	0.0000	1.0000
SYM	16	-1603.9951	3239.9902	0.0000	1.0000
JC69G	12	-1613.4438	3250.8875	0.0000	1.0000
JC69IG	13	-1613.4330	3252.8660	0.0000	1.0000
JC69I	12	-1614.6859	3253.3719	0.0000	1.0000
JC69	11	-1625.3956	3272.7913	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1568.6803	3237.4218	0.7615	0.7615
HKYI	16	-1570.0015	3240.0643	0.2032	0.9647
HKYIG	17	-1568.6778	3243.6706	0.0335	0.9981
GTRG	20	-1562.4647	3250.0059	0.0014	0.9995
GTRI	20	-1564.1644	3253.4054	0.0003	0.9998
HKY	15	-1580.4194	3254.6463	0.0001	0.9999
GTRIG	21	-1562.4673	3256.2650	0.0001	1.0000
GTR	19	-1572.7350	3264.2928	0.0000	1.0000
K2PG	13	-1594.3145	3269.9287	0.0000	1.0000
K2PI	13	-1595.7887	3272.8773	0.0000	1.0000
F81G	15	-1590.7337	3275.2748	0.0000	1.0000
K2PIG	14	-1594.2921	3276.1378	0.0000	1.0000
F81I	15	-1591.8348	3277.4771	0.0000	1.0000
F81IG	16	-1590.7336	3281.5284	0.0000	1.0000
K2P	12	-1606.8623	3288.7706	0.0000	1.0000
SYMG	17	-1592.1621	3290.6392	0.0000	1.0000
F81	14	-1602.2311	3292.0157	0.0000	1.0000
SYMI	17	-1593.7821	3293.8792	0.0000	1.0000
SYMIG	18	-1592.1462	3296.8612	0.0000	1.0000
JC69G	12	-1613.4438	3301.9335	0.0000	1.0000
JC69I	12	-1614.6859	3304.4178	0.0000	1.0000
SYM	16	-1603.9951	3308.0515	0.0000	1.0000
JC69IG	13	-1613.4330	3308.1658	0.0000	1.0000
JC69	11	-1625.3956	3319.5834	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_18800.phylip = 1 - 520;
 Partition Dummy = 1:683-tomodel-phylip/chr3_18800.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_18800.phylip = 1 - 520;
 Partition Dummy = 1:683-tomodel-phylip/chr3_18800.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_18800.phylip = 1 - 520;
 Partition Dummy = 1:683-tomodel-phylip/chr3_18800.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

