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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-1413.4924	2859.9084	0.3490	0.3490
HKYG	16	-1413.6923	2860.3083	0.2857	0.6347
HKYIG	17	-1412.7132	2860.4672	0.2639	0.8986
K2PI	13	-1419.2429	2865.1006	0.0260	0.9246
K2PG	13	-1419.4830	2865.5808	0.0205	0.9450
K2PIG	14	-1418.5772	2865.8651	0.0178	0.9628
GTRI	20	-1412.4664	2866.3686	0.0138	0.9766
GTRG	20	-1412.6646	2866.7652	0.0113	0.9879
GTRIG	21	-1411.7334	2867.0491	0.0098	0.9977
SYMI	17	-1418.3580	2871.7568	0.0009	0.9987
SYMG	17	-1418.6027	2872.2462	0.0007	0.9994
SYMIG	18	-1417.7443	2872.6538	0.0006	1.0000
F81I	15	-1437.4196	2905.6528	0.0000	1.0000
F81G	15	-1437.5464	2905.9063	0.0000	1.0000
F81IG	16	-1436.9851	2906.8939	0.0000	1.0000
HKY	15	-1439.4518	2909.7171	0.0000	1.0000
JC69I	12	-1443.0633	2910.6527	0.0000	1.0000
JC69G	12	-1443.1968	2910.9197	0.0000	1.0000
JC69IG	13	-1442.6789	2911.9726	0.0000	1.0000
K2P	12	-1445.4259	2915.3779	0.0000	1.0000
GTR	19	-1438.4356	2916.1682	0.0000	1.0000
SYM	16	-1444.5603	2922.0442	0.0000	1.0000
F81	14	-1462.5663	2953.8432	0.0000	1.0000
JC69	11	-1468.2566	2958.9576	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1413.4924	2858.9848	0.3425	0.3425
HKYG	16	-1413.6923	2859.3847	0.2804	0.6229
HKYIG	17	-1412.7132	2859.4264	0.2746	0.8976
K2PI	13	-1419.2429	2864.4858	0.0219	0.9195
GTRI	20	-1412.4664	2864.9327	0.0175	0.9370
K2PG	13	-1419.4830	2864.9660	0.0172	0.9542
K2PIG	14	-1418.5772	2865.1544	0.0157	0.9699
GTRG	20	-1412.6646	2865.3293	0.0144	0.9842
GTRIG	21	-1411.7334	2865.4669	0.0134	0.9976
SYMI	17	-1418.3580	2870.7159	0.0010	0.9986
SYMG	17	-1418.6027	2871.2054	0.0008	0.9993
SYMIG	18	-1417.7443	2871.4886	0.0007	1.0000
F81I	15	-1437.4196	2904.8392	0.0000	1.0000
F81G	15	-1437.5464	2905.0927	0.0000	1.0000
F81IG	16	-1436.9851	2905.9703	0.0000	1.0000
HKY	15	-1439.4518	2908.9036	0.0000	1.0000
JC69I	12	-1443.0633	2910.1266	0.0000	1.0000
JC69G	12	-1443.1968	2910.3936	0.0000	1.0000
JC69IG	13	-1442.6789	2911.3578	0.0000	1.0000
K2P	12	-1445.4259	2914.8518	0.0000	1.0000
GTR	19	-1438.4356	2914.8713	0.0000	1.0000
SYM	16	-1444.5603	2921.1206	0.0000	1.0000
F81	14	-1462.5663	2953.1325	0.0000	1.0000
JC69	11	-1468.2566	2958.5132	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1419.2429	2921.7752	0.5247	0.5247
K2PG	13	-1419.4830	2922.2554	0.4127	0.9374
K2PIG	14	-1418.5772	2926.8507	0.0415	0.9789
HKYI	16	-1413.4924	2929.4949	0.0111	0.9900
HKYG	16	-1413.6923	2929.8948	0.0091	0.9990
HKYIG	17	-1412.7132	2934.3434	0.0010	1.0000
SYMI	17	-1418.3580	2945.6329	0.0000	1.0000
SYMG	17	-1418.6027	2946.1223	0.0000	1.0000
SYMIG	18	-1417.7443	2950.8124	0.0000	1.0000
GTRI	20	-1412.4664	2953.0703	0.0000	1.0000
GTRG	20	-1412.6646	2953.4669	0.0000	1.0000
GTRIG	21	-1411.7334	2958.0114	0.0000	1.0000
JC69I	12	-1443.0633	2963.0091	0.0000	1.0000
JC69G	12	-1443.1968	2963.2761	0.0000	1.0000
K2P	12	-1445.4259	2967.7343	0.0000	1.0000
JC69IG	13	-1442.6789	2968.6472	0.0000	1.0000
F81I	15	-1437.4196	2970.9424	0.0000	1.0000
F81G	15	-1437.5464	2971.1959	0.0000	1.0000
HKY	15	-1439.4518	2975.0068	0.0000	1.0000
F81IG	16	-1436.9851	2976.4803	0.0000	1.0000
SYM	16	-1444.5603	2991.6306	0.0000	1.0000
GTR	19	-1438.4356	2998.6020	0.0000	1.0000
JC69	11	-1468.2566	3006.9888	0.0000	1.0000
F81	14	-1462.5663	3014.8289	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/chr7_10316.phylip = 1 - 606;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10316.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10316.phylip = 1 - 606;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10316.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 BIC model (K2PI)]

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

