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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1381.4404	2919.1801	0.6724	0.6724
GTRIG	65	-1381.0034	2921.2901	0.2341	0.9065
HKYIG	61	-1388.5106	2924.4892	0.0473	0.9538
HKYG	60	-1390.7773	2926.1184	0.0209	0.9747
SYMG	61	-1389.5805	2926.6290	0.0162	0.9909
SYMIG	62	-1389.3858	2929.1636	0.0046	0.9955
K2PIG	58	-1395.8948	2930.6030	0.0022	0.9977
K2PG	57	-1397.3248	2930.6164	0.0022	1.0000
GTRI	64	-1390.9940	2938.2874	0.0000	1.0000
HKYI	60	-1400.0081	2944.5800	0.0000	1.0000
SYMI	61	-1401.6468	2950.7616	0.0000	1.0000
K2PI	57	-1407.6743	2951.3153	0.0000	1.0000
F81G	59	-1419.9574	2981.5938	0.0000	1.0000
F81IG	60	-1419.8464	2984.2565	0.0000	1.0000
JC69G	56	-1426.3443	2985.8278	0.0000	1.0000
JC69IG	57	-1425.5371	2987.0410	0.0000	1.0000
F81I	59	-1430.5901	3002.8591	0.0000	1.0000
JC69I	56	-1436.5479	3006.2349	0.0000	1.0000
GTR	63	-1430.7775	3014.8906	0.0000	1.0000
HKY	59	-1439.4686	3020.6160	0.0000	1.0000
K2P	56	-1444.7141	3022.5673	0.0000	1.0000
SYM	60	-1440.6822	3025.9281	0.0000	1.0000
F81	58	-1467.9774	3074.7681	0.0000	1.0000
JC69	55	-1474.7295	3079.7890	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1381.4404	2890.8808	0.6249	0.6249
GTRIG	65	-1381.0034	2892.0068	0.3559	0.9808
HKYIG	61	-1388.5106	2899.0212	0.0107	0.9914
SYMG	61	-1389.5805	2901.1610	0.0037	0.9951
HKYG	60	-1390.7773	2901.5546	0.0030	0.9981
SYMIG	62	-1389.3858	2902.7717	0.0016	0.9997
K2PIG	58	-1395.8948	2907.7896	0.0001	0.9999
K2PG	57	-1397.3248	2908.6496	0.0001	1.0000
GTRI	64	-1390.9940	2909.9881	0.0000	1.0000
HKYI	60	-1400.0081	2920.0162	0.0000	1.0000
SYMI	61	-1401.6468	2925.2936	0.0000	1.0000
K2PI	57	-1407.6743	2929.3486	0.0000	1.0000
F81G	59	-1419.9574	2957.9148	0.0000	1.0000
F81IG	60	-1419.8464	2959.6927	0.0000	1.0000
JC69G	56	-1426.3443	2964.6887	0.0000	1.0000
JC69IG	57	-1425.5371	2965.0743	0.0000	1.0000
F81I	59	-1430.5901	2979.1801	0.0000	1.0000
JC69I	56	-1436.5479	2985.0958	0.0000	1.0000
GTR	63	-1430.7775	2987.5550	0.0000	1.0000
HKY	59	-1439.4686	2996.9371	0.0000	1.0000
SYM	60	-1440.6822	3001.3644	0.0000	1.0000
K2P	56	-1444.7141	3001.4282	0.0000	1.0000
F81	58	-1467.9774	3051.9548	0.0000	1.0000
JC69	55	-1474.7295	3059.4589	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1397.3248	3129.9990	0.7122	0.7122
K2PIG	58	-1395.8948	3133.0223	0.1571	0.8692
HKYG	60	-1390.7773	3134.5539	0.0730	0.9423
HKYIG	61	-1388.5106	3135.9039	0.0372	0.9794
SYMG	61	-1389.5805	3138.0437	0.0128	0.9922
GTRG	64	-1381.4404	3139.4134	0.0064	0.9986
SYMIG	62	-1389.3858	3143.5377	0.0008	0.9994
GTRIG	65	-1381.0034	3144.4227	0.0005	1.0000
K2PI	57	-1407.6743	3150.6979	0.0000	1.0000
HKYI	60	-1400.0081	3153.0156	0.0000	1.0000
GTRI	64	-1390.9940	3158.5207	0.0000	1.0000
SYMI	61	-1401.6468	3162.1762	0.0000	1.0000
JC69G	56	-1426.3443	3182.1547	0.0000	1.0000
JC69IG	57	-1425.5371	3186.4236	0.0000	1.0000
F81G	59	-1419.9574	3187.0309	0.0000	1.0000
F81IG	60	-1419.8464	3192.6920	0.0000	1.0000
JC69I	56	-1436.5479	3202.5618	0.0000	1.0000
F81I	59	-1430.5901	3208.2962	0.0000	1.0000
K2P	56	-1444.7141	3218.8942	0.0000	1.0000
HKY	59	-1439.4686	3226.0531	0.0000	1.0000
GTR	63	-1430.7775	3232.2043	0.0000	1.0000
SYM	60	-1440.6822	3234.3637	0.0000	1.0000
JC69	55	-1474.7295	3273.0417	0.0000	1.0000
F81	58	-1467.9774	3277.1875	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/5201_hg18.phylip = 1 - 359;
 Partition Dummy = 1:261-aggressive-phylip/5201_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/5201_hg18.phylip = 1 - 359;
 Partition Dummy = 1:261-aggressive-phylip/5201_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5201_hg18.phylip = 1 - 359;
 Partition Dummy = 1:261-aggressive-phylip/5201_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

