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

Input data from file "591-tomodel-phylip/chr19_372.phylip" (7 taxa, 606 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1990.7949	4023.0257	0.2965	0.2965
HKYG	16	-1995.1321	4023.1878	0.2734	0.5700
HKYI	16	-1995.8843	4024.6921	0.1289	0.6988
GTRIG	21	-1990.7992	4025.1805	0.1010	0.7998
GTRI	20	-1991.9028	4025.2415	0.0979	0.8977
HKYIG	17	-1995.1418	4025.3243	0.0940	0.9917
K2PG	13	-2003.0409	4032.6966	0.0024	0.9940
SYMG	17	-1998.9218	4032.8844	0.0021	0.9962
K2PI	13	-2003.5057	4033.6263	0.0015	0.9977
K2PIG	14	-2003.0334	4034.7774	0.0008	0.9985
SYMI	17	-1999.9430	4034.9268	0.0008	0.9993
SYMIG	18	-1998.9341	4035.0334	0.0007	1.0000
GTR	19	-2015.7077	4070.7123	0.0000	1.0000
HKY	15	-2020.6490	4072.1116	0.0000	1.0000
SYM	16	-2023.7919	4080.5074	0.0000	1.0000
K2P	12	-2028.9322	4082.3906	0.0000	1.0000
F81G	15	-2037.0073	4104.8281	0.0000	1.0000
F81IG	16	-2037.0104	4106.9445	0.0000	1.0000
F81I	15	-2038.6679	4108.1493	0.0000	1.0000
JC69G	12	-2045.5806	4115.6873	0.0000	1.0000
JC69I	12	-2046.9992	4118.5244	0.0000	1.0000
JC69IG	13	-2046.7007	4120.0163	0.0000	1.0000
F81	14	-2063.2115	4155.1336	0.0000	1.0000
JC69	11	-2070.1885	4162.8214	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1990.7949	4021.5898	0.3303	0.3303
HKYG	16	-1995.1321	4022.2642	0.2358	0.5660
GTRIG	21	-1990.7992	4023.5983	0.1210	0.6870
HKYI	16	-1995.8843	4023.7685	0.1111	0.7982
GTRI	20	-1991.9028	4023.8056	0.1091	0.9072
HKYIG	17	-1995.1418	4024.2835	0.0859	0.9931
SYMG	17	-1998.9218	4031.8436	0.0020	0.9951
K2PG	13	-2003.0409	4032.0817	0.0017	0.9968
K2PI	13	-2003.5057	4033.0115	0.0011	0.9979
SYMIG	18	-1998.9341	4033.8682	0.0007	0.9986
SYMI	17	-1999.9430	4033.8860	0.0007	0.9994
K2PIG	14	-2003.0334	4034.0667	0.0006	1.0000
GTR	19	-2015.7077	4069.4153	0.0000	1.0000
HKY	15	-2020.6490	4071.2981	0.0000	1.0000
SYM	16	-2023.7919	4079.5838	0.0000	1.0000
K2P	12	-2028.9322	4081.8645	0.0000	1.0000
F81G	15	-2037.0073	4104.0145	0.0000	1.0000
F81IG	16	-2037.0104	4106.0209	0.0000	1.0000
F81I	15	-2038.6679	4107.3357	0.0000	1.0000
JC69G	12	-2045.5806	4115.1611	0.0000	1.0000
JC69I	12	-2046.9992	4117.9983	0.0000	1.0000
JC69IG	13	-2046.7007	4119.4015	0.0000	1.0000
F81	14	-2063.2115	4154.4230	0.0000	1.0000
JC69	11	-2070.1885	4162.3770	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2003.0409	4089.3712	0.5141	0.5141
K2PI	13	-2003.5057	4090.3009	0.3230	0.8371
HKYG	16	-1995.1321	4092.7743	0.0938	0.9308
HKYI	16	-1995.8843	4094.2786	0.0442	0.9750
K2PIG	14	-2003.0334	4095.7630	0.0210	0.9961
HKYIG	17	-1995.1418	4099.2005	0.0038	0.9999
SYMG	17	-1998.9218	4106.7605	0.0001	0.9999
SYMI	17	-1999.9430	4108.8030	0.0000	1.0000
GTRG	20	-1990.7949	4109.7274	0.0000	1.0000
GTRI	20	-1991.9028	4111.9432	0.0000	1.0000
SYMIG	18	-1998.9341	4113.1920	0.0000	1.0000
GTRIG	21	-1990.7992	4116.1428	0.0000	1.0000
K2P	12	-2028.9322	4134.7470	0.0000	1.0000
HKY	15	-2020.6490	4137.4013	0.0000	1.0000
SYM	16	-2023.7919	4150.0939	0.0000	1.0000
GTR	19	-2015.7077	4153.1461	0.0000	1.0000
JC69G	12	-2045.5806	4168.0437	0.0000	1.0000
F81G	15	-2037.0073	4170.1177	0.0000	1.0000
JC69I	12	-2046.9992	4170.8809	0.0000	1.0000
F81I	15	-2038.6679	4173.4389	0.0000	1.0000
F81IG	16	-2037.0104	4176.5310	0.0000	1.0000
JC69IG	13	-2046.7007	4176.6909	0.0000	1.0000
JC69	11	-2070.1885	4210.8526	0.0000	1.0000
F81	14	-2063.2115	4216.1193	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/chr19_372.phylip = 1 - 606;
 Partition Dummy = 1:591-tomodel-phylip/chr19_372.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/chr19_372.phylip = 1 - 606;
 Partition Dummy = 1:591-tomodel-phylip/chr19_372.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr19_372.phylip = 1 - 606;
 Partition Dummy = 1:591-tomodel-phylip/chr19_372.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

