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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


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
GTRIG	65	-2227.3959	4601.7483	0.6854	0.6854
GTRG	64	-2229.4480	4603.3064	0.3145	0.9998
GTRI	64	-2237.3639	4619.1380	0.0001	1.0000
HKYIG	61	-2242.9594	4622.7503	0.0000	1.0000
HKYG	60	-2244.5396	4623.4041	0.0000	1.0000
SYMIG	62	-2243.6150	4626.5778	0.0000	1.0000
SYMG	61	-2245.4386	4627.7086	0.0000	1.0000
K2PIG	58	-2253.5182	4636.3775	0.0000	1.0000
K2PG	57	-2255.0699	4637.0037	0.0000	1.0000
HKYI	60	-2251.4194	4637.1636	0.0000	1.0000
SYMI	61	-2252.5814	4641.9942	0.0000	1.0000
K2PI	57	-2261.8312	4650.5262	0.0000	1.0000
GTR	63	-2295.2860	4732.4461	0.0000	1.0000
HKY	59	-2306.4681	4744.7643	0.0000	1.0000
SYM	60	-2308.2861	4750.8970	0.0000	1.0000
K2P	56	-2316.8270	4758.0501	0.0000	1.0000
F81IG	60	-2315.5367	4765.3982	0.0000	1.0000
F81G	59	-2316.9738	4765.7757	0.0000	1.0000
JC69IG	57	-2324.2528	4775.3694	0.0000	1.0000
JC69G	56	-2325.7250	4775.8461	0.0000	1.0000
F81I	59	-2323.4629	4778.7538	0.0000	1.0000
JC69I	56	-2332.3704	4789.1369	0.0000	1.0000
F81	58	-2376.3333	4882.0078	0.0000	1.0000
JC69	55	-2385.0508	4892.0396	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2227.3959	4584.7917	0.7411	0.7411
GTRG	64	-2229.4480	4586.8961	0.2588	0.9999
GTRI	64	-2237.3639	4602.7277	0.0001	1.0000
HKYIG	61	-2242.9594	4607.9189	0.0000	1.0000
HKYG	60	-2244.5396	4609.0793	0.0000	1.0000
SYMIG	62	-2243.6150	4611.2300	0.0000	1.0000
SYMG	61	-2245.4386	4612.8772	0.0000	1.0000
HKYI	60	-2251.4194	4622.8387	0.0000	1.0000
K2PIG	58	-2253.5182	4623.0364	0.0000	1.0000
K2PG	57	-2255.0699	4624.1399	0.0000	1.0000
SYMI	61	-2252.5814	4627.1628	0.0000	1.0000
K2PI	57	-2261.8312	4637.6624	0.0000	1.0000
GTR	63	-2295.2860	4716.5721	0.0000	1.0000
HKY	59	-2306.4681	4730.9362	0.0000	1.0000
SYM	60	-2308.2861	4736.5721	0.0000	1.0000
K2P	56	-2316.8270	4745.6540	0.0000	1.0000
F81IG	60	-2315.5367	4751.0733	0.0000	1.0000
F81G	59	-2316.9738	4751.9475	0.0000	1.0000
JC69IG	57	-2324.2528	4762.5056	0.0000	1.0000
JC69G	56	-2325.7250	4763.4500	0.0000	1.0000
F81I	59	-2323.4629	4764.9257	0.0000	1.0000
JC69I	56	-2332.3704	4776.7408	0.0000	1.0000
F81	58	-2376.3333	4868.6666	0.0000	1.0000
JC69	55	-2385.0508	4880.1016	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2229.4480	4865.2410	0.6767	0.6767
GTRIG	65	-2227.3959	4867.4858	0.2203	0.8969
HKYG	60	-2244.5396	4870.0276	0.0618	0.9587
K2PG	57	-2255.0699	4872.0408	0.0226	0.9813
HKYIG	61	-2242.9594	4873.2164	0.0125	0.9939
K2PIG	58	-2253.5182	4875.2865	0.0045	0.9983
SYMG	61	-2245.4386	4878.1747	0.0011	0.9994
SYMIG	62	-2243.6150	4880.8767	0.0003	0.9997
GTRI	64	-2237.3639	4881.0726	0.0002	0.9999
HKYI	60	-2251.4194	4883.7871	0.0001	1.0000
K2PI	57	-2261.8312	4885.5633	0.0000	1.0000
SYMI	61	-2252.5814	4892.4603	0.0000	1.0000
HKY	59	-2306.4681	4987.5354	0.0000	1.0000
K2P	56	-2316.8270	4989.2057	0.0000	1.0000
GTR	63	-2295.2860	4990.5678	0.0000	1.0000
SYM	60	-2308.2861	4997.5205	0.0000	1.0000
JC69G	56	-2325.7250	5007.0018	0.0000	1.0000
F81G	59	-2316.9738	5008.5467	0.0000	1.0000
JC69IG	57	-2324.2528	5010.4065	0.0000	1.0000
F81IG	60	-2315.5367	5012.0217	0.0000	1.0000
JC69I	56	-2332.3704	5020.2926	0.0000	1.0000
F81I	59	-2323.4629	5021.5249	0.0000	1.0000
JC69	55	-2385.0508	5119.3042	0.0000	1.0000
F81	58	-2376.3333	5120.9167	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/2097_hg18.phylip = 1 - 572;
 Partition Dummy = 1:261-aggressive-phylip/2097_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 (GTRIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2097_hg18.phylip = 1 - 572;
 Partition Dummy = 1:261-aggressive-phylip/2097_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 BIC model (GTRG)]

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


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

