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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PI


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
GTRI	20	-2627.8787	5296.8570	0.5902	0.5902
GTRIG	21	-2627.8409	5298.8927	0.2133	0.8035
GTRG	20	-2629.0746	5299.2486	0.1785	0.9820
HKYI	16	-2636.3325	5305.3733	0.0084	0.9903
HKYG	16	-2636.9401	5306.5885	0.0045	0.9949
HKYIG	17	-2636.3139	5307.4257	0.0030	0.9979
SYMI	17	-2637.8159	5310.4297	0.0007	0.9985
K2PI	13	-2642.3617	5311.1956	0.0005	0.9990
SYMG	17	-2638.3779	5311.5537	0.0004	0.9994
SYMIG	18	-2637.7770	5312.4469	0.0002	0.9996
K2PG	13	-2643.0057	5312.4835	0.0002	0.9998
K2PIG	14	-2642.3459	5313.2372	0.0002	1.0000
GTR	19	-2663.9054	5366.8043	0.0000	1.0000
HKY	15	-2670.8737	5372.3716	0.0000	1.0000
SYM	16	-2671.9208	5376.5499	0.0000	1.0000
K2P	12	-2676.9180	5378.2401	0.0000	1.0000
F81I	15	-2681.7251	5394.0744	0.0000	1.0000
F81IG	16	-2681.7365	5396.1813	0.0000	1.0000
F81G	15	-2682.9918	5396.6077	0.0000	1.0000
JC69I	12	-2686.2339	5396.8720	0.0000	1.0000
JC69IG	13	-2686.2414	5398.9550	0.0000	1.0000
JC69G	12	-2687.4168	5399.2377	0.0000	1.0000
F81	14	-2715.9164	5460.3783	0.0000	1.0000
JC69	11	-2720.5050	5463.3516	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2627.8787	5295.7575	0.5849	0.5849
GTRIG	21	-2627.8409	5297.6817	0.2235	0.8083
GTRG	20	-2629.0746	5298.1491	0.1769	0.9852
HKYI	16	-2636.3325	5304.6650	0.0068	0.9920
HKYG	16	-2636.9401	5305.8801	0.0037	0.9957
HKYIG	17	-2636.3139	5306.6278	0.0026	0.9983
SYMI	17	-2637.8159	5309.6318	0.0006	0.9988
K2PI	13	-2642.3617	5310.7235	0.0003	0.9992
SYMG	17	-2638.3779	5310.7558	0.0003	0.9995
SYMIG	18	-2637.7770	5311.5539	0.0002	0.9997
K2PG	13	-2643.0057	5312.0114	0.0002	0.9999
K2PIG	14	-2642.3459	5312.6918	0.0001	1.0000
GTR	19	-2663.9054	5365.8109	0.0000	1.0000
HKY	15	-2670.8737	5371.7475	0.0000	1.0000
SYM	16	-2671.9208	5375.8415	0.0000	1.0000
K2P	12	-2676.9180	5377.8359	0.0000	1.0000
F81I	15	-2681.7251	5393.4503	0.0000	1.0000
F81IG	16	-2681.7365	5395.4730	0.0000	1.0000
F81G	15	-2682.9918	5395.9835	0.0000	1.0000
JC69I	12	-2686.2339	5396.4678	0.0000	1.0000
JC69IG	13	-2686.2414	5398.4829	0.0000	1.0000
JC69G	12	-2687.4168	5398.8336	0.0000	1.0000
F81	14	-2715.9164	5459.8329	0.0000	1.0000
JC69	11	-2720.5050	5463.0101	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2642.3617	5371.3774	0.6282	0.6282
K2PG	13	-2643.0057	5372.6653	0.3300	0.9582
K2PIG	14	-2642.3459	5378.0114	0.0228	0.9810
HKYI	16	-2636.3325	5379.3159	0.0119	0.9928
HKYG	16	-2636.9401	5380.5311	0.0065	0.9993
HKYIG	17	-2636.3139	5385.9444	0.0004	0.9997
SYMI	17	-2637.8159	5388.9484	0.0001	0.9998
GTRI	20	-2627.8787	5389.0712	0.0001	0.9999
SYMG	17	-2638.3779	5390.0724	0.0001	1.0000
GTRG	20	-2629.0746	5391.4628	0.0000	1.0000
SYMIG	18	-2637.7770	5395.5362	0.0000	1.0000
GTRIG	21	-2627.8409	5395.6611	0.0000	1.0000
K2P	12	-2676.9180	5433.8241	0.0000	1.0000
HKY	15	-2670.8737	5441.7327	0.0000	1.0000
SYM	16	-2671.9208	5450.4925	0.0000	1.0000
JC69I	12	-2686.2339	5452.4560	0.0000	1.0000
GTR	19	-2663.9054	5454.4589	0.0000	1.0000
JC69G	12	-2687.4168	5454.8218	0.0000	1.0000
JC69IG	13	-2686.2414	5459.1368	0.0000	1.0000
F81I	15	-2681.7251	5463.4355	0.0000	1.0000
F81G	15	-2682.9918	5465.9688	0.0000	1.0000
F81IG	16	-2681.7365	5470.1239	0.0000	1.0000
JC69	11	-2720.5050	5514.3326	0.0000	1.0000
F81	14	-2715.9164	5525.1525	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_25881.phylip = 1 - 785;
 Partition Dummy = 1:591-tomodel-phylip/chr2_25881.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_25881.phylip = 1 - 785;
 Partition Dummy = 1:591-tomodel-phylip/chr2_25881.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

