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

Input data from file "261-aggressive-phylip/9035_hg18.phylip" (29 taxa, 478 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	-1913.8494	3975.8440	0.7332	0.7332
GTRIG	65	-1913.7673	3978.3599	0.2084	0.9416
HKYG	60	-1921.7992	3981.1523	0.0516	0.9932
HKYIG	61	-1923.2668	3986.7163	0.0032	0.9963
GTRI	64	-1919.4577	3987.0606	0.0027	0.9990
HKYI	60	-1925.7762	3989.1064	0.0010	1.0000
K2PG	57	-1943.0903	4015.9235	0.0000	1.0000
K2PIG	58	-1943.0482	4018.4305	0.0000	1.0000
K2PI	57	-1946.4058	4022.5544	0.0000	1.0000
SYMG	61	-1941.3725	4022.9277	0.0000	1.0000
SYMIG	62	-1941.3301	4025.4843	0.0000	1.0000
HKY	59	-1946.9625	4028.8627	0.0000	1.0000
GTR	63	-1942.0563	4029.5909	0.0000	1.0000
SYMI	61	-1944.9511	4030.0848	0.0000	1.0000
F81G	59	-1959.2683	4053.4743	0.0000	1.0000
F81IG	60	-1960.7782	4059.1103	0.0000	1.0000
K2P	56	-1965.9828	4059.1295	0.0000	1.0000
F81I	59	-1962.8675	4060.6729	0.0000	1.0000
SYM	60	-1964.6510	4066.8560	0.0000	1.0000
JC69G	56	-1973.6148	4074.3936	0.0000	1.0000
JC69IG	57	-1974.8039	4079.3507	0.0000	1.0000
JC69I	56	-1978.0982	4083.3604	0.0000	1.0000
F81	58	-1983.2579	4098.8499	0.0000	1.0000
JC69	55	-1997.3103	4119.2177	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1913.8494	3955.6987	0.7021	0.7021
GTRIG	65	-1913.7673	3957.5347	0.2804	0.9825
HKYG	60	-1921.7992	3963.5983	0.0135	0.9960
GTRI	64	-1919.4577	3966.9153	0.0026	0.9986
HKYIG	61	-1923.2668	3968.5336	0.0011	0.9997
HKYI	60	-1925.7762	3971.5524	0.0003	1.0000
K2PG	57	-1943.0903	4000.1807	0.0000	1.0000
K2PIG	58	-1943.0482	4002.0964	0.0000	1.0000
SYMG	61	-1941.3725	4004.7450	0.0000	1.0000
SYMIG	62	-1941.3301	4006.6602	0.0000	1.0000
K2PI	57	-1946.4058	4006.8115	0.0000	1.0000
GTR	63	-1942.0563	4010.1127	0.0000	1.0000
SYMI	61	-1944.9511	4011.9021	0.0000	1.0000
HKY	59	-1946.9625	4011.9249	0.0000	1.0000
F81G	59	-1959.2683	4036.5365	0.0000	1.0000
F81IG	60	-1960.7782	4041.5563	0.0000	1.0000
F81I	59	-1962.8675	4043.7351	0.0000	1.0000
K2P	56	-1965.9828	4043.9656	0.0000	1.0000
SYM	60	-1964.6510	4049.3021	0.0000	1.0000
JC69G	56	-1973.6148	4059.2297	0.0000	1.0000
JC69IG	57	-1974.8039	4063.6078	0.0000	1.0000
JC69I	56	-1978.0982	4068.1965	0.0000	1.0000
F81	58	-1983.2579	4082.5158	0.0000	1.0000
JC69	55	-1997.3103	4104.6206	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1921.7992	4213.7750	0.9594	0.9594
HKYI	60	-1925.7762	4221.7290	0.0180	0.9773
GTRG	64	-1913.8494	4222.5538	0.0119	0.9892
HKYIG	61	-1923.2668	4222.8799	0.0101	0.9994
GTRIG	65	-1913.7673	4228.5594	0.0006	1.0000
GTRI	64	-1919.4577	4233.7704	0.0000	1.0000
K2PG	57	-1943.0903	4237.8485	0.0000	1.0000
K2PIG	58	-1943.0482	4243.9338	0.0000	1.0000
K2PI	57	-1946.4058	4244.4794	0.0000	1.0000
HKY	59	-1946.9625	4257.9320	0.0000	1.0000
SYMG	61	-1941.3725	4259.0913	0.0000	1.0000
SYMIG	62	-1941.3301	4265.1761	0.0000	1.0000
SYMI	61	-1944.9511	4266.2484	0.0000	1.0000
GTR	63	-1942.0563	4272.7981	0.0000	1.0000
K2P	56	-1965.9828	4277.4638	0.0000	1.0000
F81G	59	-1959.2683	4282.5436	0.0000	1.0000
F81I	59	-1962.8675	4289.7421	0.0000	1.0000
F81IG	60	-1960.7782	4291.7330	0.0000	1.0000
JC69G	56	-1973.6148	4292.7279	0.0000	1.0000
SYM	60	-1964.6510	4299.4787	0.0000	1.0000
JC69IG	57	-1974.8039	4301.2756	0.0000	1.0000
JC69I	56	-1978.0982	4301.6947	0.0000	1.0000
F81	58	-1983.2579	4324.3532	0.0000	1.0000
JC69	55	-1997.3103	4333.9492	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/9035_hg18.phylip = 1 - 478;
 Partition Dummy = 1:261-aggressive-phylip/9035_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/9035_hg18.phylip = 1 - 478;
 Partition Dummy = 1:261-aggressive-phylip/9035_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/9035_hg18.phylip = 1 - 478;
 Partition Dummy = 1:261-aggressive-phylip/9035_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

