
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:14 04/04/2011

Input data from file "683-tomodel-phylip/chr3_5533.phylip" (7 taxa, 157 characters)

Minimum AIC  model: HKY
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	12	-329.6128	685.3922	0.3077	0.3077
HKY	15	-326.4041	686.2124	0.2042	0.5120
K2PI	13	-329.5716	687.6888	0.0976	0.6096
K2PG	13	-329.5734	687.6922	0.0974	0.7070
HKYG	16	-326.2690	688.4238	0.0676	0.7746
HKYI	16	-326.2714	688.4284	0.0674	0.8420
SYM	16	-326.6695	689.2247	0.0453	0.8873
K2PIG	14	-329.5716	690.1010	0.0292	0.9165
GTR	19	-323.5794	690.7063	0.0216	0.9381
HKYIG	17	-326.2686	690.9402	0.0192	0.9573
SYMG	17	-326.6691	691.7411	0.0129	0.9702
SYMI	17	-326.6725	691.7478	0.0128	0.9830
GTRI	20	-323.5685	693.3136	0.0059	0.9889
GTRG	20	-323.5704	693.3172	0.0059	0.9947
SYMIG	18	-326.6709	694.2984	0.0036	0.9983
GTRIG	21	-323.5682	695.9808	0.0015	0.9999
JC69	11	-339.5376	702.8959	0.0000	0.9999
F81	14	-336.3892	703.7361	0.0000	0.9999
JC69G	12	-339.5348	705.2363	0.0000	1.0000
JC69I	12	-339.5364	705.2394	0.0000	1.0000
F81I	15	-336.3813	706.1669	0.0000	1.0000
F81G	15	-336.3819	706.1681	0.0000	1.0000
JC69IG	13	-339.5362	707.6179	0.0000	1.0000
F81IG	16	-336.3799	708.6454	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-326.4041	682.8082	0.2149	0.2149
K2P	12	-329.6128	683.2255	0.1745	0.3894
HKYG	16	-326.2690	684.5381	0.0905	0.4799
HKYI	16	-326.2714	684.5427	0.0903	0.5702
K2PI	13	-329.5716	685.1433	0.0669	0.6371
K2PG	13	-329.5734	685.1467	0.0668	0.7038
GTR	19	-323.5794	685.1588	0.0664	0.7702
SYM	16	-326.6695	685.3390	0.0606	0.8308
HKYIG	17	-326.2686	686.5373	0.0333	0.8642
GTRI	20	-323.5685	687.1371	0.0247	0.8888
GTRG	20	-323.5704	687.1407	0.0246	0.9135
K2PIG	14	-329.5716	687.1432	0.0246	0.9381
SYMG	17	-326.6691	687.3383	0.0223	0.9604
SYMI	17	-326.6725	687.3449	0.0222	0.9826
GTRIG	21	-323.5682	689.1363	0.0091	0.9917
SYMIG	18	-326.6709	689.3419	0.0082	0.9999
F81	14	-336.3892	700.7784	0.0000	0.9999
JC69	11	-339.5376	701.0752	0.0000	1.0000
F81I	15	-336.3813	702.7627	0.0000	1.0000
F81G	15	-336.3819	702.7639	0.0000	1.0000
JC69G	12	-339.5348	703.0696	0.0000	1.0000
JC69I	12	-339.5364	703.0727	0.0000	1.0000
F81IG	16	-336.3799	704.7597	0.0000	1.0000
JC69IG	13	-339.5362	705.0724	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-329.6128	719.9005	0.8408	0.8408
K2PI	13	-329.5716	724.8745	0.0699	0.9107
K2PG	13	-329.5734	724.8779	0.0698	0.9805
HKY	15	-326.4041	728.6518	0.0106	0.9910
K2PIG	14	-329.5716	729.9307	0.0056	0.9966
HKYG	16	-326.2690	733.4380	0.0010	0.9976
HKYI	16	-326.2714	733.4427	0.0010	0.9986
SYM	16	-326.6695	734.2389	0.0006	0.9992
JC69	11	-339.5376	734.6939	0.0005	0.9997
HKYIG	17	-326.2686	738.4935	0.0001	0.9998
SYMG	17	-326.6691	739.2944	0.0001	0.9998
SYMI	17	-326.6725	739.3011	0.0001	0.9999
JC69G	12	-339.5348	739.7446	0.0000	0.9999
JC69I	12	-339.5364	739.7477	0.0000	1.0000
GTR	19	-323.5794	743.2275	0.0000	1.0000
F81	14	-336.3892	743.5658	0.0000	1.0000
SYMIG	18	-326.6709	744.3543	0.0000	1.0000
JC69IG	13	-339.5362	744.8036	0.0000	1.0000
GTRI	20	-323.5685	748.2620	0.0000	1.0000
GTRG	20	-323.5704	748.2657	0.0000	1.0000
F81I	15	-336.3813	748.6063	0.0000	1.0000
F81G	15	-336.3819	748.6076	0.0000	1.0000
GTRIG	21	-323.5682	753.3175	0.0000	1.0000
F81IG	16	-336.3799	753.6597	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_5533.phylip = 1 - 157;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5533.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_5533.phylip = 1 - 157;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5533.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_5533.phylip = 1 - 157;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5533.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

