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

Input data from file "591-tomodel-phylip/chr2_8593.phylip" (7 taxa, 532 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	-1897.7049	3828.4660	0.3928	0.3928
HKYI	16	-1897.8030	3828.6623	0.3561	0.7489
HKYIG	17	-1897.4250	3830.0406	0.1787	0.9276
GTRG	20	-1895.8903	3833.4245	0.0329	0.9605
GTRI	20	-1896.1522	3833.9482	0.0253	0.9859
GTRIG	21	-1895.6511	3835.1141	0.0141	1.0000
HKY	15	-1916.2034	3863.3371	0.0000	1.0000
GTR	19	-1914.4554	3868.3951	0.0000	1.0000
F81G	15	-1929.8934	3890.7171	0.0000	1.0000
F81I	15	-1929.9157	3890.7616	0.0000	1.0000
F81IG	16	-1929.6656	3892.3876	0.0000	1.0000
F81	14	-1945.8707	3920.5538	0.0000	1.0000
SYMG	17	-1955.8265	3946.8436	0.0000	1.0000
SYMI	17	-1956.5668	3948.3243	0.0000	1.0000
SYMIG	18	-1955.8097	3948.9527	0.0000	1.0000
K2PG	13	-1961.9357	3950.5741	0.0000	1.0000
K2PI	13	-1962.3062	3951.3151	0.0000	1.0000
K2PIG	14	-1961.8832	3952.5787	0.0000	1.0000
SYM	16	-1974.2639	3981.5840	0.0000	1.0000
K2P	12	-1979.7067	3984.0145	0.0000	1.0000
JC69G	12	-1983.3780	3991.3572	0.0000	1.0000
JC69I	12	-1983.6218	3991.8447	0.0000	1.0000
JC69IG	13	-1983.2789	3993.2605	0.0000	1.0000
JC69	11	-1999.7267	4021.9611	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1897.7049	3827.4097	0.3781	0.3781
HKYI	16	-1897.8030	3827.6059	0.3428	0.7209
HKYIG	17	-1897.4250	3828.8500	0.1840	0.9049
GTRG	20	-1895.8903	3831.7807	0.0425	0.9474
GTRI	20	-1896.1522	3832.3044	0.0327	0.9801
GTRIG	21	-1895.6511	3833.3023	0.0199	1.0000
HKY	15	-1916.2034	3862.4069	0.0000	1.0000
GTR	19	-1914.4554	3866.9107	0.0000	1.0000
F81G	15	-1929.8934	3889.7868	0.0000	1.0000
F81I	15	-1929.9157	3889.8314	0.0000	1.0000
F81IG	16	-1929.6656	3891.3312	0.0000	1.0000
F81	14	-1945.8707	3919.7414	0.0000	1.0000
SYMG	17	-1955.8265	3945.6530	0.0000	1.0000
SYMI	17	-1956.5668	3947.1336	0.0000	1.0000
SYMIG	18	-1955.8097	3947.6193	0.0000	1.0000
K2PG	13	-1961.9357	3949.8714	0.0000	1.0000
K2PI	13	-1962.3062	3950.6124	0.0000	1.0000
K2PIG	14	-1961.8832	3951.7664	0.0000	1.0000
SYM	16	-1974.2639	3980.5277	0.0000	1.0000
K2P	12	-1979.7067	3983.4134	0.0000	1.0000
JC69G	12	-1983.3780	3990.7560	0.0000	1.0000
JC69I	12	-1983.6218	3991.2435	0.0000	1.0000
JC69IG	13	-1983.2789	3992.5578	0.0000	1.0000
JC69	11	-1999.7267	4021.4534	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1897.7049	3895.8360	0.5092	0.5092
HKYI	16	-1897.8030	3896.0322	0.4616	0.9708
HKYIG	17	-1897.4250	3901.5529	0.0292	1.0000
GTRG	20	-1895.8903	3917.3136	0.0000	1.0000
GTRI	20	-1896.1522	3917.8373	0.0000	1.0000
GTRIG	21	-1895.6511	3923.1118	0.0000	1.0000
HKY	15	-1916.2034	3926.5565	0.0000	1.0000
GTR	19	-1914.4554	3948.1669	0.0000	1.0000
F81G	15	-1929.8934	3953.9365	0.0000	1.0000
F81I	15	-1929.9157	3953.9811	0.0000	1.0000
F81IG	16	-1929.6656	3959.7575	0.0000	1.0000
F81	14	-1945.8707	3979.6144	0.0000	1.0000
K2PG	13	-1961.9357	4005.4677	0.0000	1.0000
K2PI	13	-1962.3062	4006.2088	0.0000	1.0000
K2PIG	14	-1961.8832	4011.6394	0.0000	1.0000
SYMG	17	-1955.8265	4018.3559	0.0000	1.0000
SYMI	17	-1956.5668	4019.8365	0.0000	1.0000
SYMIG	18	-1955.8097	4024.5989	0.0000	1.0000
K2P	12	-1979.7067	4034.7331	0.0000	1.0000
JC69G	12	-1983.3780	4042.0757	0.0000	1.0000
JC69I	12	-1983.6218	4042.5632	0.0000	1.0000
JC69IG	13	-1983.2789	4048.1541	0.0000	1.0000
SYM	16	-1974.2639	4048.9540	0.0000	1.0000
JC69	11	-1999.7267	4068.4965	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/chr2_8593.phylip = 1 - 532;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8593.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/chr2_8593.phylip = 1 - 532;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8593.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/chr2_8593.phylip = 1 - 532;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8593.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

