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

Input data from file "683-tomodel-phylip/chr9_7150.phylip" (7 taxa, 1123 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	-3268.3633	6568.7266	0.5058	0.5058
HKYG	16	-3268.9049	6569.8098	0.2943	0.8002
HKYIG	17	-3269.0749	6572.1498	0.0913	0.8915
GTRI	20	-3266.7031	6573.4061	0.0487	0.9403
GTRG	20	-3266.7814	6573.5628	0.0451	0.9853
GTRIG	21	-3266.9030	6575.8060	0.0147	1.0000
HKY	15	-3310.4039	6650.8079	0.0000	1.0000
GTR	19	-3308.3072	6654.6144	0.0000	1.0000
F81I	15	-3324.2015	6678.4030	0.0000	1.0000
F81G	15	-3325.1088	6680.2176	0.0000	1.0000
F81IG	16	-3325.2581	6682.5161	0.0000	1.0000
F81	14	-3362.0565	6752.1130	0.0000	1.0000
SYMI	17	-3375.1630	6784.3260	0.0000	1.0000
SYMG	17	-3375.3201	6784.6401	0.0000	1.0000
SYMIG	18	-3374.9985	6785.9970	0.0000	1.0000
K2PI	13	-3381.9097	6789.8194	0.0000	1.0000
K2PG	13	-3382.4612	6790.9225	0.0000	1.0000
K2PIG	14	-3382.6020	6793.2041	0.0000	1.0000
SYM	16	-3416.8290	6865.6581	0.0000	1.0000
JC69I	12	-3422.3206	6868.6412	0.0000	1.0000
K2P	12	-3423.1962	6870.3925	0.0000	1.0000
JC69G	12	-3423.3259	6870.6519	0.0000	1.0000
JC69IG	13	-3422.3293	6870.6586	0.0000	1.0000
JC69	11	-3460.6986	6943.3972	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3268.3633	6649.1068	0.6264	0.6264
HKYG	16	-3268.9049	6650.1899	0.3644	0.9908
HKYIG	17	-3269.0749	6657.5537	0.0092	1.0000
GTRI	20	-3266.7031	6673.8813	0.0000	1.0000
GTRG	20	-3266.7814	6674.0380	0.0000	1.0000
GTRIG	21	-3266.9030	6681.3049	0.0000	1.0000
HKY	15	-3310.4039	6726.1643	0.0000	1.0000
GTR	19	-3308.3072	6750.0658	0.0000	1.0000
F81I	15	-3324.2015	6753.7594	0.0000	1.0000
F81G	15	-3325.1088	6755.5740	0.0000	1.0000
F81IG	16	-3325.2581	6762.8963	0.0000	1.0000
F81	14	-3362.0565	6822.4456	0.0000	1.0000
K2PI	13	-3381.9097	6855.1282	0.0000	1.0000
K2PG	13	-3382.4612	6856.2313	0.0000	1.0000
K2PIG	14	-3382.6020	6863.5367	0.0000	1.0000
SYMI	17	-3375.1630	6869.7299	0.0000	1.0000
SYMG	17	-3375.3201	6870.0440	0.0000	1.0000
SYMIG	18	-3374.9985	6876.4246	0.0000	1.0000
JC69I	12	-3422.3206	6928.9263	0.0000	1.0000
K2P	12	-3423.1962	6930.6776	0.0000	1.0000
JC69G	12	-3423.3259	6930.9370	0.0000	1.0000
JC69IG	13	-3422.3293	6935.9675	0.0000	1.0000
SYM	16	-3416.8290	6946.0382	0.0000	1.0000
JC69	11	-3460.6986	6998.6586	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/chr9_7150.phylip = 1 - 1123;
 Partition Dummy = 1:683-tomodel-phylip/chr9_7150.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/chr9_7150.phylip = 1 - 1123;
 Partition Dummy = 1:683-tomodel-phylip/chr9_7150.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/chr9_7150.phylip = 1 - 1123;
 Partition Dummy = 1:683-tomodel-phylip/chr9_7150.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

