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

Input data from file "261-aggressive-phylip/12530_hg18.phylip" (29 taxa, 709 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	-2153.1258	4447.1708	0.4211	0.4211
GTRI	64	-2153.2119	4447.3432	0.3863	0.8074
GTRIG	65	-2153.0588	4449.4613	0.1340	0.9414
HKYG	60	-2160.6803	4452.6570	0.0271	0.9685
HKYI	60	-2160.8361	4452.9686	0.0232	0.9917
HKYIG	61	-2160.6703	4455.0315	0.0083	1.0000
GTR	63	-2168.0999	4474.7021	0.0000	1.0000
HKY	59	-2174.9600	4478.8290	0.0000	1.0000
F81G	59	-2192.3064	4513.5220	0.0000	1.0000
F81I	59	-2192.3287	4513.5664	0.0000	1.0000
F81IG	60	-2192.2414	4515.7790	0.0000	1.0000
SYMG	61	-2193.4852	4520.6612	0.0000	1.0000
K2PG	57	-2198.3855	4520.9277	0.0000	1.0000
SYMI	61	-2193.7879	4521.2666	0.0000	1.0000
K2PI	57	-2198.6071	4521.3709	0.0000	1.0000
SYMIG	62	-2193.4644	4523.0217	0.0000	1.0000
K2PIG	58	-2198.3797	4523.2887	0.0000	1.0000
F81	58	-2205.7931	4538.1155	0.0000	1.0000
SYM	60	-2207.6191	4546.5344	0.0000	1.0000
K2P	56	-2212.6423	4547.0761	0.0000	1.0000
JC69G	56	-2226.1755	4574.1424	0.0000	1.0000
JC69I	56	-2226.2801	4574.3517	0.0000	1.0000
JC69IG	57	-2226.1620	4576.4806	0.0000	1.0000
JC69	55	-2239.5894	4598.6123	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2153.1258	4434.2516	0.4211	0.4211
GTRI	64	-2153.2119	4434.4239	0.3864	0.8075
GTRIG	65	-2153.0588	4436.1176	0.1657	0.9732
HKYG	60	-2160.6803	4441.3607	0.0120	0.9852
HKYI	60	-2160.8361	4441.6723	0.0103	0.9955
HKYIG	61	-2160.6703	4443.3406	0.0045	1.0000
GTR	63	-2168.0999	4462.1998	0.0000	1.0000
HKY	59	-2174.9600	4467.9199	0.0000	1.0000
F81G	59	-2192.3064	4502.6129	0.0000	1.0000
F81I	59	-2192.3287	4502.6573	0.0000	1.0000
F81IG	60	-2192.2414	4504.4827	0.0000	1.0000
SYMG	61	-2193.4852	4508.9703	0.0000	1.0000
SYMI	61	-2193.7879	4509.5758	0.0000	1.0000
K2PG	57	-2198.3855	4510.7710	0.0000	1.0000
SYMIG	62	-2193.4644	4510.9288	0.0000	1.0000
K2PI	57	-2198.6071	4511.2142	0.0000	1.0000
K2PIG	58	-2198.3797	4512.7595	0.0000	1.0000
F81	58	-2205.7931	4527.5862	0.0000	1.0000
SYM	60	-2207.6191	4535.2381	0.0000	1.0000
K2P	56	-2212.6423	4537.2847	0.0000	1.0000
JC69G	56	-2226.1755	4564.3510	0.0000	1.0000
JC69I	56	-2226.2801	4564.5603	0.0000	1.0000
JC69IG	57	-2226.1620	4566.3239	0.0000	1.0000
JC69	55	-2239.5894	4589.1789	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2160.6803	4715.1920	0.5260	0.5260
HKYI	60	-2160.8361	4715.5036	0.4501	0.9761
HKYIG	61	-2160.6703	4721.7358	0.0200	0.9961
GTRG	64	-2153.1258	4726.3383	0.0020	0.9981
GTRI	64	-2153.2119	4726.5107	0.0018	0.9999
GTRIG	65	-2153.0588	4732.7682	0.0001	1.0000
HKY	59	-2174.9600	4737.1874	0.0000	1.0000
GTR	63	-2168.0999	4749.7227	0.0000	1.0000
K2PG	57	-2198.3855	4770.9107	0.0000	1.0000
K2PI	57	-2198.6071	4771.3540	0.0000	1.0000
F81G	59	-2192.3064	4771.8804	0.0000	1.0000
F81I	59	-2192.3287	4771.9248	0.0000	1.0000
K2PIG	58	-2198.3797	4777.4631	0.0000	1.0000
F81IG	60	-2192.2414	4778.3141	0.0000	1.0000
SYMG	61	-2193.4852	4787.3655	0.0000	1.0000
SYMI	61	-2193.7879	4787.9709	0.0000	1.0000
F81	58	-2205.7931	4792.2898	0.0000	1.0000
K2P	56	-2212.6423	4792.8606	0.0000	1.0000
SYMIG	62	-2193.4644	4793.8879	0.0000	1.0000
SYM	60	-2207.6191	4809.0695	0.0000	1.0000
JC69G	56	-2226.1755	4819.9269	0.0000	1.0000
JC69I	56	-2226.2801	4820.1362	0.0000	1.0000
JC69IG	57	-2226.1620	4826.4637	0.0000	1.0000
JC69	55	-2239.5894	4840.1910	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/12530_hg18.phylip = 1 - 709;
 Partition Dummy = 1:261-aggressive-phylip/12530_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/12530_hg18.phylip = 1 - 709;
 Partition Dummy = 1:261-aggressive-phylip/12530_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/12530_hg18.phylip = 1 - 709;
 Partition Dummy = 1:261-aggressive-phylip/12530_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

