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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1505.3620	3043.4091	0.2491	0.2491
HKYG	16	-1505.3798	3043.4448	0.2447	0.4939
GTRI	20	-1501.5211	3044.1054	0.1759	0.6697
GTRG	20	-1501.5258	3044.1149	0.1751	0.8448
HKYIG	17	-1505.3615	3045.4948	0.0878	0.9326
GTRIG	21	-1501.5258	3046.2227	0.0610	0.9936
GTR	19	-1506.4353	3051.8315	0.0037	0.9973
HKY	15	-1510.9336	3052.4709	0.0027	1.0000
SYMI	17	-1518.2914	3071.3545	0.0000	1.0000
SYMG	17	-1518.2943	3071.3604	0.0000	1.0000
SYMIG	18	-1518.2904	3073.4444	0.0000	1.0000
K2PI	13	-1524.1875	3074.8317	0.0000	1.0000
K2PG	13	-1524.2523	3074.9613	0.0000	1.0000
K2PIG	14	-1524.1879	3076.9034	0.0000	1.0000
SYM	16	-1523.6355	3079.9561	0.0000	1.0000
K2P	12	-1530.1164	3084.6237	0.0000	1.0000
F81G	15	-1530.1435	3090.8908	0.0000	1.0000
F81I	15	-1530.1444	3090.8926	0.0000	1.0000
F81IG	16	-1530.1435	3092.9722	0.0000	1.0000
F81	14	-1535.2597	3099.0471	0.0000	1.0000
JC69G	12	-1548.2032	3120.7974	0.0000	1.0000
JC69I	12	-1548.2035	3120.7980	0.0000	1.0000
JC69IG	13	-1548.2029	3122.8625	0.0000	1.0000
JC69	11	-1553.3891	3129.1087	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1505.3620	3042.7240	0.2277	0.2277
HKYG	16	-1505.3798	3042.7597	0.2237	0.4514
GTRI	20	-1501.5211	3043.0421	0.1942	0.6456
GTRG	20	-1501.5258	3043.0516	0.1933	0.8389
HKYIG	17	-1505.3615	3044.7231	0.0838	0.9227
GTRIG	21	-1501.5258	3045.0516	0.0711	0.9938
GTR	19	-1506.4353	3050.8707	0.0039	0.9976
HKY	15	-1510.9336	3051.8671	0.0024	1.0000
SYMI	17	-1518.2914	3070.5827	0.0000	1.0000
SYMG	17	-1518.2943	3070.5886	0.0000	1.0000
SYMIG	18	-1518.2904	3072.5808	0.0000	1.0000
K2PI	13	-1524.1875	3074.3749	0.0000	1.0000
K2PG	13	-1524.2523	3074.5046	0.0000	1.0000
K2PIG	14	-1524.1879	3076.3758	0.0000	1.0000
SYM	16	-1523.6355	3079.2710	0.0000	1.0000
K2P	12	-1530.1164	3084.2327	0.0000	1.0000
F81G	15	-1530.1435	3090.2870	0.0000	1.0000
F81I	15	-1530.1444	3090.2888	0.0000	1.0000
F81IG	16	-1530.1435	3092.2870	0.0000	1.0000
F81	14	-1535.2597	3098.5195	0.0000	1.0000
JC69G	12	-1548.2032	3120.4065	0.0000	1.0000
JC69I	12	-1548.2035	3120.4070	0.0000	1.0000
JC69IG	13	-1548.2029	3122.4057	0.0000	1.0000
JC69	11	-1553.3891	3128.7783	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1505.3620	3117.8963	0.4703	0.4703
HKYG	16	-1505.3798	3117.9320	0.4620	0.9323
HKY	15	-1510.9336	3122.3411	0.0510	0.9833
HKYIG	17	-1505.3615	3124.5936	0.0165	0.9998
K2PI	13	-1524.1875	3135.4524	0.0001	0.9998
K2PG	13	-1524.2523	3135.5821	0.0001	0.9999
GTRI	20	-1501.5211	3137.0075	0.0000	1.0000
GTRG	20	-1501.5258	3137.0169	0.0000	1.0000
GTR	19	-1506.4353	3140.1378	0.0000	1.0000
K2P	12	-1530.1164	3140.6120	0.0000	1.0000
K2PIG	14	-1524.1879	3142.1516	0.0000	1.0000
GTRIG	21	-1501.5258	3143.7152	0.0000	1.0000
SYMI	17	-1518.2914	3150.4533	0.0000	1.0000
SYMG	17	-1518.2943	3150.4592	0.0000	1.0000
SYM	16	-1523.6355	3154.4433	0.0000	1.0000
SYMIG	18	-1518.2904	3157.1496	0.0000	1.0000
F81G	15	-1530.1435	3160.7610	0.0000	1.0000
F81I	15	-1530.1444	3160.7629	0.0000	1.0000
F81	14	-1535.2597	3164.2952	0.0000	1.0000
F81IG	16	-1530.1435	3167.4593	0.0000	1.0000
JC69G	12	-1548.2032	3176.7857	0.0000	1.0000
JC69I	12	-1548.2035	3176.7862	0.0000	1.0000
JC69	11	-1553.3891	3180.4592	0.0000	1.0000
JC69IG	13	-1548.2029	3183.4832	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24815.phylip = 1 - 811;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24815.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24815.phylip = 1 - 811;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24815.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24815.phylip = 1 - 811;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24815.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

