
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:18 08/01/2011

Input data from file "261-aggressive-phylip/10558_hg18.phylip" (29 taxa, 751 characters)

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


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	64	-1385.2695	2910.6673	0.7167	0.7167
GTRIG	65	-1385.2962	2913.1179	0.2105	0.9272
HKYG	60	-1392.5788	2915.7664	0.0560	0.9832
HKYIG	61	-1392.6055	2918.1892	0.0167	0.9998
GTR	63	-1395.2575	2928.2531	0.0001	0.9999
GTRI	64	-1395.2822	2930.6927	0.0000	1.0000
HKY	59	-1401.8075	2931.8610	0.0000	1.0000
HKYI	60	-1401.8324	2934.2734	0.0000	1.0000
K2PG	57	-1411.7298	2947.0008	0.0000	1.0000
K2PIG	58	-1411.7565	2949.4032	0.0000	1.0000
SYMG	61	-1409.4060	2951.7901	0.0000	1.0000
SYMIG	62	-1409.4323	2954.2193	0.0000	1.0000
K2P	56	-1421.7175	2964.6338	0.0000	1.0000
K2PI	57	-1421.7423	2967.0257	0.0000	1.0000
SYM	60	-1419.4863	2969.5814	0.0000	1.0000
SYMI	61	-1419.5111	2972.0004	0.0000	1.0000
F81G	59	-1427.7450	2983.7361	0.0000	1.0000
F81IG	60	-1427.7716	2986.1520	0.0000	1.0000
F81	58	-1436.5782	2999.0466	0.0000	1.0000
F81I	59	-1436.6031	3001.4523	0.0000	1.0000
JC69G	56	-1443.9756	3009.1500	0.0000	1.0000
JC69IG	57	-1444.0021	3011.5454	0.0000	1.0000
JC69	55	-1453.2163	3025.2960	0.0000	1.0000
JC69I	56	-1453.2411	3027.6810	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1385.2695	2898.5390	0.7102	0.7102
GTRIG	65	-1385.2962	2900.5923	0.2544	0.9646
HKYG	60	-1392.5788	2905.1577	0.0260	0.9906
HKYIG	61	-1392.6055	2907.2110	0.0093	0.9999
GTR	63	-1395.2575	2916.5151	0.0001	1.0000
GTRI	64	-1395.2822	2918.5645	0.0000	1.0000
HKY	59	-1401.8075	2921.6150	0.0000	1.0000
HKYI	60	-1401.8324	2923.6647	0.0000	1.0000
K2PG	57	-1411.7298	2937.4597	0.0000	1.0000
K2PIG	58	-1411.7565	2939.5130	0.0000	1.0000
SYMG	61	-1409.4060	2940.8119	0.0000	1.0000
SYMIG	62	-1409.4323	2942.8646	0.0000	1.0000
K2P	56	-1421.7175	2955.4349	0.0000	1.0000
K2PI	57	-1421.7423	2957.4845	0.0000	1.0000
SYM	60	-1419.4863	2958.9727	0.0000	1.0000
SYMI	61	-1419.5111	2961.0222	0.0000	1.0000
F81G	59	-1427.7450	2973.4901	0.0000	1.0000
F81IG	60	-1427.7716	2975.5433	0.0000	1.0000
F81	58	-1436.5782	2989.1564	0.0000	1.0000
F81I	59	-1436.6031	2991.2062	0.0000	1.0000
JC69G	56	-1443.9756	2999.9512	0.0000	1.0000
JC69IG	57	-1444.0021	3002.0043	0.0000	1.0000
JC69	55	-1453.2163	3016.4327	0.0000	1.0000
JC69I	56	-1453.2411	3018.4822	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1392.5788	3182.4420	0.9605	0.9605
HKYIG	61	-1392.6055	3189.1167	0.0341	0.9946
HKY	59	-1401.8075	3194.2779	0.0026	0.9972
GTRG	64	-1385.2695	3194.3090	0.0025	0.9997
K2PG	57	-1411.7298	3200.8798	0.0001	0.9998
HKYI	60	-1401.8324	3200.9491	0.0001	0.9999
GTRIG	65	-1385.2962	3200.9837	0.0001	1.0000
K2PIG	58	-1411.7565	3207.5545	0.0000	1.0000
GTR	63	-1395.2575	3207.6636	0.0000	1.0000
K2P	56	-1421.7175	3214.2337	0.0000	1.0000
GTRI	64	-1395.2822	3214.3344	0.0000	1.0000
K2PI	57	-1421.7423	3220.9047	0.0000	1.0000
SYMG	61	-1409.4060	3222.7176	0.0000	1.0000
SYMIG	62	-1409.4323	3229.3918	0.0000	1.0000
SYM	60	-1419.4863	3236.2570	0.0000	1.0000
SYMI	61	-1419.5111	3242.9279	0.0000	1.0000
F81G	59	-1427.7450	3246.1530	0.0000	1.0000
F81IG	60	-1427.7716	3252.8276	0.0000	1.0000
F81	58	-1436.5782	3257.1979	0.0000	1.0000
JC69G	56	-1443.9756	3258.7499	0.0000	1.0000
F81I	59	-1436.6031	3263.8692	0.0000	1.0000
JC69IG	57	-1444.0021	3265.4244	0.0000	1.0000
JC69	55	-1453.2163	3270.6100	0.0000	1.0000
JC69I	56	-1453.2411	3277.2809	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 261-aggressive-phylip/10558_hg18.phylip = 1 - 751;
 Partition Dummy = 1:261-aggressive-phylip/10558_hg18.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 261-aggressive-phylip/10558_hg18.phylip = 1 - 751;
 Partition Dummy = 1:261-aggressive-phylip/10558_hg18.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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10558_hg18.phylip = 1 - 751;
 Partition Dummy = 1:261-aggressive-phylip/10558_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

