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

Input data from file "683-tomodel-phylip/chr10_5022.phylip" (7 taxa, 699 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-2093.0339	4227.3067	0.4778	0.4778
GTRIG	21	-2092.6463	4228.6575	0.2432	0.7210
GTRG	20	-2094.4412	4230.1212	0.1170	0.8380
HKYIG	17	-2098.2541	4231.4068	0.0615	0.8995
HKYG	16	-2099.4844	4231.7665	0.0514	0.9508
HKYI	16	-2099.5286	4231.8549	0.0492	1.0000
SYMI	17	-2114.1296	4263.1579	0.0000	1.0000
SYMG	17	-2114.1620	4263.2226	0.0000	1.0000
SYMIG	18	-2114.1619	4265.3298	0.0000	1.0000
K2PI	13	-2120.1603	4266.8519	0.0000	1.0000
K2PG	13	-2120.7101	4267.9515	0.0000	1.0000
K2PIG	14	-2119.7756	4268.1653	0.0000	1.0000
GTR	19	-2123.1663	4285.4518	0.0000	1.0000
HKY	15	-2128.8632	4288.4293	0.0000	1.0000
SYM	16	-2141.4738	4315.7453	0.0000	1.0000
K2P	12	-2147.8321	4320.1191	0.0000	1.0000
F81I	15	-2147.4856	4325.6739	0.0000	1.0000
F81G	15	-2147.6931	4326.0890	0.0000	1.0000
F81IG	16	-2147.5908	4327.9793	0.0000	1.0000
JC69I	12	-2163.4534	4351.3616	0.0000	1.0000
JC69G	12	-2163.7535	4351.9619	0.0000	1.0000
JC69IG	13	-2163.9417	4354.4147	0.0000	1.0000
F81	14	-2172.8546	4374.3233	0.0000	1.0000
JC69	11	-2188.2681	4398.9205	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2093.0339	4226.0678	0.4845	0.4845
GTRIG	21	-2092.6463	4227.2926	0.2626	0.7470
GTRG	20	-2094.4412	4228.8823	0.1186	0.8656
HKYIG	17	-2098.2541	4230.5082	0.0526	0.9182
HKYG	16	-2099.4844	4230.9689	0.0418	0.9600
HKYI	16	-2099.5286	4231.0572	0.0400	1.0000
SYMI	17	-2114.1296	4262.2592	0.0000	1.0000
SYMG	17	-2114.1620	4262.3240	0.0000	1.0000
SYMIG	18	-2114.1619	4264.3239	0.0000	1.0000
K2PI	13	-2120.1603	4266.3205	0.0000	1.0000
K2PG	13	-2120.7101	4267.4201	0.0000	1.0000
K2PIG	14	-2119.7756	4267.5513	0.0000	1.0000
GTR	19	-2123.1663	4284.3325	0.0000	1.0000
HKY	15	-2128.8632	4287.7265	0.0000	1.0000
SYM	16	-2141.4738	4314.9477	0.0000	1.0000
K2P	12	-2147.8321	4319.6643	0.0000	1.0000
F81I	15	-2147.4856	4324.9712	0.0000	1.0000
F81G	15	-2147.6931	4325.3863	0.0000	1.0000
F81IG	16	-2147.5908	4327.1817	0.0000	1.0000
JC69I	12	-2163.4534	4350.9068	0.0000	1.0000
JC69G	12	-2163.7535	4351.5071	0.0000	1.0000
JC69IG	13	-2163.9417	4353.8833	0.0000	1.0000
F81	14	-2172.8546	4373.7093	0.0000	1.0000
JC69	11	-2188.2681	4398.5362	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2099.4844	4303.7633	0.4789	0.4789
HKYI	16	-2099.5286	4303.8516	0.4582	0.9372
HKYIG	17	-2098.2541	4307.8522	0.0620	0.9992
GTRI	20	-2093.0339	4317.0608	0.0006	0.9998
GTRG	20	-2094.4412	4319.8753	0.0002	1.0000
GTRIG	21	-2092.6463	4322.8353	0.0000	1.0000
K2PI	13	-2120.1603	4325.4660	0.0000	1.0000
K2PG	13	-2120.7101	4326.5656	0.0000	1.0000
K2PIG	14	-2119.7756	4331.2464	0.0000	1.0000
SYMI	17	-2114.1296	4339.6033	0.0000	1.0000
SYMG	17	-2114.1620	4339.6680	0.0000	1.0000
SYMIG	18	-2114.1619	4346.2176	0.0000	1.0000
HKY	15	-2128.8632	4355.9712	0.0000	1.0000
GTR	19	-2123.1663	4370.7759	0.0000	1.0000
K2P	12	-2147.8321	4374.2601	0.0000	1.0000
SYM	16	-2141.4738	4387.7421	0.0000	1.0000
F81I	15	-2147.4856	4393.2159	0.0000	1.0000
F81G	15	-2147.6931	4393.6310	0.0000	1.0000
F81IG	16	-2147.5908	4399.9761	0.0000	1.0000
JC69I	12	-2163.4534	4405.5026	0.0000	1.0000
JC69G	12	-2163.7535	4406.1029	0.0000	1.0000
JC69IG	13	-2163.9417	4413.0288	0.0000	1.0000
F81	14	-2172.8546	4437.4044	0.0000	1.0000
JC69	11	-2188.2681	4448.5824	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/chr10_5022.phylip = 1 - 699;
 Partition Dummy = 1:683-tomodel-phylip/chr10_5022.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/chr10_5022.phylip = 1 - 699;
 Partition Dummy = 1:683-tomodel-phylip/chr10_5022.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 (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr10_5022.phylip = 1 - 699;
 Partition Dummy = 1:683-tomodel-phylip/chr10_5022.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

