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

Input data from file "183-phylip/chr2_7982.phylip" (29 taxa, 367 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	-1063.0396	2270.0007	0.8650	0.8650
HKYI	60	-1065.6446	2275.2108	0.0639	0.9289
HKYIG	61	-1064.8740	2276.5480	0.0328	0.9616
GTRG	64	-1060.7712	2277.0920	0.0250	0.9866
K2PG	57	-1072.6240	2280.6460	0.0042	0.9908
K2PI	57	-1072.6668	2280.7317	0.0040	0.9949
K2PIG	58	-1071.9159	2282.0526	0.0021	0.9969
GTRI	64	-1063.4695	2282.4887	0.0017	0.9986
GTRIG	65	-1062.6364	2283.7777	0.0009	0.9995
SYMG	61	-1069.9436	2286.6872	0.0002	0.9997
SYMI	61	-1070.0070	2286.8141	0.0002	0.9999
SYMIG	62	-1069.2404	2288.1781	0.0001	1.0000
HKY	59	-1087.3086	2315.6792	0.0000	1.0000
GTR	63	-1085.4948	2323.6035	0.0000	1.0000
F81G	59	-1092.6102	2326.2823	0.0000	1.0000
F81I	59	-1092.6469	2326.3557	0.0000	1.0000
F81IG	60	-1092.0683	2328.0582	0.0000	1.0000
K2P	56	-1098.4662	2329.5259	0.0000	1.0000
JC69G	56	-1099.4668	2331.5271	0.0000	1.0000
JC69I	56	-1099.5044	2331.6023	0.0000	1.0000
JC69IG	57	-1098.8255	2333.0491	0.0000	1.0000
SYM	60	-1095.9101	2335.7417	0.0000	1.0000
F81	58	-1118.2564	2374.7336	0.0000	1.0000
JC69	55	-1124.9676	2379.7423	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1063.0396	2246.0792	0.7483	0.7483
GTRG	64	-1060.7712	2249.5423	0.1324	0.8807
HKYI	60	-1065.6446	2251.2892	0.0553	0.9360
HKYIG	61	-1064.8740	2251.7480	0.0440	0.9800
GTRI	64	-1063.4695	2254.9391	0.0089	0.9889
GTRIG	65	-1062.6364	2255.2728	0.0075	0.9965
K2PG	57	-1072.6240	2259.2479	0.0010	0.9975
K2PI	57	-1072.6668	2259.3337	0.0010	0.9985
K2PIG	58	-1071.9159	2259.8318	0.0008	0.9993
SYMG	61	-1069.9436	2261.8872	0.0003	0.9995
SYMI	61	-1070.0070	2262.0141	0.0003	0.9998
SYMIG	62	-1069.2404	2262.4807	0.0002	1.0000
HKY	59	-1087.3086	2292.6173	0.0000	1.0000
GTR	63	-1085.4948	2296.9896	0.0000	1.0000
F81G	59	-1092.6102	2303.2204	0.0000	1.0000
F81I	59	-1092.6469	2303.2939	0.0000	1.0000
F81IG	60	-1092.0683	2304.1367	0.0000	1.0000
K2P	56	-1098.4662	2308.9324	0.0000	1.0000
JC69G	56	-1099.4668	2310.9336	0.0000	1.0000
JC69I	56	-1099.5044	2311.0088	0.0000	1.0000
JC69IG	57	-1098.8255	2311.6510	0.0000	1.0000
SYM	60	-1095.9101	2311.8201	0.0000	1.0000
F81	58	-1118.2564	2352.5129	0.0000	1.0000
JC69	55	-1124.9676	2359.9352	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1063.0396	2480.4009	0.4806	0.4806
K2PG	57	-1072.6240	2481.8536	0.2325	0.7131
K2PI	57	-1072.6668	2481.9393	0.2227	0.9358
HKYI	60	-1065.6446	2485.6109	0.0355	0.9713
K2PIG	58	-1071.9159	2486.3428	0.0246	0.9959
HKYIG	61	-1064.8740	2489.9750	0.0040	0.9999
GTRG	64	-1060.7712	2499.4855	0.0000	0.9999
SYMG	61	-1069.9436	2500.1142	0.0000	1.0000
SYMI	61	-1070.0070	2500.2412	0.0000	1.0000
SYMIG	62	-1069.2404	2504.6131	0.0000	1.0000
GTRI	64	-1063.4695	2504.8822	0.0000	1.0000
GTRIG	65	-1062.6364	2509.1213	0.0000	1.0000
HKY	59	-1087.3086	2523.0336	0.0000	1.0000
K2P	56	-1098.4662	2527.6326	0.0000	1.0000
JC69G	56	-1099.4668	2529.6338	0.0000	1.0000
JC69I	56	-1099.5044	2529.7090	0.0000	1.0000
F81G	59	-1092.6102	2533.6368	0.0000	1.0000
F81I	59	-1092.6469	2533.7102	0.0000	1.0000
JC69IG	57	-1098.8255	2534.2566	0.0000	1.0000
F81IG	60	-1092.0683	2538.4584	0.0000	1.0000
GTR	63	-1085.4948	2543.0274	0.0000	1.0000
SYM	60	-1095.9101	2546.1418	0.0000	1.0000
JC69	55	-1124.9676	2574.7301	0.0000	1.0000
F81	58	-1118.2564	2579.0238	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/chr2_7982.phylip = 1 - 367;
 Partition Dummy = 1:183-phylip/chr2_7982.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/chr2_7982.phylip = 1 - 367;
 Partition Dummy = 1:183-phylip/chr2_7982.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/chr2_7982.phylip = 1 - 367;
 Partition Dummy = 1:183-phylip/chr2_7982.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

