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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-2054.7740	4142.2628	0.3548	0.3548
HKYG	16	-2055.3222	4143.3593	0.2051	0.5599
GTRI	20	-2051.2446	4143.5988	0.1819	0.7418
GTRG	20	-2051.5202	4144.1500	0.1381	0.8799
HKYIG	17	-2055.3223	4145.4499	0.0721	0.9520
GTRIG	21	-2051.5202	4146.2627	0.0480	1.0000
K2PI	13	-2070.3868	4167.2499	0.0000	1.0000
K2PG	13	-2070.8873	4168.2509	0.0000	1.0000
SYMI	17	-2067.7562	4170.3177	0.0000	1.0000
K2PIG	14	-2070.8874	4170.3252	0.0000	1.0000
SYMG	17	-2068.1341	4171.0735	0.0000	1.0000
SYMIG	18	-2068.1342	4173.1695	0.0000	1.0000
HKY	15	-2081.1624	4192.9547	0.0000	1.0000
GTR	19	-2078.0328	4195.0683	0.0000	1.0000
F81I	15	-2086.3975	4203.4249	0.0000	1.0000
F81G	15	-2086.8763	4204.3825	0.0000	1.0000
F81IG	16	-2086.8764	4206.4675	0.0000	1.0000
K2P	12	-2097.0201	4218.4480	0.0000	1.0000
SYM	16	-2094.5760	4221.8668	0.0000	1.0000
JC69I	12	-2099.7517	4223.9113	0.0000	1.0000
JC69G	12	-2100.2343	4224.8764	0.0000	1.0000
JC69IG	13	-2100.2344	4226.9452	0.0000	1.0000
F81	14	-2111.4570	4251.4645	0.0000	1.0000
JC69	11	-2125.1824	4272.7095	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2054.7740	4141.5480	0.3264	0.3264
GTRI	20	-2051.2446	4142.4891	0.2039	0.5303
HKYG	16	-2055.3222	4142.6444	0.1886	0.7189
GTRG	20	-2051.5202	4143.0404	0.1548	0.8737
HKYIG	17	-2055.3223	4144.6446	0.0694	0.9431
GTRIG	21	-2051.5202	4145.0405	0.0569	1.0000
K2PI	13	-2070.3868	4166.7735	0.0000	1.0000
K2PG	13	-2070.8873	4167.7745	0.0000	1.0000
SYMI	17	-2067.7562	4169.5124	0.0000	1.0000
K2PIG	14	-2070.8874	4169.7747	0.0000	1.0000
SYMG	17	-2068.1341	4170.2682	0.0000	1.0000
SYMIG	18	-2068.1342	4172.2683	0.0000	1.0000
HKY	15	-2081.1624	4192.3248	0.0000	1.0000
GTR	19	-2078.0328	4194.0656	0.0000	1.0000
F81I	15	-2086.3975	4202.7949	0.0000	1.0000
F81G	15	-2086.8763	4203.7526	0.0000	1.0000
F81IG	16	-2086.8764	4205.7527	0.0000	1.0000
K2P	12	-2097.0201	4218.0402	0.0000	1.0000
SYM	16	-2094.5760	4221.1520	0.0000	1.0000
JC69I	12	-2099.7517	4223.5035	0.0000	1.0000
JC69G	12	-2100.2343	4224.4686	0.0000	1.0000
JC69IG	13	-2100.2344	4226.4688	0.0000	1.0000
F81	14	-2111.4570	4250.9141	0.0000	1.0000
JC69	11	-2125.1824	4272.3648	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2054.7740	4216.0556	0.6232	0.6232
HKYG	16	-2055.3222	4217.1520	0.3602	0.9834
HKYIG	17	-2055.3223	4223.8090	0.0129	0.9963
K2PI	13	-2070.3868	4227.3109	0.0022	0.9985
K2PG	13	-2070.8873	4228.3119	0.0014	0.9999
K2PIG	14	-2070.8874	4234.9689	0.0000	0.9999
GTRI	20	-2051.2446	4235.6237	0.0000	1.0000
GTRG	20	-2051.5202	4236.1749	0.0000	1.0000
GTRIG	21	-2051.5202	4242.8317	0.0000	1.0000
SYMI	17	-2067.7562	4248.6768	0.0000	1.0000
SYMG	17	-2068.1341	4249.4326	0.0000	1.0000
SYMIG	18	-2068.1342	4256.0894	0.0000	1.0000
HKY	15	-2081.1624	4262.1757	0.0000	1.0000
F81I	15	-2086.3975	4272.6458	0.0000	1.0000
F81G	15	-2086.8763	4273.6035	0.0000	1.0000
K2P	12	-2097.0201	4273.9209	0.0000	1.0000
JC69I	12	-2099.7517	4279.3842	0.0000	1.0000
F81IG	16	-2086.8764	4280.2603	0.0000	1.0000
JC69G	12	-2100.2343	4280.3493	0.0000	1.0000
GTR	19	-2078.0328	4282.5434	0.0000	1.0000
JC69IG	13	-2100.2344	4287.0062	0.0000	1.0000
SYM	16	-2094.5760	4295.6596	0.0000	1.0000
F81	14	-2111.4570	4316.1082	0.0000	1.0000
JC69	11	-2125.1824	4323.5888	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_5478.phylip = 1 - 778;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5478.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_5478.phylip = 1 - 778;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5478.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_5478.phylip = 1 - 778;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5478.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

