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

Input data from file "261-aggressive-phylip/6253_hg18.phylip" (29 taxa, 370 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	-1249.4067	2654.0920	0.6598	0.6598
GTRIG	65	-1248.6471	2655.5179	0.3234	0.9833
HKYG	60	-1259.5901	2662.8695	0.0082	0.9914
HKYIG	61	-1258.5834	2663.7252	0.0053	0.9968
GTRI	64	-1254.7525	2664.7837	0.0031	0.9999
HKYI	60	-1264.3688	2672.4269	0.0001	1.0000
K2PIG	58	-1279.8381	2697.6826	0.0000	1.0000
SYMG	61	-1276.7965	2700.1514	0.0000	1.0000
SYMIG	62	-1275.8307	2701.1076	0.0000	1.0000
K2PG	57	-1284.6112	2704.4147	0.0000	1.0000
K2PI	57	-1284.9072	2705.0067	0.0000	1.0000
GTR	63	-1277.4824	2707.3178	0.0000	1.0000
SYMI	61	-1281.1902	2708.9389	0.0000	1.0000
HKY	59	-1292.3678	2725.5743	0.0000	1.0000
F81G	59	-1293.2857	2727.4101	0.0000	1.0000
F81IG	60	-1292.7126	2729.1145	0.0000	1.0000
F81I	59	-1296.9038	2734.6463	0.0000	1.0000
JC69IG	57	-1310.2386	2755.6696	0.0000	1.0000
JC69G	56	-1313.2803	2758.9567	0.0000	1.0000
K2P	56	-1314.3505	2761.0972	0.0000	1.0000
SYM	60	-1310.1423	2763.9740	0.0000	1.0000
JC69I	56	-1317.0540	2766.5042	0.0000	1.0000
F81	58	-1320.4995	2779.0054	0.0000	1.0000
JC69	55	-1341.8023	2813.2224	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1249.4067	2626.8133	0.5570	0.5570
GTRIG	65	-1248.6471	2627.2942	0.4380	0.9950
GTRI	64	-1254.7525	2637.5050	0.0027	0.9977
HKYIG	61	-1258.5834	2639.1668	0.0012	0.9988
HKYG	60	-1259.5901	2639.1802	0.0011	1.0000
HKYI	60	-1264.3688	2648.7375	0.0000	1.0000
SYMG	61	-1276.7965	2675.5929	0.0000	1.0000
SYMIG	62	-1275.8307	2675.6614	0.0000	1.0000
K2PIG	58	-1279.8381	2675.6762	0.0000	1.0000
GTR	63	-1277.4824	2680.9649	0.0000	1.0000
K2PG	57	-1284.6112	2683.2224	0.0000	1.0000
K2PI	57	-1284.9072	2683.8144	0.0000	1.0000
SYMI	61	-1281.1902	2684.3804	0.0000	1.0000
HKY	59	-1292.3678	2702.7356	0.0000	1.0000
F81G	59	-1293.2857	2704.5714	0.0000	1.0000
F81IG	60	-1292.7126	2705.4252	0.0000	1.0000
F81I	59	-1296.9038	2711.8076	0.0000	1.0000
JC69IG	57	-1310.2386	2734.4773	0.0000	1.0000
JC69G	56	-1313.2803	2738.5605	0.0000	1.0000
SYM	60	-1310.1423	2740.2847	0.0000	1.0000
K2P	56	-1314.3505	2740.7011	0.0000	1.0000
JC69I	56	-1317.0540	2746.1080	0.0000	1.0000
F81	58	-1320.4995	2756.9990	0.0000	1.0000
JC69	55	-1341.8023	2793.6045	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1259.5901	2873.9904	0.7319	0.7319
GTRG	64	-1249.4067	2877.2775	0.1415	0.8733
HKYIG	61	-1258.5834	2877.8904	0.1041	0.9775
GTRIG	65	-1248.6471	2881.6719	0.0157	0.9932
HKYI	60	-1264.3688	2883.5477	0.0062	0.9993
GTRI	64	-1254.7525	2887.9692	0.0007	1.0000
K2PIG	58	-1279.8381	2902.6594	0.0000	1.0000
K2PG	57	-1284.6112	2906.2921	0.0000	1.0000
K2PI	57	-1284.9072	2906.8840	0.0000	1.0000
SYMG	61	-1276.7965	2914.3166	0.0000	1.0000
SYMIG	62	-1275.8307	2918.2985	0.0000	1.0000
SYMI	61	-1281.1902	2923.1041	0.0000	1.0000
GTR	63	-1277.4824	2927.5155	0.0000	1.0000
HKY	59	-1292.3678	2933.6322	0.0000	1.0000
F81G	59	-1293.2857	2935.4680	0.0000	1.0000
F81IG	60	-1292.7126	2940.2354	0.0000	1.0000
F81I	59	-1296.9038	2942.7043	0.0000	1.0000
JC69IG	57	-1310.2386	2957.5470	0.0000	1.0000
JC69G	56	-1313.2803	2957.7167	0.0000	1.0000
K2P	56	-1314.3505	2959.8572	0.0000	1.0000
JC69I	56	-1317.0540	2965.2642	0.0000	1.0000
SYM	60	-1310.1423	2975.0949	0.0000	1.0000
F81	58	-1320.4995	2983.9822	0.0000	1.0000
JC69	55	-1341.8023	3008.8472	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/6253_hg18.phylip = 1 - 370;
 Partition Dummy = 1:261-aggressive-phylip/6253_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/6253_hg18.phylip = 1 - 370;
 Partition Dummy = 1:261-aggressive-phylip/6253_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/6253_hg18.phylip = 1 - 370;
 Partition Dummy = 1:261-aggressive-phylip/6253_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

