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

Input data from file "683-tomodel-phylip/chr1_22334.phylip" (7 taxa, 582 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	-1988.9629	4019.4232	0.4938	0.4938
GTRG	20	-1989.3705	4020.2383	0.3285	0.8223
GTRIG	21	-1988.9666	4021.5832	0.1677	0.9900
HKYI	16	-1997.7672	4028.4972	0.0053	0.9953
HKYG	16	-1998.3955	4029.7538	0.0028	0.9981
HKYIG	17	-1997.7778	4030.6407	0.0018	0.9999
GTR	19	-1998.6668	4036.6858	0.0001	1.0000
HKY	15	-2008.1986	4047.2453	0.0000	1.0000
SYMI	17	-2006.8016	4048.6883	0.0000	1.0000
SYMG	17	-2007.0615	4049.2081	0.0000	1.0000
SYMIG	18	-2006.8069	4050.8287	0.0000	1.0000
K2PI	13	-2012.6031	4051.8471	0.0000	1.0000
K2PG	13	-2013.2251	4053.0911	0.0000	1.0000
K2PIG	14	-2012.6146	4053.9700	0.0000	1.0000
SYM	16	-2014.6191	4062.2009	0.0000	1.0000
K2P	12	-2022.0264	4068.6011	0.0000	1.0000
F81I	15	-2033.7640	4098.3760	0.0000	1.0000
F81G	15	-2034.5293	4099.9066	0.0000	1.0000
F81IG	16	-2033.7802	4100.5233	0.0000	1.0000
JC69I	12	-2044.7563	4114.0609	0.0000	1.0000
F81	14	-2042.6849	4114.1105	0.0000	1.0000
JC69G	12	-2045.5304	4115.6091	0.0000	1.0000
JC69IG	13	-2044.7732	4116.1873	0.0000	1.0000
JC69	11	-2053.4987	4129.4606	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1988.9629	4017.9259	0.4884	0.4884
GTRG	20	-1989.3705	4018.7410	0.3249	0.8133
GTRIG	21	-1988.9666	4019.9332	0.1790	0.9923
HKYI	16	-1997.7672	4027.5343	0.0040	0.9963
HKYG	16	-1998.3955	4028.7909	0.0021	0.9985
HKYIG	17	-1997.7778	4029.5556	0.0015	0.9999
GTR	19	-1998.6668	4035.3335	0.0001	1.0000
HKY	15	-2008.1986	4046.3972	0.0000	1.0000
SYMI	17	-2006.8016	4047.6032	0.0000	1.0000
SYMG	17	-2007.0615	4048.1230	0.0000	1.0000
SYMIG	18	-2006.8069	4049.6138	0.0000	1.0000
K2PI	13	-2012.6031	4051.2063	0.0000	1.0000
K2PG	13	-2013.2251	4052.4503	0.0000	1.0000
K2PIG	14	-2012.6146	4053.2293	0.0000	1.0000
SYM	16	-2014.6191	4061.2381	0.0000	1.0000
K2P	12	-2022.0264	4068.0528	0.0000	1.0000
F81I	15	-2033.7640	4097.5279	0.0000	1.0000
F81G	15	-2034.5293	4099.0586	0.0000	1.0000
F81IG	16	-2033.7802	4099.5604	0.0000	1.0000
F81	14	-2042.6849	4113.3698	0.0000	1.0000
JC69I	12	-2044.7563	4113.5125	0.0000	1.0000
JC69G	12	-2045.5304	4115.0607	0.0000	1.0000
JC69IG	13	-2044.7732	4115.5464	0.0000	1.0000
JC69	11	-2053.4987	4128.9974	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1997.7672	4097.3978	0.6185	0.6185
HKYG	16	-1998.3955	4098.6544	0.3300	0.9485
HKYIG	17	-1997.7778	4103.7856	0.0254	0.9738
GTRI	20	-1988.9629	4105.2553	0.0122	0.9860
GTRG	20	-1989.3705	4106.0704	0.0081	0.9941
K2PI	13	-2012.6031	4107.9704	0.0031	0.9972
K2PG	13	-2013.2251	4109.2144	0.0017	0.9989
GTRIG	21	-1988.9666	4111.6291	0.0005	0.9994
HKY	15	-2008.1986	4111.8943	0.0004	0.9998
K2PIG	14	-2012.6146	4114.3599	0.0001	1.0000
GTR	19	-1998.6668	4118.2964	0.0000	1.0000
K2P	12	-2022.0264	4120.4504	0.0000	1.0000
SYMI	17	-2006.8016	4121.8332	0.0000	1.0000
SYMG	17	-2007.0615	4122.3530	0.0000	1.0000
SYMIG	18	-2006.8069	4128.2103	0.0000	1.0000
SYM	16	-2014.6191	4131.1016	0.0000	1.0000
F81I	15	-2033.7640	4163.0250	0.0000	1.0000
F81G	15	-2034.5293	4164.5556	0.0000	1.0000
JC69I	12	-2044.7563	4165.9102	0.0000	1.0000
JC69G	12	-2045.5304	4167.4584	0.0000	1.0000
F81IG	16	-2033.7802	4169.4240	0.0000	1.0000
JC69IG	13	-2044.7732	4172.3105	0.0000	1.0000
F81	14	-2042.6849	4174.5004	0.0000	1.0000
JC69	11	-2053.4987	4177.0286	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 683-tomodel-phylip/chr1_22334.phylip = 1 - 582;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22334.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 683-tomodel-phylip/chr1_22334.phylip = 1 - 582;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22334.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 683-tomodel-phylip/chr1_22334.phylip = 1 - 582;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22334.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

