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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PI


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	16	-2576.6332	5186.0347	0.3446	0.3446
HKYI	16	-2576.8923	5186.5530	0.2660	0.6106
HKYIG	17	-2576.2828	5187.4312	0.1714	0.7821
K2PI	13	-2581.6252	5189.7624	0.0534	0.8355
GTRG	20	-2574.4535	5190.1002	0.0451	0.8807
K2PG	13	-2581.8444	5190.2008	0.0429	0.9236
K2PIG	14	-2581.0355	5190.6626	0.0341	0.9577
GTRIG	21	-2574.2713	5191.8570	0.0188	0.9764
GTRI	20	-2575.3609	5191.9149	0.0182	0.9946
SYMG	17	-2580.4386	5195.7429	0.0027	0.9973
SYMI	17	-2581.0377	5196.9409	0.0015	0.9988
SYMIG	18	-2580.1880	5197.3449	0.0012	1.0000
HKY	15	-2608.6666	5248.0101	0.0000	1.0000
GTR	19	-2604.9059	5248.8899	0.0000	1.0000
K2P	12	-2612.7308	5249.8998	0.0000	1.0000
SYM	16	-2613.4286	5259.6256	0.0000	1.0000
F81I	15	-2630.9358	5292.5485	0.0000	1.0000
F81G	15	-2631.5696	5293.8162	0.0000	1.0000
F81IG	16	-2630.7765	5294.3214	0.0000	1.0000
JC69I	12	-2635.2384	5294.9150	0.0000	1.0000
JC69G	12	-2635.9623	5296.3628	0.0000	1.0000
JC69IG	13	-2635.0981	5296.7083	0.0000	1.0000
F81	14	-2662.4355	5353.4626	0.0000	1.0000
JC69	11	-2668.0535	5358.4773	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2576.6332	5185.2664	0.3395	0.3395
HKYI	16	-2576.8923	5185.7846	0.2620	0.6014
HKYIG	17	-2576.2828	5186.5656	0.1773	0.7787
GTRG	20	-2574.4535	5188.9070	0.0550	0.8337
K2PI	13	-2581.6252	5189.2504	0.0463	0.8800
K2PG	13	-2581.8444	5189.6888	0.0372	0.9172
K2PIG	14	-2581.0355	5190.0711	0.0307	0.9479
GTRIG	21	-2574.2713	5190.5426	0.0243	0.9722
GTRI	20	-2575.3609	5190.7218	0.0222	0.9944
SYMG	17	-2580.4386	5194.8773	0.0028	0.9972
SYMI	17	-2581.0377	5196.0753	0.0015	0.9987
SYMIG	18	-2580.1880	5196.3761	0.0013	1.0000
HKY	15	-2608.6666	5247.3331	0.0000	1.0000
GTR	19	-2604.9059	5247.8119	0.0000	1.0000
K2P	12	-2612.7308	5249.4616	0.0000	1.0000
SYM	16	-2613.4286	5258.8573	0.0000	1.0000
F81I	15	-2630.9358	5291.8715	0.0000	1.0000
F81G	15	-2631.5696	5293.1392	0.0000	1.0000
F81IG	16	-2630.7765	5293.5531	0.0000	1.0000
JC69I	12	-2635.2384	5294.4768	0.0000	1.0000
JC69G	12	-2635.9623	5295.9246	0.0000	1.0000
JC69IG	13	-2635.0981	5296.1963	0.0000	1.0000
F81	14	-2662.4355	5352.8711	0.0000	1.0000
JC69	11	-2668.0535	5358.1071	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2581.6252	5248.8706	0.5308	0.5308
K2PG	13	-2581.8444	5249.3090	0.4263	0.9571
K2PIG	14	-2581.0355	5254.2775	0.0356	0.9927
HKYG	16	-2576.6332	5258.6451	0.0040	0.9967
HKYI	16	-2576.8923	5259.1634	0.0031	0.9998
HKYIG	17	-2576.2828	5264.5305	0.0002	1.0000
SYMG	17	-2580.4386	5272.8422	0.0000	1.0000
SYMI	17	-2581.0377	5274.0402	0.0000	1.0000
SYMIG	18	-2580.1880	5278.9272	0.0000	1.0000
GTRG	20	-2574.4535	5280.6305	0.0000	1.0000
GTRI	20	-2575.3609	5282.4452	0.0000	1.0000
GTRIG	21	-2574.2713	5286.8522	0.0000	1.0000
K2P	12	-2612.7308	5304.4957	0.0000	1.0000
HKY	15	-2608.6666	5316.1257	0.0000	1.0000
SYM	16	-2613.4286	5332.2360	0.0000	1.0000
GTR	19	-2604.9059	5334.9492	0.0000	1.0000
JC69I	12	-2635.2384	5349.5108	0.0000	1.0000
JC69G	12	-2635.9623	5350.9586	0.0000	1.0000
JC69IG	13	-2635.0981	5355.8165	0.0000	1.0000
F81I	15	-2630.9358	5360.6641	0.0000	1.0000
F81G	15	-2631.5696	5361.9317	0.0000	1.0000
F81IG	16	-2630.7765	5366.9318	0.0000	1.0000
JC69	11	-2668.0535	5408.5549	0.0000	1.0000
F81	14	-2662.4355	5417.0775	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/chrZ_483.phylip = 1 - 725;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_483.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/chrZ_483.phylip = 1 - 725;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_483.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 (K2PI)]

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


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

