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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3757.9929	7555.9857	0.6414	0.6414
GTRIG	21	-3757.9972	7557.9944	0.2350	0.8764
GTRI	20	-3759.6397	7559.2793	0.1236	1.0000
HKYG	16	-3772.4510	7576.9021	0.0000	1.0000
HKYIG	17	-3772.4744	7578.9487	0.0000	1.0000
HKYI	16	-3773.7792	7579.5583	0.0000	1.0000
GTR	19	-3786.4561	7610.9123	0.0000	1.0000
HKY	15	-3799.7100	7629.4199	0.0000	1.0000
SYMG	17	-3828.5220	7691.0440	0.0000	1.0000
K2PG	13	-3832.6157	7691.2314	0.0000	1.0000
SYMIG	18	-3828.5220	7693.0439	0.0000	1.0000
K2PIG	14	-3832.6157	7693.2314	0.0000	1.0000
K2PI	13	-3833.9720	7693.9440	0.0000	1.0000
SYMI	17	-3830.0979	7694.1959	0.0000	1.0000
F81G	15	-3854.4933	7738.9867	0.0000	1.0000
F81IG	16	-3854.5194	7741.0389	0.0000	1.0000
SYM	16	-3854.5782	7741.1564	0.0000	1.0000
F81I	15	-3855.6842	7741.3685	0.0000	1.0000
K2P	12	-3859.8732	7743.7464	0.0000	1.0000
F81	14	-3879.6341	7787.2683	0.0000	1.0000
JC69G	12	-3901.1881	7826.3762	0.0000	1.0000
JC69IG	13	-3901.2074	7828.4147	0.0000	1.0000
JC69I	12	-3902.5020	7829.0040	0.0000	1.0000
JC69	11	-3926.7602	7875.5204	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3757.9929	7655.8277	0.4929	0.4929
HKYG	16	-3772.4510	7656.7756	0.3068	0.7998
GTRI	20	-3759.6397	7659.1212	0.0950	0.8947
HKYI	16	-3773.7792	7659.4319	0.0813	0.9760
GTRIG	21	-3757.9972	7662.8284	0.0149	0.9909
HKYIG	17	-3772.4744	7663.8144	0.0091	1.0000
HKY	15	-3799.7100	7704.3014	0.0000	1.0000
GTR	19	-3786.4561	7705.7621	0.0000	1.0000
K2PG	13	-3832.6157	7756.1286	0.0000	1.0000
K2PI	13	-3833.9720	7758.8413	0.0000	1.0000
K2PIG	14	-3832.6157	7763.1207	0.0000	1.0000
SYMG	17	-3828.5220	7775.9097	0.0000	1.0000
SYMI	17	-3830.0979	7779.0615	0.0000	1.0000
SYMIG	18	-3828.5220	7782.9017	0.0000	1.0000
K2P	12	-3859.8732	7803.6516	0.0000	1.0000
F81G	15	-3854.4933	7813.8681	0.0000	1.0000
F81I	15	-3855.6842	7816.2499	0.0000	1.0000
F81IG	16	-3854.5194	7820.9124	0.0000	1.0000
SYM	16	-3854.5782	7821.0299	0.0000	1.0000
F81	14	-3879.6341	7857.1576	0.0000	1.0000
JC69G	12	-3901.1881	7886.2814	0.0000	1.0000
JC69I	12	-3902.5020	7888.9092	0.0000	1.0000
JC69IG	13	-3901.2074	7893.3120	0.0000	1.0000
JC69	11	-3926.7602	7930.4335	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 683-tomodel-phylip/chrZ_6575.phylip = 1 - 1088;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6575.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 683-tomodel-phylip/chrZ_6575.phylip = 1 - 1088;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6575.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 (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6575.phylip = 1 - 1088;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6575.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;


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

