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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRIG


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	-5427.7019	10995.6059	1.0000	1.0000
HKYIG	61	-5450.4099	11031.7712	0.0000	1.0000
K2PIG	58	-5455.7412	11035.5531	0.0000	1.0000
GTRG	64	-5449.5898	11037.0608	0.0000	1.0000
SYMIG	62	-5451.9893	11037.2345	0.0000	1.0000
HKYG	60	-5457.4483	11043.5491	0.0000	1.0000
K2PG	57	-5482.1515	11086.0909	0.0000	1.0000
SYMG	61	-5477.6845	11086.3204	0.0000	1.0000
GTRI	64	-5480.9439	11099.7690	0.0000	1.0000
HKYI	60	-5487.8445	11104.3415	0.0000	1.0000
SYMI	61	-5498.7637	11128.4789	0.0000	1.0000
K2PI	57	-5503.7489	11129.2859	0.0000	1.0000
F81IG	60	-5604.8248	11338.3021	0.0000	1.0000
F81G	59	-5612.2158	11350.7906	0.0000	1.0000
JC69IG	57	-5617.0430	11355.8739	0.0000	1.0000
JC69G	56	-5633.1799	11385.8704	0.0000	1.0000
F81I	59	-5641.3693	11409.0975	0.0000	1.0000
JC69I	56	-5663.8569	11447.2244	0.0000	1.0000
GTR	63	-5908.9672	11953.5003	0.0000	1.0000
HKY	59	-5919.2771	11964.9131	0.0000	1.0000
SYM	60	-5923.5684	11975.7892	0.0000	1.0000
K2P	56	-5931.4904	11982.4914	0.0000	1.0000
F81	58	-6049.5913	12223.2533	0.0000	1.0000
JC69	55	-6068.8836	12255.0058	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-5427.7019	10985.4037	1.0000	1.0000
HKYIG	61	-5450.4099	11022.8197	0.0000	1.0000
GTRG	64	-5449.5898	11027.1796	0.0000	1.0000
K2PIG	58	-5455.7412	11027.4823	0.0000	1.0000
SYMIG	62	-5451.9893	11027.9786	0.0000	1.0000
HKYG	60	-5457.4483	11034.8966	0.0000	1.0000
SYMG	61	-5477.6845	11077.3690	0.0000	1.0000
K2PG	57	-5482.1515	11078.3029	0.0000	1.0000
GTRI	64	-5480.9439	11089.8877	0.0000	1.0000
HKYI	60	-5487.8445	11095.6890	0.0000	1.0000
SYMI	61	-5498.7637	11119.5274	0.0000	1.0000
K2PI	57	-5503.7489	11121.4979	0.0000	1.0000
F81IG	60	-5604.8248	11329.6496	0.0000	1.0000
F81G	59	-5612.2158	11342.4317	0.0000	1.0000
JC69IG	57	-5617.0430	11348.0859	0.0000	1.0000
JC69G	56	-5633.1799	11378.3598	0.0000	1.0000
F81I	59	-5641.3693	11400.7386	0.0000	1.0000
JC69I	56	-5663.8569	11439.7138	0.0000	1.0000
GTR	63	-5908.9672	11943.9345	0.0000	1.0000
HKY	59	-5919.2771	11956.5542	0.0000	1.0000
SYM	60	-5923.5684	11967.1367	0.0000	1.0000
K2P	56	-5931.4904	11974.9809	0.0000	1.0000
F81	58	-6049.5913	12215.1825	0.0000	1.0000
JC69	55	-6068.8836	12247.7672	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-5427.7019	11298.0630	0.9852	0.9852
K2PIG	58	-5455.7412	11306.4706	0.0147	0.9999
HKYIG	61	-5450.4099	11316.2384	0.0001	1.0000
HKYG	60	-5457.4483	11323.5052	0.0000	1.0000
SYMIG	62	-5451.9893	11326.2074	0.0000	1.0000
GTRG	64	-5449.5898	11335.0287	0.0000	1.0000
K2PG	57	-5482.1515	11352.4810	0.0000	1.0000
SYMG	61	-5477.6845	11370.7876	0.0000	1.0000
HKYI	60	-5487.8445	11384.2975	0.0000	1.0000
K2PI	57	-5503.7489	11395.6760	0.0000	1.0000
GTRI	64	-5480.9439	11397.7369	0.0000	1.0000
SYMI	61	-5498.7637	11412.9461	0.0000	1.0000
F81IG	60	-5604.8248	11618.2581	0.0000	1.0000
JC69IG	57	-5617.0430	11622.2640	0.0000	1.0000
F81G	59	-5612.2158	11626.2301	0.0000	1.0000
JC69G	56	-5633.1799	11647.7278	0.0000	1.0000
F81I	59	-5641.3693	11684.5370	0.0000	1.0000
JC69I	56	-5663.8569	11709.0818	0.0000	1.0000
HKY	59	-5919.2771	12240.3526	0.0000	1.0000
K2P	56	-5931.4904	12244.3488	0.0000	1.0000
GTR	63	-5908.9672	12246.9734	0.0000	1.0000
SYM	60	-5923.5684	12255.7453	0.0000	1.0000
F81	58	-6049.5913	12494.1708	0.0000	1.0000
JC69	55	-6068.8836	12512.3251	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/4142_hg18.phylip = 1 - 907;
 Partition Dummy = 1:261-aggressive-phylip/4142_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/4142_hg18.phylip = 1 - 907;
 Partition Dummy = 1:261-aggressive-phylip/4142_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 (GTRIG)]

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


-------------------------------------------------------------
End of Output

