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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3319.1008	6678.2016	0.6077	0.6077
GTRIG	21	-3319.0718	6680.1435	0.2301	0.8378
GTRI	20	-3320.7133	6681.4266	0.1212	0.9590
HKYG	16	-3326.3305	6684.6610	0.0240	0.9830
HKYIG	17	-3326.2630	6686.5260	0.0095	0.9925
HKYI	16	-3327.4955	6686.9909	0.0075	1.0000
GTR	19	-3339.5403	6717.0805	0.0000	1.0000
K2PG	13	-3346.5418	6719.0837	0.0000	1.0000
SYMG	17	-3342.9689	6719.9378	0.0000	1.0000
K2PIG	14	-3346.4620	6720.9240	0.0000	1.0000
K2PI	13	-3347.5513	6721.1025	0.0000	1.0000
SYMIG	18	-3342.9397	6721.8793	0.0000	1.0000
SYMI	17	-3344.3499	6722.6998	0.0000	1.0000
HKY	15	-3346.5822	6723.1645	0.0000	1.0000
SYM	16	-3361.8691	6755.7381	0.0000	1.0000
K2P	12	-3366.6633	6757.3266	0.0000	1.0000
F81G	15	-3440.8367	6911.6734	0.0000	1.0000
F81I	15	-3441.0918	6912.1837	0.0000	1.0000
F81IG	16	-3440.5960	6913.1919	0.0000	1.0000
JC69G	12	-3455.0500	6934.1000	0.0000	1.0000
JC69I	12	-3455.4844	6934.9689	0.0000	1.0000
JC69IG	13	-3454.9117	6935.8234	0.0000	1.0000
F81	14	-3457.8924	6943.7847	0.0000	1.0000
JC69	11	-3471.8262	6965.6524	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3326.3305	6760.9937	0.7405	0.7405
HKYI	16	-3327.4955	6763.3236	0.2310	0.9715
HKYIG	17	-3326.2630	6767.6294	0.0268	0.9983
GTRG	20	-3319.1008	6773.6174	0.0013	0.9996
GTRI	20	-3320.7133	6776.8424	0.0003	0.9999
GTRIG	21	-3319.0718	6780.3301	0.0000	1.0000
K2PG	13	-3346.5418	6781.1039	0.0000	1.0000
K2PI	13	-3347.5513	6783.1228	0.0000	1.0000
K2PIG	14	-3346.4620	6787.7151	0.0000	1.0000
HKY	15	-3346.5822	6794.7263	0.0000	1.0000
SYMG	17	-3342.9689	6801.0412	0.0000	1.0000
SYMI	17	-3344.3499	6803.8033	0.0000	1.0000
GTR	19	-3339.5403	6807.7255	0.0000	1.0000
SYMIG	18	-3342.9397	6807.7535	0.0000	1.0000
K2P	12	-3366.6633	6814.5761	0.0000	1.0000
SYM	16	-3361.8691	6832.0708	0.0000	1.0000
F81G	15	-3440.8367	6983.2352	0.0000	1.0000
F81I	15	-3441.0918	6983.7455	0.0000	1.0000
F81IG	16	-3440.5960	6989.5246	0.0000	1.0000
JC69G	12	-3455.0500	6991.3495	0.0000	1.0000
JC69I	12	-3455.4844	6992.2184	0.0000	1.0000
JC69IG	13	-3454.9117	6997.8437	0.0000	1.0000
F81	14	-3457.8924	7010.5758	0.0000	1.0000
JC69	11	-3471.8262	7018.1311	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5924.phylip = 1 - 872;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5924.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5924.phylip = 1 - 872;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5924.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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/chr3_5924.phylip = 1 - 872;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5924.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

