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

Input data from file "591-tomodel-phylip/chr8_199.phylip" (7 taxa, 908 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	-3205.9145	6451.8289	0.6138	0.6138
GTRIG	21	-3205.9145	6453.8290	0.2258	0.8396
HKYG	16	-3211.8603	6455.7206	0.0877	0.9272
HKYIG	17	-3211.8603	6457.7206	0.0323	0.9595
GTRI	20	-3208.8680	6457.7361	0.0320	0.9915
HKYI	16	-3214.1956	6460.3913	0.0085	1.0000
GTR	19	-3224.3005	6486.6009	0.0000	1.0000
HKY	15	-3229.0447	6488.0894	0.0000	1.0000
K2PG	13	-3283.4858	6592.9716	0.0000	1.0000
K2PIG	14	-3283.4862	6594.9724	0.0000	1.0000
SYMG	17	-3282.3349	6598.6697	0.0000	1.0000
SYMIG	18	-3282.3348	6600.6697	0.0000	1.0000
K2PI	13	-3287.3817	6600.7635	0.0000	1.0000
SYMI	17	-3284.8996	6603.7992	0.0000	1.0000
K2P	12	-3301.7415	6627.4829	0.0000	1.0000
SYM	16	-3300.3831	6632.7662	0.0000	1.0000
F81G	15	-3306.0050	6642.0100	0.0000	1.0000
F81IG	16	-3306.0051	6644.0101	0.0000	1.0000
F81I	15	-3307.6591	6645.3182	0.0000	1.0000
F81	14	-3320.9779	6669.9558	0.0000	1.0000
JC69G	12	-3360.9464	6745.8929	0.0000	1.0000
JC69IG	13	-3360.9474	6747.8948	0.0000	1.0000
JC69I	12	-3362.4653	6748.9305	0.0000	1.0000
JC69	11	-3376.4486	6774.8971	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3211.8603	6532.7005	0.8846	0.8846
HKYI	16	-3214.1956	6537.3712	0.0856	0.9702
HKYIG	17	-3211.8603	6539.5117	0.0294	0.9996
GTRG	20	-3205.9145	6548.0538	0.0004	1.0000
GTRI	20	-3208.8680	6553.9610	0.0000	1.0000
GTRIG	21	-3205.9145	6554.8651	0.0000	1.0000
HKY	15	-3229.0447	6560.2580	0.0000	1.0000
GTR	19	-3224.3005	6578.0146	0.0000	1.0000
K2PG	13	-3283.4858	6655.5178	0.0000	1.0000
K2PIG	14	-3283.4862	6662.3298	0.0000	1.0000
K2PI	13	-3287.3817	6663.3097	0.0000	1.0000
SYMG	17	-3282.3349	6680.4609	0.0000	1.0000
K2P	12	-3301.7415	6685.2179	0.0000	1.0000
SYMI	17	-3284.8996	6685.5903	0.0000	1.0000
SYMIG	18	-3282.3348	6687.2721	0.0000	1.0000
SYM	16	-3300.3831	6709.7462	0.0000	1.0000
F81G	15	-3306.0050	6714.1787	0.0000	1.0000
F81I	15	-3307.6591	6717.4869	0.0000	1.0000
F81IG	16	-3306.0051	6720.9900	0.0000	1.0000
F81	14	-3320.9779	6737.3132	0.0000	1.0000
JC69G	12	-3360.9464	6803.6278	0.0000	1.0000
JC69I	12	-3362.4653	6806.6655	0.0000	1.0000
JC69IG	13	-3360.9474	6810.4410	0.0000	1.0000
JC69	11	-3376.4486	6827.8208	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/chr8_199.phylip = 1 - 908;
 Partition Dummy = 1:591-tomodel-phylip/chr8_199.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/chr8_199.phylip = 1 - 908;
 Partition Dummy = 1:591-tomodel-phylip/chr8_199.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/chr8_199.phylip = 1 - 908;
 Partition Dummy = 1:591-tomodel-phylip/chr8_199.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

