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

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

Minimum AIC  model: HKYG
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	-2485.7514	5099.2901	0.8404	0.8404
HKYIG	61	-2486.8031	5103.6617	0.0945	0.9349
GTRG	64	-2483.7689	5104.4267	0.0644	0.9993
GTRIG	65	-2487.1396	5113.4556	0.0007	1.0000
F81G	59	-2537.7114	5200.9467	0.0000	1.0000
HKY	59	-2538.2679	5202.0597	0.0000	1.0000
GTR	63	-2534.8375	5204.2812	0.0000	1.0000
F81IG	60	-2539.8753	5207.5379	0.0000	1.0000
HKYI	60	-2543.4638	5214.7148	0.0000	1.0000
GTRI	64	-2539.5851	5216.0592	0.0000	1.0000
K2PIG	58	-2554.6653	5232.5959	0.0000	1.0000
SYMIG	62	-2552.2056	5236.7395	0.0000	1.0000
K2PG	57	-2559.0629	5239.1374	0.0000	1.0000
SYMG	61	-2556.3084	5242.6722	0.0000	1.0000
F81	58	-2591.3416	5305.9487	0.0000	1.0000
F81I	59	-2593.7228	5312.9695	0.0000	1.0000
JC69IG	57	-2599.0929	5319.1975	0.0000	1.0000
JC69G	56	-2602.9281	5324.6189	0.0000	1.0000
K2P	56	-2616.0016	5350.7659	0.0000	1.0000
K2PI	57	-2615.5863	5352.1842	0.0000	1.0000
SYM	60	-2613.5484	5354.8841	0.0000	1.0000
SYMI	61	-2613.1964	5356.4482	0.0000	1.0000
JC69	55	-2659.1126	5434.7437	0.0000	1.0000
JC69I	56	-2660.0005	5438.7637	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-2485.7514	5091.5028	0.7917	0.7917
GTRG	64	-2483.7689	5095.5378	0.1053	0.8969
HKYIG	61	-2486.8031	5095.6063	0.1017	0.9987
GTRIG	65	-2487.1396	5104.2791	0.0013	1.0000
F81G	59	-2537.7114	5193.4228	0.0000	1.0000
HKY	59	-2538.2679	5194.5357	0.0000	1.0000
GTR	63	-2534.8375	5195.6750	0.0000	1.0000
F81IG	60	-2539.8753	5199.7507	0.0000	1.0000
HKYI	60	-2543.4638	5206.9276	0.0000	1.0000
GTRI	64	-2539.5851	5207.1703	0.0000	1.0000
K2PIG	58	-2554.6653	5225.3305	0.0000	1.0000
SYMIG	62	-2552.2056	5228.4111	0.0000	1.0000
K2PG	57	-2559.0629	5232.1258	0.0000	1.0000
SYMG	61	-2556.3084	5234.6168	0.0000	1.0000
F81	58	-2591.3416	5298.6833	0.0000	1.0000
F81I	59	-2593.7228	5305.4456	0.0000	1.0000
JC69IG	57	-2599.0929	5312.1859	0.0000	1.0000
JC69G	56	-2602.9281	5317.8562	0.0000	1.0000
K2P	56	-2616.0016	5344.0032	0.0000	1.0000
K2PI	57	-2615.5863	5345.1726	0.0000	1.0000
SYM	60	-2613.5484	5347.0969	0.0000	1.0000
SYMI	61	-2613.1964	5348.3928	0.0000	1.0000
JC69	55	-2659.1126	5428.2252	0.0000	1.0000
JC69I	56	-2660.0005	5432.0010	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2485.7514	5386.0281	0.9891	0.9891
HKYIG	61	-2486.8031	5395.0403	0.0109	1.0000
GTRG	64	-2483.7689	5409.6981	0.0000	1.0000
GTRIG	65	-2487.1396	5423.3482	0.0000	1.0000
F81G	59	-2537.7114	5483.0393	0.0000	1.0000
HKY	59	-2538.2679	5484.1523	0.0000	1.0000
F81IG	60	-2539.8753	5494.2760	0.0000	1.0000
HKYI	60	-2543.4638	5501.4529	0.0000	1.0000
GTR	63	-2534.8375	5504.9266	0.0000	1.0000
K2PIG	58	-2554.6653	5510.0383	0.0000	1.0000
K2PG	57	-2559.0629	5511.9248	0.0000	1.0000
GTRI	64	-2539.5851	5521.3306	0.0000	1.0000
SYMIG	62	-2552.2056	5532.7539	0.0000	1.0000
SYMG	61	-2556.3084	5534.0509	0.0000	1.0000
F81	58	-2591.3416	5583.3910	0.0000	1.0000
JC69IG	57	-2599.0929	5591.9849	0.0000	1.0000
JC69G	56	-2602.9281	5592.7464	0.0000	1.0000
F81I	59	-2593.7228	5595.0621	0.0000	1.0000
K2P	56	-2616.0016	5618.8934	0.0000	1.0000
K2PI	57	-2615.5863	5624.9716	0.0000	1.0000
SYM	60	-2613.5484	5641.6222	0.0000	1.0000
SYMI	61	-2613.1964	5647.8269	0.0000	1.0000
JC69	55	-2659.1126	5698.2067	0.0000	1.0000
JC69I	56	-2660.0005	5706.8913	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/3387_hg18.phylip = 1 - 1001;
 Partition Dummy = 1:261-aggressive-phylip/3387_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/3387_hg18.phylip = 1 - 1001;
 Partition Dummy = 1:261-aggressive-phylip/3387_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/3387_hg18.phylip = 1 - 1001;
 Partition Dummy = 1:261-aggressive-phylip/3387_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

