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

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

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


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
GTRG	20	-2750.8132	5542.7434	0.4975	0.4975
HKYG	16	-2755.8727	5544.4650	0.2103	0.7078
GTRIG	21	-2750.6967	5544.6237	0.1943	0.9021
HKYIG	17	-2755.7580	5546.3267	0.0829	0.9850
GTRI	20	-2754.6704	5550.4578	0.0105	0.9956
HKYI	16	-2759.7302	5552.1799	0.0044	1.0000
K2PG	13	-2781.8717	5590.2230	0.0000	1.0000
K2PIG	14	-2781.8004	5592.1549	0.0000	1.0000
SYMG	17	-2780.1775	5595.1655	0.0000	1.0000
SYMIG	18	-2780.1223	5597.1517	0.0000	1.0000
K2PI	13	-2785.9168	5598.3132	0.0000	1.0000
SYMI	17	-2784.5108	5603.8322	0.0000	1.0000
HKY	15	-2793.3255	5617.2851	0.0000	1.0000
GTR	19	-2789.6487	5618.3068	0.0000	1.0000
K2P	12	-2816.3806	5657.1717	0.0000	1.0000
SYM	16	-2814.8407	5662.4010	0.0000	1.0000
F81G	15	-2843.0493	5716.7326	0.0000	1.0000
F81IG	16	-2842.9649	5718.6493	0.0000	1.0000
F81I	15	-2845.6331	5721.9002	0.0000	1.0000
JC69G	12	-2860.3154	5745.0413	0.0000	1.0000
JC69IG	13	-2860.2514	5746.9824	0.0000	1.0000
JC69I	12	-2863.0576	5750.5257	0.0000	1.0000
F81	14	-2874.2393	5777.0327	0.0000	1.0000
JC69	11	-2890.6339	5803.6147	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2750.8132	5541.6264	0.5178	0.5178
GTRIG	21	-2750.6967	5543.3934	0.2140	0.7318
HKYG	16	-2755.8727	5543.7455	0.1795	0.9112
HKYIG	17	-2755.7580	5545.5161	0.0740	0.9853
GTRI	20	-2754.6704	5549.3408	0.0109	0.9962
HKYI	16	-2759.7302	5551.4603	0.0038	1.0000
K2PG	13	-2781.8717	5589.7434	0.0000	1.0000
K2PIG	14	-2781.8004	5591.6009	0.0000	1.0000
SYMG	17	-2780.1775	5594.3549	0.0000	1.0000
SYMIG	18	-2780.1223	5596.2445	0.0000	1.0000
K2PI	13	-2785.9168	5597.8336	0.0000	1.0000
SYMI	17	-2784.5108	5603.0216	0.0000	1.0000
HKY	15	-2793.3255	5616.6511	0.0000	1.0000
GTR	19	-2789.6487	5617.2975	0.0000	1.0000
K2P	12	-2816.3806	5656.7611	0.0000	1.0000
SYM	16	-2814.8407	5661.6814	0.0000	1.0000
F81G	15	-2843.0493	5716.0985	0.0000	1.0000
F81IG	16	-2842.9649	5717.9297	0.0000	1.0000
F81I	15	-2845.6331	5721.2661	0.0000	1.0000
JC69G	12	-2860.3154	5744.6308	0.0000	1.0000
JC69IG	13	-2860.2514	5746.5029	0.0000	1.0000
JC69I	12	-2863.0576	5750.1152	0.0000	1.0000
F81	14	-2874.2393	5776.4787	0.0000	1.0000
JC69	11	-2890.6339	5803.2678	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2755.8727	5618.1499	0.9418	0.9418
HKYIG	17	-2755.7580	5624.5708	0.0380	0.9798
HKYI	16	-2759.7302	5625.8648	0.0199	0.9997
GTRG	20	-2750.8132	5634.6319	0.0002	1.0000
GTRIG	21	-2750.6967	5641.0492	0.0000	1.0000
GTRI	20	-2754.6704	5642.3464	0.0000	1.0000
K2PG	13	-2781.8717	5650.1971	0.0000	1.0000
K2PIG	14	-2781.8004	5656.7048	0.0000	1.0000
K2PI	13	-2785.9168	5658.2872	0.0000	1.0000
SYMG	17	-2780.1775	5673.4096	0.0000	1.0000
SYMIG	18	-2780.1223	5679.9495	0.0000	1.0000
SYMI	17	-2784.5108	5682.0764	0.0000	1.0000
HKY	15	-2793.3255	5686.4052	0.0000	1.0000
GTR	19	-2789.6487	5705.6528	0.0000	1.0000
K2P	12	-2816.3806	5712.5645	0.0000	1.0000
SYM	16	-2814.8407	5736.0859	0.0000	1.0000
F81G	15	-2843.0493	5785.8527	0.0000	1.0000
F81I	15	-2845.6331	5791.0203	0.0000	1.0000
F81IG	16	-2842.9649	5792.3342	0.0000	1.0000
JC69G	12	-2860.3154	5800.4342	0.0000	1.0000
JC69I	12	-2863.0576	5805.9185	0.0000	1.0000
JC69IG	13	-2860.2514	5806.9565	0.0000	1.0000
F81	14	-2874.2393	5841.5826	0.0000	1.0000
JC69	11	-2890.6339	5854.4209	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/chrZ_3271.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_3271.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/chrZ_3271.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_3271.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/chrZ_3271.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_3271.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

