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

Input data from file "683-tomodel-phylip/chr2_6794.phylip" (7 taxa, 1079 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	-2966.3215	5972.6430	0.5600	0.5600
GTRIG	21	-2966.2761	5974.5523	0.2156	0.7756
GTRI	20	-2967.9006	5975.8012	0.1155	0.8910
HKYG	16	-2972.4142	5976.8283	0.0691	0.9601
HKYIG	17	-2972.3778	5978.7555	0.0264	0.9865
HKYI	16	-2974.0437	5980.0875	0.0135	1.0000
GTR	19	-2983.4920	6004.9840	0.0000	1.0000
HKY	15	-2989.0800	6008.1600	0.0000	1.0000
SYMG	17	-2987.5331	6009.0663	0.0000	1.0000
K2PG	13	-2991.9850	6009.9701	0.0000	1.0000
SYMIG	18	-2987.4861	6010.9722	0.0000	1.0000
K2PIG	14	-2991.9468	6011.8936	0.0000	1.0000
SYMI	17	-2989.3315	6012.6631	0.0000	1.0000
K2PI	13	-2993.7669	6013.5338	0.0000	1.0000
SYM	16	-3004.0158	6040.0315	0.0000	1.0000
K2P	12	-3008.3810	6040.7621	0.0000	1.0000
F81G	15	-3050.1346	6130.2692	0.0000	1.0000
F81IG	16	-3050.1213	6132.2425	0.0000	1.0000
F81I	15	-3051.4607	6132.9214	0.0000	1.0000
JC69G	12	-3065.8491	6155.6982	0.0000	1.0000
F81	14	-3064.6906	6157.3811	0.0000	1.0000
JC69IG	13	-3065.8347	6157.6694	0.0000	1.0000
JC69I	12	-3067.2358	6158.4715	0.0000	1.0000
JC69	11	-3080.2569	6182.5138	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2972.4142	6056.5690	0.8142	0.8142
HKYI	16	-2974.0437	6059.8281	0.1596	0.9738
HKYIG	17	-2972.3778	6063.4800	0.0257	0.9995
GTRG	20	-2966.3215	6072.3188	0.0003	0.9998
K2PG	13	-2991.9850	6074.7594	0.0001	0.9999
GTRI	20	-2967.9006	6075.4770	0.0001	1.0000
K2PI	13	-2993.7669	6078.3231	0.0000	1.0000
GTRIG	21	-2966.2761	6079.2118	0.0000	1.0000
K2PIG	14	-2991.9468	6081.6666	0.0000	1.0000
HKY	15	-2989.0800	6082.9168	0.0000	1.0000
SYMG	17	-2987.5331	6093.7907	0.0000	1.0000
SYMI	17	-2989.3315	6097.3875	0.0000	1.0000
GTR	19	-2983.4920	6099.6760	0.0000	1.0000
K2P	12	-3008.3810	6100.5675	0.0000	1.0000
SYMIG	18	-2987.4861	6100.6804	0.0000	1.0000
SYM	16	-3004.0158	6119.7722	0.0000	1.0000
F81G	15	-3050.1346	6205.0260	0.0000	1.0000
F81I	15	-3051.4607	6207.6783	0.0000	1.0000
F81IG	16	-3050.1213	6211.9831	0.0000	1.0000
JC69G	12	-3065.8491	6215.5037	0.0000	1.0000
JC69I	12	-3067.2358	6218.2770	0.0000	1.0000
JC69IG	13	-3065.8347	6222.4586	0.0000	1.0000
F81	14	-3064.6906	6227.1542	0.0000	1.0000
JC69	11	-3080.2569	6237.3354	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/chr2_6794.phylip = 1 - 1079;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6794.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/chr2_6794.phylip = 1 - 1079;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6794.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/chr2_6794.phylip = 1 - 1079;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6794.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

