
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:40 02/17/2011

Input data from file "183-phylip/chr11_4735.phylip" (29 taxa, 232 characters)

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-978.0078	2118.8227	0.3957	0.3957
K2PG	57	-983.9256	2119.8511	0.2366	0.6324
K2PI	57	-984.9327	2121.8654	0.0864	0.7188
GTRG	64	-972.0966	2122.0137	0.0803	0.7991
HKYI	60	-979.6798	2122.1665	0.0744	0.8734
HKYIG	61	-977.9930	2122.4801	0.0636	0.9370
K2PIG	58	-983.8712	2123.3032	0.0421	0.9791
GTRIG	65	-972.0889	2125.8646	0.0117	0.9908
GTRI	64	-974.5570	2126.9344	0.0069	0.9977
SYMG	61	-981.6952	2129.8845	0.0016	0.9992
SYMI	61	-982.8145	2132.1232	0.0005	0.9997
SYMIG	62	-981.6597	2133.5443	0.0003	1.0000
HKY	59	-1000.4615	2160.0859	0.0000	1.0000
K2P	56	-1006.9095	2162.2991	0.0000	1.0000
GTR	63	-994.3246	2162.6492	0.0000	1.0000
SYM	60	-1003.8233	2170.4536	0.0000	1.0000
JC69G	56	-1025.0284	2198.5368	0.0000	1.0000
JC69I	56	-1025.8969	2200.2739	0.0000	1.0000
JC69IG	57	-1024.9856	2201.9713	0.0000	1.0000
F81G	59	-1022.2866	2203.7359	0.0000	1.0000
F81I	59	-1023.4601	2206.0831	0.0000	1.0000
F81IG	60	-1022.2717	2207.3505	0.0000	1.0000
JC69	55	-1047.0736	2239.1473	0.0000	1.0000
F81	58	-1044.6343	2244.8292	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-972.0966	2072.1933	0.5871	0.5871
GTRIG	65	-972.0889	2074.1778	0.2176	0.8047
HKYG	60	-978.0078	2076.0157	0.0868	0.8915
GTRI	64	-974.5570	2077.1140	0.0501	0.9417
HKYIG	61	-977.9930	2077.9860	0.0324	0.9741
HKYI	60	-979.6798	2079.3595	0.0163	0.9904
K2PG	57	-983.9256	2081.8511	0.0047	0.9951
K2PIG	58	-983.8712	2083.7425	0.0018	0.9969
K2PI	57	-984.9327	2083.8654	0.0017	0.9986
SYMG	61	-981.6952	2085.3904	0.0008	0.9994
SYMIG	62	-981.6597	2087.3195	0.0003	0.9997
SYMI	61	-982.8145	2087.6291	0.0003	1.0000
GTR	63	-994.3246	2114.6492	0.0000	1.0000
HKY	59	-1000.4615	2118.9231	0.0000	1.0000
K2P	56	-1006.9095	2125.8191	0.0000	1.0000
SYM	60	-1003.8233	2127.6466	0.0000	1.0000
JC69G	56	-1025.0284	2162.0568	0.0000	1.0000
F81G	59	-1022.2866	2162.5731	0.0000	1.0000
JC69I	56	-1025.8969	2163.7939	0.0000	1.0000
JC69IG	57	-1024.9856	2163.9713	0.0000	1.0000
F81IG	60	-1022.2717	2164.5434	0.0000	1.0000
F81I	59	-1023.4601	2164.9203	0.0000	1.0000
JC69	55	-1047.0736	2204.1473	0.0000	1.0000
F81	58	-1044.6343	2205.2685	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-983.9256	2278.3151	0.6379	0.6379
K2PI	57	-984.9327	2280.3295	0.2330	0.8709
HKYG	60	-978.0078	2282.8199	0.0671	0.9380
K2PIG	58	-983.8712	2283.6532	0.0442	0.9822
HKYI	60	-979.6798	2286.1638	0.0126	0.9948
HKYIG	61	-977.9930	2288.2370	0.0045	0.9993
GTRG	64	-972.0966	2292.7845	0.0005	0.9998
SYMG	61	-981.6952	2295.6413	0.0001	0.9999
GTRI	64	-974.5570	2297.7052	0.0000	0.9999
SYMI	61	-982.8145	2297.8800	0.0000	1.0000
GTRIG	65	-972.0889	2298.2158	0.0000	1.0000
SYMIG	62	-981.6597	2301.0172	0.0000	1.0000
K2P	56	-1006.9095	2318.8364	0.0000	1.0000
HKY	59	-1000.4615	2322.2806	0.0000	1.0000
GTR	63	-994.3246	2331.7937	0.0000	1.0000
SYM	60	-1003.8233	2334.4509	0.0000	1.0000
JC69G	56	-1025.0284	2355.0741	0.0000	1.0000
JC69I	56	-1025.8969	2356.8112	0.0000	1.0000
JC69IG	57	-1024.9856	2360.4353	0.0000	1.0000
F81G	59	-1022.2866	2365.9306	0.0000	1.0000
F81I	59	-1023.4601	2368.2778	0.0000	1.0000
F81IG	60	-1022.2717	2371.3477	0.0000	1.0000
JC69	55	-1047.0736	2393.7178	0.0000	1.0000
F81	58	-1044.6343	2405.1793	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 183-phylip/chr11_4735.phylip = 1 - 232;
 Partition Dummy = 1:183-phylip/chr11_4735.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr11_4735.phylip = 1 - 232;
 Partition Dummy = 1:183-phylip/chr11_4735.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr11_4735.phylip = 1 - 232;
 Partition Dummy = 1:183-phylip/chr11_4735.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

