
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:15 04/04/2011

Input data from file "683-tomodel-phylip/chr5_11202.phylip" (7 taxa, 676 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	-1956.9904	3955.2633	0.5215	0.5215
GTRI	20	-1957.6880	3956.6584	0.2596	0.7811
GTRIG	21	-1956.9830	3957.3788	0.1811	0.9622
HKYG	16	-1964.7132	3962.2519	0.0158	0.9780
HKYI	16	-1964.7187	3962.2629	0.0158	0.9937
HKYIG	17	-1964.5891	3964.1082	0.0063	1.0000
GTR	19	-1971.3783	3981.9152	0.0000	1.0000
HKY	15	-1978.6388	3988.0048	0.0000	1.0000
F81I	15	-2017.7194	4066.1660	0.0000	1.0000
F81G	15	-2018.4117	4067.5508	0.0000	1.0000
F81IG	16	-2017.7248	4068.2750	0.0000	1.0000
F81	14	-2030.5283	4089.6920	0.0000	1.0000
K2PI	13	-2040.1899	4106.9296	0.0000	1.0000
K2PG	13	-2041.0350	4108.6199	0.0000	1.0000
K2PIG	14	-2040.1995	4109.0343	0.0000	1.0000
SYMI	17	-2038.7571	4112.4444	0.0000	1.0000
SYMG	17	-2039.5967	4114.1235	0.0000	1.0000
SYMIG	18	-2038.7665	4114.5740	0.0000	1.0000
K2P	12	-2055.5947	4135.6599	0.0000	1.0000
SYM	16	-2054.4562	4141.7380	0.0000	1.0000
JC69I	12	-2080.7845	4186.0396	0.0000	1.0000
JC69G	12	-2081.8283	4188.1271	0.0000	1.0000
JC69IG	13	-2080.8004	4188.1506	0.0000	1.0000
JC69	11	-2094.6747	4211.7470	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1956.9904	3953.9809	0.5190	0.5190
GTRI	20	-1957.6880	3955.3760	0.2584	0.7774
GTRIG	21	-1956.9830	3955.9660	0.1924	0.9698
HKYG	16	-1964.7132	3961.4264	0.0125	0.9823
HKYI	16	-1964.7187	3961.4374	0.0125	0.9948
HKYIG	17	-1964.5891	3963.1781	0.0052	1.0000
GTR	19	-1971.3783	3980.7567	0.0000	1.0000
HKY	15	-1978.6388	3987.2775	0.0000	1.0000
F81I	15	-2017.7194	4065.4387	0.0000	1.0000
F81G	15	-2018.4117	4066.8235	0.0000	1.0000
F81IG	16	-2017.7248	4067.4495	0.0000	1.0000
F81	14	-2030.5283	4089.0566	0.0000	1.0000
K2PI	13	-2040.1899	4106.3798	0.0000	1.0000
K2PG	13	-2041.0350	4108.0701	0.0000	1.0000
K2PIG	14	-2040.1995	4108.3989	0.0000	1.0000
SYMI	17	-2038.7571	4111.5143	0.0000	1.0000
SYMG	17	-2039.5967	4113.1934	0.0000	1.0000
SYMIG	18	-2038.7665	4113.5329	0.0000	1.0000
K2P	12	-2055.5947	4135.1893	0.0000	1.0000
SYM	16	-2054.4562	4140.9125	0.0000	1.0000
JC69I	12	-2080.7845	4185.5690	0.0000	1.0000
JC69IG	13	-2080.8004	4187.6008	0.0000	1.0000
JC69G	12	-2081.8283	4187.6566	0.0000	1.0000
JC69	11	-2094.6747	4211.3494	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1964.7132	4033.6855	0.4888	0.4888
HKYI	16	-1964.7187	4033.6964	0.4862	0.9750
HKYIG	17	-1964.5891	4039.9534	0.0213	0.9963
GTRG	20	-1956.9904	4044.3047	0.0024	0.9987
GTRI	20	-1957.6880	4045.6998	0.0012	0.9999
GTRIG	21	-1956.9830	4050.8061	0.0001	1.0000
HKY	15	-1978.6388	4055.0204	0.0000	1.0000
GTR	19	-1971.3783	4066.5643	0.0000	1.0000
F81I	15	-2017.7194	4133.1816	0.0000	1.0000
F81G	15	-2018.4117	4134.5664	0.0000	1.0000
F81IG	16	-2017.7248	4139.7086	0.0000	1.0000
F81	14	-2030.5283	4152.2833	0.0000	1.0000
K2PI	13	-2040.1899	4165.0903	0.0000	1.0000
K2PG	13	-2041.0350	4166.7806	0.0000	1.0000
K2PIG	14	-2040.1995	4171.6256	0.0000	1.0000
SYMI	17	-2038.7571	4188.2896	0.0000	1.0000
K2P	12	-2055.5947	4189.3837	0.0000	1.0000
SYMG	17	-2039.5967	4189.9687	0.0000	1.0000
SYMIG	18	-2038.7665	4194.8244	0.0000	1.0000
SYM	16	-2054.4562	4213.1715	0.0000	1.0000
JC69I	12	-2080.7845	4239.7633	0.0000	1.0000
JC69G	12	-2081.8283	4241.8509	0.0000	1.0000
JC69IG	13	-2080.8004	4246.3113	0.0000	1.0000
JC69	11	-2094.6747	4261.0275	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 683-tomodel-phylip/chr5_11202.phylip = 1 - 676;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11202.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 683-tomodel-phylip/chr5_11202.phylip = 1 - 676;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11202.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 683-tomodel-phylip/chr5_11202.phylip = 1 - 676;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11202.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

