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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-1186.2655	2508.4440	0.3883	0.3883
HKYG	60	-1186.5528	2509.0187	0.2913	0.6795
HKYIG	61	-1185.6171	2509.7134	0.2058	0.8854
GTRI	64	-1183.1524	2512.5504	0.0498	0.9352
GTRG	64	-1183.3803	2513.0061	0.0397	0.9748
GTRIG	65	-1182.5562	2513.9695	0.0245	0.9994
F81I	59	-1194.7827	2522.9233	0.0003	0.9996
F81G	59	-1195.0976	2523.5531	0.0002	0.9998
F81IG	60	-1194.0317	2523.9765	0.0002	1.0000
HKY	59	-1203.4766	2540.3112	0.0000	1.0000
GTR	63	-1200.5130	2544.6715	0.0000	1.0000
F81	58	-1211.9990	2554.8119	0.0000	1.0000
K2PI	57	-1216.7751	2561.8309	0.0000	1.0000
K2PG	57	-1217.2359	2562.7525	0.0000	1.0000
K2PIG	58	-1216.1163	2563.0464	0.0000	1.0000
SYMI	61	-1215.8671	2570.2135	0.0000	1.0000
SYMG	61	-1216.3507	2571.1806	0.0000	1.0000
SYMIG	62	-1215.3402	2571.7372	0.0000	1.0000
JC69I	56	-1223.6337	2573.0259	0.0000	1.0000
JC69G	56	-1224.0462	2573.8510	0.0000	1.0000
JC69IG	57	-1222.8664	2574.0136	0.0000	1.0000
K2P	56	-1234.7036	2595.1658	0.0000	1.0000
SYM	60	-1233.8570	2603.6270	0.0000	1.0000
JC69	55	-1241.3252	2605.8977	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1186.2655	2492.5310	0.2882	0.2882
HKYG	60	-1186.5528	2493.1057	0.2162	0.5044
HKYIG	61	-1185.6171	2493.2341	0.2027	0.7071
GTRI	64	-1183.1524	2494.3048	0.1187	0.8258
GTRG	64	-1183.3803	2494.7605	0.0945	0.9203
GTRIG	65	-1182.5562	2495.1124	0.0793	0.9996
F81I	59	-1194.7827	2507.5654	0.0002	0.9998
F81IG	60	-1194.0317	2508.0634	0.0001	0.9999
F81G	59	-1195.0976	2508.1951	0.0001	1.0000
HKY	59	-1203.4766	2524.9533	0.0000	1.0000
GTR	63	-1200.5130	2527.0260	0.0000	1.0000
F81	58	-1211.9990	2539.9980	0.0000	1.0000
K2PI	57	-1216.7751	2547.5501	0.0000	1.0000
K2PIG	58	-1216.1163	2548.2325	0.0000	1.0000
K2PG	57	-1217.2359	2548.4717	0.0000	1.0000
SYMI	61	-1215.8671	2553.7342	0.0000	1.0000
SYMIG	62	-1215.3402	2554.6805	0.0000	1.0000
SYMG	61	-1216.3507	2554.7013	0.0000	1.0000
JC69I	56	-1223.6337	2559.2673	0.0000	1.0000
JC69IG	57	-1222.8664	2559.7329	0.0000	1.0000
JC69G	56	-1224.0462	2560.0924	0.0000	1.0000
K2P	56	-1234.7036	2581.4072	0.0000	1.0000
SYM	60	-1233.8570	2587.7140	0.0000	1.0000
JC69	55	-1241.3252	2592.6504	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1186.2655	2747.8760	0.5427	0.5427
HKYG	60	-1186.5528	2748.4507	0.4072	0.9499
HKYIG	61	-1185.6171	2752.8349	0.0455	0.9954
F81I	59	-1194.7827	2758.6547	0.0025	0.9979
F81G	59	-1195.0976	2759.2844	0.0018	0.9997
F81IG	60	-1194.0317	2763.4084	0.0002	0.9999
GTRI	64	-1183.1524	2766.6728	0.0000	1.0000
GTRG	64	-1183.3803	2767.1285	0.0000	1.0000
GTRIG	65	-1182.5562	2771.7361	0.0000	1.0000
HKY	59	-1203.4766	2776.0425	0.0000	1.0000
F81	58	-1211.9990	2786.8315	0.0000	1.0000
K2PI	57	-1216.7751	2790.1279	0.0000	1.0000
K2PG	57	-1217.2359	2791.0495	0.0000	1.0000
K2PIG	58	-1216.1163	2795.0660	0.0000	1.0000
GTR	63	-1200.5130	2795.1382	0.0000	1.0000
JC69I	56	-1223.6337	2797.5893	0.0000	1.0000
JC69G	56	-1224.0462	2798.4144	0.0000	1.0000
JC69IG	57	-1222.8664	2802.3106	0.0000	1.0000
SYMI	61	-1215.8671	2813.3350	0.0000	1.0000
SYMG	61	-1216.3507	2814.3021	0.0000	1.0000
SYMIG	62	-1215.3402	2818.5370	0.0000	1.0000
K2P	56	-1234.7036	2819.7292	0.0000	1.0000
JC69	55	-1241.3252	2826.7167	0.0000	1.0000
SYM	60	-1233.8570	2843.0590	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6050_hg18.phylip = 1 - 521;
 Partition Dummy = 1:261-aggressive-phylip/6050_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6050_hg18.phylip = 1 - 521;
 Partition Dummy = 1:261-aggressive-phylip/6050_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6050_hg18.phylip = 1 - 521;
 Partition Dummy = 1:261-aggressive-phylip/6050_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

