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

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

Minimum AIC  model: HKYG
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	-1058.2188	2267.4545	0.3237	0.3237
K2PG	57	-1062.9331	2267.5314	0.3115	0.6352
HKYI	60	-1059.3221	2269.6612	0.1074	0.7426
K2PI	57	-1064.2770	2270.2192	0.0812	0.8238
K2PIG	58	-1062.9372	2270.6308	0.0661	0.8900
HKYIG	61	-1058.2264	2270.6401	0.0658	0.9558
GTRG	64	-1054.6620	2273.1861	0.0184	0.9742
SYMG	61	-1060.0425	2274.2723	0.0107	0.9849
GTRI	64	-1055.7361	2275.3343	0.0063	0.9912
GTRIG	65	-1054.6654	2276.4736	0.0036	0.9948
SYMI	61	-1061.3719	2276.9310	0.0028	0.9976
SYMIG	62	-1060.0462	2277.4770	0.0022	0.9998
HKY	59	-1067.1938	2282.2610	0.0002	1.0000
K2P	56	-1074.0245	2286.6490	0.0000	1.0000
GTR	63	-1064.4898	2289.5890	0.0000	1.0000
SYM	60	-1069.3663	2289.7496	0.0000	1.0000
JC69G	56	-1082.5123	2303.6246	0.0000	1.0000
F81G	59	-1078.2584	2304.3902	0.0000	1.0000
JC69IG	57	-1082.5131	2306.6915	0.0000	1.0000
JC69I	56	-1084.1212	2306.8423	0.0000	1.0000
F81I	59	-1079.4868	2306.8470	0.0000	1.0000
F81IG	60	-1079.5116	2310.0402	0.0000	1.0000
F81	58	-1086.6367	2318.0297	0.0000	1.0000
JC69	55	-1091.9022	2319.3645	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1058.2188	2236.4376	0.3141	0.3141
GTRG	64	-1054.6620	2237.3240	0.2016	0.5157
HKYIG	61	-1058.2264	2238.4528	0.1147	0.6304
HKYI	60	-1059.3221	2238.6442	0.1042	0.7346
GTRIG	65	-1054.6654	2239.3308	0.0739	0.8085
GTRI	64	-1055.7361	2239.4723	0.0689	0.8774
K2PG	57	-1062.9331	2239.8661	0.0566	0.9340
K2PIG	58	-1062.9372	2241.8745	0.0207	0.9547
SYMG	61	-1060.0425	2242.0850	0.0187	0.9733
K2PI	57	-1064.2770	2242.5539	0.0148	0.9881
SYMIG	62	-1060.0462	2244.0923	0.0068	0.9949
SYMI	61	-1061.3719	2244.7437	0.0049	0.9999
HKY	59	-1067.1938	2252.3876	0.0001	1.0000
GTR	63	-1064.4898	2254.9796	0.0000	1.0000
SYM	60	-1069.3663	2258.7326	0.0000	1.0000
K2P	56	-1074.0245	2260.0490	0.0000	1.0000
F81G	59	-1078.2584	2274.5168	0.0000	1.0000
F81I	59	-1079.4868	2276.9736	0.0000	1.0000
JC69G	56	-1082.5123	2277.0246	0.0000	1.0000
F81IG	60	-1079.5116	2279.0233	0.0000	1.0000
JC69IG	57	-1082.5131	2279.0263	0.0000	1.0000
JC69I	56	-1084.1212	2280.2423	0.0000	1.0000
F81	58	-1086.6367	2289.2734	0.0000	1.0000
JC69	55	-1091.9022	2293.8043	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1062.9331	2450.4088	0.7410	0.7410
K2PI	57	-1064.2770	2453.0967	0.1933	0.9343
K2PIG	58	-1062.9372	2456.1110	0.0428	0.9771
HKYG	60	-1058.2188	2458.0615	0.0161	0.9933
HKYI	60	-1059.3221	2460.2681	0.0054	0.9986
HKYIG	61	-1058.2264	2463.7705	0.0009	0.9996
K2P	56	-1074.0245	2466.8980	0.0002	0.9998
SYMG	61	-1060.0425	2467.4027	0.0002	0.9999
SYMI	61	-1061.3719	2470.0614	0.0000	0.9999
HKY	59	-1067.1938	2470.3178	0.0000	1.0000
SYMIG	62	-1060.0462	2473.1037	0.0000	1.0000
GTRG	64	-1054.6620	2473.7229	0.0000	1.0000
GTRI	64	-1055.7361	2475.8711	0.0000	1.0000
GTRIG	65	-1054.6654	2479.4234	0.0000	1.0000
SYM	60	-1069.3663	2480.3566	0.0000	1.0000
JC69G	56	-1082.5123	2483.8736	0.0000	1.0000
JC69I	56	-1084.1212	2487.0913	0.0000	1.0000
GTR	63	-1064.4898	2487.6847	0.0000	1.0000
JC69IG	57	-1082.5131	2489.5690	0.0000	1.0000
F81G	59	-1078.2584	2492.4470	0.0000	1.0000
F81I	59	-1079.4868	2494.9038	0.0000	1.0000
JC69	55	-1091.9022	2496.9596	0.0000	1.0000
F81IG	60	-1079.5116	2500.6472	0.0000	1.0000
F81	58	-1086.6367	2503.5098	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/chr3_9696.phylip = 1 - 297;
 Partition Dummy = 1:183-phylip/chr3_9696.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/chr3_9696.phylip = 1 - 297;
 Partition Dummy = 1:183-phylip/chr3_9696.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/chr3_9696.phylip = 1 - 297;
 Partition Dummy = 1:183-phylip/chr3_9696.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

