
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:43 02/17/2011

Input data from file "183-phylip/chr2_7957.phylip" (29 taxa, 670 characters)

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


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	64	-2264.4985	4670.7491	0.5279	0.5279
GTRIG	65	-2264.3994	4673.0040	0.1710	0.6988
SYMG	61	-2269.5344	4673.5095	0.1328	0.8316
HKYG	60	-2271.4486	4674.9169	0.0657	0.8973
K2PG	57	-2275.2031	4675.2102	0.0567	0.9540
SYMIG	62	-2269.9252	4676.7202	0.0267	0.9807
HKYIG	61	-2272.0503	4678.5413	0.0107	0.9914
K2PIG	58	-2275.9300	4679.0612	0.0083	0.9997
GTRI	64	-2272.2707	4686.2934	0.0002	0.9999
SYMI	61	-2277.1054	4688.6515	0.0001	1.0000
HKYI	60	-2279.6967	4691.4132	0.0000	1.0000
K2PI	57	-2283.4104	4691.6248	0.0000	1.0000
JC69G	56	-2318.0255	4758.4653	0.0000	1.0000
JC69IG	57	-2317.7982	4760.4002	0.0000	1.0000
F81G	59	-2315.4622	4760.5310	0.0000	1.0000
F81IG	60	-2315.1030	4762.2258	0.0000	1.0000
JC69I	56	-2324.4548	4771.3239	0.0000	1.0000
F81I	59	-2322.0170	4773.6405	0.0000	1.0000
K2P	56	-2325.8106	4774.0356	0.0000	1.0000
HKY	59	-2322.7369	4775.0803	0.0000	1.0000
SYM	60	-2322.5723	4777.1643	0.0000	1.0000
GTR	63	-2319.3156	4777.9382	0.0000	1.0000
JC69	55	-2366.6215	4853.2755	0.0000	1.0000
F81	58	-2364.2513	4855.7038	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2264.4985	4656.9970	0.6018	0.6018
GTRIG	65	-2264.3994	4658.7987	0.2445	0.8463
SYMG	61	-2269.5344	4661.0687	0.0786	0.9249
HKYG	60	-2271.4486	4662.8972	0.0315	0.9563
SYMIG	62	-2269.9252	4663.8503	0.0196	0.9759
K2PG	57	-2275.2031	4664.4063	0.0148	0.9907
HKYIG	61	-2272.0503	4666.1006	0.0063	0.9971
K2PIG	58	-2275.9300	4667.8599	0.0026	0.9997
GTRI	64	-2272.2707	4672.5414	0.0003	0.9999
SYMI	61	-2277.1054	4676.2107	0.0000	1.0000
HKYI	60	-2279.6967	4679.3935	0.0000	1.0000
K2PI	57	-2283.4104	4680.8209	0.0000	1.0000
JC69G	56	-2318.0255	4748.0509	0.0000	1.0000
F81G	59	-2315.4622	4748.9244	0.0000	1.0000
JC69IG	57	-2317.7982	4749.5963	0.0000	1.0000
F81IG	60	-2315.1030	4750.2061	0.0000	1.0000
JC69I	56	-2324.4548	4760.9095	0.0000	1.0000
F81I	59	-2322.0170	4762.0340	0.0000	1.0000
HKY	59	-2322.7369	4763.4738	0.0000	1.0000
K2P	56	-2325.8106	4763.6213	0.0000	1.0000
GTR	63	-2319.3156	4764.6312	0.0000	1.0000
SYM	60	-2322.5723	4765.1446	0.0000	1.0000
JC69	55	-2366.6215	4843.2430	0.0000	1.0000
F81	58	-2364.2513	4844.5025	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2275.2031	4921.3211	0.9783	0.9783
K2PIG	58	-2275.9300	4929.2820	0.0183	0.9966
HKYG	60	-2271.4486	4933.3339	0.0024	0.9990
SYMG	61	-2269.5344	4936.0127	0.0006	0.9997
K2PI	57	-2283.4104	4937.7357	0.0003	0.9999
HKYIG	61	-2272.0503	4941.0445	0.0001	1.0000
SYMIG	62	-2269.9252	4943.3015	0.0000	1.0000
GTRG	64	-2264.4985	4945.4628	0.0000	1.0000
HKYI	60	-2279.6967	4949.8301	0.0000	1.0000
SYMI	61	-2277.1054	4951.1547	0.0000	1.0000
GTRIG	65	-2264.3994	4951.7718	0.0000	1.0000
GTRI	64	-2272.2707	4961.0072	0.0000	1.0000
JC69G	56	-2318.0255	5000.4585	0.0000	1.0000
JC69IG	57	-2317.7982	5006.5111	0.0000	1.0000
JC69I	56	-2324.4548	5013.3171	0.0000	1.0000
F81G	59	-2315.4622	5014.8538	0.0000	1.0000
K2P	56	-2325.8106	5016.0288	0.0000	1.0000
F81IG	60	-2315.1030	5020.6427	0.0000	1.0000
F81I	59	-2322.0170	5027.9634	0.0000	1.0000
HKY	59	-2322.7369	5029.4031	0.0000	1.0000
SYM	60	-2322.5723	5035.5812	0.0000	1.0000
GTR	63	-2319.3156	5048.5897	0.0000	1.0000
JC69	55	-2366.6215	5091.1432	0.0000	1.0000
F81	58	-2364.2513	5105.9246	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 183-phylip/chr2_7957.phylip = 1 - 670;
 Partition Dummy = 1:183-phylip/chr2_7957.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 183-phylip/chr2_7957.phylip = 1 - 670;
 Partition Dummy = 1:183-phylip/chr2_7957.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_7957.phylip = 1 - 670;
 Partition Dummy = 1:183-phylip/chr2_7957.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

