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

Input data from file "261-aggressive-phylip/1054_hg18.phylip" (29 taxa, 1225 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	-8602.9986	17343.4001	0.9188	0.9188
GTRG	64	-8606.5402	17348.2529	0.0812	1.0000
HKYG	60	-8636.6076	17399.5038	0.0000	1.0000
HKYIG	61	-8638.0593	17404.6224	0.0000	1.0000
GTRI	64	-8651.3571	17437.8867	0.0000	1.0000
HKYI	60	-8678.6208	17483.5302	0.0000	1.0000
SYMG	61	-8720.6687	17569.8412	0.0000	1.0000
SYMIG	62	-8724.7409	17580.2047	0.0000	1.0000
K2PG	57	-8732.8262	17585.3182	0.0000	1.0000
K2PIG	58	-8732.1833	17586.2361	0.0000	1.0000
SYMI	61	-8768.0726	17664.6491	0.0000	1.0000
K2PI	57	-8779.8430	17679.3519	0.0000	1.0000
F81G	59	-8833.9735	17792.0243	0.0000	1.0000
F81IG	60	-8836.9616	17800.2119	0.0000	1.0000
F81I	59	-8875.9411	17875.9595	0.0000	1.0000
JC69G	56	-8888.5868	17894.6394	0.0000	1.0000
JC69IG	57	-8887.9727	17895.6112	0.0000	1.0000
GTR	63	-8885.2951	17903.5360	0.0000	1.0000
HKY	59	-8901.9033	17927.8839	0.0000	1.0000
JC69I	56	-8933.6199	17984.7056	0.0000	1.0000
SYM	60	-8974.6668	18075.6222	0.0000	1.0000
K2P	56	-8984.2306	18085.9270	0.0000	1.0000
F81	58	-9082.3218	18286.5132	0.0000	1.0000
JC69	55	-9133.6453	18382.5601	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-8602.9986	17335.9971	0.9270	0.9270
GTRG	64	-8606.5402	17341.0805	0.0730	1.0000
HKYG	60	-8636.6076	17393.2152	0.0000	1.0000
HKYIG	61	-8638.0593	17398.1185	0.0000	1.0000
GTRI	64	-8651.3571	17430.7143	0.0000	1.0000
HKYI	60	-8678.6208	17477.2415	0.0000	1.0000
SYMG	61	-8720.6687	17563.3373	0.0000	1.0000
SYMIG	62	-8724.7409	17573.4818	0.0000	1.0000
K2PG	57	-8732.8262	17579.6524	0.0000	1.0000
K2PIG	58	-8732.1833	17580.3665	0.0000	1.0000
SYMI	61	-8768.0726	17658.1453	0.0000	1.0000
K2PI	57	-8779.8430	17673.6861	0.0000	1.0000
F81G	59	-8833.9735	17785.9471	0.0000	1.0000
F81IG	60	-8836.9616	17793.9232	0.0000	1.0000
F81I	59	-8875.9411	17869.8822	0.0000	1.0000
JC69G	56	-8888.5868	17889.1736	0.0000	1.0000
JC69IG	57	-8887.9727	17889.9454	0.0000	1.0000
GTR	63	-8885.2951	17896.5903	0.0000	1.0000
HKY	59	-8901.9033	17921.8066	0.0000	1.0000
JC69I	56	-8933.6199	17979.2399	0.0000	1.0000
SYM	60	-8974.6668	18069.3336	0.0000	1.0000
K2P	56	-8984.2306	18080.4612	0.0000	1.0000
F81	58	-9082.3218	18280.6436	0.0000	1.0000
JC69	55	-9133.6453	18377.2907	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-8606.5402	17668.1650	0.5034	0.5034
GTRIG	65	-8602.9986	17668.1924	0.4966	1.0000
HKYG	60	-8636.6076	17699.8569	0.0000	1.0000
HKYIG	61	-8638.0593	17709.8710	0.0000	1.0000
GTRI	64	-8651.3571	17757.7988	0.0000	1.0000
HKYI	60	-8678.6208	17783.8833	0.0000	1.0000
K2PG	57	-8732.8262	17870.9621	0.0000	1.0000
SYMG	61	-8720.6687	17875.0898	0.0000	1.0000
K2PIG	58	-8732.1833	17876.7869	0.0000	1.0000
SYMIG	62	-8724.7409	17890.3449	0.0000	1.0000
K2PI	57	-8779.8430	17964.9957	0.0000	1.0000
SYMI	61	-8768.0726	17969.8977	0.0000	1.0000
F81G	59	-8833.9735	18087.4782	0.0000	1.0000
F81IG	60	-8836.9616	18100.5650	0.0000	1.0000
F81I	59	-8875.9411	18171.4133	0.0000	1.0000
JC69G	56	-8888.5868	18175.3726	0.0000	1.0000
JC69IG	57	-8887.9727	18181.2551	0.0000	1.0000
GTR	63	-8885.2951	18218.5641	0.0000	1.0000
HKY	59	-8901.9033	18223.3377	0.0000	1.0000
JC69I	56	-8933.6199	18265.4388	0.0000	1.0000
K2P	56	-8984.2306	18366.6602	0.0000	1.0000
SYM	60	-8974.6668	18375.9753	0.0000	1.0000
F81	58	-9082.3218	18577.0640	0.0000	1.0000
JC69	55	-9133.6453	18658.3789	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/1054_hg18.phylip = 1 - 1225;
 Partition Dummy = 1:261-aggressive-phylip/1054_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/1054_hg18.phylip = 1 - 1225;
 Partition Dummy = 1:261-aggressive-phylip/1054_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/1054_hg18.phylip = 1 - 1225;
 Partition Dummy = 1:261-aggressive-phylip/1054_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

