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

Input data from file "683-tomodel-phylip/chr21_1554.phylip" (7 taxa, 1221 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4963.0563	9966.1126	0.3562	0.3562
GTRIG	21	-4962.0991	9966.1982	0.3413	0.6975
GTRG	20	-4963.2234	9966.4469	0.3014	0.9988
HKYIG	17	-4972.8018	9979.6036	0.0004	0.9992
HKYG	16	-4973.8805	9979.7610	0.0004	0.9996
HKYI	16	-4973.9406	9979.8812	0.0004	1.0000
SYMG	17	-4985.8747	10005.7494	0.0000	1.0000
SYMIG	18	-4985.1882	10006.3764	0.0000	1.0000
SYMI	17	-4986.7697	10007.5393	0.0000	1.0000
K2PI	13	-5006.0436	10038.0872	0.0000	1.0000
K2PIG	14	-5005.1393	10038.2786	0.0000	1.0000
K2PG	13	-5006.2901	10038.5803	0.0000	1.0000
GTR	19	-5042.6367	10123.2733	0.0000	1.0000
HKY	15	-5053.2462	10136.4923	0.0000	1.0000
SYM	16	-5064.4986	10160.9972	0.0000	1.0000
K2P	12	-5083.0306	10190.0611	0.0000	1.0000
F81I	15	-5115.5514	10261.1029	0.0000	1.0000
F81IG	16	-5115.3856	10262.7712	0.0000	1.0000
F81G	15	-5117.6045	10265.2090	0.0000	1.0000
JC69I	12	-5142.3186	10308.6372	0.0000	1.0000
JC69IG	13	-5142.1675	10310.3350	0.0000	1.0000
JC69G	12	-5144.3829	10312.7658	0.0000	1.0000
F81	14	-5186.2455	10400.4911	0.0000	1.0000
JC69	11	-5211.8509	10445.7018	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4973.8805	10061.4798	0.4784	0.4784
HKYI	16	-4973.9406	10061.6000	0.4504	0.9288
HKYIG	17	-4972.8018	10066.4298	0.0403	0.9691
GTRI	20	-4963.0563	10068.2611	0.0161	0.9852
GTRG	20	-4963.2234	10068.5954	0.0136	0.9988
GTRIG	21	-4962.0991	10073.4541	0.0012	1.0000
SYMG	17	-4985.8747	10092.5757	0.0000	1.0000
SYMI	17	-4986.7697	10094.3656	0.0000	1.0000
SYMIG	18	-4985.1882	10098.3101	0.0000	1.0000
K2PI	13	-5006.0436	10104.4837	0.0000	1.0000
K2PG	13	-5006.2901	10104.9768	0.0000	1.0000
K2PIG	14	-5005.1393	10109.7826	0.0000	1.0000
HKY	15	-5053.2462	10213.1037	0.0000	1.0000
GTR	19	-5042.6367	10220.3144	0.0000	1.0000
SYM	16	-5064.4986	10242.7160	0.0000	1.0000
K2P	12	-5083.0306	10251.3502	0.0000	1.0000
F81I	15	-5115.5514	10337.7143	0.0000	1.0000
F81G	15	-5117.6045	10341.8204	0.0000	1.0000
F81IG	16	-5115.3856	10344.4900	0.0000	1.0000
JC69I	12	-5142.3186	10369.9263	0.0000	1.0000
JC69G	12	-5144.3829	10374.0549	0.0000	1.0000
JC69IG	13	-5142.1675	10376.7315	0.0000	1.0000
F81	14	-5186.2455	10471.9950	0.0000	1.0000
JC69	11	-5211.8509	10501.8835	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 683-tomodel-phylip/chr21_1554.phylip = 1 - 1221;
 Partition Dummy = 1:683-tomodel-phylip/chr21_1554.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 683-tomodel-phylip/chr21_1554.phylip = 1 - 1221;
 Partition Dummy = 1:683-tomodel-phylip/chr21_1554.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 (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr21_1554.phylip = 1 - 1221;
 Partition Dummy = 1:683-tomodel-phylip/chr21_1554.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

