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

Input data from file "683-tomodel-phylip/chr1_24640.phylip" (7 taxa, 1042 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	-2360.6438	4753.2877	0.2574	0.2574
HKYI	16	-2360.6541	4753.3081	0.2548	0.5122
GTRI	20	-2357.0741	4754.1481	0.1674	0.6797
GTRG	20	-2357.0883	4754.1766	0.1651	0.8447
HKYIG	17	-2360.6446	4755.2892	0.0946	0.9393
GTRIG	21	-2357.0893	4756.1785	0.0607	1.0000
SYMG	17	-2374.1337	4782.2673	0.0000	1.0000
SYMI	17	-2374.1559	4782.3119	0.0000	1.0000
SYMIG	18	-2374.1319	4784.2638	0.0000	1.0000
K2PG	13	-2379.4701	4784.9402	0.0000	1.0000
K2PI	13	-2379.4775	4784.9551	0.0000	1.0000
K2PIG	14	-2379.4701	4786.9402	0.0000	1.0000
HKY	15	-2383.5466	4797.0932	0.0000	1.0000
GTR	19	-2380.1689	4798.3379	0.0000	1.0000
F81G	15	-2396.4001	4822.8001	0.0000	1.0000
F81I	15	-2396.4040	4822.8080	0.0000	1.0000
F81IG	16	-2396.4030	4824.8060	0.0000	1.0000
SYM	16	-2397.4194	4826.8389	0.0000	1.0000
K2P	12	-2402.2974	4828.5948	0.0000	1.0000
JC69G	12	-2414.2998	4852.5995	0.0000	1.0000
JC69I	12	-2414.3070	4852.6141	0.0000	1.0000
JC69IG	13	-2414.3020	4854.6039	0.0000	1.0000
F81	14	-2418.1048	4864.2096	0.0000	1.0000
JC69	11	-2436.0280	4894.0559	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2360.6438	4832.4700	0.4947	0.4947
HKYI	16	-2360.6541	4832.4905	0.4897	0.9844
HKYIG	17	-2360.6446	4839.4205	0.0153	0.9997
K2PG	13	-2379.4701	4849.2759	0.0001	0.9999
K2PI	13	-2379.4775	4849.2907	0.0001	1.0000
GTRI	20	-2357.0741	4853.1261	0.0000	1.0000
GTRG	20	-2357.0883	4853.1545	0.0000	1.0000
K2PIG	14	-2379.4701	4856.2247	0.0000	1.0000
GTRIG	21	-2357.0893	4860.1054	0.0000	1.0000
SYMG	17	-2374.1337	4866.3986	0.0000	1.0000
SYMI	17	-2374.1559	4866.4431	0.0000	1.0000
HKY	15	-2383.5466	4871.3267	0.0000	1.0000
SYMIG	18	-2374.1319	4873.3439	0.0000	1.0000
K2P	12	-2402.2974	4887.9815	0.0000	1.0000
GTR	19	-2380.1689	4892.3669	0.0000	1.0000
F81G	15	-2396.4001	4897.0336	0.0000	1.0000
F81I	15	-2396.4040	4897.0415	0.0000	1.0000
F81IG	16	-2396.4030	4903.9883	0.0000	1.0000
SYM	16	-2397.4194	4906.0212	0.0000	1.0000
JC69G	12	-2414.2998	4911.9863	0.0000	1.0000
JC69I	12	-2414.3070	4912.0009	0.0000	1.0000
JC69IG	13	-2414.3020	4918.9396	0.0000	1.0000
F81	14	-2418.1048	4933.4942	0.0000	1.0000
JC69	11	-2436.0280	4948.4938	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/chr1_24640.phylip = 1 - 1042;
 Partition Dummy = 1:683-tomodel-phylip/chr1_24640.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/chr1_24640.phylip = 1 - 1042;
 Partition Dummy = 1:683-tomodel-phylip/chr1_24640.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/chr1_24640.phylip = 1 - 1042;
 Partition Dummy = 1:683-tomodel-phylip/chr1_24640.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

