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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2415.9755	4864.7957	0.5592	0.5592
HKYIG	17	-2415.8525	4866.6568	0.2205	0.7798
GTRG	20	-2413.3610	4868.0346	0.1107	0.8905
HKYI	16	-2418.2906	4869.4260	0.0552	0.9457
GTRIG	21	-2413.2447	4869.9355	0.0428	0.9885
GTRI	20	-2415.6280	4872.5685	0.0115	1.0000
HKY	15	-2449.7313	4930.2068	0.0000	1.0000
GTR	19	-2447.2459	4933.6775	0.0000	1.0000
SYMG	17	-2460.4854	4955.9226	0.0000	1.0000
K2PG	13	-2465.6811	4957.9248	0.0000	1.0000
SYMIG	18	-2460.4402	4957.9459	0.0000	1.0000
K2PIG	14	-2465.6301	4959.9105	0.0000	1.0000
SYMI	17	-2463.7185	4962.3887	0.0000	1.0000
K2PI	13	-2468.3794	4963.3214	0.0000	1.0000
F81G	15	-2467.8522	4966.4485	0.0000	1.0000
F81IG	16	-2467.8197	4968.4842	0.0000	1.0000
F81I	15	-2469.6031	4969.9504	0.0000	1.0000
K2P	12	-2498.1426	5020.7667	0.0000	1.0000
SYM	16	-2494.0766	5020.9980	0.0000	1.0000
F81	14	-2497.6294	5023.9090	0.0000	1.0000
JC69G	12	-2507.7341	5039.9497	0.0000	1.0000
JC69IG	13	-2507.7202	5042.0030	0.0000	1.0000
JC69I	12	-2509.7460	5043.9734	0.0000	1.0000
JC69	11	-2537.3809	5097.1686	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2415.9755	4863.9510	0.5279	0.5279
HKYIG	17	-2415.8525	4865.7050	0.2196	0.7475
GTRG	20	-2413.3610	4866.7221	0.1321	0.8796
GTRIG	21	-2413.2447	4868.4895	0.0546	0.9342
HKYI	16	-2418.2906	4868.5813	0.0521	0.9863
GTRI	20	-2415.6280	4871.2560	0.0137	1.0000
HKY	15	-2449.7313	4929.4626	0.0000	1.0000
GTR	19	-2447.2459	4932.4919	0.0000	1.0000
SYMG	17	-2460.4854	4954.9708	0.0000	1.0000
SYMIG	18	-2460.4402	4956.8804	0.0000	1.0000
K2PG	13	-2465.6811	4957.3622	0.0000	1.0000
K2PIG	14	-2465.6301	4959.2603	0.0000	1.0000
SYMI	17	-2463.7185	4961.4369	0.0000	1.0000
K2PI	13	-2468.3794	4962.7588	0.0000	1.0000
F81G	15	-2467.8522	4965.7044	0.0000	1.0000
F81IG	16	-2467.8197	4967.6395	0.0000	1.0000
F81I	15	-2469.6031	4969.2062	0.0000	1.0000
SYM	16	-2494.0766	5020.1532	0.0000	1.0000
K2P	12	-2498.1426	5020.2852	0.0000	1.0000
F81	14	-2497.6294	5023.2589	0.0000	1.0000
JC69G	12	-2507.7341	5039.4683	0.0000	1.0000
JC69IG	13	-2507.7202	5041.4404	0.0000	1.0000
JC69I	12	-2509.7460	5043.4920	0.0000	1.0000
JC69	11	-2537.3809	5096.7618	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2415.9755	4935.8511	0.8751	0.8751
HKYI	16	-2418.2906	4940.4814	0.0864	0.9615
HKYIG	17	-2415.8525	4942.0989	0.0385	1.0000
GTRG	20	-2413.3610	4956.5972	0.0000	1.0000
GTRI	20	-2415.6280	4961.1311	0.0000	1.0000
GTRIG	21	-2413.2447	4962.8583	0.0000	1.0000
HKY	15	-2449.7313	4996.8689	0.0000	1.0000
K2PG	13	-2465.6811	5015.7810	0.0000	1.0000
GTR	19	-2447.2459	5017.8732	0.0000	1.0000
K2PI	13	-2468.3794	5021.1776	0.0000	1.0000
K2PIG	14	-2465.6301	5022.1729	0.0000	1.0000
SYMG	17	-2460.4854	5031.3646	0.0000	1.0000
F81G	15	-2467.8522	5033.1107	0.0000	1.0000
F81I	15	-2469.6031	5036.6125	0.0000	1.0000
SYMIG	18	-2460.4402	5037.7680	0.0000	1.0000
SYMI	17	-2463.7185	5037.8307	0.0000	1.0000
F81IG	16	-2467.8197	5039.5395	0.0000	1.0000
K2P	12	-2498.1426	5074.2102	0.0000	1.0000
F81	14	-2497.6294	5086.1714	0.0000	1.0000
SYM	16	-2494.0766	5092.0533	0.0000	1.0000
JC69G	12	-2507.7341	5093.3933	0.0000	1.0000
JC69I	12	-2509.7460	5097.4170	0.0000	1.0000
JC69IG	13	-2507.7202	5099.8592	0.0000	1.0000
JC69	11	-2537.3809	5146.1931	0.0000	1.0000



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

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

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


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

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


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

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

