
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:37 04/04/2011

Input data from file "591-tomodel-phylip/chr2_21409.phylip" (7 taxa, 638 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	20	-1945.7577	3932.8768	0.5347	0.5347
GTRI	20	-1946.4379	3934.2373	0.2708	0.8055
GTRIG	21	-1945.7404	3934.9809	0.1867	0.9923
HKYG	16	-1954.8584	3942.5928	0.0042	0.9964
HKYI	16	-1955.5260	3943.9279	0.0021	0.9985
HKYIG	17	-1954.8434	3944.6739	0.0015	1.0000
GTR	19	-1962.4657	3964.1612	0.0000	1.0000
HKY	15	-1971.9893	3974.7503	0.0000	1.0000
SYMG	17	-1973.1778	3981.3428	0.0000	1.0000
SYMI	17	-1973.6671	3982.3213	0.0000	1.0000
SYMIG	18	-1973.1627	3983.4304	0.0000	1.0000
K2PG	13	-1980.0801	3986.7434	0.0000	1.0000
K2PI	13	-1980.6090	3987.8014	0.0000	1.0000
K2PIG	14	-1980.0630	3988.8003	0.0000	1.0000
SYM	16	-1991.6712	4016.2183	0.0000	1.0000
K2P	12	-1998.5085	4021.5162	0.0000	1.0000
F81G	15	-2003.1230	4037.0176	0.0000	1.0000
F81I	15	-2003.2721	4037.3159	0.0000	1.0000
F81IG	16	-2003.1194	4039.1149	0.0000	1.0000
F81	14	-2018.5404	4065.7551	0.0000	1.0000
JC69G	12	-2024.6995	4073.8982	0.0000	1.0000
JC69I	12	-2024.9057	4074.3106	0.0000	1.0000
JC69IG	13	-2024.6931	4075.9695	0.0000	1.0000
JC69	11	-2041.2763	4104.9743	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1945.7577	3931.5153	0.5285	0.5285
GTRI	20	-1946.4379	3932.8758	0.2677	0.7961
GTRIG	21	-1945.7404	3933.4809	0.1978	0.9939
HKYG	16	-1954.8584	3941.7168	0.0032	0.9971
HKYI	16	-1955.5260	3943.0519	0.0017	0.9988
HKYIG	17	-1954.8434	3943.6868	0.0012	1.0000
GTR	19	-1962.4657	3962.9314	0.0000	1.0000
HKY	15	-1971.9893	3973.9786	0.0000	1.0000
SYMG	17	-1973.1778	3980.3557	0.0000	1.0000
SYMI	17	-1973.6671	3981.3342	0.0000	1.0000
SYMIG	18	-1973.1627	3982.3254	0.0000	1.0000
K2PG	13	-1980.0801	3986.1601	0.0000	1.0000
K2PI	13	-1980.6090	3987.2181	0.0000	1.0000
K2PIG	14	-1980.0630	3988.1261	0.0000	1.0000
SYM	16	-1991.6712	4015.3423	0.0000	1.0000
K2P	12	-1998.5085	4021.0170	0.0000	1.0000
F81G	15	-2003.1230	4036.2459	0.0000	1.0000
F81I	15	-2003.2721	4036.5442	0.0000	1.0000
F81IG	16	-2003.1194	4038.2389	0.0000	1.0000
F81	14	-2018.5404	4065.0809	0.0000	1.0000
JC69G	12	-2024.6995	4073.3990	0.0000	1.0000
JC69I	12	-2024.9057	4073.8114	0.0000	1.0000
JC69IG	13	-2024.6931	4075.3862	0.0000	1.0000
JC69	11	-2041.2763	4104.5526	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1954.8584	4013.0503	0.6300	0.6300
HKYI	16	-1955.5260	4014.3853	0.3232	0.9532
HKYIG	17	-1954.8434	4019.4785	0.0253	0.9785
GTRG	20	-1945.7577	4020.6821	0.0139	0.9924
GTRI	20	-1946.4379	4022.0426	0.0070	0.9994
GTRIG	21	-1945.7404	4027.1060	0.0006	1.0000
HKY	15	-1971.9893	4040.8536	0.0000	1.0000
K2PG	13	-1980.0801	4044.1185	0.0000	1.0000
K2PI	13	-1980.6090	4045.1765	0.0000	1.0000
GTR	19	-1962.4657	4047.6398	0.0000	1.0000
K2PIG	14	-1980.0630	4050.5428	0.0000	1.0000
SYMG	17	-1973.1778	4056.1474	0.0000	1.0000
SYMI	17	-1973.6671	4057.1260	0.0000	1.0000
SYMIG	18	-1973.1627	4062.5755	0.0000	1.0000
K2P	12	-1998.5085	4074.5170	0.0000	1.0000
SYM	16	-1991.6712	4086.6758	0.0000	1.0000
F81G	15	-2003.1230	4103.1210	0.0000	1.0000
F81I	15	-2003.2721	4103.4193	0.0000	1.0000
F81IG	16	-2003.1194	4109.5723	0.0000	1.0000
JC69G	12	-2024.6995	4126.8991	0.0000	1.0000
JC69I	12	-2024.9057	4127.3114	0.0000	1.0000
F81	14	-2018.5404	4127.4976	0.0000	1.0000
JC69IG	13	-2024.6931	4133.3446	0.0000	1.0000
JC69	11	-2041.2763	4153.5943	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 591-tomodel-phylip/chr2_21409.phylip = 1 - 638;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21409.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 591-tomodel-phylip/chr2_21409.phylip = 1 - 638;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21409.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 591-tomodel-phylip/chr2_21409.phylip = 1 - 638;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21409.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

