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

Input data from file "683-tomodel-phylip/chrZ_8050.phylip" (7 taxa, 664 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-2309.8418	4660.9899	0.6208	0.6208
GTRIG	21	-2309.8417	4663.1227	0.2137	0.8345
GTRI	20	-2311.1694	4663.6452	0.1646	0.9990
HKYG	16	-2321.4589	4675.7586	0.0004	0.9994
HKYI	16	-2322.0304	4676.9016	0.0002	0.9996
HKYIG	17	-2321.4726	4677.8926	0.0001	0.9998
GTR	19	-2319.4313	4678.0427	0.0001	0.9999
SYMG	17	-2322.4685	4679.8844	0.0000	0.9999
SYMI	17	-2323.3404	4681.6281	0.0000	1.0000
SYMIG	18	-2322.4684	4681.9973	0.0000	1.0000
K2PG	13	-2328.6653	4683.8906	0.0000	1.0000
K2PI	13	-2329.3073	4685.1746	0.0000	1.0000
K2PIG	14	-2328.6716	4685.9902	0.0000	1.0000
HKY	15	-2330.8820	4692.5047	0.0000	1.0000
SYM	16	-2332.4956	4697.8320	0.0000	1.0000
K2P	12	-2338.7725	4702.0242	0.0000	1.0000
F81G	15	-2404.3964	4839.5336	0.0000	1.0000
F81I	15	-2404.9360	4840.6128	0.0000	1.0000
F81IG	16	-2404.3968	4841.6345	0.0000	1.0000
JC69G	12	-2411.0349	4846.5490	0.0000	1.0000
JC69I	12	-2411.5192	4847.5176	0.0000	1.0000
JC69IG	13	-2411.0418	4848.6436	0.0000	1.0000
F81	14	-2411.1946	4851.0363	0.0000	1.0000
JC69	11	-2418.2395	4858.8839	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2309.8418	4659.6835	0.6119	0.6119
GTRIG	21	-2309.8417	4661.6834	0.2251	0.8370
GTRI	20	-2311.1694	4662.3388	0.1622	0.9992
HKYG	16	-2321.4589	4674.9178	0.0003	0.9995
HKYI	16	-2322.0304	4676.0608	0.0002	0.9997
GTR	19	-2319.4313	4676.8626	0.0001	0.9998
HKYIG	17	-2321.4726	4676.9452	0.0001	0.9999
SYMG	17	-2322.4685	4678.9371	0.0000	1.0000
SYMI	17	-2323.3404	4680.6808	0.0000	1.0000
SYMIG	18	-2322.4684	4680.9369	0.0000	1.0000
K2PG	13	-2328.6653	4683.3306	0.0000	1.0000
K2PI	13	-2329.3073	4684.6146	0.0000	1.0000
K2PIG	14	-2328.6716	4685.3431	0.0000	1.0000
HKY	15	-2330.8820	4691.7640	0.0000	1.0000
SYM	16	-2332.4956	4696.9912	0.0000	1.0000
K2P	12	-2338.7725	4701.5450	0.0000	1.0000
F81G	15	-2404.3964	4838.7929	0.0000	1.0000
F81I	15	-2404.9360	4839.8720	0.0000	1.0000
F81IG	16	-2404.3968	4840.7937	0.0000	1.0000
JC69G	12	-2411.0349	4846.0698	0.0000	1.0000
JC69I	12	-2411.5192	4847.0384	0.0000	1.0000
JC69IG	13	-2411.0418	4848.0836	0.0000	1.0000
F81	14	-2411.1946	4850.3892	0.0000	1.0000
JC69	11	-2418.2395	4858.4790	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2328.6653	4741.8083	0.5815	0.5815
K2PI	13	-2329.3073	4743.0923	0.3060	0.8874
HKYG	16	-2321.4589	4746.8903	0.0458	0.9332
HKYI	16	-2322.0304	4748.0333	0.0259	0.9591
K2PIG	14	-2328.6716	4748.3191	0.0224	0.9815
GTRG	20	-2309.8418	4749.6492	0.0115	0.9931
GTRI	20	-2311.1694	4752.3045	0.0031	0.9961
HKYIG	17	-2321.4726	4753.4160	0.0018	0.9979
SYMG	17	-2322.4685	4755.4079	0.0006	0.9985
K2P	12	-2338.7725	4755.5243	0.0006	0.9991
GTRIG	21	-2309.8417	4756.1474	0.0004	0.9996
SYMI	17	-2323.3404	4757.1516	0.0003	0.9999
HKY	15	-2330.8820	4759.2382	0.0001	1.0000
SYMIG	18	-2322.4684	4761.9059	0.0000	1.0000
GTR	19	-2319.4313	4762.3300	0.0000	1.0000
SYM	16	-2332.4956	4768.9637	0.0000	1.0000
JC69G	12	-2411.0349	4900.0492	0.0000	1.0000
JC69I	12	-2411.5192	4901.0178	0.0000	1.0000
F81G	15	-2404.3964	4906.2671	0.0000	1.0000
JC69IG	13	-2411.0418	4906.5613	0.0000	1.0000
F81I	15	-2404.9360	4907.3463	0.0000	1.0000
JC69	11	-2418.2395	4907.9601	0.0000	1.0000
F81IG	16	-2404.3968	4912.7662	0.0000	1.0000
F81	14	-2411.1946	4913.3652	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 683-tomodel-phylip/chrZ_8050.phylip = 1 - 664;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_8050.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 683-tomodel-phylip/chrZ_8050.phylip = 1 - 664;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_8050.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_8050.phylip = 1 - 664;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_8050.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

