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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2042.3159	4117.5106	0.3765	0.3765
GTRG	20	-2038.5907	4118.5473	0.2242	0.6006
HKYIG	17	-2042.1962	4119.3827	0.1476	0.7483
GTRI	20	-2039.3630	4120.0919	0.1036	0.8518
GTRIG	21	-2038.5097	4120.5244	0.0834	0.9352
HKYI	16	-2044.0761	4121.0310	0.0648	1.0000
HKY	15	-2086.5107	4203.7956	0.0000	1.0000
GTR	19	-2083.9266	4207.0869	0.0000	1.0000
K2PG	13	-2091.5746	4209.7345	0.0000	1.0000
K2PIG	14	-2091.2787	4211.2337	0.0000	1.0000
SYMG	17	-2089.6353	4214.2608	0.0000	1.0000
K2PI	13	-2094.1023	4214.7898	0.0000	1.0000
SYMIG	18	-2089.3553	4215.8192	0.0000	1.0000
SYMI	17	-2092.0597	4219.1098	0.0000	1.0000
F81G	15	-2096.2036	4223.1815	0.0000	1.0000
F81IG	16	-2096.0403	4224.9594	0.0000	1.0000
F81I	15	-2097.9100	4226.5942	0.0000	1.0000
K2P	12	-2135.9600	4296.4208	0.0000	1.0000
JC69G	12	-2136.5260	4297.5527	0.0000	1.0000
JC69IG	13	-2136.3037	4299.1926	0.0000	1.0000
F81	14	-2135.6843	4300.0449	0.0000	1.0000
JC69I	12	-2138.5666	4301.6341	0.0000	1.0000
SYM	16	-2134.8218	4302.5224	0.0000	1.0000
JC69	11	-2176.9950	4376.4131	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2042.3159	4116.6317	0.3333	0.3333
GTRG	20	-2038.5907	4117.1814	0.2532	0.5865
HKYIG	17	-2042.1962	4118.3924	0.1382	0.7247
GTRI	20	-2039.3630	4118.7260	0.1170	0.8417
GTRIG	21	-2038.5097	4119.0195	0.1010	0.9427
HKYI	16	-2044.0761	4120.1522	0.0573	1.0000
HKY	15	-2086.5107	4203.0214	0.0000	1.0000
GTR	19	-2083.9266	4205.8531	0.0000	1.0000
K2PG	13	-2091.5746	4209.1493	0.0000	1.0000
K2PIG	14	-2091.2787	4210.5574	0.0000	1.0000
SYMG	17	-2089.6353	4213.2705	0.0000	1.0000
K2PI	13	-2094.1023	4214.2046	0.0000	1.0000
SYMIG	18	-2089.3553	4214.7106	0.0000	1.0000
SYMI	17	-2092.0597	4218.1195	0.0000	1.0000
F81G	15	-2096.2036	4222.4073	0.0000	1.0000
F81IG	16	-2096.0403	4224.0806	0.0000	1.0000
F81I	15	-2097.9100	4225.8200	0.0000	1.0000
K2P	12	-2135.9600	4295.9200	0.0000	1.0000
JC69G	12	-2136.5260	4297.0519	0.0000	1.0000
JC69IG	13	-2136.3037	4298.6074	0.0000	1.0000
F81	14	-2135.6843	4299.3686	0.0000	1.0000
JC69I	12	-2138.5666	4301.1333	0.0000	1.0000
SYM	16	-2134.8218	4301.6435	0.0000	1.0000
JC69	11	-2176.9950	4375.9900	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2042.3159	4187.9149	0.8218	0.8218
HKYI	16	-2044.0761	4191.4354	0.1414	0.9631
HKYIG	17	-2042.1962	4194.1308	0.0367	0.9999
GTRG	20	-2038.5907	4206.2854	0.0001	1.0000
GTRI	20	-2039.3630	4207.8300	0.0000	1.0000
GTRIG	21	-2038.5097	4212.5786	0.0000	1.0000
K2PG	13	-2091.5746	4267.0669	0.0000	1.0000
HKY	15	-2086.5107	4269.8493	0.0000	1.0000
K2PI	13	-2094.1023	4272.1222	0.0000	1.0000
K2PIG	14	-2091.2787	4272.9301	0.0000	1.0000
SYMG	17	-2089.6353	4289.0089	0.0000	1.0000
F81G	15	-2096.2036	4289.2353	0.0000	1.0000
GTR	19	-2083.9266	4290.5019	0.0000	1.0000
F81I	15	-2097.9100	4292.6480	0.0000	1.0000
SYMI	17	-2092.0597	4293.8579	0.0000	1.0000
SYMIG	18	-2089.3553	4294.9042	0.0000	1.0000
F81IG	16	-2096.0403	4295.3638	0.0000	1.0000
K2P	12	-2135.9600	4349.3824	0.0000	1.0000
JC69G	12	-2136.5260	4350.5143	0.0000	1.0000
JC69I	12	-2138.5666	4354.5957	0.0000	1.0000
JC69IG	13	-2136.3037	4356.5250	0.0000	1.0000
F81	14	-2135.6843	4361.7413	0.0000	1.0000
SYM	16	-2134.8218	4372.9267	0.0000	1.0000
JC69	11	-2176.9950	4424.9972	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_2737.phylip = 1 - 636;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2737.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_2737.phylip = 1 - 636;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2737.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_2737.phylip = 1 - 636;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2737.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

