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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYIG


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	-1398.5525	2937.6898	0.9240	0.9240
GTRIG	65	-1395.9917	2943.2737	0.0566	0.9807
HKYI	60	-1404.2317	2946.4047	0.0118	0.9925
HKYG	60	-1404.8015	2947.5442	0.0067	0.9992
GTRI	64	-1402.1533	2952.9007	0.0005	0.9997
GTRG	64	-1402.4937	2953.5814	0.0003	1.0000
F81IG	60	-1422.5552	2983.0516	0.0000	1.0000
F81I	59	-1427.6163	2990.5430	0.0000	1.0000
K2PIG	58	-1429.3622	2991.4172	0.0000	1.0000
F81G	59	-1428.3695	2992.0496	0.0000	1.0000
K2PI	57	-1434.6665	2999.4206	0.0000	1.0000
SYMIG	62	-1428.3482	2999.9378	0.0000	1.0000
K2PG	57	-1435.3203	3000.7281	0.0000	1.0000
SYMI	61	-1433.2116	3007.0081	0.0000	1.0000
SYMG	61	-1434.1519	3008.8886	0.0000	1.0000
JC69IG	57	-1450.1114	3030.3105	0.0000	1.0000
JC69I	56	-1455.0951	3037.6853	0.0000	1.0000
JC69G	56	-1455.7771	3039.0493	0.0000	1.0000
HKY	59	-1461.7402	3058.7910	0.0000	1.0000
GTR	63	-1459.7107	3065.3325	0.0000	1.0000
F81	58	-1485.0969	3102.8865	0.0000	1.0000
K2P	56	-1493.2781	3114.0513	0.0000	1.0000
SYM	60	-1492.4416	3122.8243	0.0000	1.0000
JC69	55	-1513.2890	3151.4932	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1398.5525	2919.1050	0.7975	0.7975
GTRIG	65	-1395.9917	2921.9834	0.1891	0.9866
HKYI	60	-1404.2317	2928.4635	0.0074	0.9940
HKYG	60	-1404.8015	2929.6031	0.0042	0.9981
GTRI	64	-1402.1533	2932.3066	0.0011	0.9992
GTRG	64	-1402.4937	2932.9874	0.0008	1.0000
F81IG	60	-1422.5552	2965.1105	0.0000	1.0000
F81I	59	-1427.6163	2973.2325	0.0000	1.0000
K2PIG	58	-1429.3622	2974.7245	0.0000	1.0000
F81G	59	-1428.3695	2974.7391	0.0000	1.0000
SYMIG	62	-1428.3482	2980.6964	0.0000	1.0000
K2PI	57	-1434.6665	2983.3330	0.0000	1.0000
K2PG	57	-1435.3203	2984.6405	0.0000	1.0000
SYMI	61	-1433.2116	2988.4233	0.0000	1.0000
SYMG	61	-1434.1519	2990.3039	0.0000	1.0000
JC69IG	57	-1450.1114	3014.2229	0.0000	1.0000
JC69I	56	-1455.0951	3022.1902	0.0000	1.0000
JC69G	56	-1455.7771	3023.5541	0.0000	1.0000
HKY	59	-1461.7402	3041.4805	0.0000	1.0000
GTR	63	-1459.7107	3045.4213	0.0000	1.0000
F81	58	-1485.0969	3086.1938	0.0000	1.0000
K2P	56	-1493.2781	3098.5561	0.0000	1.0000
SYM	60	-1492.4416	3104.8832	0.0000	1.0000
JC69	55	-1513.2890	3136.5779	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-1398.5525	3172.2918	0.8961	0.8961
HKYI	60	-1404.2317	3177.4997	0.0663	0.9624
HKYG	60	-1404.8015	3178.6392	0.0375	0.9999
GTRIG	65	-1395.9917	3191.7725	0.0001	1.0000
GTRI	64	-1402.1533	3197.9452	0.0000	1.0000
GTRG	64	-1402.4937	3198.6259	0.0000	1.0000
F81IG	60	-1422.5552	3214.1466	0.0000	1.0000
K2PIG	58	-1429.3622	3215.4594	0.0000	1.0000
F81I	59	-1427.6163	3218.1181	0.0000	1.0000
F81G	59	-1428.3695	3219.6246	0.0000	1.0000
K2PI	57	-1434.6665	3219.9174	0.0000	1.0000
K2PG	57	-1435.3203	3221.2249	0.0000	1.0000
SYMIG	62	-1428.3482	3238.0338	0.0000	1.0000
SYMI	61	-1433.2116	3241.6101	0.0000	1.0000
SYMG	61	-1434.1519	3243.4906	0.0000	1.0000
JC69IG	57	-1450.1114	3250.8072	0.0000	1.0000
JC69I	56	-1455.0951	3254.6239	0.0000	1.0000
JC69G	56	-1455.7771	3255.9879	0.0000	1.0000
HKY	59	-1461.7402	3286.3660	0.0000	1.0000
GTR	63	-1459.7107	3306.9093	0.0000	1.0000
F81	58	-1485.0969	3326.9288	0.0000	1.0000
K2P	56	-1493.2781	3330.9899	0.0000	1.0000
SYM	60	-1492.4416	3353.9193	0.0000	1.0000
JC69	55	-1513.2890	3364.8611	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1615_hg18.phylip = 1 - 469;
 Partition Dummy = 1:261-aggressive-phylip/1615_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 AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1615_hg18.phylip = 1 - 469;
 Partition Dummy = 1:261-aggressive-phylip/1615_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 (HKYIG)]

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


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

