
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:38 04/04/2011

Input data from file "591-tomodel-phylip/chr2_27968.phylip" (7 taxa, 494 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	16	-1157.3765	2347.8935	0.3510	0.3510
HKYI	16	-1157.7585	2348.6575	0.2395	0.5905
HKYIG	17	-1157.2001	2349.6860	0.1432	0.7338
GTRG	20	-1154.0694	2349.9146	0.1278	0.8615
GTRI	20	-1154.4826	2350.7410	0.0845	0.9460
GTRIG	21	-1153.8633	2351.6843	0.0527	0.9988
HKY	15	-1164.7407	2360.4856	0.0006	0.9994
GTR	19	-1161.4793	2362.5619	0.0002	0.9997
F81G	15	-1166.5246	2364.0535	0.0001	0.9998
F81I	15	-1166.9172	2364.8387	0.0001	0.9998
K2PG	13	-1169.3836	2365.5256	0.0001	0.9999
F81IG	16	-1166.3971	2365.9347	0.0000	0.9999
K2PI	13	-1169.8002	2366.3587	0.0000	1.0000
K2PIG	14	-1169.2464	2367.3697	0.0000	1.0000
SYMG	17	-1167.3325	2369.9507	0.0000	1.0000
SYMI	17	-1167.7839	2370.8536	0.0000	1.0000
SYMIG	18	-1167.1985	2371.8370	0.0000	1.0000
F81	14	-1173.1629	2375.2026	0.0000	1.0000
K2P	12	-1176.0970	2376.8426	0.0000	1.0000
JC69G	12	-1177.3573	2379.3632	0.0000	1.0000
JC69I	12	-1177.7692	2380.1870	0.0000	1.0000
SYM	16	-1173.8014	2380.7433	0.0000	1.0000
JC69IG	13	-1177.2479	2381.2542	0.0000	1.0000
JC69	11	-1183.6908	2389.9292	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1157.3765	2346.7530	0.3143	0.3143
HKYI	16	-1157.7585	2347.5171	0.2145	0.5288
GTRG	20	-1154.0694	2348.1387	0.1572	0.6860
HKYIG	17	-1157.2001	2348.4003	0.1379	0.8239
GTRI	20	-1154.4826	2348.9651	0.1040	0.9279
GTRIG	21	-1153.8633	2349.7266	0.0711	0.9989
HKY	15	-1164.7407	2359.4814	0.0005	0.9995
GTR	19	-1161.4793	2360.9585	0.0003	0.9997
F81G	15	-1166.5246	2363.0493	0.0001	0.9998
F81I	15	-1166.9172	2363.8345	0.0001	0.9999
K2PG	13	-1169.3836	2364.7673	0.0000	0.9999
F81IG	16	-1166.3971	2364.7942	0.0000	0.9999
K2PI	13	-1169.8002	2365.6004	0.0000	1.0000
K2PIG	14	-1169.2464	2366.4929	0.0000	1.0000
SYMG	17	-1167.3325	2368.6649	0.0000	1.0000
SYMI	17	-1167.7839	2369.5679	0.0000	1.0000
SYMIG	18	-1167.1985	2370.3970	0.0000	1.0000
F81	14	-1173.1629	2374.3258	0.0000	1.0000
K2P	12	-1176.0970	2376.1939	0.0000	1.0000
JC69G	12	-1177.3573	2378.7145	0.0000	1.0000
JC69I	12	-1177.7692	2379.5384	0.0000	1.0000
SYM	16	-1173.8014	2379.6028	0.0000	1.0000
JC69IG	13	-1177.2479	2380.4958	0.0000	1.0000
JC69	11	-1183.6908	2389.3815	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1157.3765	2413.9936	0.5342	0.5342
HKYI	16	-1157.7585	2414.7576	0.3646	0.8988
K2PG	13	-1169.3836	2419.4003	0.0358	0.9346
HKYIG	17	-1157.2001	2419.8434	0.0287	0.9633
K2PI	13	-1169.8002	2420.2333	0.0236	0.9869
HKY	15	-1164.7407	2422.5195	0.0075	0.9944
K2PIG	14	-1169.2464	2425.3284	0.0018	0.9962
F81G	15	-1166.5246	2426.0873	0.0013	0.9975
K2P	12	-1176.0970	2426.6243	0.0010	0.9985
F81I	15	-1166.9172	2426.8725	0.0009	0.9993
JC69G	12	-1177.3573	2429.1449	0.0003	0.9996
JC69I	12	-1177.7692	2429.9688	0.0002	0.9998
F81IG	16	-1166.3971	2432.0348	0.0001	0.9998
GTRG	20	-1154.0694	2432.1895	0.0001	0.9999
GTRI	20	-1154.4826	2433.0158	0.0000	0.9999
F81	14	-1173.1629	2433.1613	0.0000	1.0000
JC69IG	13	-1177.2479	2435.1288	0.0000	1.0000
JC69	11	-1183.6908	2435.6094	0.0000	1.0000
GTRIG	21	-1153.8633	2437.9799	0.0000	1.0000
SYMG	17	-1167.3325	2440.1080	0.0000	1.0000
GTR	19	-1161.4793	2440.8067	0.0000	1.0000
SYMI	17	-1167.7839	2441.0110	0.0000	1.0000
SYMIG	18	-1167.1985	2446.0426	0.0000	1.0000
SYM	16	-1173.8014	2446.8434	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 591-tomodel-phylip/chr2_27968.phylip = 1 - 494;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27968.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 591-tomodel-phylip/chr2_27968.phylip = 1 - 494;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27968.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 591-tomodel-phylip/chr2_27968.phylip = 1 - 494;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27968.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

