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

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

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


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	-2405.1218	4851.3517	0.6144	0.6144
GTRIG	21	-2405.0539	4853.3285	0.2287	0.8431
GTRI	20	-2406.4868	4854.0818	0.1569	1.0000
SYMG	17	-2418.6771	4872.1584	0.0000	1.0000
SYMIG	18	-2418.5894	4874.0787	0.0000	1.0000
SYMI	17	-2420.1699	4875.1441	0.0000	1.0000
HKYG	16	-2426.5038	4885.7215	0.0000	1.0000
HKYIG	17	-2426.4061	4887.6164	0.0000	1.0000
HKYI	16	-2428.0054	4888.7248	0.0000	1.0000
K2PG	13	-2434.3078	4895.0914	0.0000	1.0000
K2PIG	14	-2434.2115	4896.9727	0.0000	1.0000
K2PI	13	-2435.8801	4898.2361	0.0000	1.0000
GTR	19	-2435.1363	4909.2740	0.0000	1.0000
SYM	16	-2445.9642	4924.6423	0.0000	1.0000
HKY	15	-2454.4763	4939.5817	0.0000	1.0000
F81G	15	-2455.3542	4941.3374	0.0000	1.0000
F81IG	16	-2455.2994	4943.3127	0.0000	1.0000
F81I	15	-2456.6120	4943.8532	0.0000	1.0000
JC69G	12	-2461.1183	4946.6439	0.0000	1.0000
K2P	12	-2461.3265	4947.0603	0.0000	1.0000
JC69IG	13	-2461.0603	4948.5964	0.0000	1.0000
JC69I	12	-2462.4157	4949.2387	0.0000	1.0000
F81	14	-2481.4356	4991.4210	0.0000	1.0000
JC69	11	-2486.9206	4996.1853	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2405.1218	4850.2435	0.6064	0.6064
GTRIG	21	-2405.0539	4852.1079	0.2387	0.8451
GTRI	20	-2406.4868	4852.9736	0.1549	1.0000
SYMG	17	-2418.6771	4871.3542	0.0000	1.0000
SYMIG	18	-2418.5894	4873.1787	0.0000	1.0000
SYMI	17	-2420.1699	4874.3399	0.0000	1.0000
HKYG	16	-2426.5038	4885.0076	0.0000	1.0000
HKYIG	17	-2426.4061	4886.8122	0.0000	1.0000
HKYI	16	-2428.0054	4888.0109	0.0000	1.0000
K2PG	13	-2434.3078	4894.6156	0.0000	1.0000
K2PIG	14	-2434.2115	4896.4229	0.0000	1.0000
K2PI	13	-2435.8801	4897.7602	0.0000	1.0000
GTR	19	-2435.1363	4908.2727	0.0000	1.0000
SYM	16	-2445.9642	4923.9283	0.0000	1.0000
HKY	15	-2454.4763	4938.9526	0.0000	1.0000
F81G	15	-2455.3542	4940.7083	0.0000	1.0000
F81IG	16	-2455.2994	4942.5988	0.0000	1.0000
F81I	15	-2456.6120	4943.2241	0.0000	1.0000
JC69G	12	-2461.1183	4946.2365	0.0000	1.0000
K2P	12	-2461.3265	4946.6529	0.0000	1.0000
JC69IG	13	-2461.0603	4948.1206	0.0000	1.0000
JC69I	12	-2462.4157	4948.8314	0.0000	1.0000
F81	14	-2481.4356	4990.8712	0.0000	1.0000
JC69	11	-2486.9206	4995.8411	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2405.1218	4943.4037	0.7501	0.7501
GTRI	20	-2406.4868	4946.1338	0.1915	0.9416
GTRIG	21	-2405.0539	4949.9261	0.0288	0.9704
SYMG	17	-2418.6771	4950.5404	0.0212	0.9915
SYMI	17	-2420.1699	4953.5260	0.0048	0.9963
K2PG	13	-2434.3078	4955.1697	0.0021	0.9984
SYMIG	18	-2418.5894	4957.0229	0.0008	0.9992
K2PI	13	-2435.8801	4958.3144	0.0004	0.9996
HKYG	16	-2426.5038	4959.5358	0.0002	0.9999
K2PIG	14	-2434.2115	4961.6351	0.0001	0.9999
HKYI	16	-2428.0054	4962.5390	0.0001	1.0000
HKYIG	17	-2426.4061	4965.9984	0.0000	1.0000
GTR	19	-2435.1363	4996.7749	0.0000	1.0000
SYM	16	-2445.9642	4998.4565	0.0000	1.0000
JC69G	12	-2461.1183	5002.1327	0.0000	1.0000
K2P	12	-2461.3265	5002.5491	0.0000	1.0000
JC69I	12	-2462.4157	5004.7275	0.0000	1.0000
JC69IG	13	-2461.0603	5008.6747	0.0000	1.0000
HKY	15	-2454.4763	5008.8228	0.0000	1.0000
F81G	15	-2455.3542	5010.5785	0.0000	1.0000
F81I	15	-2456.6120	5013.0942	0.0000	1.0000
F81IG	16	-2455.2994	5017.1270	0.0000	1.0000
JC69	11	-2486.9206	5047.0793	0.0000	1.0000
F81	14	-2481.4356	5056.0834	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/chr2_4311.phylip = 1 - 779;
 Partition Dummy = 1:591-tomodel-phylip/chr2_4311.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/chr2_4311.phylip = 1 - 779;
 Partition Dummy = 1:591-tomodel-phylip/chr2_4311.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_4311.phylip = 1 - 779;
 Partition Dummy = 1:591-tomodel-phylip/chr2_4311.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;


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

