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

Input data from file "591-tomodel-phylip/chr5_10260.phylip" (7 taxa, 704 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	-2184.7417	4402.2752	0.5330	0.5330
HKYIG	17	-2184.5473	4403.9867	0.2265	0.7595
HKYI	16	-2186.0060	4404.8039	0.1505	0.9100
GTRG	20	-2182.8322	4406.8942	0.0529	0.9630
GTRIG	21	-2182.6884	4408.7316	0.0211	0.9841
GTRI	20	-2184.0337	4409.2973	0.0159	1.0000
K2PG	13	-2206.2212	4438.9700	0.0000	1.0000
K2PIG	14	-2206.0349	4440.6794	0.0000	1.0000
K2PI	13	-2207.5224	4441.5723	0.0000	1.0000
SYMG	17	-2204.1854	4443.2628	0.0000	1.0000
SYMIG	18	-2204.0005	4444.9996	0.0000	1.0000
SYMI	17	-2205.5352	4445.9625	0.0000	1.0000
HKY	15	-2224.6818	4480.0612	0.0000	1.0000
GTR	19	-2223.0894	4485.2899	0.0000	1.0000
F81G	15	-2231.7375	4494.1726	0.0000	1.0000
F81IG	16	-2231.5337	4495.8593	0.0000	1.0000
F81I	15	-2233.2544	4497.2065	0.0000	1.0000
K2P	12	-2244.3131	4513.0777	0.0000	1.0000
SYM	16	-2241.9199	4516.6316	0.0000	1.0000
JC69G	12	-2249.6881	4523.8277	0.0000	1.0000
JC69IG	13	-2249.5200	4525.5675	0.0000	1.0000
JC69I	12	-2251.1437	4526.7389	0.0000	1.0000
F81	14	-2268.8366	4566.2828	0.0000	1.0000
JC69	11	-2285.7193	4593.8201	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2184.7417	4401.4834	0.5148	0.5148
HKYIG	17	-2184.5473	4403.0946	0.2300	0.7448
HKYI	16	-2186.0060	4404.0120	0.1454	0.8902
GTRG	20	-2182.8322	4405.6644	0.0636	0.9538
GTRIG	21	-2182.6884	4407.3768	0.0270	0.9809
GTRI	20	-2184.0337	4408.0674	0.0191	1.0000
K2PG	13	-2206.2212	4438.4425	0.0000	1.0000
K2PIG	14	-2206.0349	4440.0698	0.0000	1.0000
K2PI	13	-2207.5224	4441.0448	0.0000	1.0000
SYMG	17	-2204.1854	4442.3707	0.0000	1.0000
SYMIG	18	-2204.0005	4444.0010	0.0000	1.0000
SYMI	17	-2205.5352	4445.0704	0.0000	1.0000
HKY	15	-2224.6818	4479.3635	0.0000	1.0000
GTR	19	-2223.0894	4484.1788	0.0000	1.0000
F81G	15	-2231.7375	4493.4750	0.0000	1.0000
F81IG	16	-2231.5337	4495.0674	0.0000	1.0000
F81I	15	-2233.2544	4496.5088	0.0000	1.0000
K2P	12	-2244.3131	4512.6261	0.0000	1.0000
SYM	16	-2241.9199	4515.8397	0.0000	1.0000
JC69G	12	-2249.6881	4523.3762	0.0000	1.0000
JC69IG	13	-2249.5200	4525.0400	0.0000	1.0000
JC69I	12	-2251.1437	4526.2874	0.0000	1.0000
F81	14	-2268.8366	4565.6733	0.0000	1.0000
JC69	11	-2285.7193	4593.4386	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2184.7417	4474.3918	0.7529	0.7529
HKYI	16	-2186.0060	4476.9205	0.2126	0.9655
HKYIG	17	-2184.5473	4480.5598	0.0345	1.0000
GTRG	20	-2182.8322	4496.7999	0.0000	1.0000
K2PG	13	-2206.2212	4497.6806	0.0000	1.0000
GTRI	20	-2184.0337	4499.2030	0.0000	1.0000
K2PI	13	-2207.5224	4500.2829	0.0000	1.0000
GTRIG	21	-2182.6884	4503.0691	0.0000	1.0000
K2PIG	14	-2206.0349	4503.8647	0.0000	1.0000
SYMG	17	-2204.1854	4519.8360	0.0000	1.0000
SYMI	17	-2205.5352	4522.5356	0.0000	1.0000
SYMIG	18	-2204.0005	4526.0231	0.0000	1.0000
HKY	15	-2224.6818	4547.7152	0.0000	1.0000
F81G	15	-2231.7375	4561.8266	0.0000	1.0000
F81I	15	-2233.2544	4564.8605	0.0000	1.0000
K2P	12	-2244.3131	4567.3075	0.0000	1.0000
F81IG	16	-2231.5337	4567.9759	0.0000	1.0000
GTR	19	-2223.0894	4570.7576	0.0000	1.0000
JC69G	12	-2249.6881	4578.0575	0.0000	1.0000
JC69I	12	-2251.1437	4580.9687	0.0000	1.0000
JC69IG	13	-2249.5200	4584.2781	0.0000	1.0000
SYM	16	-2241.9199	4588.7482	0.0000	1.0000
F81	14	-2268.8366	4629.4682	0.0000	1.0000
JC69	11	-2285.7193	4643.5631	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/chr5_10260.phylip = 1 - 704;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10260.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/chr5_10260.phylip = 1 - 704;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10260.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/chr5_10260.phylip = 1 - 704;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10260.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

