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

Input data from file "591-tomodel-phylip/chr7_10361.phylip" (7 taxa, 931 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	-2516.8793	5065.7586	0.3585	0.3585
HKYI	16	-2517.1647	5066.3294	0.2695	0.6280
HKYIG	17	-2516.8724	5067.7448	0.1328	0.7608
GTRG	20	-2513.9525	5067.9050	0.1226	0.8833
GTRI	20	-2514.5052	5069.0104	0.0705	0.9539
GTRIG	21	-2513.9297	5069.8593	0.0461	1.0000
HKY	15	-2534.5336	5099.0672	0.0000	1.0000
GTR	19	-2531.7422	5101.4845	0.0000	1.0000
K2PG	13	-2560.6070	5147.2141	0.0000	1.0000
K2PI	13	-2560.6217	5147.2433	0.0000	1.0000
SYMG	17	-2557.3929	5148.7858	0.0000	1.0000
SYMI	17	-2557.5089	5149.0178	0.0000	1.0000
K2PIG	14	-2560.6236	5149.2472	0.0000	1.0000
SYMIG	18	-2557.4078	5150.8157	0.0000	1.0000
K2P	12	-2578.9352	5181.8703	0.0000	1.0000
SYM	16	-2575.5939	5183.1878	0.0000	1.0000
F81I	15	-2579.5988	5189.1975	0.0000	1.0000
F81G	15	-2580.4321	5190.8642	0.0000	1.0000
F81IG	16	-2580.4626	5192.9253	0.0000	1.0000
F81	14	-2596.1771	5220.3542	0.0000	1.0000
JC69I	12	-2617.6416	5259.2831	0.0000	1.0000
JC69IG	13	-2617.6658	5261.3316	0.0000	1.0000
JC69G	12	-2618.6713	5261.3425	0.0000	1.0000
JC69	11	-2634.7187	5291.4373	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2516.8793	5143.1388	0.5603	0.5603
HKYI	16	-2517.1647	5143.7095	0.4212	0.9815
HKYIG	17	-2516.8724	5149.9612	0.0185	1.0000
GTRG	20	-2513.9525	5164.6301	0.0000	1.0000
GTRI	20	-2514.5052	5165.7356	0.0000	1.0000
GTRIG	21	-2513.9297	5171.4208	0.0000	1.0000
HKY	15	-2534.5336	5171.6110	0.0000	1.0000
GTR	19	-2531.7422	5193.3734	0.0000	1.0000
K2PG	13	-2560.6070	5210.0854	0.0000	1.0000
K2PI	13	-2560.6217	5210.1147	0.0000	1.0000
K2PIG	14	-2560.6236	5216.9549	0.0000	1.0000
SYMG	17	-2557.3929	5231.0022	0.0000	1.0000
SYMI	17	-2557.5089	5231.2342	0.0000	1.0000
SYMIG	18	-2557.4078	5237.8683	0.0000	1.0000
K2P	12	-2578.9352	5239.9054	0.0000	1.0000
SYM	16	-2575.5939	5260.5679	0.0000	1.0000
F81I	15	-2579.5988	5261.7414	0.0000	1.0000
F81G	15	-2580.4321	5263.4081	0.0000	1.0000
F81IG	16	-2580.4626	5270.3054	0.0000	1.0000
F81	14	-2596.1771	5288.0618	0.0000	1.0000
JC69I	12	-2617.6416	5317.3183	0.0000	1.0000
JC69G	12	-2618.6713	5319.3776	0.0000	1.0000
JC69IG	13	-2617.6658	5324.2030	0.0000	1.0000
JC69	11	-2634.7187	5344.6362	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 591-tomodel-phylip/chr7_10361.phylip = 1 - 931;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10361.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 591-tomodel-phylip/chr7_10361.phylip = 1 - 931;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10361.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 591-tomodel-phylip/chr7_10361.phylip = 1 - 931;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10361.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

