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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-1948.4758	3938.2828	0.4999	0.4999
GTRG	20	-1948.7260	3938.7833	0.3892	0.8890
GTRIG	21	-1948.9727	3941.4121	0.1046	0.9936
HKYI	16	-1957.8712	3948.5991	0.0029	0.9965
HKYG	16	-1957.9658	3948.7883	0.0026	0.9991
HKYIG	17	-1957.9658	3950.8970	0.0009	1.0000
K2PG	13	-1971.0646	3968.6997	0.0000	1.0000
K2PI	13	-1971.1228	3968.8161	0.0000	1.0000
SYMI	17	-1967.4725	3969.9104	0.0000	1.0000
SYMG	17	-1967.5240	3970.0133	0.0000	1.0000
K2PIG	14	-1971.0665	3970.7924	0.0000	1.0000
SYMIG	18	-1967.5240	3972.1286	0.0000	1.0000
F81I	15	-1979.9772	3990.7092	0.0000	1.0000
F81G	15	-1980.2370	3991.2287	0.0000	1.0000
F81IG	16	-1980.2371	3993.3308	0.0000	1.0000
JC69I	12	-1991.4147	4007.3177	0.0000	1.0000
JC69G	12	-1991.5293	4007.5469	0.0000	1.0000
JC69IG	13	-1991.5294	4009.6292	0.0000	1.0000
GTR	19	-1986.4982	4012.1990	0.0000	1.0000
HKY	15	-1993.3373	4017.4294	0.0000	1.0000
K2P	12	-2007.1701	4038.8284	0.0000	1.0000
SYM	16	-2003.0543	4038.9654	0.0000	1.0000
F81	14	-2017.1806	4063.0205	0.0000	1.0000
JC69	11	-2026.1671	4074.7467	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1948.4758	3936.9516	0.4969	0.4969
GTRG	20	-1948.7260	3937.4521	0.3869	0.8837
GTRIG	21	-1948.9727	3939.9454	0.1112	0.9949
HKYI	16	-1957.8712	3947.7424	0.0023	0.9972
HKYG	16	-1957.9658	3947.9316	0.0021	0.9992
HKYIG	17	-1957.9658	3949.9317	0.0008	1.0000
K2PG	13	-1971.0646	3968.1292	0.0000	1.0000
K2PI	13	-1971.1228	3968.2456	0.0000	1.0000
SYMI	17	-1967.4725	3968.9451	0.0000	1.0000
SYMG	17	-1967.5240	3969.0480	0.0000	1.0000
K2PIG	14	-1971.0665	3970.1331	0.0000	1.0000
SYMIG	18	-1967.5240	3971.0481	0.0000	1.0000
F81I	15	-1979.9772	3989.9545	0.0000	1.0000
F81G	15	-1980.2370	3990.4740	0.0000	1.0000
F81IG	16	-1980.2371	3992.4741	0.0000	1.0000
JC69I	12	-1991.4147	4006.8294	0.0000	1.0000
JC69G	12	-1991.5293	4007.0586	0.0000	1.0000
JC69IG	13	-1991.5294	4009.0587	0.0000	1.0000
GTR	19	-1986.4982	4010.9965	0.0000	1.0000
HKY	15	-1993.3373	4016.6747	0.0000	1.0000
SYM	16	-2003.0543	4038.1087	0.0000	1.0000
K2P	12	-2007.1701	4038.3402	0.0000	1.0000
F81	14	-2017.1806	4062.3612	0.0000	1.0000
JC69	11	-2026.1671	4074.3342	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1957.8712	4019.4231	0.4859	0.4859
HKYG	16	-1957.9658	4019.6123	0.4421	0.9280
HKYIG	17	-1957.9658	4026.0924	0.0173	0.9453
K2PG	13	-1971.0646	4026.3697	0.0151	0.9604
K2PI	13	-1971.1228	4026.4861	0.0142	0.9746
GTRI	20	-1948.4758	4026.5525	0.0138	0.9884
GTRG	20	-1948.7260	4027.0530	0.0107	0.9991
K2PIG	14	-1971.0665	4032.8537	0.0006	0.9997
GTRIG	21	-1948.9727	4034.0264	0.0003	1.0000
SYMI	17	-1967.4725	4045.1059	0.0000	1.0000
SYMG	17	-1967.5240	4045.2088	0.0000	1.0000
SYMIG	18	-1967.5240	4051.6889	0.0000	1.0000
F81I	15	-1979.9772	4057.1552	0.0000	1.0000
F81G	15	-1980.2370	4057.6747	0.0000	1.0000
JC69I	12	-1991.4147	4060.5899	0.0000	1.0000
JC69G	12	-1991.5293	4060.8192	0.0000	1.0000
F81IG	16	-1980.2371	4064.1549	0.0000	1.0000
JC69IG	13	-1991.5294	4067.2993	0.0000	1.0000
HKY	15	-1993.3373	4083.8753	0.0000	1.0000
K2P	12	-2007.1701	4092.1007	0.0000	1.0000
GTR	19	-1986.4982	4096.1173	0.0000	1.0000
SYM	16	-2003.0543	4109.7894	0.0000	1.0000
JC69	11	-2026.1671	4123.6147	0.0000	1.0000
F81	14	-2017.1806	4125.0818	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr19_346.phylip = 1 - 652;
 Partition Dummy = 1:591-tomodel-phylip/chr19_346.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr19_346.phylip = 1 - 652;
 Partition Dummy = 1:591-tomodel-phylip/chr19_346.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr19_346.phylip = 1 - 652;
 Partition Dummy = 1:591-tomodel-phylip/chr19_346.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

