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

Input data from file "261-aggressive-phylip/5745_hg18.phylip" (29 taxa, 328 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	64	-1297.9175	2755.4700	0.9892	0.9892
GTRIG	65	-1301.0337	2764.8154	0.0092	0.9984
HKYG	60	-1310.9134	2769.2426	0.0010	0.9994
GTRI	64	-1305.9515	2771.5380	0.0003	0.9997
HKYIG	61	-1310.8401	2772.1163	0.0002	1.0000
HKYI	60	-1314.7738	2776.9634	0.0000	1.0000
K2PG	57	-1325.8289	2790.1467	0.0000	1.0000
SYMG	61	-1320.9679	2792.3719	0.0000	1.0000
K2PIG	58	-1325.7378	2792.9180	0.0000	1.0000
SYMIG	62	-1320.9319	2795.3431	0.0000	1.0000
K2PI	57	-1329.0149	2796.5187	0.0000	1.0000
SYMI	61	-1324.1044	2798.6448	0.0000	1.0000
GTR	63	-1321.6968	2799.9390	0.0000	1.0000
HKY	59	-1330.4071	2805.2322	0.0000	1.0000
K2P	56	-1345.7917	2827.1405	0.0000	1.0000
SYM	60	-1339.9790	2827.3737	0.0000	1.0000
F81G	59	-1358.3311	2861.0801	0.0000	1.0000
F81IG	60	-1358.1967	2863.8092	0.0000	1.0000
F81I	59	-1361.7978	2868.0136	0.0000	1.0000
JC69G	56	-1368.3684	2872.2939	0.0000	1.0000
JC69IG	57	-1368.2486	2874.9861	0.0000	1.0000
JC69I	56	-1371.5785	2878.7143	0.0000	1.0000
F81	58	-1377.3958	2896.2339	0.0000	1.0000
JC69	55	-1387.6119	2907.8709	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1297.9175	2723.8350	0.9835	0.9835
GTRIG	65	-1301.0337	2732.0673	0.0160	0.9995
GTRI	64	-1305.9515	2739.9031	0.0003	0.9998
HKYG	60	-1310.9134	2741.8269	0.0001	0.9999
HKYIG	61	-1310.8401	2743.6802	0.0000	1.0000
HKYI	60	-1314.7738	2749.5477	0.0000	1.0000
SYMG	61	-1320.9679	2763.9358	0.0000	1.0000
K2PG	57	-1325.8289	2765.6578	0.0000	1.0000
SYMIG	62	-1320.9319	2765.8639	0.0000	1.0000
K2PIG	58	-1325.7378	2767.4756	0.0000	1.0000
GTR	63	-1321.6968	2769.3935	0.0000	1.0000
SYMI	61	-1324.1044	2770.2087	0.0000	1.0000
K2PI	57	-1329.0149	2772.0298	0.0000	1.0000
HKY	59	-1330.4071	2778.8143	0.0000	1.0000
SYM	60	-1339.9790	2799.9580	0.0000	1.0000
K2P	56	-1345.7917	2803.5833	0.0000	1.0000
F81G	59	-1358.3311	2834.6622	0.0000	1.0000
F81IG	60	-1358.1967	2836.3934	0.0000	1.0000
F81I	59	-1361.7978	2841.5956	0.0000	1.0000
JC69G	56	-1368.3684	2848.7367	0.0000	1.0000
JC69IG	57	-1368.2486	2850.4973	0.0000	1.0000
JC69I	56	-1371.5785	2855.1571	0.0000	1.0000
F81	58	-1377.3958	2870.7916	0.0000	1.0000
JC69	55	-1387.6119	2885.2238	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1297.9175	2966.5879	0.7892	0.7892
HKYG	60	-1310.9134	2969.4077	0.1927	0.9819
HKYIG	61	-1310.8401	2975.0540	0.0114	0.9933
HKYI	60	-1314.7738	2977.1285	0.0041	0.9974
GTRIG	65	-1301.0337	2978.6132	0.0019	0.9993
K2PG	57	-1325.8289	2981.8596	0.0004	0.9997
GTRI	64	-1305.9515	2982.6559	0.0003	1.0000
K2PIG	58	-1325.7378	2987.4704	0.0000	1.0000
K2PI	57	-1329.0149	2988.2316	0.0000	1.0000
SYMG	61	-1320.9679	2995.3096	0.0000	1.0000
SYMIG	62	-1320.9319	3001.0307	0.0000	1.0000
SYMI	61	-1324.1044	3001.5826	0.0000	1.0000
HKY	59	-1330.4071	3002.6021	0.0000	1.0000
GTR	63	-1321.6968	3008.3534	0.0000	1.0000
K2P	56	-1345.7917	3015.9921	0.0000	1.0000
SYM	60	-1339.9790	3027.5388	0.0000	1.0000
F81G	59	-1358.3311	3058.4500	0.0000	1.0000
JC69G	56	-1368.3684	3061.1455	0.0000	1.0000
F81IG	60	-1358.1967	3063.9742	0.0000	1.0000
F81I	59	-1361.7978	3065.3834	0.0000	1.0000
JC69IG	57	-1368.2486	3066.6990	0.0000	1.0000
JC69I	56	-1371.5785	3067.5658	0.0000	1.0000
F81	58	-1377.3958	3090.7863	0.0000	1.0000
JC69	55	-1387.6119	3093.8396	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/5745_hg18.phylip = 1 - 328;
 Partition Dummy = 1:261-aggressive-phylip/5745_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/5745_hg18.phylip = 1 - 328;
 Partition Dummy = 1:261-aggressive-phylip/5745_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5745_hg18.phylip = 1 - 328;
 Partition Dummy = 1:261-aggressive-phylip/5745_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;


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

