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

Input data from file "591-tomodel-phylip/chr9_3633.phylip" (7 taxa, 788 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	-2177.3044	4395.7040	0.4316	0.4316
GTRI	20	-2177.6852	4396.4657	0.2949	0.7264
GTRIG	21	-2177.3132	4397.8327	0.1489	0.8753
HKYG	16	-2183.4621	4399.6298	0.0606	0.9359
HKYI	16	-2183.8114	4400.3284	0.0427	0.9786
HKYIG	17	-2183.4700	4401.7347	0.0212	0.9998
GTR	19	-2186.2756	4411.5408	0.0002	1.0000
HKY	15	-2191.9590	4414.5398	0.0000	1.0000
K2PG	13	-2211.8277	4450.1258	0.0000	1.0000
K2PI	13	-2212.2401	4450.9504	0.0000	1.0000
K2PIG	14	-2211.8341	4452.2116	0.0000	1.0000
SYMG	17	-2209.0649	4452.9245	0.0000	1.0000
SYMI	17	-2209.4903	4453.7754	0.0000	1.0000
SYMIG	18	-2209.0679	4455.0252	0.0000	1.0000
K2P	12	-2220.1694	4464.7413	0.0000	1.0000
SYM	16	-2217.3898	4467.4852	0.0000	1.0000
F81G	15	-2237.4428	4505.5073	0.0000	1.0000
F81I	15	-2237.7741	4506.1699	0.0000	1.0000
F81IG	16	-2237.4519	4507.6094	0.0000	1.0000
F81	14	-2244.9746	4518.4926	0.0000	1.0000
JC69G	12	-2260.3173	4545.0372	0.0000	1.0000
JC69I	12	-2260.6806	4545.7638	0.0000	1.0000
JC69IG	13	-2260.3261	4547.1225	0.0000	1.0000
JC69	11	-2267.8992	4558.1386	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2177.3044	4394.6088	0.4371	0.4371
GTRI	20	-2177.6852	4395.3705	0.2987	0.7358
GTRIG	21	-2177.3132	4396.6264	0.1594	0.8952
HKYG	16	-2183.4621	4398.9242	0.0505	0.9457
HKYI	16	-2183.8114	4399.6228	0.0356	0.9814
HKYIG	17	-2183.4700	4400.9399	0.0184	0.9998
GTR	19	-2186.2756	4410.5512	0.0002	1.0000
HKY	15	-2191.9590	4413.9180	0.0000	1.0000
K2PG	13	-2211.8277	4449.6555	0.0000	1.0000
K2PI	13	-2212.2401	4450.4801	0.0000	1.0000
K2PIG	14	-2211.8341	4451.6683	0.0000	1.0000
SYMG	17	-2209.0649	4452.1297	0.0000	1.0000
SYMI	17	-2209.4903	4452.9806	0.0000	1.0000
SYMIG	18	-2209.0679	4454.1357	0.0000	1.0000
K2P	12	-2220.1694	4464.3387	0.0000	1.0000
SYM	16	-2217.3898	4466.7796	0.0000	1.0000
F81G	15	-2237.4428	4504.8856	0.0000	1.0000
F81I	15	-2237.7741	4505.5482	0.0000	1.0000
F81IG	16	-2237.4519	4506.9038	0.0000	1.0000
F81	14	-2244.9746	4517.9493	0.0000	1.0000
JC69G	12	-2260.3173	4544.6346	0.0000	1.0000
JC69I	12	-2260.6806	4545.3612	0.0000	1.0000
JC69IG	13	-2260.3261	4546.6522	0.0000	1.0000
JC69	11	-2267.8992	4557.7984	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2183.4621	4473.6361	0.5722	0.5722
HKYI	16	-2183.8114	4474.3348	0.4035	0.9757
HKYIG	17	-2183.4700	4480.3214	0.0202	0.9960
HKY	15	-2191.9590	4483.9605	0.0033	0.9993
GTRG	20	-2177.3044	4487.9987	0.0004	0.9997
GTRI	20	-2177.6852	4488.7605	0.0003	1.0000
GTRIG	21	-2177.3132	4494.6859	0.0000	1.0000
GTR	19	-2186.2756	4499.2717	0.0000	1.0000
K2PG	13	-2211.8277	4510.3590	0.0000	1.0000
K2PI	13	-2212.2401	4511.1836	0.0000	1.0000
K2PIG	14	-2211.8341	4517.0412	0.0000	1.0000
K2P	12	-2220.1694	4520.3727	0.0000	1.0000
SYMG	17	-2209.0649	4531.5112	0.0000	1.0000
SYMI	17	-2209.4903	4532.3621	0.0000	1.0000
SYMIG	18	-2209.0679	4538.1867	0.0000	1.0000
SYM	16	-2217.3898	4541.4915	0.0000	1.0000
F81G	15	-2237.4428	4574.9281	0.0000	1.0000
F81I	15	-2237.7741	4575.5907	0.0000	1.0000
F81IG	16	-2237.4519	4581.6158	0.0000	1.0000
F81	14	-2244.9746	4583.3223	0.0000	1.0000
JC69G	12	-2260.3173	4600.6686	0.0000	1.0000
JC69I	12	-2260.6806	4601.3952	0.0000	1.0000
JC69IG	13	-2260.3261	4607.3557	0.0000	1.0000
JC69	11	-2267.8992	4609.1629	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/chr9_3633.phylip = 1 - 788;
 Partition Dummy = 1:591-tomodel-phylip/chr9_3633.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/chr9_3633.phylip = 1 - 788;
 Partition Dummy = 1:591-tomodel-phylip/chr9_3633.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/chr9_3633.phylip = 1 - 788;
 Partition Dummy = 1:591-tomodel-phylip/chr9_3633.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

