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

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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYIG
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
HKYIG	61	-1277.9668	2701.7949	0.5640	0.5640
GTRIG	65	-1272.9098	2703.2317	0.2750	0.8390
HKYG	60	-1280.8072	2704.6332	0.1364	0.9754
GTRG	64	-1277.1858	2708.8684	0.0164	0.9918
HKYI	60	-1283.8092	2710.6372	0.0068	0.9986
GTRI	64	-1279.8631	2714.2230	0.0011	0.9997
K2PIG	58	-1290.3103	2718.0081	0.0002	0.9999
K2PG	57	-1292.6533	2719.9047	0.0001	1.0000
SYMIG	62	-1287.4933	2723.7081	0.0000	1.0000
K2PI	57	-1295.2447	2725.0875	0.0000	1.0000
SYMG	61	-1290.6794	2727.2200	0.0000	1.0000
SYMI	61	-1292.6313	2731.1238	0.0000	1.0000
F81IG	60	-1318.7033	2780.4255	0.0000	1.0000
F81G	59	-1321.4127	2783.0198	0.0000	1.0000
F81I	59	-1323.3878	2786.9699	0.0000	1.0000
JC69IG	57	-1327.2922	2789.1825	0.0000	1.0000
JC69G	56	-1329.7979	2791.4218	0.0000	1.0000
HKY	59	-1327.3439	2794.8821	0.0000	1.0000
JC69I	56	-1331.6715	2795.1690	0.0000	1.0000
GTR	63	-1324.6029	2800.8058	0.0000	1.0000
K2P	56	-1338.7519	2809.3298	0.0000	1.0000
SYM	60	-1337.0745	2817.1679	0.0000	1.0000
F81	58	-1366.0919	2869.5712	0.0000	1.0000
JC69	55	-1374.4269	2877.9251	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1272.9098	2675.8196	0.6917	0.6917
HKYIG	61	-1277.9668	2677.9337	0.2403	0.9320
HKYG	60	-1280.8072	2681.6144	0.0382	0.9701
GTRG	64	-1277.1858	2682.3716	0.0261	0.9963
HKYI	60	-1283.8092	2687.6183	0.0019	0.9982
GTRI	64	-1279.8631	2687.7262	0.0018	1.0000
K2PIG	58	-1290.3103	2696.6206	0.0000	1.0000
SYMIG	62	-1287.4933	2698.9866	0.0000	1.0000
K2PG	57	-1292.6533	2699.3065	0.0000	1.0000
SYMG	61	-1290.6794	2703.3588	0.0000	1.0000
K2PI	57	-1295.2447	2704.4894	0.0000	1.0000
SYMI	61	-1292.6313	2707.2626	0.0000	1.0000
F81IG	60	-1318.7033	2757.4067	0.0000	1.0000
F81G	59	-1321.4127	2760.8255	0.0000	1.0000
F81I	59	-1323.3878	2764.7755	0.0000	1.0000
JC69IG	57	-1327.2922	2768.5844	0.0000	1.0000
JC69G	56	-1329.7979	2771.5957	0.0000	1.0000
HKY	59	-1327.3439	2772.6878	0.0000	1.0000
GTR	63	-1324.6029	2775.2058	0.0000	1.0000
JC69I	56	-1331.6715	2775.3429	0.0000	1.0000
K2P	56	-1338.7519	2789.5038	0.0000	1.0000
SYM	60	-1337.0745	2794.1490	0.0000	1.0000
F81	58	-1366.0919	2848.1837	0.0000	1.0000
JC69	55	-1374.4269	2858.8539	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1280.8072	2917.8665	0.4961	0.4961
HKYIG	61	-1277.9668	2918.1234	0.4364	0.9325
K2PG	57	-1292.6533	2923.7461	0.0262	0.9587
HKYI	60	-1283.8092	2923.8705	0.0247	0.9834
K2PIG	58	-1290.3103	2924.9977	0.0140	0.9974
K2PI	57	-1295.2447	2928.9289	0.0020	0.9994
GTRIG	65	-1272.9098	2931.7594	0.0005	0.9999
GTRG	64	-1277.1858	2934.3739	0.0001	1.0000
GTRI	64	-1279.8631	2939.7285	0.0000	1.0000
SYMIG	62	-1287.4933	2943.1139	0.0000	1.0000
SYMG	61	-1290.6794	2943.5485	0.0000	1.0000
SYMI	61	-1292.6313	2947.4523	0.0000	1.0000
JC69G	56	-1329.7979	2992.0977	0.0000	1.0000
JC69IG	57	-1327.2922	2993.0239	0.0000	1.0000
F81G	59	-1321.4127	2993.1401	0.0000	1.0000
F81IG	60	-1318.7033	2993.6588	0.0000	1.0000
JC69I	56	-1331.6715	2995.8449	0.0000	1.0000
F81I	59	-1323.3878	2997.0902	0.0000	1.0000
HKY	59	-1327.3439	3005.0024	0.0000	1.0000
K2P	56	-1338.7519	3010.0058	0.0000	1.0000
GTR	63	-1324.6029	3023.2706	0.0000	1.0000
SYM	60	-1337.0745	3030.4012	0.0000	1.0000
JC69	55	-1374.4269	3075.4184	0.0000	1.0000
F81	58	-1366.0919	3076.5608	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5590_hg18.phylip = 1 - 379;
 Partition Dummy = 1:261-aggressive-phylip/5590_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5590_hg18.phylip = 1 - 379;
 Partition Dummy = 1:261-aggressive-phylip/5590_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

