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

Input data from file "261-aggressive-phylip/7639_hg18.phylip" (29 taxa, 657 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	-2296.5271	4735.1083	0.7166	0.7166
GTRIG	65	-2296.4008	4737.3193	0.2372	0.9539
GTRI	64	-2299.4811	4741.0163	0.0374	0.9913
HKYG	60	-2306.1602	4744.6024	0.0062	0.9975
HKYIG	61	-2306.0418	4746.7962	0.0021	0.9995
HKYI	60	-2308.7781	4749.8381	0.0005	1.0000
GTR	63	-2315.5306	4770.6598	0.0000	1.0000
HKY	59	-2324.3194	4778.4981	0.0000	1.0000
K2PG	57	-2345.0201	4815.0785	0.0000	1.0000
K2PIG	58	-2344.9113	4817.2675	0.0000	1.0000
K2PI	57	-2347.5679	4820.1742	0.0000	1.0000
SYMG	61	-2343.0463	4820.8053	0.0000	1.0000
SYMIG	62	-2342.9103	4822.9721	0.0000	1.0000
SYMI	61	-2345.5742	4825.8610	0.0000	1.0000
F81G	59	-2355.3187	4840.4967	0.0000	1.0000
F81IG	60	-2355.2018	4842.6856	0.0000	1.0000
F81I	59	-2357.8212	4845.5017	0.0000	1.0000
K2P	56	-2364.2777	4851.1954	0.0000	1.0000
SYM	60	-2361.8733	4856.0285	0.0000	1.0000
F81	58	-2372.9912	4873.4272	0.0000	1.0000
JC69G	56	-2387.4489	4897.5377	0.0000	1.0000
JC69IG	57	-2387.3238	4899.6861	0.0000	1.0000
JC69I	56	-2389.6918	4902.0236	0.0000	1.0000
JC69	55	-2405.6597	4931.5689	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2296.5271	4721.0542	0.6780	0.6780
GTRIG	65	-2296.4008	4722.8015	0.2830	0.9610
GTRI	64	-2299.4811	4726.9622	0.0353	0.9964
HKYG	60	-2306.1602	4732.3205	0.0024	0.9988
HKYIG	61	-2306.0418	4734.0836	0.0010	0.9998
HKYI	60	-2308.7781	4737.5563	0.0002	1.0000
GTR	63	-2315.5306	4757.0611	0.0000	1.0000
HKY	59	-2324.3194	4766.6388	0.0000	1.0000
K2PG	57	-2345.0201	4804.0401	0.0000	1.0000
K2PIG	58	-2344.9113	4805.8227	0.0000	1.0000
SYMG	61	-2343.0463	4808.0927	0.0000	1.0000
K2PI	57	-2347.5679	4809.1358	0.0000	1.0000
SYMIG	62	-2342.9103	4809.8206	0.0000	1.0000
SYMI	61	-2345.5742	4813.1484	0.0000	1.0000
F81G	59	-2355.3187	4828.6374	0.0000	1.0000
F81IG	60	-2355.2018	4830.4037	0.0000	1.0000
F81I	59	-2357.8212	4833.6424	0.0000	1.0000
K2P	56	-2364.2777	4840.5554	0.0000	1.0000
SYM	60	-2361.8733	4843.7466	0.0000	1.0000
F81	58	-2372.9912	4861.9824	0.0000	1.0000
JC69G	56	-2387.4489	4886.8977	0.0000	1.0000
JC69IG	57	-2387.3238	4888.6477	0.0000	1.0000
JC69I	56	-2389.6918	4891.3836	0.0000	1.0000
JC69	55	-2405.6597	4921.3194	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2306.1602	5001.5815	0.8653	0.8653
HKYI	60	-2308.7781	5006.8173	0.0631	0.9285
HKYIG	61	-2306.0418	5007.8323	0.0380	0.9665
GTRG	64	-2296.5271	5008.2660	0.0306	0.9971
GTRI	64	-2299.4811	5014.1740	0.0016	0.9986
GTRIG	65	-2296.4008	5014.5010	0.0014	1.0000
HKY	59	-2324.3194	5031.4122	0.0000	1.0000
GTR	63	-2315.5306	5039.7852	0.0000	1.0000
K2PG	57	-2345.0201	5059.8381	0.0000	1.0000
K2PI	57	-2347.5679	5064.9338	0.0000	1.0000
K2PIG	58	-2344.9113	5066.1083	0.0000	1.0000
SYMG	61	-2343.0463	5081.8414	0.0000	1.0000
SYMI	61	-2345.5742	5086.8971	0.0000	1.0000
SYMIG	62	-2342.9103	5088.0570	0.0000	1.0000
K2P	56	-2364.2777	5091.8657	0.0000	1.0000
F81G	59	-2355.3187	5093.4107	0.0000	1.0000
F81I	59	-2357.8212	5098.4157	0.0000	1.0000
F81IG	60	-2355.2018	5099.6647	0.0000	1.0000
SYM	60	-2361.8733	5113.0077	0.0000	1.0000
F81	58	-2372.9912	5122.2681	0.0000	1.0000
JC69G	56	-2387.4489	5138.2080	0.0000	1.0000
JC69I	56	-2389.6918	5142.6939	0.0000	1.0000
JC69IG	57	-2387.3238	5144.4456	0.0000	1.0000
JC69	55	-2405.6597	5168.1420	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/7639_hg18.phylip = 1 - 657;
 Partition Dummy = 1:261-aggressive-phylip/7639_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/7639_hg18.phylip = 1 - 657;
 Partition Dummy = 1:261-aggressive-phylip/7639_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/7639_hg18.phylip = 1 - 657;
 Partition Dummy = 1:261-aggressive-phylip/7639_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

