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

Input data from file "683-tomodel-phylip/chr2_23113.phylip" (7 taxa, 678 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	-2088.4773	4218.2331	0.5141	0.5141
HKYG	16	-2093.1419	4219.1067	0.3321	0.8462
GTRIG	21	-2089.6876	4222.7838	0.0528	0.8990
GTRI	20	-2090.7692	4222.8169	0.0520	0.9510
HKYI	16	-2095.3187	4223.4605	0.0377	0.9886
HKYIG	17	-2095.4639	4225.8551	0.0114	1.0000
HKY	15	-2110.7914	4252.3079	0.0000	1.0000
GTR	19	-2106.8650	4252.8851	0.0000	1.0000
F81G	15	-2143.5234	4317.7720	0.0000	1.0000
F81IG	16	-2143.5145	4319.8520	0.0000	1.0000
K2PG	13	-2146.9337	4320.4156	0.0000	1.0000
F81I	15	-2145.3355	4321.3962	0.0000	1.0000
SYMG	17	-2145.9832	4326.8937	0.0000	1.0000
K2PIG	14	-2149.1443	4326.9221	0.0000	1.0000
K2PI	13	-2151.4502	4329.4486	0.0000	1.0000
SYMIG	18	-2148.2150	4333.4679	0.0000	1.0000
SYMI	17	-2150.3860	4335.6993	0.0000	1.0000
F81	14	-2160.0219	4348.6773	0.0000	1.0000
K2P	12	-2164.9467	4354.3627	0.0000	1.0000
SYM	16	-2164.1255	4361.0740	0.0000	1.0000
JC69G	12	-2186.7202	4397.9097	0.0000	1.0000
JC69IG	13	-2186.7072	4399.9625	0.0000	1.0000
JC69I	12	-2188.6330	4401.7353	0.0000	1.0000
JC69	11	-2203.5660	4429.5284	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2088.4773	4216.9545	0.5549	0.5549
HKYG	16	-2093.1419	4218.2837	0.2855	0.8404
GTRIG	21	-2089.6876	4221.3752	0.0609	0.9012
GTRI	20	-2090.7692	4221.5384	0.0561	0.9573
HKYI	16	-2095.3187	4222.6375	0.0324	0.9897
HKYIG	17	-2095.4639	4224.9279	0.0103	1.0000
HKY	15	-2110.7914	4251.5828	0.0000	1.0000
GTR	19	-2106.8650	4251.7300	0.0000	1.0000
F81G	15	-2143.5234	4317.0469	0.0000	1.0000
F81IG	16	-2143.5145	4319.0290	0.0000	1.0000
K2PG	13	-2146.9337	4319.8674	0.0000	1.0000
F81I	15	-2145.3355	4320.6711	0.0000	1.0000
SYMG	17	-2145.9832	4325.9664	0.0000	1.0000
K2PIG	14	-2149.1443	4326.2886	0.0000	1.0000
K2PI	13	-2151.4502	4328.9004	0.0000	1.0000
SYMIG	18	-2148.2150	4332.4300	0.0000	1.0000
SYMI	17	-2150.3860	4334.7720	0.0000	1.0000
F81	14	-2160.0219	4348.0438	0.0000	1.0000
K2P	12	-2164.9467	4353.8935	0.0000	1.0000
SYM	16	-2164.1255	4360.2510	0.0000	1.0000
JC69G	12	-2186.7202	4397.4405	0.0000	1.0000
JC69IG	13	-2186.7072	4399.4143	0.0000	1.0000
JC69I	12	-2188.6330	4401.2661	0.0000	1.0000
JC69	11	-2203.5660	4429.1320	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2093.1419	4290.5901	0.8949	0.8949
HKYI	16	-2095.3187	4294.9438	0.1015	0.9964
HKYIG	17	-2095.4639	4301.7534	0.0034	0.9998
GTRG	20	-2088.4773	4307.3375	0.0002	1.0000
GTRI	20	-2090.7692	4311.9213	0.0000	1.0000
GTRIG	21	-2089.6876	4316.2773	0.0000	1.0000
HKY	15	-2110.7914	4319.3700	0.0000	1.0000
GTR	19	-2106.8650	4337.5938	0.0000	1.0000
K2PG	13	-2146.9337	4378.6163	0.0000	1.0000
F81G	15	-2143.5234	4384.8341	0.0000	1.0000
K2PI	13	-2151.4502	4387.6493	0.0000	1.0000
F81I	15	-2145.3355	4388.4583	0.0000	1.0000
K2PIG	14	-2149.1443	4389.5566	0.0000	1.0000
F81IG	16	-2143.5145	4391.3354	0.0000	1.0000
SYMG	17	-2145.9832	4402.7919	0.0000	1.0000
K2P	12	-2164.9467	4408.1232	0.0000	1.0000
F81	14	-2160.0219	4411.3119	0.0000	1.0000
SYMI	17	-2150.3860	4411.5975	0.0000	1.0000
SYMIG	18	-2148.2150	4413.7746	0.0000	1.0000
SYM	16	-2164.1255	4432.5574	0.0000	1.0000
JC69G	12	-2186.7202	4451.6703	0.0000	1.0000
JC69I	12	-2188.6330	4455.4958	0.0000	1.0000
JC69IG	13	-2186.7072	4458.1632	0.0000	1.0000
JC69	11	-2203.5660	4478.8426	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 683-tomodel-phylip/chr2_23113.phylip = 1 - 678;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23113.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 683-tomodel-phylip/chr2_23113.phylip = 1 - 678;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23113.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 683-tomodel-phylip/chr2_23113.phylip = 1 - 678;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23113.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

