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

Input data from file "261-aggressive-phylip/11690_hg18.phylip" (29 taxa, 312 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	-1009.1337	2167.4307	0.4889	0.4889
HKYG	60	-1009.3484	2167.8601	0.3945	0.8834
HKYIG	61	-1009.1254	2170.5069	0.1050	0.9884
GTRG	64	-1007.5225	2176.7291	0.0047	0.9931
GTRI	64	-1007.5489	2176.7821	0.0046	0.9977
K2PG	57	-1019.7921	2179.6158	0.0011	0.9988
GTRIG	65	-1007.5308	2179.9396	0.0009	0.9997
K2PI	57	-1021.7130	2183.4575	0.0002	0.9999
HKY	59	-1019.3442	2184.7837	0.0001	1.0000
K2PIG	58	-1021.7684	2186.5882	0.0000	1.0000
SYMG	61	-1018.8797	2190.0154	0.0000	1.0000
SYMI	61	-1021.1527	2194.5613	0.0000	1.0000
GTR	63	-1018.2361	2194.9884	0.0000	1.0000
K2P	56	-1029.0472	2195.1296	0.0000	1.0000
SYMIG	62	-1020.9240	2197.2216	0.0000	1.0000
SYM	60	-1028.0024	2205.1682	0.0000	1.0000
F81I	59	-1030.3716	2206.8385	0.0000	1.0000
F81G	59	-1030.6326	2207.3603	0.0000	1.0000
F81IG	60	-1030.3740	2209.9113	0.0000	1.0000
JC69G	56	-1037.7938	2212.6229	0.0000	1.0000
JC69I	56	-1039.7613	2216.5579	0.0000	1.0000
JC69IG	57	-1039.7533	2219.5381	0.0000	1.0000
F81	58	-1038.4598	2219.9709	0.0000	1.0000
JC69	55	-1046.9342	2227.9308	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1009.1337	2138.2674	0.4179	0.4179
HKYG	60	-1009.3484	2138.6967	0.3372	0.7550
HKYIG	61	-1009.1254	2140.2509	0.1550	0.9100
GTRG	64	-1007.5225	2143.0449	0.0383	0.9484
GTRI	64	-1007.5489	2143.0978	0.0373	0.9857
GTRIG	65	-1007.5308	2145.0616	0.0140	0.9997
K2PG	57	-1019.7921	2153.5843	0.0002	0.9999
HKY	59	-1019.3442	2156.6885	0.0000	0.9999
K2PI	57	-1021.7130	2157.4260	0.0000	1.0000
K2PIG	58	-1021.7684	2159.5369	0.0000	1.0000
SYMG	61	-1018.8797	2159.7594	0.0000	1.0000
GTR	63	-1018.2361	2162.4723	0.0000	1.0000
SYMI	61	-1021.1527	2164.3053	0.0000	1.0000
SYMIG	62	-1020.9240	2165.8481	0.0000	1.0000
K2P	56	-1029.0472	2170.0943	0.0000	1.0000
SYM	60	-1028.0024	2176.0048	0.0000	1.0000
F81I	59	-1030.3716	2178.7433	0.0000	1.0000
F81G	59	-1030.6326	2179.2651	0.0000	1.0000
F81IG	60	-1030.3740	2180.7480	0.0000	1.0000
JC69G	56	-1037.7938	2187.5876	0.0000	1.0000
JC69I	56	-1039.7613	2191.5226	0.0000	1.0000
F81	58	-1038.4598	2192.9195	0.0000	1.0000
JC69IG	57	-1039.7533	2193.5066	0.0000	1.0000
JC69	55	-1046.9342	2203.8683	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1009.1337	2362.8476	0.4964	0.4964
HKYG	60	-1009.3484	2363.2769	0.4005	0.8970
K2PG	57	-1019.7921	2366.9354	0.0643	0.9613
HKYIG	61	-1009.1254	2368.5741	0.0283	0.9896
K2PI	57	-1021.7130	2370.7771	0.0094	0.9990
K2PIG	58	-1021.7684	2376.6310	0.0005	0.9995
HKY	59	-1019.3442	2377.5257	0.0003	0.9998
K2P	56	-1029.0472	2379.7025	0.0001	0.9999
GTRG	64	-1007.5225	2382.5971	0.0000	1.0000
GTRI	64	-1007.5489	2382.6500	0.0000	1.0000
SYMG	61	-1018.8797	2388.0826	0.0000	1.0000
GTRIG	65	-1007.5308	2388.3568	0.0000	1.0000
SYMI	61	-1021.1527	2392.6285	0.0000	1.0000
JC69G	56	-1037.7938	2397.1957	0.0000	1.0000
SYMIG	62	-1020.9240	2397.9143	0.0000	1.0000
GTR	63	-1018.2361	2398.2815	0.0000	1.0000
F81I	59	-1030.3716	2399.5805	0.0000	1.0000
F81G	59	-1030.6326	2400.1023	0.0000	1.0000
SYM	60	-1028.0024	2400.5850	0.0000	1.0000
JC69I	56	-1039.7613	2401.1308	0.0000	1.0000
F81IG	60	-1030.3740	2405.3282	0.0000	1.0000
JC69IG	57	-1039.7533	2406.8578	0.0000	1.0000
JC69	55	-1046.9342	2409.7335	0.0000	1.0000
F81	58	-1038.4598	2410.0137	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/11690_hg18.phylip = 1 - 312;
 Partition Dummy = 1:261-aggressive-phylip/11690_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/11690_hg18.phylip = 1 - 312;
 Partition Dummy = 1:261-aggressive-phylip/11690_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/11690_hg18.phylip = 1 - 312;
 Partition Dummy = 1:261-aggressive-phylip/11690_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

