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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-2157.5741	4454.8787	0.7428	0.7428
GTRG	64	-2153.4022	4457.4747	0.2029	0.9457
GTRIG	65	-2153.4025	4460.2475	0.0507	0.9964
HKYIG	61	-2161.6169	4465.6770	0.0034	0.9998
K2PG	57	-2169.6533	4470.9858	0.0002	1.0000
K2PIG	58	-2173.5414	4481.4314	0.0000	1.0000
GTRI	64	-2166.7623	4484.1949	0.0000	1.0000
SYMG	61	-2172.0148	4486.4729	0.0000	1.0000
HKYI	60	-2173.5176	4486.7656	0.0000	1.0000
SYMIG	62	-2171.9893	4489.1493	0.0000	1.0000
K2PI	57	-2185.5541	4502.7874	0.0000	1.0000
SYMI	61	-2185.8074	4514.0581	0.0000	1.0000
GTR	63	-2199.7148	4547.3427	0.0000	1.0000
HKY	59	-2205.3753	4547.7828	0.0000	1.0000
K2P	56	-2216.9937	4563.0114	0.0000	1.0000
SYM	60	-2220.7764	4581.2832	0.0000	1.0000
F81G	59	-2223.7743	4584.5808	0.0000	1.0000
JC69G	56	-2230.0381	4589.1002	0.0000	1.0000
JC69IG	57	-2230.0401	4591.7594	0.0000	1.0000
F81IG	60	-2227.1925	4594.1155	0.0000	1.0000
F81I	59	-2237.3771	4611.7864	0.0000	1.0000
JC69I	56	-2244.1120	4617.2481	0.0000	1.0000
F81	58	-2268.6630	4671.6744	0.0000	1.0000
JC69	55	-2274.9862	4676.3555	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2153.4022	4434.8044	0.4514	0.4514
HKYG	60	-2157.5741	4435.1482	0.3801	0.8315
GTRIG	65	-2153.4025	4436.8049	0.1660	0.9975
HKYIG	61	-2161.6169	4445.2338	0.0025	1.0000
K2PG	57	-2169.6533	4453.3066	0.0000	1.0000
GTRI	64	-2166.7623	4461.5246	0.0000	1.0000
K2PIG	58	-2173.5414	4463.0828	0.0000	1.0000
SYMG	61	-2172.0148	4466.0296	0.0000	1.0000
HKYI	60	-2173.5176	4467.0352	0.0000	1.0000
SYMIG	62	-2171.9893	4467.9785	0.0000	1.0000
K2PI	57	-2185.5541	4485.1083	0.0000	1.0000
SYMI	61	-2185.8074	4493.6149	0.0000	1.0000
GTR	63	-2199.7148	4525.4297	0.0000	1.0000
HKY	59	-2205.3753	4528.7505	0.0000	1.0000
K2P	56	-2216.9937	4545.9874	0.0000	1.0000
SYM	60	-2220.7764	4561.5527	0.0000	1.0000
F81G	59	-2223.7743	4565.5486	0.0000	1.0000
JC69G	56	-2230.0381	4572.0762	0.0000	1.0000
JC69IG	57	-2230.0401	4574.0803	0.0000	1.0000
F81IG	60	-2227.1925	4574.3850	0.0000	1.0000
F81I	59	-2237.3771	4592.7542	0.0000	1.0000
JC69I	56	-2244.1120	4600.2241	0.0000	1.0000
F81	58	-2268.6630	4653.3259	0.0000	1.0000
JC69	55	-2274.9862	4659.9725	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2157.5741	4679.2538	0.9504	0.9504
K2PG	57	-2169.6533	4685.2069	0.0484	0.9988
HKYIG	61	-2161.6169	4693.4077	0.0008	0.9996
GTRG	64	-2153.4022	4695.1837	0.0003	0.9999
K2PIG	58	-2173.5414	4699.0515	0.0000	1.0000
GTRIG	65	-2153.4025	4701.2526	0.0000	1.0000
HKYI	60	-2173.5176	4711.1407	0.0000	1.0000
SYMG	61	-2172.0148	4714.2036	0.0000	1.0000
K2PI	57	-2185.5541	4717.0086	0.0000	1.0000
SYMIG	62	-2171.9893	4720.2209	0.0000	1.0000
GTRI	64	-2166.7623	4721.9038	0.0000	1.0000
SYMI	61	-2185.8074	4741.7888	0.0000	1.0000
HKY	59	-2205.3753	4768.7876	0.0000	1.0000
K2P	56	-2216.9937	4773.8192	0.0000	1.0000
GTR	63	-2199.7148	4781.7405	0.0000	1.0000
JC69G	56	-2230.0381	4799.9081	0.0000	1.0000
F81G	59	-2223.7743	4805.5857	0.0000	1.0000
SYM	60	-2220.7764	4805.6582	0.0000	1.0000
JC69IG	57	-2230.0401	4805.9805	0.0000	1.0000
F81IG	60	-2227.1925	4818.4906	0.0000	1.0000
JC69I	56	-2244.1120	4828.0559	0.0000	1.0000
F81I	59	-2237.3771	4832.7913	0.0000	1.0000
JC69	55	-2274.9862	4883.7359	0.0000	1.0000
F81	58	-2268.6630	4889.2946	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/7192_hg18.phylip = 1 - 432;
 Partition Dummy = 1:261-aggressive-phylip/7192_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 (HKYG)]

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


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

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

