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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-4292.0011	8721.3610	0.6397	0.6397
GTRIG	65	-4291.4234	8722.5090	0.3603	1.0000
GTRI	64	-4302.2548	8741.8684	0.0000	1.0000
HKYG	60	-4328.7212	8785.6394	0.0000	1.0000
HKYIG	61	-4328.0847	8786.6493	0.0000	1.0000
HKYI	60	-4335.8643	8799.9257	0.0000	1.0000
SYMG	61	-4372.9037	8876.2872	0.0000	1.0000
SYMIG	62	-4372.2976	8877.3628	0.0000	1.0000
SYMI	61	-4380.7868	8892.0534	0.0000	1.0000
K2PG	57	-4387.3295	8896.0384	0.0000	1.0000
K2PIG	58	-4386.6280	8896.9029	0.0000	1.0000
K2PI	57	-4395.8142	8913.0079	0.0000	1.0000
GTR	63	-4402.7851	8940.6308	0.0000	1.0000
HKY	59	-4431.4414	8988.8022	0.0000	1.0000
F81G	59	-4435.2501	8996.4197	0.0000	1.0000
F81IG	60	-4434.5764	8997.3498	0.0000	1.0000
F81I	59	-4441.6961	9009.3116	0.0000	1.0000
JC69G	56	-4474.8384	9068.7939	0.0000	1.0000
JC69IG	57	-4474.1838	9069.7470	0.0000	1.0000
SYM	60	-4473.9406	9076.0783	0.0000	1.0000
JC69I	56	-4481.6508	9082.4187	0.0000	1.0000
K2P	56	-4490.8846	9100.8863	0.0000	1.0000
F81	58	-4533.5174	9190.6816	0.0000	1.0000
JC69	55	-4575.0884	9267.0366	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-4292.0011	8712.0021	0.6040	0.6040
GTRIG	65	-4291.4234	8712.8468	0.3960	1.0000
GTRI	64	-4302.2548	8732.5095	0.0000	1.0000
HKYG	60	-4328.7212	8777.4424	0.0000	1.0000
HKYIG	61	-4328.0847	8778.1695	0.0000	1.0000
HKYI	60	-4335.8643	8791.7286	0.0000	1.0000
SYMG	61	-4372.9037	8867.8073	0.0000	1.0000
SYMIG	62	-4372.2976	8868.5951	0.0000	1.0000
SYMI	61	-4380.7868	8883.5736	0.0000	1.0000
K2PG	57	-4387.3295	8888.6589	0.0000	1.0000
K2PIG	58	-4386.6280	8889.2560	0.0000	1.0000
K2PI	57	-4395.8142	8905.6284	0.0000	1.0000
GTR	63	-4402.7851	8931.5701	0.0000	1.0000
HKY	59	-4431.4414	8980.8827	0.0000	1.0000
F81G	59	-4435.2501	8988.5003	0.0000	1.0000
F81IG	60	-4434.5764	8989.1527	0.0000	1.0000
F81I	59	-4441.6961	9001.3921	0.0000	1.0000
JC69G	56	-4474.8384	9061.6769	0.0000	1.0000
JC69IG	57	-4474.1838	9062.3676	0.0000	1.0000
SYM	60	-4473.9406	9067.8812	0.0000	1.0000
JC69I	56	-4481.6508	9075.3017	0.0000	1.0000
K2P	56	-4490.8846	9093.7693	0.0000	1.0000
F81	58	-4533.5174	9183.0347	0.0000	1.0000
JC69	55	-4575.0884	9260.1769	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-4292.0011	9023.0846	0.9454	0.9454
GTRIG	65	-4291.4234	9028.7899	0.0545	1.0000
GTRI	64	-4302.2548	9043.5920	0.0000	1.0000
HKYG	60	-4328.7212	9069.0822	0.0000	1.0000
HKYIG	61	-4328.0847	9074.6699	0.0000	1.0000
HKYI	60	-4335.8643	9083.3684	0.0000	1.0000
SYMG	61	-4372.9037	9164.3078	0.0000	1.0000
K2PG	57	-4387.3295	9165.7168	0.0000	1.0000
SYMIG	62	-4372.2976	9169.9562	0.0000	1.0000
K2PIG	58	-4386.6280	9171.1745	0.0000	1.0000
SYMI	61	-4380.7868	9180.0740	0.0000	1.0000
K2PI	57	-4395.8142	9182.6862	0.0000	1.0000
GTR	63	-4402.7851	9237.7919	0.0000	1.0000
HKY	59	-4431.4414	9267.6619	0.0000	1.0000
F81G	59	-4435.2501	9275.2794	0.0000	1.0000
F81IG	60	-4434.5764	9280.7926	0.0000	1.0000
F81I	59	-4441.6961	9288.1713	0.0000	1.0000
JC69G	56	-4474.8384	9333.8740	0.0000	1.0000
JC69IG	57	-4474.1838	9339.4254	0.0000	1.0000
JC69I	56	-4481.6508	9347.4988	0.0000	1.0000
SYM	60	-4473.9406	9359.5210	0.0000	1.0000
K2P	56	-4490.8846	9365.9664	0.0000	1.0000
F81	58	-4533.5174	9464.9532	0.0000	1.0000
JC69	55	-4575.0884	9527.5134	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/2361_hg18.phylip = 1 - 954;
 Partition Dummy = 1:261-aggressive-phylip/2361_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/2361_hg18.phylip = 1 - 954;
 Partition Dummy = 1:261-aggressive-phylip/2361_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 (GTRG)]

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

