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

Input data from file "261-aggressive-phylip/9836_hg18.phylip" (29 taxa, 549 characters)

Minimum AIC  model: GTRIG
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	60	-2162.5013	4460.0027	0.4383	0.4383
HKYIG	61	-2161.8181	4461.1680	0.2447	0.6830
GTRIG	65	-2156.7402	4461.2443	0.2356	0.9186
GTRG	64	-2159.0902	4463.3705	0.0814	1.0000
HKYI	60	-2171.9859	4478.9717	0.0000	1.0000
GTRI	64	-2169.3765	4483.9431	0.0000	1.0000
SYMG	61	-2181.9834	4501.4986	0.0000	1.0000
K2PIG	58	-2186.5091	4502.9856	0.0000	1.0000
K2PG	57	-2189.0478	4505.5619	0.0000	1.0000
SYMIG	62	-2185.5658	4511.2056	0.0000	1.0000
K2PI	57	-2192.9119	4513.2903	0.0000	1.0000
SYMI	61	-2195.3538	4528.2395	0.0000	1.0000
HKY	59	-2227.2104	4586.8994	0.0000	1.0000
F81G	59	-2228.1724	4588.8234	0.0000	1.0000
GTR	63	-2223.5088	4589.6444	0.0000	1.0000
F81IG	60	-2228.2313	4591.4626	0.0000	1.0000
F81I	59	-2238.6321	4609.7426	0.0000	1.0000
JC69IG	57	-2241.9022	4611.2708	0.0000	1.0000
JC69G	56	-2249.1053	4623.1861	0.0000	1.0000
JC69I	56	-2255.1188	4635.2133	0.0000	1.0000
K2P	56	-2255.3066	4635.5888	0.0000	1.0000
SYM	60	-2254.3553	4643.7107	0.0000	1.0000
F81	58	-2290.1595	4710.2864	0.0000	1.0000
JC69	55	-2307.3047	4737.1044	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2156.7402	4443.4803	0.4839	0.4839
HKYG	60	-2162.5013	4445.0027	0.2260	0.7099
HKYIG	61	-2161.8181	4445.6362	0.1647	0.8746
GTRG	64	-2159.0902	4446.1805	0.1254	1.0000
HKYI	60	-2171.9859	4463.9717	0.0000	1.0000
GTRI	64	-2169.3765	4466.7530	0.0000	1.0000
SYMG	61	-2181.9834	4485.9668	0.0000	1.0000
K2PIG	58	-2186.5091	4489.0183	0.0000	1.0000
K2PG	57	-2189.0478	4492.0955	0.0000	1.0000
SYMIG	62	-2185.5658	4495.1315	0.0000	1.0000
K2PI	57	-2192.9119	4499.8239	0.0000	1.0000
SYMI	61	-2195.3538	4512.7077	0.0000	1.0000
HKY	59	-2227.2104	4572.4208	0.0000	1.0000
GTR	63	-2223.5088	4573.0176	0.0000	1.0000
F81G	59	-2228.1724	4574.3449	0.0000	1.0000
F81IG	60	-2228.2313	4576.4626	0.0000	1.0000
F81I	59	-2238.6321	4595.2641	0.0000	1.0000
JC69IG	57	-2241.9022	4597.8044	0.0000	1.0000
JC69G	56	-2249.1053	4610.2105	0.0000	1.0000
JC69I	56	-2255.1188	4622.2377	0.0000	1.0000
K2P	56	-2255.3066	4622.6132	0.0000	1.0000
SYM	60	-2254.3553	4628.7107	0.0000	1.0000
F81	58	-2290.1595	4696.3191	0.0000	1.0000
JC69	55	-2307.3047	4724.6094	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2162.5013	4703.4886	0.9219	0.9219
HKYIG	61	-2161.8181	4708.4302	0.0779	0.9998
GTRG	64	-2159.0902	4721.8988	0.0001	0.9999
HKYI	60	-2171.9859	4722.4576	0.0001	1.0000
GTRIG	65	-2156.7402	4723.5067	0.0000	1.0000
K2PG	57	-2189.0478	4737.6572	0.0000	1.0000
K2PIG	58	-2186.5091	4738.8880	0.0000	1.0000
GTRI	64	-2169.3765	4742.4713	0.0000	1.0000
K2PI	57	-2192.9119	4745.3855	0.0000	1.0000
SYMG	61	-2181.9834	4748.7608	0.0000	1.0000
SYMIG	62	-2185.5658	4762.2336	0.0000	1.0000
SYMI	61	-2195.3538	4775.5017	0.0000	1.0000
HKY	59	-2227.2104	4826.5986	0.0000	1.0000
F81G	59	-2228.1724	4828.5227	0.0000	1.0000
F81IG	60	-2228.2313	4834.9485	0.0000	1.0000
JC69IG	57	-2241.9022	4843.3660	0.0000	1.0000
GTR	63	-2223.5088	4844.4278	0.0000	1.0000
F81I	59	-2238.6321	4849.4419	0.0000	1.0000
JC69G	56	-2249.1053	4851.4641	0.0000	1.0000
JC69I	56	-2255.1188	4863.4912	0.0000	1.0000
K2P	56	-2255.3066	4863.8667	0.0000	1.0000
SYM	60	-2254.3553	4887.1966	0.0000	1.0000
F81	58	-2290.1595	4946.1888	0.0000	1.0000
JC69	55	-2307.3047	4961.5549	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9836_hg18.phylip = 1 - 549;
 Partition Dummy = 1:261-aggressive-phylip/9836_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9836_hg18.phylip = 1 - 549;
 Partition Dummy = 1:261-aggressive-phylip/9836_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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9836_hg18.phylip = 1 - 549;
 Partition Dummy = 1:261-aggressive-phylip/9836_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;


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

