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

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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYIG
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
HKYIG	61	-2129.0871	4393.7542	0.4742	0.4742
HKYG	60	-2130.9309	4394.9801	0.2569	0.7311
GTRIG	65	-2125.0113	4395.5381	0.1944	0.9255
GTRG	64	-2127.2371	4397.4923	0.0732	0.9987
HKYI	60	-2136.3384	4405.7950	0.0012	0.9998
GTRI	64	-2133.2074	4409.4328	0.0002	1.0000
K2PIG	58	-2150.5820	4429.3854	0.0000	1.0000
K2PG	57	-2152.3721	4430.5302	0.0000	1.0000
SYMIG	62	-2149.7223	4437.4950	0.0000	1.0000
SYMG	61	-2151.3070	4438.1938	0.0000	1.0000
K2PI	57	-2157.2345	4440.2552	0.0000	1.0000
SYMI	61	-2156.0723	4447.7244	0.0000	1.0000
HKY	59	-2182.5024	4495.6704	0.0000	1.0000
GTR	63	-2178.3661	4497.2619	0.0000	1.0000
F81IG	60	-2192.7312	4518.5807	0.0000	1.0000
F81G	59	-2194.0401	4518.7458	0.0000	1.0000
K2P	56	-2200.4763	4524.3120	0.0000	1.0000
F81I	59	-2198.7561	4528.1777	0.0000	1.0000
SYM	60	-2199.1453	4531.4088	0.0000	1.0000
JC69IG	57	-2211.1072	4548.0004	0.0000	1.0000
JC69G	56	-2212.3279	4548.0153	0.0000	1.0000
JC69I	56	-2216.6820	4556.7235	0.0000	1.0000
F81	58	-2239.2759	4606.7732	0.0000	1.0000
JC69	55	-2256.9921	4634.9256	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2125.0113	4380.0227	0.3814	0.3814
HKYIG	61	-2129.0871	4380.1743	0.3536	0.7350
HKYG	60	-2130.9309	4381.8618	0.1521	0.8871
GTRG	64	-2127.2371	4382.4742	0.1120	0.9990
HKYI	60	-2136.3384	4392.6767	0.0007	0.9997
GTRI	64	-2133.2074	4394.4147	0.0003	1.0000
K2PIG	58	-2150.5820	4417.1640	0.0000	1.0000
K2PG	57	-2152.3721	4418.7441	0.0000	1.0000
SYMIG	62	-2149.7223	4423.4446	0.0000	1.0000
SYMG	61	-2151.3070	4424.6139	0.0000	1.0000
K2PI	57	-2157.2345	4428.4691	0.0000	1.0000
SYMI	61	-2156.0723	4434.1445	0.0000	1.0000
GTR	63	-2178.3661	4482.7321	0.0000	1.0000
HKY	59	-2182.5024	4483.0049	0.0000	1.0000
F81IG	60	-2192.7312	4505.4625	0.0000	1.0000
F81G	59	-2194.0401	4506.0803	0.0000	1.0000
K2P	56	-2200.4763	4512.9526	0.0000	1.0000
F81I	59	-2198.7561	4515.5122	0.0000	1.0000
SYM	60	-2199.1453	4518.2905	0.0000	1.0000
JC69IG	57	-2211.1072	4536.2143	0.0000	1.0000
JC69G	56	-2212.3279	4536.6559	0.0000	1.0000
JC69I	56	-2216.6820	4545.3641	0.0000	1.0000
F81	58	-2239.2759	4594.5518	0.0000	1.0000
JC69	55	-2256.9921	4623.9843	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2130.9309	4647.5482	0.7945	0.7945
HKYIG	61	-2129.0871	4650.2887	0.2018	0.9963
HKYI	60	-2136.3384	4658.3631	0.0036	0.9999
GTRG	64	-2127.2371	4665.8729	0.0001	1.0000
GTRIG	65	-2125.0113	4667.8495	0.0000	1.0000
K2PG	57	-2152.3721	4671.1461	0.0000	1.0000
K2PIG	58	-2150.5820	4673.9941	0.0000	1.0000
GTRI	64	-2133.2074	4677.8135	0.0000	1.0000
K2PI	57	-2157.2345	4680.8711	0.0000	1.0000
SYMG	61	-2151.3070	4694.7284	0.0000	1.0000
SYMIG	62	-2149.7223	4697.9872	0.0000	1.0000
SYMI	61	-2156.0723	4704.2589	0.0000	1.0000
HKY	59	-2182.5024	4744.2631	0.0000	1.0000
K2P	56	-2200.4763	4760.9265	0.0000	1.0000
GTR	63	-2178.3661	4761.7028	0.0000	1.0000
F81G	59	-2194.0401	4767.3385	0.0000	1.0000
F81IG	60	-2192.7312	4771.1488	0.0000	1.0000
F81I	59	-2198.7561	4776.7704	0.0000	1.0000
SYM	60	-2199.1453	4783.9768	0.0000	1.0000
JC69G	56	-2212.3279	4784.6298	0.0000	1.0000
JC69IG	57	-2211.1072	4788.6163	0.0000	1.0000
JC69I	56	-2216.6820	4793.3380	0.0000	1.0000
F81	58	-2239.2759	4851.3819	0.0000	1.0000
JC69	55	-2256.9921	4867.5300	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/4303_hg18.phylip = 1 - 619;
 Partition Dummy = 1:261-aggressive-phylip/4303_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 (HKYIG)]

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


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

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

