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

Input data from file "683-tomodel-phylip/chr2_24773.phylip" (7 taxa, 1149 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	-3009.9840	6059.9680	0.4495	0.4495
GTRIG	21	-3009.7949	6061.5898	0.1998	0.6493
GTRI	20	-3011.0660	6062.1319	0.1524	0.8016
HKYG	16	-3015.3516	6062.7033	0.1145	0.9161
HKYIG	17	-3015.1466	6064.2933	0.0517	0.9678
HKYI	16	-3016.6208	6065.2415	0.0322	1.0000
K2PG	13	-3035.1716	6096.3431	0.0000	1.0000
K2PIG	14	-3035.0075	6098.0150	0.0000	1.0000
K2PI	13	-3036.4301	6098.8602	0.0000	1.0000
SYMG	17	-3032.5486	6099.0972	0.0000	1.0000
GTR	19	-3031.2615	6100.5230	0.0000	1.0000
SYMIG	18	-3032.4086	6100.8172	0.0000	1.0000
SYMI	17	-3033.6398	6101.2795	0.0000	1.0000
HKY	15	-3036.1520	6102.3041	0.0000	1.0000
K2P	12	-3054.6160	6133.2320	0.0000	1.0000
SYM	16	-3051.6741	6135.3483	0.0000	1.0000
F81G	15	-3096.8993	6223.7985	0.0000	1.0000
F81IG	16	-3096.7924	6225.5848	0.0000	1.0000
F81I	15	-3098.1008	6226.2016	0.0000	1.0000
JC69G	12	-3112.4679	6248.9359	0.0000	1.0000
JC69IG	13	-3112.3670	6250.7339	0.0000	1.0000
JC69I	12	-3113.6907	6251.3814	0.0000	1.0000
F81	14	-3114.8326	6257.6652	0.0000	1.0000
JC69	11	-3130.0604	6282.1208	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3015.3516	6143.4496	0.7589	0.7589
HKYI	16	-3016.6208	6145.9879	0.2133	0.9723
HKYIG	17	-3015.1466	6150.0863	0.0275	0.9997
GTRG	20	-3009.9840	6160.9010	0.0001	0.9999
K2PG	13	-3035.1716	6161.9496	0.0001	0.9999
GTRI	20	-3011.0660	6163.0648	0.0000	1.0000
K2PI	13	-3036.4301	6164.4666	0.0000	1.0000
GTRIG	21	-3009.7949	6167.5694	0.0000	1.0000
K2PIG	14	-3035.0075	6168.6681	0.0000	1.0000
HKY	15	-3036.1520	6178.0038	0.0000	1.0000
SYMG	17	-3032.5486	6184.8902	0.0000	1.0000
SYMI	17	-3033.6398	6187.0725	0.0000	1.0000
SYMIG	18	-3032.4086	6191.6569	0.0000	1.0000
K2P	12	-3054.6160	6193.7917	0.0000	1.0000
GTR	19	-3031.2615	6196.4093	0.0000	1.0000
SYM	16	-3051.6741	6216.0947	0.0000	1.0000
F81G	15	-3096.8993	6299.4982	0.0000	1.0000
F81I	15	-3098.1008	6301.9013	0.0000	1.0000
F81IG	16	-3096.7924	6306.3312	0.0000	1.0000
JC69G	12	-3112.4679	6309.4956	0.0000	1.0000
JC69I	12	-3113.6907	6311.9411	0.0000	1.0000
JC69IG	13	-3112.3670	6316.3404	0.0000	1.0000
F81	14	-3114.8326	6328.3183	0.0000	1.0000
JC69	11	-3130.0604	6337.6339	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_24773.phylip = 1 - 1149;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24773.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_24773.phylip = 1 - 1149;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24773.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_24773.phylip = 1 - 1149;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24773.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

