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

Input data from file "683-tomodel-phylip/chr5_7136.phylip" (7 taxa, 682 characters)

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-1175.7782	2382.2771	0.4196	0.4196
HKYG	16	-1175.4036	2383.6252	0.2138	0.6334
HKYI	16	-1175.4258	2383.6696	0.2091	0.8425
HKYIG	17	-1175.4040	2385.7297	0.0747	0.9172
GTR	19	-1173.9284	2387.0049	0.0395	0.9566
GTRG	20	-1173.6200	2388.5107	0.0186	0.9752
GTRI	20	-1173.6304	2388.5316	0.0184	0.9936
GTRIG	21	-1173.6204	2390.6408	0.0064	1.0000
F81	14	-1197.6691	2423.9678	0.0000	1.0000
F81G	15	-1197.4011	2425.5229	0.0000	1.0000
F81I	15	-1197.4060	2425.5328	0.0000	1.0000
F81IG	16	-1197.4011	2427.6202	0.0000	1.0000
K2P	12	-1228.3112	2481.0888	0.0000	1.0000
K2PG	13	-1228.0090	2482.5629	0.0000	1.0000
K2PI	13	-1228.0160	2482.5769	0.0000	1.0000
K2PIG	14	-1228.0091	2484.6479	0.0000	1.0000
SYM	16	-1227.4712	2487.7605	0.0000	1.0000
SYMG	17	-1227.1817	2489.2851	0.0000	1.0000
SYMI	17	-1227.1827	2489.2872	0.0000	1.0000
SYMIG	18	-1227.1813	2491.3942	0.0000	1.0000
JC69	11	-1249.8982	2522.1904	0.0000	1.0000
JC69G	12	-1249.6242	2523.7147	0.0000	1.0000
JC69I	12	-1249.6260	2523.7183	0.0000	1.0000
JC69IG	13	-1249.6238	2525.7926	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-1175.7782	2381.5564	0.3986	0.3986
HKYG	16	-1175.4036	2382.8071	0.2133	0.6118
HKYI	16	-1175.4258	2382.8515	0.2086	0.8204
HKYIG	17	-1175.4040	2384.8080	0.0784	0.8988
GTR	19	-1173.9284	2385.8568	0.0464	0.9452
GTRG	20	-1173.6200	2387.2399	0.0232	0.9684
GTRI	20	-1173.6304	2387.2608	0.0230	0.9915
GTRIG	21	-1173.6204	2389.2408	0.0085	1.0000
F81	14	-1197.6691	2423.3382	0.0000	1.0000
F81G	15	-1197.4011	2424.8022	0.0000	1.0000
F81I	15	-1197.4060	2424.8121	0.0000	1.0000
F81IG	16	-1197.4011	2426.8021	0.0000	1.0000
K2P	12	-1228.3112	2480.6225	0.0000	1.0000
K2PG	13	-1228.0090	2482.0180	0.0000	1.0000
K2PI	13	-1228.0160	2482.0319	0.0000	1.0000
K2PIG	14	-1228.0091	2484.0183	0.0000	1.0000
SYM	16	-1227.4712	2486.9425	0.0000	1.0000
SYMG	17	-1227.1817	2488.3634	0.0000	1.0000
SYMI	17	-1227.1827	2488.3655	0.0000	1.0000
SYMIG	18	-1227.1813	2490.3626	0.0000	1.0000
JC69	11	-1249.8982	2521.7964	0.0000	1.0000
JC69G	12	-1249.6242	2523.2484	0.0000	1.0000
JC69I	12	-1249.6260	2523.2519	0.0000	1.0000
JC69IG	13	-1249.6238	2525.2476	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1175.7782	2449.4318	0.8990	0.8990
HKYG	16	-1175.4036	2455.2076	0.0501	0.9491
HKYI	16	-1175.4258	2455.2520	0.0490	0.9981
HKYIG	17	-1175.4040	2461.7335	0.0019	1.0000
GTR	19	-1173.9284	2471.8324	0.0000	1.0000
GTRG	20	-1173.6200	2477.7405	0.0000	1.0000
GTRI	20	-1173.6304	2477.7614	0.0000	1.0000
GTRIG	21	-1173.6204	2484.2664	0.0000	1.0000
F81	14	-1197.6691	2486.6886	0.0000	1.0000
F81G	15	-1197.4011	2492.6776	0.0000	1.0000
F81I	15	-1197.4060	2492.6875	0.0000	1.0000
F81IG	16	-1197.4011	2499.2026	0.0000	1.0000
K2P	12	-1228.3112	2534.9228	0.0000	1.0000
K2PG	13	-1228.0090	2540.8434	0.0000	1.0000
K2PI	13	-1228.0160	2540.8573	0.0000	1.0000
K2PIG	14	-1228.0091	2547.3687	0.0000	1.0000
SYM	16	-1227.4712	2559.3430	0.0000	1.0000
SYMG	17	-1227.1817	2565.2889	0.0000	1.0000
SYMI	17	-1227.1827	2565.2910	0.0000	1.0000
JC69	11	-1249.8982	2571.5717	0.0000	1.0000
SYMIG	18	-1227.1813	2571.8131	0.0000	1.0000
JC69G	12	-1249.6242	2577.5487	0.0000	1.0000
JC69I	12	-1249.6260	2577.5523	0.0000	1.0000
JC69IG	13	-1249.6238	2584.0730	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_7136.phylip = 1 - 682;
 Partition Dummy = 1:683-tomodel-phylip/chr5_7136.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


[Mrbayes block for the best AICc model (HKY)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_7136.phylip = 1 - 682;
 Partition Dummy = 1:683-tomodel-phylip/chr5_7136.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_7136.phylip = 1 - 682;
 Partition Dummy = 1:683-tomodel-phylip/chr5_7136.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

