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

Input data from file "591-tomodel-phylip/chr2_22017.phylip" (7 taxa, 1198 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	-4778.7831	9597.5662	0.6042	0.6042
GTRIG	21	-4778.6564	9599.3127	0.2523	0.8565
HKYG	16	-4784.7292	9601.4585	0.0863	0.9428
HKYIG	17	-4784.5613	9603.1227	0.0376	0.9804
GTRI	20	-4782.4697	9604.9393	0.0151	0.9955
HKYI	16	-4787.6867	9607.3734	0.0045	1.0000
K2PG	13	-4803.1472	9632.2943	0.0000	1.0000
K2PIG	14	-4802.9956	9633.9912	0.0000	1.0000
SYMG	17	-4800.9758	9635.9517	0.0000	1.0000
SYMIG	18	-4800.8836	9637.7671	0.0000	1.0000
K2PI	13	-4805.9491	9637.8983	0.0000	1.0000
SYMI	17	-4804.0946	9642.1892	0.0000	1.0000
GTR	19	-4826.8798	9691.7596	0.0000	1.0000
HKY	15	-4831.1008	9692.2015	0.0000	1.0000
K2P	12	-4848.0025	9720.0049	0.0000	1.0000
SYM	16	-4846.8100	9725.6200	0.0000	1.0000
F81G	15	-4862.3631	9754.7263	0.0000	1.0000
F81IG	16	-4862.1636	9756.3271	0.0000	1.0000
F81I	15	-4864.0709	9758.1418	0.0000	1.0000
JC69G	12	-4875.5396	9775.0793	0.0000	1.0000
JC69IG	13	-4875.3525	9776.7050	0.0000	1.0000
JC69I	12	-4877.2492	9778.4985	0.0000	1.0000
F81	14	-4904.1520	9836.3040	0.0000	1.0000
JC69	11	-4916.8799	9855.7598	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4784.7292	9682.8730	0.9201	0.9201
HKYI	16	-4787.6867	9688.7880	0.0478	0.9679
HKYIG	17	-4784.5613	9689.6256	0.0314	0.9993
K2PG	13	-4803.1472	9698.4436	0.0004	0.9997
GTRG	20	-4778.7831	9699.3344	0.0002	0.9999
K2PI	13	-4805.9491	9704.0476	0.0000	1.0000
K2PIG	14	-4802.9956	9705.2289	0.0000	1.0000
GTRIG	21	-4778.6564	9706.1693	0.0000	1.0000
GTRI	20	-4782.4697	9706.7075	0.0000	1.0000
SYMG	17	-4800.9758	9722.4546	0.0000	1.0000
SYMI	17	-4804.0946	9728.6921	0.0000	1.0000
SYMIG	18	-4800.8836	9729.3585	0.0000	1.0000
HKY	15	-4831.1008	9768.5277	0.0000	1.0000
K2P	12	-4848.0025	9781.0658	0.0000	1.0000
GTR	19	-4826.8798	9788.4393	0.0000	1.0000
SYM	16	-4846.8100	9807.0346	0.0000	1.0000
F81G	15	-4862.3631	9831.0524	0.0000	1.0000
F81I	15	-4864.0709	9834.4679	0.0000	1.0000
JC69G	12	-4875.5396	9836.1402	0.0000	1.0000
F81IG	16	-4862.1636	9837.7416	0.0000	1.0000
JC69I	12	-4877.2492	9839.5594	0.0000	1.0000
JC69IG	13	-4875.3525	9842.8543	0.0000	1.0000
F81	14	-4904.1520	9907.5418	0.0000	1.0000
JC69	11	-4916.8799	9911.7323	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 591-tomodel-phylip/chr2_22017.phylip = 1 - 1198;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22017.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 591-tomodel-phylip/chr2_22017.phylip = 1 - 1198;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22017.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 591-tomodel-phylip/chr2_22017.phylip = 1 - 1198;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22017.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

