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

Input data from file "591-tomodel-phylip/chr8_7440.phylip" (7 taxa, 759 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	-2416.5376	4874.2134	0.6166	0.6166
GTRIG	21	-2416.5475	4876.3488	0.2120	0.8286
GTRI	20	-2417.8180	4876.7743	0.1714	1.0000
HKYG	16	-2431.8639	4896.4609	0.0000	1.0000
HKYIG	17	-2431.8683	4898.5626	0.0000	1.0000
HKYI	16	-2433.3058	4899.3448	0.0000	1.0000
GTR	19	-2434.2548	4907.5380	0.0000	1.0000
SYMG	17	-2438.3245	4911.4749	0.0000	1.0000
SYMIG	18	-2438.3322	4913.5887	0.0000	1.0000
SYMI	17	-2440.0287	4914.8834	0.0000	1.0000
K2PG	13	-2450.0408	4926.5701	0.0000	1.0000
HKY	15	-2448.4878	4927.6216	0.0000	1.0000
K2PIG	14	-2450.0406	4928.6458	0.0000	1.0000
K2PI	13	-2451.8131	4930.1148	0.0000	1.0000
SYM	16	-2454.1303	4940.9937	0.0000	1.0000
K2P	12	-2465.6088	4955.6359	0.0000	1.0000
F81G	15	-2476.8845	4984.4149	0.0000	1.0000
F81IG	16	-2476.8844	4986.5020	0.0000	1.0000
F81I	15	-2478.8938	4988.4337	0.0000	1.0000
JC69G	12	-2491.0637	5006.5456	0.0000	1.0000
JC69IG	13	-2491.0637	5008.6159	0.0000	1.0000
F81	14	-2490.5844	5009.7332	0.0000	1.0000
JC69I	12	-2493.1806	5010.7795	0.0000	1.0000
JC69	11	-2504.7321	5031.8176	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2416.5376	4873.0752	0.6089	0.6089
GTRIG	21	-2416.5475	4875.0951	0.2218	0.8308
GTRI	20	-2417.8180	4875.6361	0.1692	1.0000
HKYG	16	-2431.8639	4895.7278	0.0000	1.0000
HKYIG	17	-2431.8683	4897.7367	0.0000	1.0000
HKYI	16	-2433.3058	4898.6116	0.0000	1.0000
GTR	19	-2434.2548	4906.5096	0.0000	1.0000
SYMG	17	-2438.3245	4910.6490	0.0000	1.0000
SYMIG	18	-2438.3322	4912.6644	0.0000	1.0000
SYMI	17	-2440.0287	4914.0575	0.0000	1.0000
K2PG	13	-2450.0408	4926.0815	0.0000	1.0000
HKY	15	-2448.4878	4926.9756	0.0000	1.0000
K2PIG	14	-2450.0406	4928.0813	0.0000	1.0000
K2PI	13	-2451.8131	4929.6262	0.0000	1.0000
SYM	16	-2454.1303	4940.2605	0.0000	1.0000
K2P	12	-2465.6088	4955.2176	0.0000	1.0000
F81G	15	-2476.8845	4983.7689	0.0000	1.0000
F81IG	16	-2476.8844	4985.7688	0.0000	1.0000
F81I	15	-2478.8938	4987.7877	0.0000	1.0000
JC69G	12	-2491.0637	5006.1273	0.0000	1.0000
JC69IG	13	-2491.0637	5008.1273	0.0000	1.0000
F81	14	-2490.5844	5009.1687	0.0000	1.0000
JC69I	12	-2493.1806	5010.3613	0.0000	1.0000
JC69	11	-2504.7321	5031.4642	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2416.5376	4965.7152	0.6776	0.6776
GTRI	20	-2417.8180	4968.2761	0.1883	0.8660
HKYG	16	-2431.8639	4969.8398	0.0862	0.9521
GTRIG	21	-2416.5475	4972.3671	0.0244	0.9765
HKYI	16	-2433.3058	4972.7236	0.0204	0.9969
HKYIG	17	-2431.8683	4976.4807	0.0031	1.0000
K2PG	13	-2450.0408	4986.2975	0.0000	1.0000
SYMG	17	-2438.3245	4989.3931	0.0000	1.0000
K2PI	13	-2451.8131	4989.8422	0.0000	1.0000
SYMI	17	-2440.0287	4992.8015	0.0000	1.0000
K2PIG	14	-2450.0406	4992.9293	0.0000	1.0000
GTR	19	-2434.2548	4994.5176	0.0000	1.0000
SYMIG	18	-2438.3322	4996.0404	0.0000	1.0000
HKY	15	-2448.4878	4996.4556	0.0000	1.0000
K2P	12	-2465.6088	5010.8016	0.0000	1.0000
SYM	16	-2454.1303	5014.3726	0.0000	1.0000
F81G	15	-2476.8845	5053.2489	0.0000	1.0000
F81I	15	-2478.8938	5057.2677	0.0000	1.0000
F81IG	16	-2476.8844	5059.8809	0.0000	1.0000
JC69G	12	-2491.0637	5061.7114	0.0000	1.0000
JC69I	12	-2493.1806	5065.9453	0.0000	1.0000
JC69IG	13	-2491.0637	5068.3434	0.0000	1.0000
F81	14	-2490.5844	5074.0167	0.0000	1.0000
JC69	11	-2504.7321	5082.4162	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/chr8_7440.phylip = 1 - 759;
 Partition Dummy = 1:591-tomodel-phylip/chr8_7440.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/chr8_7440.phylip = 1 - 759;
 Partition Dummy = 1:591-tomodel-phylip/chr8_7440.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/chr8_7440.phylip = 1 - 759;
 Partition Dummy = 1:591-tomodel-phylip/chr8_7440.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

