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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-4314.0804	8668.1607	0.6061	0.6061
GTRIG	21	-4314.0496	8670.0993	0.2299	0.8360
GTRI	20	-4315.4063	8670.8126	0.1609	0.9969
HKYG	16	-4323.8045	8679.6091	0.0020	0.9989
HKYIG	17	-4323.8048	8681.6096	0.0007	0.9996
HKYI	16	-4325.5052	8683.0104	0.0004	1.0000
SYMG	17	-4341.0265	8716.0530	0.0000	1.0000
GTR	19	-4339.6397	8717.2795	0.0000	1.0000
K2PG	13	-4345.6470	8717.2940	0.0000	1.0000
SYMIG	18	-4341.0097	8718.0194	0.0000	1.0000
SYMI	17	-4342.4933	8718.9866	0.0000	1.0000
K2PIG	14	-4345.6492	8719.2983	0.0000	1.0000
K2PI	13	-4347.4233	8720.8465	0.0000	1.0000
HKY	15	-4348.2577	8726.5154	0.0000	1.0000
SYM	16	-4364.4730	8760.9460	0.0000	1.0000
K2P	12	-4369.3886	8762.7772	0.0000	1.0000
F81G	15	-4410.8317	8851.6634	0.0000	1.0000
F81I	15	-4411.7366	8853.4732	0.0000	1.0000
F81IG	16	-4410.7858	8853.5716	0.0000	1.0000
JC69G	12	-4427.1549	8878.3098	0.0000	1.0000
JC69IG	13	-4427.1135	8880.2269	0.0000	1.0000
JC69I	12	-4428.1629	8880.3258	0.0000	1.0000
F81	14	-4431.9172	8891.8344	0.0000	1.0000
JC69	11	-4448.2048	8918.4096	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4323.8045	8760.6862	0.8147	0.8147
HKYI	16	-4325.5052	8764.0875	0.1487	0.9634
HKYIG	17	-4323.8048	8767.7540	0.0238	0.9872
GTRG	20	-4314.0804	8769.5071	0.0099	0.9971
GTRI	20	-4315.4063	8772.1590	0.0026	0.9997
GTRIG	21	-4314.0496	8776.5130	0.0003	1.0000
K2PG	13	-4345.6470	8783.1692	0.0000	1.0000
K2PI	13	-4347.4233	8786.7217	0.0000	1.0000
K2PIG	14	-4345.6492	8790.2408	0.0000	1.0000
SYMG	17	-4341.0265	8802.1974	0.0000	1.0000
HKY	15	-4348.2577	8802.5252	0.0000	1.0000
SYMI	17	-4342.4933	8805.1311	0.0000	1.0000
SYMIG	18	-4341.0097	8809.2312	0.0000	1.0000
GTR	19	-4339.6397	8813.5585	0.0000	1.0000
K2P	12	-4369.3886	8823.5850	0.0000	1.0000
SYM	16	-4364.4730	8842.0231	0.0000	1.0000
F81G	15	-4410.8317	8927.6732	0.0000	1.0000
F81I	15	-4411.7366	8929.4830	0.0000	1.0000
F81IG	16	-4410.7858	8934.6487	0.0000	1.0000
JC69G	12	-4427.1549	8939.1177	0.0000	1.0000
JC69I	12	-4428.1629	8941.1337	0.0000	1.0000
JC69IG	13	-4427.1135	8946.1021	0.0000	1.0000
F81	14	-4431.9172	8962.7769	0.0000	1.0000
JC69	11	-4448.2048	8974.1501	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 683-tomodel-phylip/chr2_19463.phylip = 1 - 1173;
 Partition Dummy = 1:683-tomodel-phylip/chr2_19463.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 683-tomodel-phylip/chr2_19463.phylip = 1 - 1173;
 Partition Dummy = 1:683-tomodel-phylip/chr2_19463.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 (HKYG)]

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

