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

Input data from file "683-tomodel-phylip/chrZ_4747.phylip" (7 taxa, 1056 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	-3823.9700	7687.9401	0.4801	0.4801
GTRI	20	-3824.3728	7688.7456	0.3209	0.8011
GTRIG	21	-3823.8519	7689.7039	0.1988	0.9998
HKYG	16	-3836.5824	7705.1648	0.0001	0.9999
HKYI	16	-3837.2949	7706.5898	0.0000	1.0000
HKYIG	17	-3836.5506	7707.1012	0.0000	1.0000
GTR	19	-3840.5040	7719.0080	0.0000	1.0000
HKY	15	-3853.4001	7736.8002	0.0000	1.0000
SYMG	17	-3870.2431	7774.4863	0.0000	1.0000
SYMI	17	-3870.3348	7774.6697	0.0000	1.0000
SYMIG	18	-3870.0212	7776.0425	0.0000	1.0000
K2PG	13	-3878.0895	7782.1790	0.0000	1.0000
K2PI	13	-3878.6053	7783.2107	0.0000	1.0000
K2PIG	14	-3878.0070	7784.0141	0.0000	1.0000
SYM	16	-3887.1637	7806.3273	0.0000	1.0000
K2P	12	-3895.7330	7815.4660	0.0000	1.0000
F81I	15	-3974.0355	7978.0710	0.0000	1.0000
F81G	15	-3974.8447	7979.6893	0.0000	1.0000
F81IG	16	-3974.0366	7980.0733	0.0000	1.0000
F81	14	-3988.6719	8005.3437	0.0000	1.0000
JC69I	12	-4004.0478	8032.0955	0.0000	1.0000
JC69G	12	-4004.8597	8033.7195	0.0000	1.0000
JC69IG	13	-4004.0466	8034.0933	0.0000	1.0000
JC69	11	-4019.2662	8060.5325	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3836.5824	7784.5607	0.5049	0.5049
HKYI	16	-3837.2949	7785.9857	0.2476	0.7525
GTRG	20	-3823.9700	7787.1849	0.1359	0.8884
GTRI	20	-3824.3728	7787.9905	0.0909	0.9793
HKYIG	17	-3836.5506	7791.4594	0.0160	0.9953
GTRIG	21	-3823.8519	7793.9110	0.0047	1.0000
HKY	15	-3853.4001	7811.2338	0.0000	1.0000
GTR	19	-3840.5040	7813.2906	0.0000	1.0000
K2PG	13	-3878.0895	7846.6882	0.0000	1.0000
K2PI	13	-3878.6053	7847.7198	0.0000	1.0000
K2PIG	14	-3878.0070	7853.4855	0.0000	1.0000
SYMG	17	-3870.2431	7858.8444	0.0000	1.0000
SYMI	17	-3870.3348	7859.0278	0.0000	1.0000
SYMIG	18	-3870.0212	7865.3629	0.0000	1.0000
K2P	12	-3895.7330	7875.0130	0.0000	1.0000
SYM	16	-3887.1637	7885.7232	0.0000	1.0000
F81I	15	-3974.0355	8052.5047	0.0000	1.0000
F81G	15	-3974.8447	8054.1230	0.0000	1.0000
F81IG	16	-3974.0366	8059.4692	0.0000	1.0000
F81	14	-3988.6719	8074.8151	0.0000	1.0000
JC69I	12	-4004.0478	8091.6424	0.0000	1.0000
JC69G	12	-4004.8597	8093.2664	0.0000	1.0000
JC69IG	13	-4004.0466	8098.6024	0.0000	1.0000
JC69	11	-4019.2662	8115.1171	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_4747.phylip = 1 - 1056;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_4747.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_4747.phylip = 1 - 1056;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_4747.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 683-tomodel-phylip/chrZ_4747.phylip = 1 - 1056;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_4747.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

