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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2870.8415	5773.6830	0.4126	0.4126
HKYI	16	-2871.1833	5774.3666	0.2932	0.7058
HKYIG	17	-2870.8452	5775.6904	0.1512	0.8570
GTRG	20	-2868.6140	5777.2280	0.0701	0.9271
GTRI	20	-2869.0119	5778.0239	0.0471	0.9742
GTRIG	21	-2868.6145	5779.2289	0.0258	1.0000
HKY	15	-2882.7972	5795.5943	0.0000	1.0000
GTR	19	-2880.3542	5798.7085	0.0000	1.0000
K2PG	13	-2910.3898	5846.7795	0.0000	1.0000
K2PI	13	-2910.7613	5847.5225	0.0000	1.0000
K2PIG	14	-2910.3814	5848.7628	0.0000	1.0000
SYMG	17	-2909.4763	5852.9527	0.0000	1.0000
SYMI	17	-2909.8852	5853.7705	0.0000	1.0000
SYMIG	18	-2909.4623	5854.9246	0.0000	1.0000
K2P	12	-2924.4087	5872.8174	0.0000	1.0000
SYM	16	-2923.5100	5879.0200	0.0000	1.0000
F81G	15	-2930.5561	5891.1122	0.0000	1.0000
F81I	15	-2930.8412	5891.6824	0.0000	1.0000
F81IG	16	-2930.5644	5893.1289	0.0000	1.0000
F81	14	-2942.3437	5912.6875	0.0000	1.0000
JC69G	12	-2966.7722	5957.5444	0.0000	1.0000
JC69I	12	-2967.1067	5958.2135	0.0000	1.0000
JC69IG	13	-2966.7736	5959.5472	0.0000	1.0000
JC69	11	-2979.6594	5981.3187	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2870.8415	5854.0773	0.5746	0.5746
HKYI	16	-2871.1833	5854.7610	0.4082	0.9828
HKYIG	17	-2870.8452	5861.1095	0.0171	0.9999
HKY	15	-2882.7972	5870.9641	0.0001	1.0000
GTRG	20	-2868.6140	5877.7210	0.0000	1.0000
GTRI	20	-2869.0119	5878.5169	0.0000	1.0000
GTRIG	21	-2868.6145	5884.7466	0.0000	1.0000
GTR	19	-2880.3542	5894.1768	0.0000	1.0000
K2PG	13	-2910.3898	5912.1000	0.0000	1.0000
K2PI	13	-2910.7613	5912.8430	0.0000	1.0000
K2PIG	14	-2910.3814	5919.1078	0.0000	1.0000
K2P	12	-2924.4087	5933.1132	0.0000	1.0000
SYMG	17	-2909.4763	5938.3717	0.0000	1.0000
SYMI	17	-2909.8852	5939.1895	0.0000	1.0000
SYMIG	18	-2909.4623	5945.3683	0.0000	1.0000
SYM	16	-2923.5100	5959.4144	0.0000	1.0000
F81G	15	-2930.5561	5966.4819	0.0000	1.0000
F81I	15	-2930.8412	5967.0522	0.0000	1.0000
F81IG	16	-2930.5644	5973.5233	0.0000	1.0000
F81	14	-2942.3437	5983.0326	0.0000	1.0000
JC69G	12	-2966.7722	6017.8402	0.0000	1.0000
JC69I	12	-2967.1067	6018.5093	0.0000	1.0000
JC69IG	13	-2966.7736	6024.8677	0.0000	1.0000
JC69	11	-2979.6594	6036.5899	0.0000	1.0000



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

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

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


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

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


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

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

