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

Input data from file "261-aggressive-phylip/3383_hg18.phylip" (29 taxa, 744 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	-1895.9646	3922.6467	0.4335	0.4335
GTRG	64	-1891.3601	3922.9735	0.3681	0.8016
HKYIG	61	-1895.8828	3924.8565	0.1436	0.9452
GTRIG	65	-1892.0632	3926.7812	0.0548	1.0000
K2PG	57	-1912.5104	3948.6594	0.0000	1.0000
K2PIG	58	-1912.5146	3951.0204	0.0000	1.0000
SYMG	61	-1910.1320	3953.3549	0.0000	1.0000
SYMIG	62	-1910.1406	3955.7526	0.0000	1.0000
F81G	59	-1917.1512	3962.6533	0.0000	1.0000
F81IG	60	-1917.1592	3965.0358	0.0000	1.0000
JC69G	56	-1931.6968	3984.6861	0.0000	1.0000
JC69IG	57	-1931.7051	3987.0487	0.0000	1.0000
GTR	63	-1926.7413	3991.3414	0.0000	1.0000
HKY	59	-1931.6190	3991.5888	0.0000	1.0000
GTRI	64	-1926.7461	3993.7455	0.0000	1.0000
HKYI	60	-1931.6176	3993.9527	0.0000	1.0000
K2P	56	-1948.0730	4017.4387	0.0000	1.0000
K2PI	57	-1948.0706	4019.7797	0.0000	1.0000
SYM	60	-1945.1354	4020.9883	0.0000	1.0000
SYMI	61	-1945.1403	4023.3716	0.0000	1.0000
F81	58	-1952.3069	4030.6050	0.0000	1.0000
F81I	59	-1952.3108	4032.9724	0.0000	1.0000
JC69	55	-1966.8598	4052.6731	0.0000	1.0000
JC69I	56	-1966.8633	4055.0192	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1891.3601	3910.7202	0.5137	0.5137
HKYG	60	-1895.9646	3911.9293	0.2807	0.7944
HKYIG	61	-1895.8828	3913.7656	0.1121	0.9064
GTRIG	65	-1892.0632	3914.1264	0.0936	1.0000
K2PG	57	-1912.5104	3939.0209	0.0000	1.0000
K2PIG	58	-1912.5146	3941.0291	0.0000	1.0000
SYMG	61	-1910.1320	3942.2640	0.0000	1.0000
SYMIG	62	-1910.1406	3944.2812	0.0000	1.0000
F81G	59	-1917.1512	3952.3025	0.0000	1.0000
F81IG	60	-1917.1592	3954.3184	0.0000	1.0000
JC69G	56	-1931.6968	3975.3935	0.0000	1.0000
JC69IG	57	-1931.7051	3977.4102	0.0000	1.0000
GTR	63	-1926.7413	3979.4826	0.0000	1.0000
HKY	59	-1931.6190	3981.2379	0.0000	1.0000
GTRI	64	-1926.7461	3981.4921	0.0000	1.0000
HKYI	60	-1931.6176	3983.2353	0.0000	1.0000
K2P	56	-1948.0730	4008.1461	0.0000	1.0000
K2PI	57	-1948.0706	4010.1412	0.0000	1.0000
SYM	60	-1945.1354	4010.2709	0.0000	1.0000
SYMI	61	-1945.1403	4012.2807	0.0000	1.0000
F81	58	-1952.3069	4020.6137	0.0000	1.0000
F81I	59	-1952.3108	4022.6215	0.0000	1.0000
JC69	55	-1966.8598	4043.7196	0.0000	1.0000
JC69I	56	-1966.8633	4045.7266	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1895.9646	4188.6517	0.9603	0.9603
HKYIG	61	-1895.8828	4195.1001	0.0382	0.9985
K2PG	57	-1912.5104	4201.9072	0.0013	0.9998
GTRG	64	-1891.3601	4205.8908	0.0002	1.0000
K2PIG	58	-1912.5146	4208.5275	0.0000	1.0000
GTRIG	65	-1892.0632	4213.9090	0.0000	1.0000
SYMG	61	-1910.1320	4223.5985	0.0000	1.0000
F81G	59	-1917.1512	4224.4129	0.0000	1.0000
SYMIG	62	-1910.1406	4230.2278	0.0000	1.0000
F81IG	60	-1917.1592	4231.0409	0.0000	1.0000
JC69G	56	-1931.6968	4233.6678	0.0000	1.0000
JC69IG	57	-1931.7051	4240.2966	0.0000	1.0000
HKY	59	-1931.6190	4253.3483	0.0000	1.0000
HKYI	60	-1931.6176	4259.9577	0.0000	1.0000
K2P	56	-1948.0730	4266.4204	0.0000	1.0000
GTR	63	-1926.7413	4270.0411	0.0000	1.0000
K2PI	57	-1948.0706	4273.0276	0.0000	1.0000
GTRI	64	-1926.7461	4276.6628	0.0000	1.0000
SYM	60	-1945.1354	4286.9933	0.0000	1.0000
F81	58	-1952.3069	4288.1121	0.0000	1.0000
SYMI	61	-1945.1403	4293.6152	0.0000	1.0000
F81I	59	-1952.3108	4294.7319	0.0000	1.0000
JC69	55	-1966.8598	4297.3818	0.0000	1.0000
JC69I	56	-1966.8633	4304.0009	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/3383_hg18.phylip = 1 - 744;
 Partition Dummy = 1:261-aggressive-phylip/3383_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/3383_hg18.phylip = 1 - 744;
 Partition Dummy = 1:261-aggressive-phylip/3383_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/3383_hg18.phylip = 1 - 744;
 Partition Dummy = 1:261-aggressive-phylip/3383_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

