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

Input data from file "591-tomodel-phylip/chr6_3170.phylip" (7 taxa, 671 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	-1986.2492	4013.7907	0.6277	0.6277
GTRIG	21	-1985.9134	4015.2504	0.3025	0.9303
GTRI	20	-1988.4526	4018.1975	0.0693	0.9996
HKYG	16	-1998.3993	4029.6304	0.0002	0.9998
HKYIG	17	-1997.8338	4030.6048	0.0001	1.0000
HKYI	16	-2000.4956	4033.8230	0.0000	1.0000
SYMG	17	-2001.0047	4036.9467	0.0000	1.0000
SYMIG	18	-2000.4821	4038.0133	0.0000	1.0000
SYMI	17	-2003.3377	4041.6126	0.0000	1.0000
K2PG	13	-2011.6909	4049.9359	0.0000	1.0000
K2PIG	14	-2011.0225	4050.6853	0.0000	1.0000
K2PI	13	-2013.8048	4054.1637	0.0000	1.0000
GTR	19	-2026.2326	4091.6327	0.0000	1.0000
HKY	15	-2040.4454	4111.6237	0.0000	1.0000
SYM	16	-2041.7010	4116.2338	0.0000	1.0000
F81G	15	-2044.5884	4119.9097	0.0000	1.0000
F81IG	16	-2044.2875	4121.4068	0.0000	1.0000
F81I	15	-2046.1178	4122.9684	0.0000	1.0000
K2P	12	-2054.8144	4134.1030	0.0000	1.0000
JC69G	12	-2059.6432	4143.7606	0.0000	1.0000
JC69IG	13	-2059.2849	4145.1239	0.0000	1.0000
JC69I	12	-2061.3341	4147.1424	0.0000	1.0000
F81	14	-2084.5534	4197.7471	0.0000	1.0000
JC69	11	-2098.8798	4220.1603	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1986.2492	4012.4984	0.6151	0.6151
GTRIG	21	-1985.9134	4013.8267	0.3166	0.9318
GTRI	20	-1988.4526	4016.9052	0.0679	0.9997
HKYG	16	-1998.3993	4028.7986	0.0002	0.9999
HKYIG	17	-1997.8338	4029.6675	0.0001	1.0000
HKYI	16	-2000.4956	4032.9912	0.0000	1.0000
SYMG	17	-2001.0047	4036.0094	0.0000	1.0000
SYMIG	18	-2000.4821	4036.9643	0.0000	1.0000
SYMI	17	-2003.3377	4040.6754	0.0000	1.0000
K2PG	13	-2011.6909	4049.3818	0.0000	1.0000
K2PIG	14	-2011.0225	4050.0451	0.0000	1.0000
K2PI	13	-2013.8048	4053.6097	0.0000	1.0000
GTR	19	-2026.2326	4090.4653	0.0000	1.0000
HKY	15	-2040.4454	4110.8909	0.0000	1.0000
SYM	16	-2041.7010	4115.4020	0.0000	1.0000
F81G	15	-2044.5884	4119.1768	0.0000	1.0000
F81IG	16	-2044.2875	4120.5750	0.0000	1.0000
F81I	15	-2046.1178	4122.2356	0.0000	1.0000
K2P	12	-2054.8144	4133.6289	0.0000	1.0000
JC69G	12	-2059.6432	4143.2864	0.0000	1.0000
JC69IG	13	-2059.2849	4144.5698	0.0000	1.0000
JC69I	12	-2061.3341	4146.6682	0.0000	1.0000
F81	14	-2084.5534	4197.1068	0.0000	1.0000
JC69	11	-2098.8798	4219.7597	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1998.3993	4100.9389	0.5819	0.5819
GTRG	20	-1986.2492	4102.6738	0.2444	0.8264
HKYI	16	-2000.4956	4105.1315	0.0715	0.8979
HKYIG	17	-1997.8338	4106.3166	0.0395	0.9374
GTRI	20	-1988.4526	4107.0806	0.0270	0.9644
K2PG	13	-2011.6909	4107.9958	0.0171	0.9815
GTRIG	21	-1985.9134	4108.5109	0.0132	0.9947
K2PI	13	-2013.8048	4112.2237	0.0021	0.9968
SYMG	17	-2001.0047	4112.6585	0.0017	0.9984
K2PIG	14	-2011.0225	4113.1678	0.0013	0.9997
SYMI	17	-2003.3377	4117.3244	0.0002	0.9999
SYMIG	18	-2000.4821	4118.1221	0.0001	1.0000
GTR	19	-2026.2326	4176.1319	0.0000	1.0000
HKY	15	-2040.4454	4178.5224	0.0000	1.0000
F81G	15	-2044.5884	4186.8084	0.0000	1.0000
SYM	16	-2041.7010	4187.5423	0.0000	1.0000
K2P	12	-2054.8144	4187.7341	0.0000	1.0000
F81I	15	-2046.1178	4189.8671	0.0000	1.0000
F81IG	16	-2044.2875	4192.7153	0.0000	1.0000
JC69G	12	-2059.6432	4197.3917	0.0000	1.0000
JC69I	12	-2061.3341	4200.7734	0.0000	1.0000
JC69IG	13	-2059.2849	4203.1838	0.0000	1.0000
F81	14	-2084.5534	4260.2296	0.0000	1.0000
JC69	11	-2098.8798	4269.3561	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/chr6_3170.phylip = 1 - 671;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3170.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/chr6_3170.phylip = 1 - 671;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3170.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/chr6_3170.phylip = 1 - 671;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3170.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

