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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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	-1315.1724	2776.2105	0.7028	0.7028
HKYIG	61	-1314.9878	2778.7983	0.1927	0.8956
GTRG	64	-1311.4448	2780.7104	0.0741	0.9696
GTRIG	65	-1310.8739	2782.6111	0.0286	0.9983
HKYI	60	-1321.2214	2788.3085	0.0017	0.9999
GTRI	64	-1318.6988	2795.2185	0.0001	1.0000
K2PG	57	-1337.2433	2811.6055	0.0000	1.0000
K2PIG	58	-1337.9376	2815.8893	0.0000	1.0000
SYMG	61	-1336.8575	2822.5378	0.0000	1.0000
K2PI	57	-1343.7200	2824.5590	0.0000	1.0000
SYMIG	62	-1337.5822	2826.9652	0.0000	1.0000
SYMI	61	-1343.5154	2835.8535	0.0000	1.0000
HKY	59	-1359.4517	2861.8329	0.0000	1.0000
GTR	63	-1356.1005	2867.0010	0.0000	1.0000
K2P	56	-1385.7912	2905.8263	0.0000	1.0000
F81G	59	-1383.4671	2909.8638	0.0000	1.0000
SYM	60	-1383.3985	2912.6627	0.0000	1.0000
F81IG	60	-1383.9296	2913.7248	0.0000	1.0000
F81I	59	-1391.0865	2925.1025	0.0000	1.0000
JC69IG	57	-1396.7136	2930.5461	0.0000	1.0000
JC69G	56	-1398.8196	2931.8832	0.0000	1.0000
JC69I	56	-1403.4888	2941.2215	0.0000	1.0000
F81	58	-1430.7604	3001.5349	0.0000	1.0000
JC69	55	-1444.4677	3020.3242	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1315.1724	2750.3448	0.3700	0.3700
GTRG	64	-1311.4448	2750.8896	0.2818	0.6518
GTRIG	65	-1310.8739	2751.7478	0.1835	0.8352
HKYIG	61	-1314.9878	2751.9756	0.1637	0.9989
HKYI	60	-1321.2214	2762.4428	0.0009	0.9998
GTRI	64	-1318.6988	2765.3977	0.0002	1.0000
K2PG	57	-1337.2433	2788.4866	0.0000	1.0000
K2PIG	58	-1337.9376	2791.8752	0.0000	1.0000
SYMG	61	-1336.8575	2795.7151	0.0000	1.0000
SYMIG	62	-1337.5822	2799.1645	0.0000	1.0000
K2PI	57	-1343.7200	2801.4401	0.0000	1.0000
SYMI	61	-1343.5154	2809.0308	0.0000	1.0000
HKY	59	-1359.4517	2836.9034	0.0000	1.0000
GTR	63	-1356.1005	2838.2010	0.0000	1.0000
K2P	56	-1385.7912	2883.5824	0.0000	1.0000
F81G	59	-1383.4671	2884.9342	0.0000	1.0000
SYM	60	-1383.3985	2886.7970	0.0000	1.0000
F81IG	60	-1383.9296	2887.8591	0.0000	1.0000
F81I	59	-1391.0865	2900.1729	0.0000	1.0000
JC69IG	57	-1396.7136	2907.4272	0.0000	1.0000
JC69G	56	-1398.8196	2909.6393	0.0000	1.0000
JC69I	56	-1403.4888	2918.9776	0.0000	1.0000
F81	58	-1430.7604	2977.5208	0.0000	1.0000
JC69	55	-1444.4677	2998.9353	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1315.1724	2980.7833	0.9367	0.9367
HKYIG	61	-1314.9878	2986.2547	0.0607	0.9974
HKYI	60	-1321.2214	2992.8813	0.0022	0.9996
GTRG	64	-1311.4448	2996.6907	0.0003	1.0000
GTRIG	65	-1310.8739	3001.3895	0.0000	1.0000
K2PG	57	-1337.2433	3007.4032	0.0000	1.0000
GTRI	64	-1318.6988	3011.1988	0.0000	1.0000
K2PIG	58	-1337.9376	3014.6325	0.0000	1.0000
K2PI	57	-1343.7200	3020.3567	0.0000	1.0000
SYMG	61	-1336.8575	3029.9942	0.0000	1.0000
SYMIG	62	-1337.5822	3037.2842	0.0000	1.0000
SYMI	61	-1343.5154	3043.3100	0.0000	1.0000
HKY	59	-1359.4517	3063.5012	0.0000	1.0000
GTR	63	-1356.1005	3080.1614	0.0000	1.0000
K2P	56	-1385.7912	3098.6583	0.0000	1.0000
F81G	59	-1383.4671	3111.5321	0.0000	1.0000
SYM	60	-1383.3985	3117.2355	0.0000	1.0000
F81IG	60	-1383.9296	3118.2976	0.0000	1.0000
JC69G	56	-1398.8196	3124.7152	0.0000	1.0000
JC69IG	57	-1396.7136	3126.3438	0.0000	1.0000
F81I	59	-1391.0865	3126.7708	0.0000	1.0000
JC69I	56	-1403.4888	3134.0536	0.0000	1.0000
F81	58	-1430.7604	3200.2781	0.0000	1.0000
JC69	55	-1444.4677	3210.1706	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr13_3036.phylip = 1 - 344;
 Partition Dummy = 1:183-phylip/chr13_3036.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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr13_3036.phylip = 1 - 344;
 Partition Dummy = 1:183-phylip/chr13_3036.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 (HKYG)]

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


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

