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

Input data from file "683-tomodel-phylip/chr5_1674.phylip" (7 taxa, 1132 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	-3222.9132	6477.8264	0.4111	0.4111
HKYG	16	-3223.0327	6478.0655	0.3647	0.7758
HKYIG	17	-3222.9076	6479.8151	0.1521	0.9279
GTRI	20	-3221.2881	6482.5762	0.0382	0.9661
GTRG	20	-3221.7457	6483.4914	0.0242	0.9903
GTRIG	21	-3221.6598	6485.3196	0.0097	1.0000
HKY	15	-3247.8631	6525.7262	0.0000	1.0000
GTR	19	-3246.7318	6531.4636	0.0000	1.0000
K2PG	13	-3266.0042	6558.0083	0.0000	1.0000
K2PI	13	-3266.0358	6558.0716	0.0000	1.0000
K2PIG	14	-3265.9641	6559.9283	0.0000	1.0000
SYMG	17	-3264.0907	6562.1814	0.0000	1.0000
SYMI	17	-3264.2524	6562.5048	0.0000	1.0000
SYMIG	18	-3264.0210	6564.0419	0.0000	1.0000
K2P	12	-3290.2964	6604.5927	0.0000	1.0000
SYM	16	-3287.9900	6607.9801	0.0000	1.0000
F81I	15	-3289.0570	6608.1141	0.0000	1.0000
F81G	15	-3289.8074	6609.6148	0.0000	1.0000
F81IG	16	-3289.0761	6610.1522	0.0000	1.0000
F81	14	-3311.1793	6650.3587	0.0000	1.0000
JC69I	12	-3325.9450	6675.8901	0.0000	1.0000
JC69G	12	-3326.4522	6676.9044	0.0000	1.0000
JC69IG	13	-3327.2027	6680.4054	0.0000	1.0000
JC69	11	-3347.8258	6717.6517	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3222.9132	6558.3343	0.5216	0.5216
HKYG	16	-3223.0327	6558.5733	0.4628	0.9844
HKYIG	17	-3222.9076	6565.3547	0.0156	1.0000
GTRI	20	-3221.2881	6583.2110	0.0000	1.0000
GTRG	20	-3221.7457	6584.1263	0.0000	1.0000
GTRIG	21	-3221.6598	6590.9861	0.0000	1.0000
HKY	15	-3247.8631	6601.2024	0.0000	1.0000
K2PG	13	-3266.0042	6623.4210	0.0000	1.0000
K2PI	13	-3266.0358	6623.4842	0.0000	1.0000
GTR	19	-3246.7318	6627.0667	0.0000	1.0000
K2PIG	14	-3265.9641	6630.3727	0.0000	1.0000
SYMG	17	-3264.0907	6647.7210	0.0000	1.0000
SYMI	17	-3264.2524	6648.0444	0.0000	1.0000
SYMIG	18	-3264.0210	6654.6133	0.0000	1.0000
K2P	12	-3290.2964	6664.9736	0.0000	1.0000
F81I	15	-3289.0570	6683.5902	0.0000	1.0000
F81G	15	-3289.8074	6685.0909	0.0000	1.0000
SYM	16	-3287.9900	6688.4880	0.0000	1.0000
F81IG	16	-3289.0761	6690.6601	0.0000	1.0000
F81	14	-3311.1793	6720.8031	0.0000	1.0000
JC69I	12	-3325.9450	6736.2710	0.0000	1.0000
JC69G	12	-3326.4522	6737.2853	0.0000	1.0000
JC69IG	13	-3327.2027	6745.8181	0.0000	1.0000
JC69	11	-3347.8258	6773.0008	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/chr5_1674.phylip = 1 - 1132;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1674.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/chr5_1674.phylip = 1 - 1132;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1674.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/chr5_1674.phylip = 1 - 1132;
 Partition Dummy = 1:683-tomodel-phylip/chr5_1674.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

