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

Input data from file "261-aggressive-phylip/6656_hg18.phylip" (29 taxa, 923 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	-5680.8558	11499.4085	0.9595	0.9595
GTRIG	65	-5682.9377	11505.8870	0.0376	0.9971
HKYG	60	-5691.3213	11511.1345	0.0027	0.9998
HKYIG	61	-5693.1296	11517.0444	0.0001	1.0000
SYMG	61	-5695.0344	11520.8539	0.0000	1.0000
K2PG	57	-5702.1299	11525.9037	0.0000	1.0000
SYMIG	62	-5706.1645	11545.4127	0.0000	1.0000
K2PIG	58	-5711.6775	11547.2762	0.0000	1.0000
GTRI	64	-5726.7887	11591.2745	0.0000	1.0000
HKYI	60	-5740.6825	11609.8570	0.0000	1.0000
SYMI	61	-5739.9125	11610.6101	0.0000	1.0000
K2PI	57	-5749.1540	11619.9520	0.0000	1.0000
F81IG	60	-5915.3748	11959.2415	0.0000	1.0000
SYM	60	-5915.8627	11960.2173	0.0000	1.0000
JC69G	56	-5920.4448	11960.2614	0.0000	1.0000
F81G	59	-5917.3949	11960.9938	0.0000	1.0000
JC69IG	57	-5920.9832	11963.6104	0.0000	1.0000
GTR	63	-5926.7820	11988.9517	0.0000	1.0000
HKY	59	-5933.5026	11993.2092	0.0000	1.0000
K2P	56	-5943.0717	12005.5152	0.0000	1.0000
F81I	59	-5951.3356	12028.8751	0.0000	1.0000
JC69I	56	-5961.8622	12043.0961	0.0000	1.0000
F81	58	-6144.8821	12413.6854	0.0000	1.0000
JC69	55	-6149.8750	12416.8551	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-5680.8558	11489.7116	0.9546	0.9546
GTRIG	65	-5682.9377	11495.8753	0.0438	0.9984
HKYG	60	-5691.3213	11502.6427	0.0015	0.9999
HKYIG	61	-5693.1296	11508.2593	0.0001	1.0000
SYMG	61	-5695.0344	11512.0687	0.0000	1.0000
K2PG	57	-5702.1299	11518.2597	0.0000	1.0000
SYMIG	62	-5706.1645	11536.3290	0.0000	1.0000
K2PIG	58	-5711.6775	11539.3549	0.0000	1.0000
GTRI	64	-5726.7887	11581.5775	0.0000	1.0000
HKYI	60	-5740.6825	11601.3651	0.0000	1.0000
SYMI	61	-5739.9125	11601.8250	0.0000	1.0000
K2PI	57	-5749.1540	11612.3080	0.0000	1.0000
F81IG	60	-5915.3748	11950.7496	0.0000	1.0000
SYM	60	-5915.8627	11951.7255	0.0000	1.0000
F81G	59	-5917.3949	11952.7899	0.0000	1.0000
JC69G	56	-5920.4448	11952.8896	0.0000	1.0000
JC69IG	57	-5920.9832	11955.9665	0.0000	1.0000
GTR	63	-5926.7820	11979.5640	0.0000	1.0000
HKY	59	-5933.5026	11985.0052	0.0000	1.0000
K2P	56	-5943.0717	11998.1434	0.0000	1.0000
F81I	59	-5951.3356	12020.6711	0.0000	1.0000
JC69I	56	-5961.8622	12035.7243	0.0000	1.0000
F81	58	-6144.8821	12405.7641	0.0000	1.0000
JC69	55	-6149.8750	12409.7501	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-5691.3213	11792.3004	0.6193	0.6193
K2PG	57	-5702.1299	11793.4346	0.3513	0.9706
GTRG	64	-5680.8558	11798.6799	0.0255	0.9961
HKYIG	61	-5693.1296	11802.7446	0.0033	0.9994
SYMG	61	-5695.0344	11806.5541	0.0005	0.9999
GTRIG	65	-5682.9377	11809.6712	0.0001	1.0000
K2PIG	58	-5711.6775	11819.3574	0.0000	1.0000
SYMIG	62	-5706.1645	11835.6420	0.0000	1.0000
K2PI	57	-5749.1540	11887.4829	0.0000	1.0000
GTRI	64	-5726.7887	11890.5458	0.0000	1.0000
HKYI	60	-5740.6825	11891.0228	0.0000	1.0000
SYMI	61	-5739.9125	11896.3103	0.0000	1.0000
JC69G	56	-5920.4448	12223.2368	0.0000	1.0000
JC69IG	57	-5920.9832	12231.1413	0.0000	1.0000
F81G	59	-5917.3949	12237.6200	0.0000	1.0000
F81IG	60	-5915.3748	12240.4074	0.0000	1.0000
SYM	60	-5915.8627	12241.3832	0.0000	1.0000
K2P	56	-5943.0717	12268.4906	0.0000	1.0000
HKY	59	-5933.5026	12269.8354	0.0000	1.0000
GTR	63	-5926.7820	12283.7047	0.0000	1.0000
F81I	59	-5951.3356	12305.5013	0.0000	1.0000
JC69I	56	-5961.8622	12306.0716	0.0000	1.0000
JC69	55	-6149.8750	12675.2697	0.0000	1.0000
F81	58	-6144.8821	12685.7666	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/6656_hg18.phylip = 1 - 923;
 Partition Dummy = 1:261-aggressive-phylip/6656_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/6656_hg18.phylip = 1 - 923;
 Partition Dummy = 1:261-aggressive-phylip/6656_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/6656_hg18.phylip = 1 - 923;
 Partition Dummy = 1:261-aggressive-phylip/6656_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

