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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2348.1446	4728.2892	0.4229	0.4229
HKYG	16	-2348.3014	4728.6028	0.3615	0.7844
HKYIG	17	-2348.3090	4730.6180	0.1320	0.9164
GTRI	20	-2346.5581	4733.1162	0.0379	0.9543
GTRG	20	-2346.7323	4733.4646	0.0318	0.9861
GTRIG	21	-2346.5604	4735.1208	0.0139	1.0000
HKY	15	-2359.1583	4748.3166	0.0000	1.0000
GTR	19	-2357.1680	4752.3359	0.0000	1.0000
SYMG	17	-2370.3663	4774.7326	0.0000	1.0000
SYMI	17	-2370.4751	4774.9501	0.0000	1.0000
SYMIG	18	-2370.3789	4776.7577	0.0000	1.0000
K2PG	13	-2375.5132	4777.0265	0.0000	1.0000
K2PI	13	-2375.5717	4777.1433	0.0000	1.0000
K2PIG	14	-2375.5259	4779.0518	0.0000	1.0000
SYM	16	-2380.2496	4792.4992	0.0000	1.0000
F81G	15	-2381.5490	4793.0979	0.0000	1.0000
F81I	15	-2381.7384	4793.4767	0.0000	1.0000
K2P	12	-2385.3557	4794.7114	0.0000	1.0000
F81IG	16	-2381.5575	4795.1150	0.0000	1.0000
F81	14	-2391.9169	4811.8337	0.0000	1.0000
JC69G	12	-2407.0026	4838.0052	0.0000	1.0000
JC69I	12	-2407.2866	4838.5732	0.0000	1.0000
JC69IG	13	-2407.0127	4840.0254	0.0000	1.0000
JC69	11	-2416.4144	4854.8288	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2348.1446	4804.7497	0.5308	0.5308
HKYG	16	-2348.3014	4805.0634	0.4537	0.9845
HKYIG	17	-2348.3090	4811.8573	0.0152	0.9997
HKY	15	-2359.1583	4819.9984	0.0003	1.0000
GTRI	20	-2346.5581	4828.6919	0.0000	1.0000
GTRG	20	-2346.7323	4829.0403	0.0000	1.0000
GTRIG	21	-2346.5604	4835.4753	0.0000	1.0000
K2PG	13	-2375.5132	4839.1507	0.0000	1.0000
K2PI	13	-2375.5717	4839.2675	0.0000	1.0000
GTR	19	-2357.1680	4843.1329	0.0000	1.0000
K2PIG	14	-2375.5259	4845.9548	0.0000	1.0000
K2P	12	-2385.3557	4852.0569	0.0000	1.0000
SYMG	17	-2370.3663	4855.9719	0.0000	1.0000
SYMI	17	-2370.4751	4856.1895	0.0000	1.0000
SYMIG	18	-2370.3789	4862.7759	0.0000	1.0000
F81G	15	-2381.5490	4864.7797	0.0000	1.0000
F81I	15	-2381.7384	4865.1585	0.0000	1.0000
SYM	16	-2380.2496	4868.9597	0.0000	1.0000
F81IG	16	-2381.5575	4871.5756	0.0000	1.0000
F81	14	-2391.9169	4878.7367	0.0000	1.0000
JC69G	12	-2407.0026	4895.3506	0.0000	1.0000
JC69I	12	-2407.2866	4895.9186	0.0000	1.0000
JC69IG	13	-2407.0127	4902.1496	0.0000	1.0000
JC69	11	-2416.4144	4907.3955	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 591-tomodel-phylip/chr1_14993.phylip = 1 - 879;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14993.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 591-tomodel-phylip/chr1_14993.phylip = 1 - 879;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14993.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 591-tomodel-phylip/chr1_14993.phylip = 1 - 879;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14993.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

