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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3421.7374	6883.4748	0.7216	0.7216
GTRIG	21	-3421.7191	6885.4381	0.2704	0.9920
GTRI	20	-3426.2414	6892.4827	0.0080	1.0000
HKYG	16	-3438.3012	6908.6024	0.0000	1.0000
HKYIG	17	-3438.2883	6910.5767	0.0000	1.0000
HKYI	16	-3442.0945	6916.1891	0.0000	1.0000
SYMG	17	-3454.3124	6942.6248	0.0000	1.0000
SYMIG	18	-3454.2971	6944.5943	0.0000	1.0000
K2PG	13	-3461.6061	6949.2121	0.0000	1.0000
SYMI	17	-3458.2648	6950.5295	0.0000	1.0000
K2PIG	14	-3461.5945	6951.1891	0.0000	1.0000
GTR	19	-3456.9888	6951.9776	0.0000	1.0000
K2PI	13	-3465.2339	6956.4679	0.0000	1.0000
HKY	15	-3471.1566	6972.3133	0.0000	1.0000
SYM	16	-3487.4256	7006.8512	0.0000	1.0000
K2P	12	-3494.6181	7013.2362	0.0000	1.0000
F81G	15	-3520.0177	7070.0355	0.0000	1.0000
F81IG	16	-3520.0212	7072.0425	0.0000	1.0000
F81I	15	-3523.0049	7076.0098	0.0000	1.0000
JC69G	12	-3537.7845	7099.5690	0.0000	1.0000
JC69IG	13	-3537.7829	7101.5657	0.0000	1.0000
JC69I	12	-3540.7348	7105.4696	0.0000	1.0000
F81	14	-3550.2938	7128.5876	0.0000	1.0000
JC69	11	-3568.2519	7158.5039	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3421.7374	6982.0064	0.8984	0.8984
HKYG	16	-3438.3012	6987.4277	0.0597	0.9582
GTRIG	21	-3421.7191	6988.8962	0.0287	0.9868
GTRI	20	-3426.2414	6991.0143	0.0099	0.9968
HKYIG	17	-3438.2883	6994.3285	0.0019	0.9987
HKYI	16	-3442.0945	6995.0143	0.0013	1.0000
K2PG	13	-3461.6061	7013.2576	0.0000	1.0000
K2PIG	14	-3461.5945	7020.1612	0.0000	1.0000
K2PI	13	-3465.2339	7020.5134	0.0000	1.0000
SYMG	17	-3454.3124	7026.3766	0.0000	1.0000
SYMIG	18	-3454.2971	7033.2727	0.0000	1.0000
SYMI	17	-3458.2648	7034.2813	0.0000	1.0000
GTR	19	-3456.9888	7045.5825	0.0000	1.0000
HKY	15	-3471.1566	7046.2120	0.0000	1.0000
K2P	12	-3494.6181	7072.3551	0.0000	1.0000
SYM	16	-3487.4256	7085.6764	0.0000	1.0000
F81G	15	-3520.0177	7143.9341	0.0000	1.0000
F81I	15	-3523.0049	7149.9085	0.0000	1.0000
F81IG	16	-3520.0212	7150.8677	0.0000	1.0000
JC69G	12	-3537.7845	7158.6880	0.0000	1.0000
JC69I	12	-3540.7348	7164.5885	0.0000	1.0000
JC69IG	13	-3537.7829	7165.6112	0.0000	1.0000
F81	14	-3550.2938	7197.5597	0.0000	1.0000
JC69	11	-3568.2519	7212.6962	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/chr4_10251.phylip = 1 - 1019;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10251.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/chr4_10251.phylip = 1 - 1019;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10251.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_10251.phylip = 1 - 1019;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10251.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;


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

