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

Input data from file "261-aggressive-phylip/1250_hg18.phylip" (29 taxa, 465 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	60	-1730.2014	3598.5217	0.3127	0.3127
K2PG	57	-1734.3705	3598.9867	0.2479	0.5606
HKYIG	61	-1729.3151	3599.3994	0.2017	0.7623
K2PIG	58	-1733.5849	3600.0269	0.1473	0.9096
GTRG	64	-1727.2098	3603.2195	0.0299	0.9395
GTRIG	65	-1726.1893	3603.8824	0.0214	0.9609
K2PI	57	-1737.1422	3604.5301	0.0155	0.9764
HKYI	60	-1733.3597	3604.8383	0.0133	0.9897
SYMG	61	-1732.8923	3606.5538	0.0056	0.9953
SYMIG	62	-1732.0585	3607.5499	0.0034	0.9988
GTRI	64	-1730.7014	3610.2027	0.0009	0.9997
SYMI	61	-1735.7709	3612.3110	0.0003	1.0000
HKY	59	-1767.4447	3670.3710	0.0000	1.0000
K2P	56	-1772.3986	3672.4443	0.0000	1.0000
JC69G	56	-1774.6713	3676.9897	0.0000	1.0000
GTR	63	-1765.5138	3677.1373	0.0000	1.0000
JC69IG	57	-1773.9335	3678.1126	0.0000	1.0000
F81G	59	-1771.8679	3679.2172	0.0000	1.0000
SYM	60	-1770.9793	3680.0773	0.0000	1.0000
F81IG	60	-1771.1018	3680.3224	0.0000	1.0000
JC69I	56	-1777.4126	3682.4722	0.0000	1.0000
F81I	59	-1774.7455	3684.9726	0.0000	1.0000
F81	58	-1807.8871	3748.6314	0.0000	1.0000
JC69	55	-1811.8576	3748.7764	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1730.2014	3580.4029	0.3020	0.3020
HKYIG	61	-1729.3151	3580.6301	0.2695	0.5715
GTRIG	65	-1726.1893	3582.3787	0.1124	0.6839
GTRG	64	-1727.2098	3582.4195	0.1102	0.7941
K2PG	57	-1734.3705	3582.7410	0.0938	0.8879
K2PIG	58	-1733.5849	3583.1697	0.0757	0.9636
HKYI	60	-1733.3597	3586.7195	0.0128	0.9764
SYMG	61	-1732.8923	3587.7845	0.0075	0.9840
SYMIG	62	-1732.0585	3588.1171	0.0064	0.9904
K2PI	57	-1737.1422	3588.2844	0.0059	0.9962
GTRI	64	-1730.7014	3589.4027	0.0034	0.9996
SYMI	61	-1735.7709	3593.5418	0.0004	1.0000
HKY	59	-1767.4447	3652.8895	0.0000	1.0000
K2P	56	-1772.3986	3656.7973	0.0000	1.0000
GTR	63	-1765.5138	3657.0275	0.0000	1.0000
JC69G	56	-1774.6713	3661.3427	0.0000	1.0000
F81G	59	-1771.8679	3661.7357	0.0000	1.0000
JC69IG	57	-1773.9335	3661.8669	0.0000	1.0000
SYM	60	-1770.9793	3661.9585	0.0000	1.0000
F81IG	60	-1771.1018	3662.2036	0.0000	1.0000
JC69I	56	-1777.4126	3666.8252	0.0000	1.0000
F81I	59	-1774.7455	3667.4911	0.0000	1.0000
F81	58	-1807.8871	3731.7743	0.0000	1.0000
JC69	55	-1811.8576	3733.7152	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1734.3705	3818.8372	0.8534	0.8534
K2PIG	58	-1733.5849	3823.4079	0.0868	0.9402
K2PI	57	-1737.1422	3824.3805	0.0534	0.9936
HKYG	60	-1730.2014	3828.9251	0.0055	0.9991
HKYIG	61	-1729.3151	3833.2944	0.0006	0.9997
HKYI	60	-1733.3597	3835.2417	0.0002	1.0000
SYMG	61	-1732.8923	3840.4488	0.0000	1.0000
SYMIG	62	-1732.0585	3844.9234	0.0000	1.0000
SYMI	61	-1735.7709	3846.2061	0.0000	1.0000
GTRG	64	-1727.2098	3847.5099	0.0000	1.0000
GTRIG	65	-1726.1893	3851.6111	0.0000	1.0000
GTRI	64	-1730.7014	3854.4931	0.0000	1.0000
K2P	56	-1772.3986	3888.7514	0.0000	1.0000
JC69G	56	-1774.6713	3893.2968	0.0000	1.0000
HKY	59	-1767.4447	3897.2697	0.0000	1.0000
JC69IG	57	-1773.9335	3897.9631	0.0000	1.0000
JC69I	56	-1777.4126	3898.7793	0.0000	1.0000
F81G	59	-1771.8679	3906.1159	0.0000	1.0000
SYM	60	-1770.9793	3910.4808	0.0000	1.0000
F81IG	60	-1771.1018	3910.7258	0.0000	1.0000
F81I	59	-1774.7455	3911.8713	0.0000	1.0000
GTR	63	-1765.5138	3917.9759	0.0000	1.0000
JC69	55	-1811.8576	3961.5273	0.0000	1.0000
F81	58	-1807.8871	3972.0124	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 261-aggressive-phylip/1250_hg18.phylip = 1 - 465;
 Partition Dummy = 1:261-aggressive-phylip/1250_hg18.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 261-aggressive-phylip/1250_hg18.phylip = 1 - 465;
 Partition Dummy = 1:261-aggressive-phylip/1250_hg18.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 261-aggressive-phylip/1250_hg18.phylip = 1 - 465;
 Partition Dummy = 1:261-aggressive-phylip/1250_hg18.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

