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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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	-1375.0016	2791.3537	0.3872	0.3872
GTRIG	21	-1374.0959	2791.6798	0.3290	0.7162
GTRG	20	-1375.3557	2792.0618	0.2718	0.9880
HKYI	16	-1383.4658	2799.8006	0.0057	0.9936
HKYG	16	-1384.0429	2800.9547	0.0032	0.9968
HKYIG	17	-1382.9925	2800.9641	0.0032	1.0000
GTR	19	-1401.7697	2842.7594	0.0000	1.0000
HKY	15	-1409.6208	2850.0071	0.0000	1.0000
SYMI	17	-1416.1512	2867.2816	0.0000	1.0000
SYMG	17	-1416.9953	2868.9698	0.0000	1.0000
F81I	15	-1419.5357	2869.8369	0.0000	1.0000
F81G	15	-1419.8389	2870.4433	0.0000	1.0000
F81IG	16	-1419.0757	2871.0204	0.0000	1.0000
SYMIG	18	-1417.0814	2871.2590	0.0000	1.0000
K2PI	13	-1426.8644	2880.3074	0.0000	1.0000
K2PG	13	-1427.8031	2882.1850	0.0000	1.0000
K2PIG	14	-1426.8171	2882.3029	0.0000	1.0000
F81	14	-1443.5509	2915.7706	0.0000	1.0000
SYM	16	-1444.4473	2921.7636	0.0000	1.0000
K2P	12	-1454.3690	2933.2333	0.0000	1.0000
JC69I	12	-1459.8091	2944.1133	0.0000	1.0000
JC69G	12	-1460.1656	2944.8264	0.0000	1.0000
JC69IG	13	-1459.5956	2945.7699	0.0000	1.0000
JC69	11	-1484.4545	2991.3273	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1375.0016	2790.0032	0.3793	0.3793
GTRIG	21	-1374.0959	2790.1919	0.3451	0.7244
GTRG	20	-1375.3557	2790.7113	0.2662	0.9906
HKYI	16	-1383.4658	2798.9316	0.0044	0.9950
HKYIG	17	-1382.9925	2799.9849	0.0026	0.9975
HKYG	16	-1384.0429	2800.0857	0.0025	1.0000
GTR	19	-1401.7697	2841.5395	0.0000	1.0000
HKY	15	-1409.6208	2849.2416	0.0000	1.0000
SYMI	17	-1416.1512	2866.3024	0.0000	1.0000
SYMG	17	-1416.9953	2867.9906	0.0000	1.0000
F81I	15	-1419.5357	2869.0714	0.0000	1.0000
F81G	15	-1419.8389	2869.6777	0.0000	1.0000
F81IG	16	-1419.0757	2870.1514	0.0000	1.0000
SYMIG	18	-1417.0814	2870.1628	0.0000	1.0000
K2PI	13	-1426.8644	2879.7287	0.0000	1.0000
K2PG	13	-1427.8031	2881.6063	0.0000	1.0000
K2PIG	14	-1426.8171	2881.6341	0.0000	1.0000
F81	14	-1443.5509	2915.1019	0.0000	1.0000
SYM	16	-1444.4473	2920.8946	0.0000	1.0000
K2P	12	-1454.3690	2932.7380	0.0000	1.0000
JC69I	12	-1459.8091	2943.6181	0.0000	1.0000
JC69G	12	-1460.1656	2944.3311	0.0000	1.0000
JC69IG	13	-1459.5956	2945.1912	0.0000	1.0000
JC69	11	-1484.4545	2990.9090	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1383.4658	2870.3899	0.6077	0.6077
HKYG	16	-1384.0429	2871.5440	0.3413	0.9490
HKYIG	17	-1382.9925	2875.9094	0.0385	0.9875
GTRI	20	-1375.0016	2879.3261	0.0070	0.9944
GTRG	20	-1375.3557	2880.0342	0.0049	0.9993
GTRIG	21	-1374.0959	2883.9809	0.0007	1.0000
HKY	15	-1409.6208	2916.2338	0.0000	1.0000
GTR	19	-1401.7697	2926.3962	0.0000	1.0000
F81I	15	-1419.5357	2936.0635	0.0000	1.0000
F81G	15	-1419.8389	2936.6699	0.0000	1.0000
K2PI	13	-1426.8644	2937.7886	0.0000	1.0000
K2PG	13	-1427.8031	2939.6662	0.0000	1.0000
F81IG	16	-1419.0757	2941.6098	0.0000	1.0000
SYMI	17	-1416.1512	2942.2269	0.0000	1.0000
SYMG	17	-1416.9953	2943.9151	0.0000	1.0000
K2PIG	14	-1426.8171	2944.1601	0.0000	1.0000
SYMIG	18	-1417.0814	2950.5534	0.0000	1.0000
F81	14	-1443.5509	2977.6279	0.0000	1.0000
K2P	12	-1454.3690	2986.3318	0.0000	1.0000
SYM	16	-1444.4473	2992.3529	0.0000	1.0000
JC69I	12	-1459.8091	2997.2118	0.0000	1.0000
JC69G	12	-1460.1656	2997.9249	0.0000	1.0000
JC69IG	13	-1459.5956	3003.2510	0.0000	1.0000
JC69	11	-1484.4545	3040.0366	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/chr11_3256.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3256.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/chr11_3256.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3256.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 (HKYI)]

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


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

