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

Input data from file "261-aggressive-phylip/9158_hg18.phylip" (29 taxa, 962 characters)

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	64	-3977.4770	8092.2294	0.5705	0.5705
GTRIG	65	-3976.6133	8092.8025	0.4283	0.9988
SYMG	61	-3988.0840	8106.5724	0.0004	0.9993
SYMIG	62	-3987.3439	8107.3774	0.0003	0.9995
HKYG	60	-3989.8172	8107.7588	0.0002	0.9998
HKYIG	61	-3989.0093	8108.4231	0.0002	1.0000
K2PG	57	-3995.6360	8112.5861	0.0000	1.0000
K2PIG	58	-3994.8660	8113.3112	0.0000	1.0000
GTRI	64	-4001.7584	8140.7922	0.0000	1.0000
SYMI	61	-4009.9519	8150.3083	0.0000	1.0000
HKYI	60	-4011.8409	8151.8060	0.0000	1.0000
K2PI	57	-4016.8292	8154.9725	0.0000	1.0000
SYM	60	-4084.9649	8298.0541	0.0000	1.0000
GTR	63	-4083.4942	8301.9683	0.0000	1.0000
K2P	56	-4092.1230	8303.3001	0.0000	1.0000
HKY	59	-4093.6009	8313.0510	0.0000	1.0000
JC69G	56	-4109.6450	8338.3442	0.0000	1.0000
F81G	59	-4106.5486	8338.9465	0.0000	1.0000
JC69IG	57	-4109.0428	8339.3998	0.0000	1.0000
F81IG	60	-4105.9383	8340.0009	0.0000	1.0000
JC69I	56	-4129.8395	8378.7331	0.0000	1.0000
F81I	59	-4126.9267	8379.7027	0.0000	1.0000
F81	58	-4206.4256	8536.4304	0.0000	1.0000
JC69	55	-4209.9521	8536.7034	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3977.4770	8082.9540	0.5336	0.5336
GTRIG	65	-3976.6133	8083.2266	0.4656	0.9993
SYMG	61	-3988.0840	8098.1680	0.0003	0.9995
SYMIG	62	-3987.3439	8098.6878	0.0002	0.9998
HKYG	60	-3989.8172	8099.6345	0.0001	0.9999
HKYIG	61	-3989.0093	8100.0186	0.0001	1.0000
K2PG	57	-3995.6360	8105.2720	0.0000	1.0000
K2PIG	58	-3994.8660	8105.7320	0.0000	1.0000
GTRI	64	-4001.7584	8131.5168	0.0000	1.0000
SYMI	61	-4009.9519	8141.9038	0.0000	1.0000
HKYI	60	-4011.8409	8143.6817	0.0000	1.0000
K2PI	57	-4016.8292	8147.6584	0.0000	1.0000
SYM	60	-4084.9649	8289.9298	0.0000	1.0000
GTR	63	-4083.4942	8292.9883	0.0000	1.0000
K2P	56	-4092.1230	8296.2460	0.0000	1.0000
HKY	59	-4093.6009	8305.2017	0.0000	1.0000
F81G	59	-4106.5486	8331.0973	0.0000	1.0000
JC69G	56	-4109.6450	8331.2900	0.0000	1.0000
F81IG	60	-4105.9383	8331.8766	0.0000	1.0000
JC69IG	57	-4109.0428	8332.0856	0.0000	1.0000
JC69I	56	-4129.8395	8371.6789	0.0000	1.0000
F81I	59	-4126.9267	8371.8535	0.0000	1.0000
F81	58	-4206.4256	8528.8512	0.0000	1.0000
JC69	55	-4209.9521	8529.9042	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3995.6360	8382.8058	0.9200	0.9200
K2PIG	58	-3994.8660	8388.1348	0.0641	0.9841
HKYG	60	-3989.8172	8391.7753	0.0104	0.9945
GTRG	64	-3977.4770	8394.5709	0.0026	0.9970
SYMG	61	-3988.0840	8395.1779	0.0019	0.9989
HKYIG	61	-3989.0093	8397.0285	0.0008	0.9997
GTRIG	65	-3976.6133	8399.7126	0.0002	0.9999
SYMIG	62	-3987.3439	8400.5667	0.0001	1.0000
K2PI	57	-4016.8292	8425.1922	0.0000	1.0000
HKYI	60	-4011.8409	8435.8226	0.0000	1.0000
SYMI	61	-4009.9519	8438.9137	0.0000	1.0000
GTRI	64	-4001.7584	8443.1337	0.0000	1.0000
K2P	56	-4092.1230	8568.9108	0.0000	1.0000
SYM	60	-4084.9649	8582.0707	0.0000	1.0000
HKY	59	-4093.6009	8592.4736	0.0000	1.0000
GTR	63	-4083.4942	8599.7362	0.0000	1.0000
JC69G	56	-4109.6450	8603.9548	0.0000	1.0000
JC69IG	57	-4109.0428	8609.6194	0.0000	1.0000
F81G	59	-4106.5486	8618.3691	0.0000	1.0000
F81IG	60	-4105.9383	8624.0174	0.0000	1.0000
JC69I	56	-4129.8395	8644.3437	0.0000	1.0000
F81I	59	-4126.9267	8659.1253	0.0000	1.0000
JC69	55	-4209.9521	8797.7000	0.0000	1.0000
F81	58	-4206.4256	8811.2541	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 261-aggressive-phylip/9158_hg18.phylip = 1 - 962;
 Partition Dummy = 1:261-aggressive-phylip/9158_hg18.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 261-aggressive-phylip/9158_hg18.phylip = 1 - 962;
 Partition Dummy = 1:261-aggressive-phylip/9158_hg18.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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9158_hg18.phylip = 1 - 962;
 Partition Dummy = 1:261-aggressive-phylip/9158_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

