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

Input data from file "261-aggressive-phylip/973_hg18.phylip" (29 taxa, 613 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	-2176.5221	4496.2266	0.6466	0.6466
GTRIG	65	-2176.0291	4497.7437	0.3028	0.9494
HKYG	60	-2184.4339	4502.1287	0.0338	0.9832
HKYIG	61	-2184.0830	4503.8937	0.0140	0.9972
GTRI	64	-2182.0611	4507.3047	0.0025	0.9997
HKYI	60	-2189.2137	4511.6883	0.0003	1.0000
K2PG	57	-2196.7322	4519.3779	0.0000	1.0000
SYMG	61	-2192.7082	4521.1442	0.0000	1.0000
K2PIG	58	-2196.4080	4521.1698	0.0000	1.0000
SYMIG	62	-2192.3386	4522.8809	0.0000	1.0000
K2PI	57	-2201.3807	4528.6749	0.0000	1.0000
SYMI	61	-2197.7270	4531.1818	0.0000	1.0000
F81IG	60	-2213.2563	4559.7734	0.0000	1.0000
F81G	59	-2216.1219	4563.0468	0.0000	1.0000
F81I	59	-2220.8757	4572.5542	0.0000	1.0000
JC69G	56	-2227.4661	4578.4142	0.0000	1.0000
GTR	63	-2219.9260	4580.5405	0.0000	1.0000
HKY	59	-2225.1063	4581.0155	0.0000	1.0000
JC69IG	57	-2227.6513	4581.2161	0.0000	1.0000
JC69I	56	-2229.8087	4583.0994	0.0000	1.0000
K2P	56	-2236.7857	4597.0534	0.0000	1.0000
SYM	60	-2233.4929	4600.2467	0.0000	1.0000
F81	58	-2255.6137	4639.5812	0.0000	1.0000
JC69	55	-2266.7002	4654.4597	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2176.5221	4481.0441	0.6109	0.6109
GTRIG	65	-2176.0291	4482.0581	0.3680	0.9789
HKYG	60	-2184.4339	4488.8678	0.0122	0.9911
HKYIG	61	-2184.0830	4490.1659	0.0064	0.9975
GTRI	64	-2182.0611	4492.1222	0.0024	0.9999
HKYI	60	-2189.2137	4498.4275	0.0001	1.0000
SYMG	61	-2192.7082	4507.4164	0.0000	1.0000
K2PG	57	-2196.7322	4507.4644	0.0000	1.0000
SYMIG	62	-2192.3386	4508.6773	0.0000	1.0000
K2PIG	58	-2196.4080	4508.8160	0.0000	1.0000
K2PI	57	-2201.3807	4516.7614	0.0000	1.0000
SYMI	61	-2197.7270	4517.4541	0.0000	1.0000
F81IG	60	-2213.2563	4546.5125	0.0000	1.0000
F81G	59	-2216.1219	4550.2439	0.0000	1.0000
F81I	59	-2220.8757	4559.7513	0.0000	1.0000
GTR	63	-2219.9260	4565.8519	0.0000	1.0000
JC69G	56	-2227.4661	4566.9322	0.0000	1.0000
HKY	59	-2225.1063	4568.2126	0.0000	1.0000
JC69IG	57	-2227.6513	4569.3025	0.0000	1.0000
JC69I	56	-2229.8087	4571.6174	0.0000	1.0000
K2P	56	-2236.7857	4585.5714	0.0000	1.0000
SYM	60	-2233.4929	4586.9858	0.0000	1.0000
F81	58	-2255.6137	4627.2274	0.0000	1.0000
JC69	55	-2266.7002	4643.4005	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2184.4339	4753.9697	0.8716	0.8716
K2PG	57	-2196.7322	4759.3112	0.0603	0.9319
HKYIG	61	-2184.0830	4759.6862	0.0500	0.9819
HKYI	60	-2189.2137	4763.5294	0.0073	0.9893
GTRG	64	-2176.5221	4763.8195	0.0063	0.9956
K2PIG	58	-2196.4080	4765.0811	0.0034	0.9990
K2PI	57	-2201.3807	4768.6082	0.0006	0.9995
GTRIG	65	-2176.0291	4769.2518	0.0004	1.0000
GTRI	64	-2182.0611	4774.8975	0.0000	1.0000
SYMG	61	-2192.7082	4776.9367	0.0000	1.0000
SYMIG	62	-2192.3386	4782.6159	0.0000	1.0000
SYMI	61	-2197.7270	4786.9743	0.0000	1.0000
F81G	59	-2216.1219	4810.9274	0.0000	1.0000
F81IG	60	-2213.2563	4811.6144	0.0000	1.0000
JC69G	56	-2227.4661	4814.3606	0.0000	1.0000
JC69I	56	-2229.8087	4819.0458	0.0000	1.0000
F81I	59	-2220.8757	4820.4349	0.0000	1.0000
JC69IG	57	-2227.6513	4821.1493	0.0000	1.0000
HKY	59	-2225.1063	4828.8962	0.0000	1.0000
K2P	56	-2236.7857	4832.9998	0.0000	1.0000
GTR	63	-2219.9260	4844.2089	0.0000	1.0000
SYM	60	-2233.4929	4852.0877	0.0000	1.0000
F81	58	-2255.6137	4883.4926	0.0000	1.0000
JC69	55	-2266.7002	4886.4106	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/973_hg18.phylip = 1 - 613;
 Partition Dummy = 1:261-aggressive-phylip/973_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/973_hg18.phylip = 1 - 613;
 Partition Dummy = 1:261-aggressive-phylip/973_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/973_hg18.phylip = 1 - 613;
 Partition Dummy = 1:261-aggressive-phylip/973_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

