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

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

Minimum AIC  model: HKYI
Minimum AICc model: K2PI
Minimum BIC  model: JC69


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
K2PI	13	-303.2381	635.2130	0.1465	0.1465
K2PG	13	-303.2382	635.2131	0.1465	0.2930
K2P	12	-305.0321	636.3925	0.0812	0.3742
HKYI	16	-300.2212	636.6270	0.0722	0.4465
HKYG	16	-300.2224	636.6293	0.0722	0.5186
JC69	11	-306.4463	636.8482	0.0647	0.5833
JC69I	12	-305.4685	637.2654	0.0525	0.6358
JC69G	12	-305.4903	637.3089	0.0514	0.6872
K2PIG	14	-303.2381	637.6580	0.0431	0.7304
HKY	15	-302.0800	637.8240	0.0397	0.7701
F81	14	-303.4577	638.0972	0.0346	0.8047
F81I	15	-302.4809	638.6260	0.0266	0.8313
F81G	15	-302.5045	638.6731	0.0260	0.8573
SYM	16	-301.3720	638.9286	0.0229	0.8801
HKYIG	17	-300.2211	639.1864	0.0201	0.9002
SYMI	17	-300.2808	639.3057	0.0189	0.9191
SYMG	17	-300.3250	639.3942	0.0181	0.9373
JC69IG	13	-305.5034	639.7436	0.0152	0.9525
GTR	19	-298.0990	640.1822	0.0122	0.9647
GTRI	20	-296.9117	640.4900	0.0105	0.9751
GTRG	20	-296.9567	640.5801	0.0100	0.9852
F81IG	16	-302.5074	641.1994	0.0073	0.9925
SYMIG	18	-300.3246	641.9930	0.0049	0.9974
GTRIG	21	-296.9565	643.3050	0.0026	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-300.2212	632.4423	0.0983	0.0983
HKYG	16	-300.2224	632.4447	0.0982	0.1966
K2PI	13	-303.2381	632.4761	0.0967	0.2932
K2PG	13	-303.2382	632.4763	0.0967	0.3899
GTRI	20	-296.9117	633.8233	0.0493	0.4392
GTRG	20	-296.9567	633.9135	0.0471	0.4863
K2P	12	-305.0321	634.0642	0.0437	0.5301
HKY	15	-302.0800	634.1599	0.0417	0.5717
GTR	19	-298.0990	634.1979	0.0409	0.6126
HKYIG	17	-300.2211	634.4422	0.0362	0.6488
K2PIG	14	-303.2381	634.4761	0.0356	0.6843
SYMI	17	-300.2808	634.5615	0.0341	0.7184
SYMG	17	-300.3250	634.6500	0.0326	0.7510
SYM	16	-301.3720	634.7439	0.0311	0.7821
JC69	11	-306.4463	634.8927	0.0289	0.8110
F81	14	-303.4577	634.9154	0.0286	0.8396
JC69I	12	-305.4685	634.9370	0.0282	0.8678
F81I	15	-302.4809	634.9619	0.0279	0.8957
JC69G	12	-305.4903	634.9805	0.0276	0.9234
F81G	15	-302.5045	635.0090	0.0273	0.9506
GTRIG	21	-296.9565	635.9130	0.0173	0.9680
SYMIG	18	-300.3246	636.6493	0.0120	0.9800
JC69IG	13	-305.5034	637.0067	0.0100	0.9900
F81IG	16	-302.5074	637.0147	0.0100	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
JC69	11	-306.4463	667.7874	0.4617	0.4617
K2P	12	-305.0321	669.9494	0.1567	0.6184
JC69I	12	-305.4685	670.8222	0.1013	0.7196
JC69G	12	-305.4903	670.8657	0.0991	0.8187
K2PI	13	-303.2381	671.3518	0.0777	0.8964
K2PG	13	-303.2382	671.3519	0.0777	0.9741
JC69IG	13	-305.5034	675.8823	0.0081	0.9822
K2PIG	14	-303.2381	676.3422	0.0064	0.9886
F81	14	-303.4577	676.7814	0.0051	0.9937
HKY	15	-302.0800	679.0164	0.0017	0.9954
F81I	15	-302.4809	679.8183	0.0011	0.9965
F81G	15	-302.5045	679.8655	0.0011	0.9976
HKYI	16	-300.2212	680.2893	0.0009	0.9985
HKYG	16	-300.2224	680.2916	0.0009	0.9994
SYM	16	-301.3720	682.5909	0.0003	0.9997
F81IG	16	-302.5074	684.8617	0.0001	0.9998
HKYIG	17	-300.2211	685.2796	0.0001	0.9999
SYMI	17	-300.2808	685.3989	0.0001	0.9999
SYMG	17	-300.3250	685.4873	0.0001	1.0000
SYMIG	18	-300.3246	690.4770	0.0000	1.0000
GTR	19	-298.0990	691.0161	0.0000	1.0000
GTRI	20	-296.9117	693.6320	0.0000	1.0000
GTRG	20	-296.9567	693.7221	0.0000	1.0000
GTRIG	21	-296.9565	698.7121	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/chr8_1548.phylip = 1 - 147;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1548.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 (K2PI)]

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


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_1548.phylip = 1 - 147;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1548.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal);
END;


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

