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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PIG


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	-1998.2820	4144.8074	0.9609	0.9609
GTRIG	65	-2000.2426	4151.4121	0.0354	0.9963
GTRI	64	-2003.8752	4155.9937	0.0036	0.9998
K2PIG	58	-2015.7712	4163.9548	0.0001	0.9999
SYMI	61	-2012.0166	4164.3037	0.0001	1.0000
SYMIG	62	-2011.4461	4165.8075	0.0000	1.0000
SYMG	61	-2015.7830	4171.8365	0.0000	1.0000
HKYIG	61	-2016.7760	4173.8225	0.0000	1.0000
K2PG	57	-2023.5987	4177.0155	0.0000	1.0000
K2PI	57	-2023.9571	4177.7325	0.0000	1.0000
HKYG	60	-2021.3943	4180.4271	0.0000	1.0000
HKYI	60	-2025.3988	4188.4362	0.0000	1.0000
JC69IG	57	-2078.2828	4286.3838	0.0000	1.0000
F81IG	60	-2075.7833	4289.2052	0.0000	1.0000
F81G	59	-2079.1607	4293.3406	0.0000	1.0000
JC69G	56	-2084.2920	4295.8203	0.0000	1.0000
JC69I	56	-2084.5278	4296.2920	0.0000	1.0000
F81I	59	-2081.6684	4298.3561	0.0000	1.0000
SYM	60	-2135.2884	4408.2153	0.0000	1.0000
GTR	63	-2132.1680	4409.9088	0.0000	1.0000
K2P	56	-2146.2945	4419.8252	0.0000	1.0000
HKY	59	-2149.3483	4433.7158	0.0000	1.0000
JC69	55	-2204.9824	4534.6314	0.0000	1.0000
F81	58	-2212.5109	4557.4343	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1998.2820	4124.5641	0.9474	0.9474
GTRIG	65	-2000.2426	4130.4853	0.0491	0.9964
GTRI	64	-2003.8752	4135.7504	0.0035	1.0000
SYMI	61	-2012.0166	4146.0332	0.0000	1.0000
SYMIG	62	-2011.4461	4146.8922	0.0000	1.0000
K2PIG	58	-2015.7712	4147.5423	0.0000	1.0000
SYMG	61	-2015.7830	4153.5660	0.0000	1.0000
HKYIG	61	-2016.7760	4155.5520	0.0000	1.0000
K2PG	57	-2023.5987	4161.1973	0.0000	1.0000
K2PI	57	-2023.9571	4161.9143	0.0000	1.0000
HKYG	60	-2021.3943	4162.7886	0.0000	1.0000
HKYI	60	-2025.3988	4170.7977	0.0000	1.0000
JC69IG	57	-2078.2828	4270.5656	0.0000	1.0000
F81IG	60	-2075.7833	4271.5667	0.0000	1.0000
F81G	59	-2079.1607	4276.3214	0.0000	1.0000
JC69G	56	-2084.2920	4280.5840	0.0000	1.0000
JC69I	56	-2084.5278	4281.0557	0.0000	1.0000
F81I	59	-2081.6684	4281.3369	0.0000	1.0000
GTR	63	-2132.1680	4390.3360	0.0000	1.0000
SYM	60	-2135.2884	4390.5768	0.0000	1.0000
K2P	56	-2146.2945	4404.5889	0.0000	1.0000
HKY	59	-2149.3483	4416.6965	0.0000	1.0000
JC69	55	-2204.9824	4519.9647	0.0000	1.0000
F81	58	-2212.5109	4541.0218	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-2015.7712	4389.1366	0.7203	0.7203
GTRG	64	-1998.2820	4391.1508	0.2631	0.9834
K2PG	57	-2023.5987	4398.6261	0.0063	0.9896
K2PI	57	-2023.9571	4399.3431	0.0044	0.9940
SYMI	61	-2012.0166	4400.1236	0.0030	0.9970
GTRIG	65	-2000.2426	4401.2374	0.0017	0.9987
GTRI	64	-2003.8752	4402.3372	0.0010	0.9997
SYMIG	62	-2011.4461	4405.1481	0.0002	0.9999
SYMG	61	-2015.7830	4407.6564	0.0001	1.0000
HKYIG	61	-2016.7760	4409.6425	0.0000	1.0000
HKYG	60	-2021.3943	4412.7137	0.0000	1.0000
HKYI	60	-2025.3988	4420.7228	0.0000	1.0000
JC69IG	57	-2078.2828	4507.9945	0.0000	1.0000
JC69G	56	-2084.2920	4513.8474	0.0000	1.0000
JC69I	56	-2084.5278	4514.3191	0.0000	1.0000
F81IG	60	-2075.7833	4521.4917	0.0000	1.0000
F81G	59	-2079.1607	4522.0810	0.0000	1.0000
F81I	59	-2081.6684	4527.0965	0.0000	1.0000
K2P	56	-2146.2945	4637.8523	0.0000	1.0000
SYM	60	-2135.2884	4640.5019	0.0000	1.0000
GTR	63	-2132.1680	4652.7573	0.0000	1.0000
HKY	59	-2149.3483	4662.4562	0.0000	1.0000
JC69	55	-2204.9824	4749.0627	0.0000	1.0000
F81	58	-2212.5109	4782.6160	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/8718_hg18.phylip = 1 - 476;
 Partition Dummy = 1:261-aggressive-phylip/8718_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/8718_hg18.phylip = 1 - 476;
 Partition Dummy = 1:261-aggressive-phylip/8718_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 (K2PIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8718_hg18.phylip = 1 - 476;
 Partition Dummy = 1:261-aggressive-phylip/8718_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

