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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	20	-1392.0641	2825.4986	0.4445	0.4445
GTRI	20	-1392.2527	2825.8758	0.3681	0.8126
GTRIG	21	-1391.8582	2827.2261	0.1874	1.0000
HKYG	16	-1407.3291	2847.5398	0.0000	1.0000
HKYI	16	-1407.4624	2847.8064	0.0000	1.0000
HKYIG	17	-1407.1308	2849.2551	0.0000	1.0000
GTR	19	-1406.9347	2853.1071	0.0000	1.0000
HKY	15	-1421.5270	2873.8308	0.0000	1.0000
F81I	15	-1423.7979	2878.3724	0.0000	1.0000
F81IG	16	-1423.5332	2879.9482	0.0000	1.0000
F81G	15	-1425.9546	2882.6860	0.0000	1.0000
F81	14	-1436.7609	2902.2003	0.0000	1.0000
SYMIG	18	-1434.4904	2906.0930	0.0000	1.0000
SYMG	17	-1436.0015	2906.9966	0.0000	1.0000
SYMI	17	-1436.0337	2907.0609	0.0000	1.0000
K2PI	13	-1446.3927	2919.3724	0.0000	1.0000
K2PIG	14	-1446.1328	2920.9441	0.0000	1.0000
K2PG	13	-1447.5426	2921.6723	0.0000	1.0000
SYM	16	-1449.7006	2932.2829	0.0000	1.0000
JC69I	12	-1459.6724	2943.8473	0.0000	1.0000
JC69IG	13	-1459.3334	2945.2539	0.0000	1.0000
JC69G	12	-1460.6728	2945.8481	0.0000	1.0000
K2P	12	-1460.8804	2946.2632	0.0000	1.0000
JC69	11	-1473.0109	2968.4462	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1392.0641	2824.1283	0.4386	0.4386
GTRI	20	-1392.2527	2824.5055	0.3632	0.8017
GTRIG	21	-1391.8582	2825.7163	0.1982	1.0000
HKYG	16	-1407.3291	2846.6582	0.0000	1.0000
HKYI	16	-1407.4624	2846.9247	0.0000	1.0000
HKYIG	17	-1407.1308	2848.2616	0.0000	1.0000
GTR	19	-1406.9347	2851.8694	0.0000	1.0000
HKY	15	-1421.5270	2873.0541	0.0000	1.0000
F81I	15	-1423.7979	2877.5957	0.0000	1.0000
F81IG	16	-1423.5332	2879.0665	0.0000	1.0000
F81G	15	-1425.9546	2881.9093	0.0000	1.0000
F81	14	-1436.7609	2901.5218	0.0000	1.0000
SYMIG	18	-1434.4904	2904.9808	0.0000	1.0000
SYMG	17	-1436.0015	2906.0031	0.0000	1.0000
SYMI	17	-1436.0337	2906.0674	0.0000	1.0000
K2PI	13	-1446.3927	2918.7853	0.0000	1.0000
K2PIG	14	-1446.1328	2920.2656	0.0000	1.0000
K2PG	13	-1447.5426	2921.0852	0.0000	1.0000
SYM	16	-1449.7006	2931.4012	0.0000	1.0000
JC69I	12	-1459.6724	2943.3449	0.0000	1.0000
JC69IG	13	-1459.3334	2944.6668	0.0000	1.0000
JC69G	12	-1460.6728	2945.3457	0.0000	1.0000
K2P	12	-1460.8804	2945.7608	0.0000	1.0000
JC69	11	-1473.0109	2968.0218	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-1392.0641	2913.1693	0.4858	0.4858
GTRI	20	-1392.2527	2913.5465	0.4023	0.8881
HKYG	16	-1407.3291	2917.8909	0.0458	0.9340
HKYI	16	-1407.4624	2918.1575	0.0401	0.9741
GTRIG	21	-1391.8582	2919.2093	0.0237	0.9978
HKYIG	17	-1407.1308	2923.9465	0.0022	1.0000
GTR	19	-1406.9347	2936.4583	0.0000	1.0000
HKY	15	-1421.5270	2939.8348	0.0000	1.0000
F81I	15	-1423.7979	2944.3765	0.0000	1.0000
F81G	15	-1425.9546	2948.6900	0.0000	1.0000
F81IG	16	-1423.5332	2950.2993	0.0000	1.0000
F81	14	-1436.7609	2963.8505	0.0000	1.0000
K2PI	13	-1446.3927	2976.6620	0.0000	1.0000
K2PG	13	-1447.5426	2978.9618	0.0000	1.0000
SYMG	17	-1436.0015	2981.6879	0.0000	1.0000
SYMI	17	-1436.0337	2981.7523	0.0000	1.0000
K2PIG	14	-1446.1328	2982.5942	0.0000	1.0000
SYMIG	18	-1434.4904	2985.1177	0.0000	1.0000
JC69I	12	-1459.6724	2996.7695	0.0000	1.0000
JC69G	12	-1460.6728	2998.7702	0.0000	1.0000
K2P	12	-1460.8804	2999.1853	0.0000	1.0000
JC69IG	13	-1459.3334	3002.5434	0.0000	1.0000
SYM	16	-1449.7006	3002.6340	0.0000	1.0000
JC69	11	-1473.0109	3016.9943	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3428.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3428.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3428.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3428.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3428.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3428.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

