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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2089.3225	4211.3292	0.4589	0.4589
HKYI	16	-2089.7911	4212.2665	0.2872	0.7461
HKYIG	17	-2089.0832	4212.9371	0.2054	0.9515
GTRG	20	-2088.1100	4217.2819	0.0234	0.9749
GTRI	20	-2088.6060	4218.2740	0.0142	0.9891
GTRIG	21	-2087.8431	4218.8557	0.0107	0.9998
K2PG	13	-2100.9478	4228.3517	0.0001	0.9999
K2PI	13	-2101.4281	4229.3124	0.0001	0.9999
K2PIG	14	-2100.7029	4229.9328	0.0000	1.0000
SYMG	17	-2098.8867	4232.5442	0.0000	1.0000
SYMI	17	-2099.5543	4233.8794	0.0000	1.0000
SYMIG	18	-2098.5145	4233.8914	0.0000	1.0000
HKY	15	-2116.0572	4262.7174	0.0000	1.0000
GTR	19	-2114.8110	4268.5816	0.0000	1.0000
K2P	12	-2127.0362	4278.4629	0.0000	1.0000
F81G	15	-2124.5371	4279.6772	0.0000	1.0000
F81I	15	-2124.9932	4280.5893	0.0000	1.0000
F81IG	16	-2124.3669	4281.4180	0.0000	1.0000
SYM	16	-2124.5412	4281.7666	0.0000	1.0000
JC69G	12	-2134.4948	4293.3801	0.0000	1.0000
JC69I	12	-2134.9479	4294.2863	0.0000	1.0000
JC69IG	13	-2134.3290	4295.1141	0.0000	1.0000
F81	14	-2148.3780	4325.2829	0.0000	1.0000
JC69	11	-2157.9529	4338.2359	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2089.3225	4210.6449	0.4500	0.4500
HKYI	16	-2089.7911	4211.5822	0.2816	0.7316
HKYIG	17	-2089.0832	4212.1663	0.2103	0.9419
GTRG	20	-2088.1100	4216.2199	0.0277	0.9696
GTRI	20	-2088.6060	4217.2121	0.0169	0.9865
GTRIG	21	-2087.8431	4217.6861	0.0133	0.9998
K2PG	13	-2100.9478	4227.8955	0.0001	0.9999
K2PI	13	-2101.4281	4228.8562	0.0000	0.9999
K2PIG	14	-2100.7029	4229.4058	0.0000	1.0000
SYMG	17	-2098.8867	4231.7735	0.0000	1.0000
SYMIG	18	-2098.5145	4233.0289	0.0000	1.0000
SYMI	17	-2099.5543	4233.1086	0.0000	1.0000
HKY	15	-2116.0572	4262.1144	0.0000	1.0000
GTR	19	-2114.8110	4267.6220	0.0000	1.0000
K2P	12	-2127.0362	4278.0724	0.0000	1.0000
F81G	15	-2124.5371	4279.0742	0.0000	1.0000
F81I	15	-2124.9932	4279.9863	0.0000	1.0000
F81IG	16	-2124.3669	4280.7338	0.0000	1.0000
SYM	16	-2124.5412	4281.0823	0.0000	1.0000
JC69G	12	-2134.4948	4292.9896	0.0000	1.0000
JC69I	12	-2134.9479	4293.8958	0.0000	1.0000
JC69IG	13	-2134.3290	4294.6579	0.0000	1.0000
F81	14	-2148.3780	4324.7559	0.0000	1.0000
JC69	11	-2157.9529	4337.9059	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2089.3225	4285.8369	0.4964	0.4964
HKYI	16	-2089.7911	4286.7742	0.3107	0.8071
K2PG	13	-2100.9478	4288.9890	0.1027	0.9098
K2PI	13	-2101.4281	4289.9497	0.0635	0.9733
HKYIG	17	-2089.0832	4292.0578	0.0221	0.9954
K2PIG	14	-2100.7029	4295.1988	0.0046	1.0000
GTRG	20	-2088.1100	4310.2099	0.0000	1.0000
GTRI	20	-2088.6060	4311.2021	0.0000	1.0000
SYMG	17	-2098.8867	4311.6650	0.0000	1.0000
SYMI	17	-2099.5543	4313.0001	0.0000	1.0000
GTRIG	21	-2087.8431	4316.3756	0.0000	1.0000
SYMIG	18	-2098.5145	4317.6199	0.0000	1.0000
HKY	15	-2116.0572	4332.6069	0.0000	1.0000
K2P	12	-2127.0362	4334.4664	0.0000	1.0000
JC69G	12	-2134.4948	4349.3836	0.0000	1.0000
F81G	15	-2124.5371	4349.5667	0.0000	1.0000
JC69I	12	-2134.9479	4350.2898	0.0000	1.0000
F81I	15	-2124.9932	4350.4788	0.0000	1.0000
JC69IG	13	-2134.3290	4355.7514	0.0000	1.0000
F81IG	16	-2124.3669	4355.9258	0.0000	1.0000
SYM	16	-2124.5412	4356.2743	0.0000	1.0000
GTR	19	-2114.8110	4356.9125	0.0000	1.0000
JC69	11	-2157.9529	4389.6004	0.0000	1.0000
F81	14	-2148.3780	4390.5489	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_8694.phylip = 1 - 812;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_8694.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_8694.phylip = 1 - 812;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_8694.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/chrZ_8694.phylip = 1 - 812;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_8694.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

