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

Input data from file "591-tomodel-phylip/chr2_24799.phylip" (7 taxa, 551 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	-1983.9708	4009.5265	0.4724	0.4724
GTRI	20	-1984.2096	4010.0041	0.3720	0.8444
GTRIG	21	-1984.0028	4011.7523	0.1552	0.9996
HKYG	16	-1995.9602	4024.9390	0.0002	0.9998
HKYI	16	-1996.5893	4026.1973	0.0001	0.9999
HKYIG	17	-1995.9600	4027.0683	0.0001	1.0000
SYMG	17	-2000.7467	4036.6415	0.0000	1.0000
SYMI	17	-2001.6908	4038.5298	0.0000	1.0000
SYMIG	18	-2000.7543	4038.7943	0.0000	1.0000
K2PG	13	-2009.2021	4045.0819	0.0000	1.0000
K2PI	13	-2010.2393	4047.1564	0.0000	1.0000
K2PIG	14	-2009.2043	4047.1921	0.0000	1.0000
GTR	19	-2009.7532	4058.9378	0.0000	1.0000
HKY	15	-2019.9869	4070.8709	0.0000	1.0000
SYM	16	-2022.5363	4078.0913	0.0000	1.0000
K2P	12	-2030.9351	4086.4501	0.0000	1.0000
F81I	15	-2072.1935	4175.2843	0.0000	1.0000
F81G	15	-2072.5237	4175.9446	0.0000	1.0000
F81IG	16	-2072.7508	4178.5202	0.0000	1.0000
JC69I	12	-2080.6517	4185.8833	0.0000	1.0000
JC69G	12	-2080.6936	4185.9671	0.0000	1.0000
JC69IG	13	-2080.6091	4187.8960	0.0000	1.0000
F81	14	-2090.7030	4210.1896	0.0000	1.0000
JC69	11	-2098.2487	4218.9872	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1983.9708	4007.9416	0.4663	0.4663
GTRI	20	-1984.2096	4008.4192	0.3672	0.8336
GTRIG	21	-1984.0028	4010.0056	0.1661	0.9997
HKYG	16	-1995.9602	4023.9203	0.0002	0.9999
HKYI	16	-1996.5893	4025.1786	0.0001	0.9999
HKYIG	17	-1995.9600	4025.9201	0.0001	1.0000
SYMG	17	-2000.7467	4035.4933	0.0000	1.0000
SYMI	17	-2001.6908	4037.3815	0.0000	1.0000
SYMIG	18	-2000.7543	4037.5086	0.0000	1.0000
K2PG	13	-2009.2021	4044.4041	0.0000	1.0000
K2PIG	14	-2009.2043	4046.4085	0.0000	1.0000
K2PI	13	-2010.2393	4046.4785	0.0000	1.0000
GTR	19	-2009.7532	4057.5065	0.0000	1.0000
HKY	15	-2019.9869	4069.9737	0.0000	1.0000
SYM	16	-2022.5363	4077.0726	0.0000	1.0000
K2P	12	-2030.9351	4085.8702	0.0000	1.0000
F81I	15	-2072.1935	4174.3871	0.0000	1.0000
F81G	15	-2072.5237	4175.0474	0.0000	1.0000
F81IG	16	-2072.7508	4177.5015	0.0000	1.0000
JC69I	12	-2080.6517	4185.3033	0.0000	1.0000
JC69G	12	-2080.6936	4185.3872	0.0000	1.0000
JC69IG	13	-2080.6091	4187.2182	0.0000	1.0000
F81	14	-2090.7030	4209.4061	0.0000	1.0000
JC69	11	-2098.2487	4218.4974	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1995.9602	4092.9081	0.3879	0.3879
HKYI	16	-1996.5893	4094.1664	0.2067	0.5946
GTRG	20	-1983.9708	4094.1763	0.2057	0.8003
GTRI	20	-1984.2096	4094.6539	0.1620	0.9624
HKYIG	17	-1995.9600	4099.2196	0.0165	0.9789
K2PG	13	-2009.2021	4100.4567	0.0089	0.9878
GTRIG	21	-1984.0028	4100.5520	0.0085	0.9963
K2PI	13	-2010.2393	4102.5311	0.0032	0.9994
K2PIG	14	-2009.2043	4106.7728	0.0004	0.9998
SYMG	17	-2000.7467	4108.7928	0.0001	0.9999
SYMI	17	-2001.6908	4110.6810	0.0001	1.0000
SYMIG	18	-2000.7543	4115.1198	0.0000	1.0000
HKY	15	-2019.9869	4134.6497	0.0000	1.0000
K2P	12	-2030.9351	4137.6110	0.0000	1.0000
GTR	19	-2009.7532	4139.4295	0.0000	1.0000
SYM	16	-2022.5363	4146.0603	0.0000	1.0000
JC69I	12	-2080.6517	4237.0442	0.0000	1.0000
JC69G	12	-2080.6936	4237.1280	0.0000	1.0000
F81I	15	-2072.1935	4239.0631	0.0000	1.0000
F81G	15	-2072.5237	4239.7235	0.0000	1.0000
JC69IG	13	-2080.6091	4243.2708	0.0000	1.0000
F81IG	16	-2072.7508	4246.4893	0.0000	1.0000
JC69	11	-2098.2487	4265.9265	0.0000	1.0000
F81	14	-2090.7030	4269.7703	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_24799.phylip = 1 - 551;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24799.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_24799.phylip = 1 - 551;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24799.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_24799.phylip = 1 - 551;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24799.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

