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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4000.0985	8040.1970	0.4312	0.4312
GTRIG	21	-3999.4919	8040.9839	0.2909	0.7221
GTRG	20	-4000.5392	8041.0783	0.2775	0.9996
HKYI	16	-4011.3490	8054.6979	0.0003	0.9999
HKYG	16	-4013.2935	8058.5870	0.0000	1.0000
HKYIG	17	-4012.8578	8059.7157	0.0000	1.0000
GTR	19	-4034.0282	8106.0564	0.0000	1.0000
HKY	15	-4044.3883	8118.7766	0.0000	1.0000
K2PI	13	-4076.8396	8179.6791	0.0000	1.0000
SYMG	17	-4073.3408	8180.6816	0.0000	1.0000
K2PG	13	-4078.1541	8182.3082	0.0000	1.0000
SYMIG	18	-4073.3432	8182.6863	0.0000	1.0000
K2PIG	14	-4078.1314	8184.2627	0.0000	1.0000
SYMI	17	-4077.7064	8189.4127	0.0000	1.0000
K2P	12	-4108.6094	8241.2189	0.0000	1.0000
SYM	16	-4106.7399	8245.4798	0.0000	1.0000
F81I	15	-4120.6822	8271.3644	0.0000	1.0000
F81G	15	-4121.1426	8272.2852	0.0000	1.0000
F81IG	16	-4120.5247	8273.0493	0.0000	1.0000
F81	14	-4149.5067	8327.0133	0.0000	1.0000
JC69G	12	-4167.4385	8358.8770	0.0000	1.0000
JC69IG	13	-4167.2073	8360.4146	0.0000	1.0000
JC69I	12	-4169.9077	8363.8153	0.0000	1.0000
JC69	11	-4196.8461	8415.6923	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-4011.3490	8134.2896	0.7910	0.7910
HKYG	16	-4013.2935	8138.1787	0.1132	0.9042
GTRI	20	-4000.0985	8139.6866	0.0532	0.9574
GTRG	20	-4000.5392	8140.5679	0.0343	0.9917
HKYIG	17	-4012.8578	8144.2818	0.0054	0.9970
GTRIG	21	-3999.4919	8145.4479	0.0030	1.0000
HKY	15	-4044.3883	8193.3938	0.0000	1.0000
GTR	19	-4034.0282	8200.5715	0.0000	1.0000
K2PI	13	-4076.8396	8244.3473	0.0000	1.0000
K2PG	13	-4078.1541	8246.9765	0.0000	1.0000
K2PIG	14	-4078.1314	8253.9054	0.0000	1.0000
SYMG	17	-4073.3408	8265.2477	0.0000	1.0000
SYMIG	18	-4073.3432	8272.2269	0.0000	1.0000
SYMI	17	-4077.7064	8273.9789	0.0000	1.0000
K2P	12	-4108.6094	8300.9126	0.0000	1.0000
SYM	16	-4106.7399	8325.0715	0.0000	1.0000
F81I	15	-4120.6822	8345.9816	0.0000	1.0000
F81G	15	-4121.1426	8346.9024	0.0000	1.0000
F81IG	16	-4120.5247	8352.6410	0.0000	1.0000
F81	14	-4149.5067	8396.6560	0.0000	1.0000
JC69G	12	-4167.4385	8418.5707	0.0000	1.0000
JC69I	12	-4169.9077	8423.5091	0.0000	1.0000
JC69IG	13	-4167.2073	8425.0828	0.0000	1.0000
JC69	11	-4196.8461	8470.4115	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 591-tomodel-phylip/chr2_24822.phylip = 1 - 1069;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24822.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 591-tomodel-phylip/chr2_24822.phylip = 1 - 1069;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24822.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24822.phylip = 1 - 1069;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24822.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

