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

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

Minimum AIC  model: GTRIG
Minimum AICc model: SYMG
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
SYMG	61	-4278.3618	8687.7069	0.5617	0.5617
GTRIG	65	-4274.1064	8688.4515	0.3871	0.9487
GTRG	64	-4277.2932	8692.5030	0.0511	0.9998
HKYG	60	-4288.3838	8705.4508	0.0001	0.9999
K2PG	57	-4291.9647	8705.7450	0.0001	0.9999
HKYIG	61	-4288.0270	8707.0374	0.0000	1.0000
SYMIG	62	-4287.8004	8708.8896	0.0000	1.0000
K2PIG	58	-4293.6288	8711.3570	0.0000	1.0000
GTRI	64	-4353.6058	8845.1281	0.0000	1.0000
SYMI	61	-4361.1788	8853.3409	0.0000	1.0000
HKYI	60	-4364.3667	8857.4167	0.0000	1.0000
K2PI	57	-4369.1637	8860.1430	0.0000	1.0000
F81G	59	-4389.4420	8905.2726	0.0000	1.0000
F81IG	60	-4393.3208	8915.3248	0.0000	1.0000
JC69G	56	-4402.7156	8924.9683	0.0000	1.0000
JC69IG	57	-4402.1425	8926.1006	0.0000	1.0000
JC69I	56	-4466.5344	9052.6060	0.0000	1.0000
F81I	59	-4464.9007	9056.1900	0.0000	1.0000
GTR	63	-4502.0325	9139.6649	0.0000	1.0000
SYM	60	-4515.4641	9159.6114	0.0000	1.0000
HKY	59	-4530.0774	9186.5433	0.0000	1.0000
K2P	56	-4535.1761	9189.8895	0.0000	1.0000
F81	58	-4632.1672	9388.4338	0.0000	1.0000
JC69	55	-4638.4277	9394.1195	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4274.1064	8678.2128	0.5299	0.5299
SYMG	61	-4278.3618	8678.7236	0.4105	0.9404
GTRG	64	-4277.2932	8682.5865	0.0595	0.9999
HKYG	60	-4288.3838	8696.7675	0.0000	0.9999
K2PG	57	-4291.9647	8697.9294	0.0000	1.0000
HKYIG	61	-4288.0270	8698.0541	0.0000	1.0000
SYMIG	62	-4287.8004	8699.6007	0.0000	1.0000
K2PIG	58	-4293.6288	8703.2576	0.0000	1.0000
GTRI	64	-4353.6058	8835.2115	0.0000	1.0000
SYMI	61	-4361.1788	8844.3576	0.0000	1.0000
HKYI	60	-4364.3667	8848.7334	0.0000	1.0000
K2PI	57	-4369.1637	8852.3274	0.0000	1.0000
F81G	59	-4389.4420	8896.8840	0.0000	1.0000
F81IG	60	-4393.3208	8906.6415	0.0000	1.0000
JC69G	56	-4402.7156	8917.4311	0.0000	1.0000
JC69IG	57	-4402.1425	8918.2850	0.0000	1.0000
JC69I	56	-4466.5344	9045.0688	0.0000	1.0000
F81I	59	-4464.9007	9047.8013	0.0000	1.0000
GTR	63	-4502.0325	9130.0649	0.0000	1.0000
SYM	60	-4515.4641	9150.9281	0.0000	1.0000
HKY	59	-4530.0774	9178.1547	0.0000	1.0000
K2P	56	-4535.1761	9182.3523	0.0000	1.0000
F81	58	-4632.1672	9380.3344	0.0000	1.0000
JC69	55	-4638.4277	9386.8553	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-4291.9647	8971.9187	0.5007	0.5007
SYMG	61	-4278.3618	8971.9402	0.4954	0.9961
K2PIG	58	-4293.6288	8982.0537	0.0032	0.9992
HKYG	60	-4288.3838	8985.1773	0.0007	0.9999
GTRG	64	-4277.2932	8990.2236	0.0001	0.9999
GTRIG	65	-4274.1064	8990.6567	0.0000	1.0000
HKYIG	61	-4288.0270	8991.2707	0.0000	1.0000
SYMIG	62	-4287.8004	8997.6241	0.0000	1.0000
K2PI	57	-4369.1637	9126.3167	0.0000	1.0000
HKYI	60	-4364.3667	9137.1432	0.0000	1.0000
SYMI	61	-4361.1788	9137.5742	0.0000	1.0000
GTRI	64	-4353.6058	9142.8486	0.0000	1.0000
F81G	59	-4389.4420	9180.4869	0.0000	1.0000
JC69G	56	-4402.7156	9186.6136	0.0000	1.0000
JC69IG	57	-4402.1425	9192.2743	0.0000	1.0000
F81IG	60	-4393.3208	9195.0513	0.0000	1.0000
JC69I	56	-4466.5344	9314.2512	0.0000	1.0000
F81I	59	-4464.9007	9331.4043	0.0000	1.0000
GTR	63	-4502.0325	9432.8952	0.0000	1.0000
SYM	60	-4515.4641	9439.3379	0.0000	1.0000
K2P	56	-4535.1761	9451.5347	0.0000	1.0000
HKY	59	-4530.0774	9461.7577	0.0000	1.0000
JC69	55	-4638.4277	9651.2309	0.0000	1.0000
F81	58	-4632.1672	9659.1305	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10269_hg18.phylip = 1 - 904;
 Partition Dummy = 1:261-aggressive-phylip/10269_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10269_hg18.phylip = 1 - 904;
 Partition Dummy = 1:261-aggressive-phylip/10269_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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10269_hg18.phylip = 1 - 904;
 Partition Dummy = 1:261-aggressive-phylip/10269_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

