
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:59 08/01/2011

Input data from file "261-aggressive-phylip/9670_hg18.phylip" (29 taxa, 435 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	64	-1718.6415	3587.7696	0.3554	0.3554
GTRG	64	-1718.6985	3587.8836	0.3357	0.6911
GTRIG	65	-1718.0087	3589.2695	0.1679	0.8590
HKYI	60	-1725.5805	3590.7333	0.0808	0.9397
HKYG	60	-1726.4172	3592.4066	0.0350	0.9747
HKYIG	61	-1725.3880	3593.0547	0.0253	1.0000
K2PI	57	-1738.8174	3609.1732	0.0000	1.0000
K2PIG	58	-1738.7792	3611.7604	0.0000	1.0000
K2PG	57	-1740.2246	3611.9877	0.0000	1.0000
SYMI	61	-1735.2323	3612.7435	0.0000	1.0000
SYMG	61	-1736.4928	3615.2645	0.0000	1.0000
SYMIG	62	-1735.1647	3615.3294	0.0000	1.0000
GTR	63	-1753.1907	3654.1173	0.0000	1.0000
HKY	59	-1761.6263	3660.1326	0.0000	1.0000
K2P	56	-1777.1404	3683.1696	0.0000	1.0000
F81I	59	-1773.4191	3683.7181	0.0000	1.0000
SYM	60	-1772.8745	3685.3213	0.0000	1.0000
F81G	59	-1774.3439	3685.5678	0.0000	1.0000
F81IG	60	-1773.2732	3686.1187	0.0000	1.0000
JC69I	56	-1784.2084	3697.3057	0.0000	1.0000
JC69G	56	-1785.2107	3699.3104	0.0000	1.0000
JC69IG	57	-1784.0983	3699.7350	0.0000	1.0000
F81	58	-1810.7623	3755.7267	0.0000	1.0000
JC69	55	-1822.8823	3772.0178	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1718.6415	3565.2831	0.3654	0.3654
GTRG	64	-1718.6985	3565.3971	0.3452	0.7106
GTRIG	65	-1718.0087	3566.0175	0.2531	0.9637
HKYI	60	-1725.5805	3571.1611	0.0193	0.9830
HKYIG	61	-1725.3880	3572.7759	0.0086	0.9916
HKYG	60	-1726.4172	3572.8344	0.0084	1.0000
K2PI	57	-1738.8174	3591.6348	0.0000	1.0000
SYMI	61	-1735.2323	3592.4646	0.0000	1.0000
K2PIG	58	-1738.7792	3593.5583	0.0000	1.0000
SYMIG	62	-1735.1647	3594.3294	0.0000	1.0000
K2PG	57	-1740.2246	3594.4493	0.0000	1.0000
SYMG	61	-1736.4928	3594.9857	0.0000	1.0000
GTR	63	-1753.1907	3632.3815	0.0000	1.0000
HKY	59	-1761.6263	3641.2526	0.0000	1.0000
F81I	59	-1773.4191	3664.8381	0.0000	1.0000
SYM	60	-1772.8745	3665.7491	0.0000	1.0000
K2P	56	-1777.1404	3666.2807	0.0000	1.0000
F81IG	60	-1773.2732	3666.5465	0.0000	1.0000
F81G	59	-1774.3439	3666.6878	0.0000	1.0000
JC69I	56	-1784.2084	3680.4168	0.0000	1.0000
JC69IG	57	-1784.0983	3682.1965	0.0000	1.0000
JC69G	56	-1785.2107	3682.4215	0.0000	1.0000
F81	58	-1810.7623	3737.5245	0.0000	1.0000
JC69	55	-1822.8823	3755.7645	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1725.5805	3815.6818	0.6565	0.6565
HKYG	60	-1726.4172	3817.3551	0.2843	0.9408
HKYIG	61	-1725.3880	3821.3720	0.0382	0.9790
K2PI	57	-1738.8174	3823.9295	0.0106	0.9896
GTRI	64	-1718.6415	3826.1052	0.0036	0.9932
GTRG	64	-1718.6985	3826.2192	0.0034	0.9965
K2PG	57	-1740.2246	3826.7440	0.0026	0.9991
K2PIG	58	-1738.7792	3829.9284	0.0005	0.9997
GTRIG	65	-1718.0087	3830.9150	0.0003	1.0000
SYMI	61	-1735.2323	3841.0607	0.0000	1.0000
SYMG	61	-1736.4928	3843.5818	0.0000	1.0000
SYMIG	62	-1735.1647	3847.0008	0.0000	1.0000
HKY	59	-1761.6263	3881.6980	0.0000	1.0000
GTR	63	-1753.1907	3889.1283	0.0000	1.0000
K2P	56	-1777.1404	3894.5001	0.0000	1.0000
F81I	59	-1773.4191	3905.2835	0.0000	1.0000
F81G	59	-1774.3439	3907.1332	0.0000	1.0000
JC69I	56	-1784.2084	3908.6362	0.0000	1.0000
SYM	60	-1772.8745	3910.2698	0.0000	1.0000
JC69G	56	-1785.2107	3910.6409	0.0000	1.0000
F81IG	60	-1773.2732	3911.0673	0.0000	1.0000
JC69IG	57	-1784.0983	3914.4913	0.0000	1.0000
F81	58	-1810.7623	3973.8946	0.0000	1.0000
JC69	55	-1822.8823	3979.9086	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 261-aggressive-phylip/9670_hg18.phylip = 1 - 435;
 Partition Dummy = 1:261-aggressive-phylip/9670_hg18.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 261-aggressive-phylip/9670_hg18.phylip = 1 - 435;
 Partition Dummy = 1:261-aggressive-phylip/9670_hg18.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 261-aggressive-phylip/9670_hg18.phylip = 1 - 435;
 Partition Dummy = 1:261-aggressive-phylip/9670_hg18.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

