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

Input data from file "591-tomodel-phylip/chr2_11806.phylip" (7 taxa, 887 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2934.6988	5909.3976	0.5368	0.5368
GTRIG	21	-2934.2974	5910.5947	0.2950	0.8318
GTRG	20	-2936.4927	5912.9854	0.0893	0.9211
HKYI	16	-2941.4707	5914.9413	0.0336	0.9547
HKYG	16	-2942.0482	5916.0965	0.0188	0.9735
HKYIG	17	-2941.6277	5917.2554	0.0106	0.9841
SYMIG	18	-2941.1819	5918.3638	0.0061	0.9901
SYMI	17	-2942.2670	5918.5341	0.0056	0.9957
SYMG	17	-2942.6532	5919.3064	0.0038	0.9995
K2PI	13	-2949.4784	5924.9567	0.0002	0.9997
K2PG	13	-2949.4954	5924.9908	0.0002	0.9999
K2PIG	14	-2949.4956	5926.9912	0.0001	1.0000
F81I	15	-3010.7652	6051.5305	0.0000	1.0000
F81G	15	-3012.2130	6054.4259	0.0000	1.0000
F81IG	16	-3012.2135	6056.4271	0.0000	1.0000
JC69I	12	-3019.2005	6062.4010	0.0000	1.0000
JC69G	12	-3020.4813	6064.9626	0.0000	1.0000
JC69IG	13	-3020.4819	6066.9637	0.0000	1.0000
GTR	19	-3033.5276	6105.0553	0.0000	1.0000
SYM	16	-3043.8952	6119.7903	0.0000	1.0000
HKY	15	-3045.4000	6120.8001	0.0000	1.0000
K2P	12	-3054.9115	6133.8230	0.0000	1.0000
F81	14	-3108.5811	6245.1622	0.0000	1.0000
JC69	11	-3116.6885	6255.3770	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2949.4784	5987.1987	0.4545	0.4545
K2PG	13	-2949.4954	5987.2327	0.4468	0.9014
HKYI	16	-2941.4707	5991.5469	0.0517	0.9530
HKYG	16	-2942.0482	5992.7020	0.0290	0.9821
K2PIG	14	-2949.4956	5994.0210	0.0150	0.9971
HKYIG	17	-2941.6277	5998.6488	0.0015	0.9985
SYMI	17	-2942.2670	5999.9274	0.0008	0.9993
SYMG	17	-2942.6532	6000.6998	0.0005	0.9999
SYMIG	18	-2941.1819	6004.5450	0.0001	0.9999
GTRI	20	-2934.6988	6005.1545	0.0001	1.0000
GTRG	20	-2936.4927	6008.7423	0.0000	1.0000
GTRIG	21	-2934.2974	6011.1395	0.0000	1.0000
JC69I	12	-3019.2005	6119.8551	0.0000	1.0000
JC69G	12	-3020.4813	6122.4167	0.0000	1.0000
F81I	15	-3010.7652	6123.3482	0.0000	1.0000
F81G	15	-3012.2130	6126.2436	0.0000	1.0000
JC69IG	13	-3020.4819	6129.2057	0.0000	1.0000
F81IG	16	-3012.2135	6133.0326	0.0000	1.0000
K2P	12	-3054.9115	6191.2772	0.0000	1.0000
HKY	15	-3045.4000	6192.6178	0.0000	1.0000
GTR	19	-3033.5276	6196.0243	0.0000	1.0000
SYM	16	-3043.8952	6196.3958	0.0000	1.0000
JC69	11	-3116.6885	6308.0433	0.0000	1.0000
F81	14	-3108.5811	6312.1920	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_11806.phylip = 1 - 887;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11806.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_11806.phylip = 1 - 887;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11806.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (K2PI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_11806.phylip = 1 - 887;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11806.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

