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

Input data from file "261-aggressive-phylip/5811_hg18.phylip" (29 taxa, 633 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	64	-2126.3483	4395.3444	0.8375	0.8375
GTRIG	65	-2127.7352	4400.6027	0.0604	0.8980
HKYG	60	-2134.2742	4401.3456	0.0417	0.9396
GTRI	64	-2129.4447	4401.5373	0.0379	0.9775
HKYIG	61	-2133.7425	4402.7319	0.0208	0.9983
HKYI	60	-2137.4946	4407.7864	0.0017	1.0000
GTR	63	-2159.6390	4459.4503	0.0000	1.0000
HKY	59	-2168.7831	4467.9222	0.0000	1.0000
K2PG	57	-2173.7090	4472.9171	0.0000	1.0000
K2PIG	58	-2173.0958	4474.1149	0.0000	1.0000
K2PI	57	-2177.4385	4480.3762	0.0000	1.0000
SYMG	61	-2173.4336	4482.1141	0.0000	1.0000
SYMIG	62	-2172.8141	4483.3334	0.0000	1.0000
SYMI	61	-2177.0161	4489.2791	0.0000	1.0000
F81G	59	-2203.1772	4536.7105	0.0000	1.0000
F81IG	60	-2202.8057	4538.4087	0.0000	1.0000
K2P	56	-2208.0909	4539.2650	0.0000	1.0000
F81I	59	-2206.2886	4542.9333	0.0000	1.0000
SYM	60	-2207.5018	4547.8009	0.0000	1.0000
JC69G	56	-2234.2896	4591.6626	0.0000	1.0000
JC69IG	57	-2233.8312	4593.1615	0.0000	1.0000
JC69I	56	-2237.7831	4598.6496	0.0000	1.0000
F81	58	-2236.2816	4600.4866	0.0000	1.0000
JC69	55	-2268.0143	4656.7046	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2126.3483	4380.6966	0.8547	0.8547
GTRIG	65	-2127.7352	4385.4704	0.0786	0.9333
GTRI	64	-2129.4447	4386.8894	0.0386	0.9719
HKYG	60	-2134.2742	4388.5484	0.0169	0.9888
HKYIG	61	-2133.7425	4389.4850	0.0106	0.9993
HKYI	60	-2137.4946	4394.9892	0.0007	1.0000
GTR	63	-2159.6390	4445.2780	0.0000	1.0000
HKY	59	-2168.7831	4455.5662	0.0000	1.0000
K2PG	57	-2173.7090	4461.4180	0.0000	1.0000
K2PIG	58	-2173.0958	4462.1916	0.0000	1.0000
SYMG	61	-2173.4336	4468.8672	0.0000	1.0000
K2PI	57	-2177.4385	4468.8771	0.0000	1.0000
SYMIG	62	-2172.8141	4469.6282	0.0000	1.0000
SYMI	61	-2177.0161	4476.0321	0.0000	1.0000
F81G	59	-2203.1772	4524.3545	0.0000	1.0000
F81IG	60	-2202.8057	4525.6115	0.0000	1.0000
K2P	56	-2208.0909	4528.1817	0.0000	1.0000
F81I	59	-2206.2886	4530.5772	0.0000	1.0000
SYM	60	-2207.5018	4535.0037	0.0000	1.0000
JC69G	56	-2234.2896	4580.5793	0.0000	1.0000
JC69IG	57	-2233.8312	4581.6624	0.0000	1.0000
JC69I	56	-2237.7831	4587.5663	0.0000	1.0000
F81	58	-2236.2816	4588.5632	0.0000	1.0000
JC69	55	-2268.0143	4646.0287	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2134.2742	4655.5766	0.8970	0.8970
HKYIG	61	-2133.7425	4660.9637	0.0607	0.9576
HKYI	60	-2137.4946	4662.0175	0.0358	0.9935
GTRG	64	-2126.3483	4665.5267	0.0062	0.9997
GTRI	64	-2129.4447	4671.7195	0.0003	0.9999
GTRIG	65	-2127.7352	4674.7510	0.0001	1.0000
K2PG	57	-2173.7090	4715.0948	0.0000	1.0000
HKY	59	-2168.7831	4718.1440	0.0000	1.0000
K2PIG	58	-2173.0958	4720.3189	0.0000	1.0000
K2PI	57	-2177.4385	4722.5539	0.0000	1.0000
GTR	63	-2159.6390	4725.6577	0.0000	1.0000
SYMG	61	-2173.4336	4740.3459	0.0000	1.0000
SYMIG	62	-2172.8141	4745.5573	0.0000	1.0000
SYMI	61	-2177.0161	4747.5108	0.0000	1.0000
K2P	56	-2208.0909	4777.4080	0.0000	1.0000
F81G	59	-2203.1772	4786.9322	0.0000	1.0000
F81IG	60	-2202.8057	4792.6397	0.0000	1.0000
F81I	59	-2206.2886	4793.1550	0.0000	1.0000
SYM	60	-2207.5018	4802.0319	0.0000	1.0000
JC69G	56	-2234.2896	4829.8056	0.0000	1.0000
JC69IG	57	-2233.8312	4835.3392	0.0000	1.0000
JC69I	56	-2237.7831	4836.7926	0.0000	1.0000
F81	58	-2236.2816	4846.6905	0.0000	1.0000
JC69	55	-2268.0143	4890.8046	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 261-aggressive-phylip/5811_hg18.phylip = 1 - 633;
 Partition Dummy = 1:261-aggressive-phylip/5811_hg18.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 261-aggressive-phylip/5811_hg18.phylip = 1 - 633;
 Partition Dummy = 1:261-aggressive-phylip/5811_hg18.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 261-aggressive-phylip/5811_hg18.phylip = 1 - 633;
 Partition Dummy = 1:261-aggressive-phylip/5811_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

