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

Input data from file "591-tomodel-phylip/chr12_5730.phylip" (7 taxa, 444 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	-1123.4677	2280.2094	0.4595	0.4595
HKYI	16	-1124.1644	2281.6027	0.2290	0.6885
HKYIG	17	-1123.0865	2281.6096	0.2282	0.9167
GTRG	20	-1121.4909	2284.9675	0.0426	0.9593
GTRI	20	-1122.2274	2286.4407	0.0204	0.9796
GTRIG	21	-1121.1270	2286.4436	0.0204	1.0000
HKY	15	-1138.6546	2308.4307	0.0000	1.0000
GTR	19	-1136.4752	2312.7429	0.0000	1.0000
F81G	15	-1145.5759	2322.2734	0.0000	1.0000
F81I	15	-1146.1919	2323.5054	0.0000	1.0000
F81IG	16	-1145.2848	2323.8435	0.0000	1.0000
K2PG	13	-1152.5323	2331.9112	0.0000	1.0000
K2PI	13	-1153.0759	2332.9983	0.0000	1.0000
K2PIG	14	-1152.2326	2333.4442	0.0000	1.0000
SYMG	17	-1151.9316	2339.2999	0.0000	1.0000
SYMI	17	-1152.4411	2340.3189	0.0000	1.0000
SYMIG	18	-1151.6328	2340.8751	0.0000	1.0000
F81	14	-1158.8413	2346.6616	0.0000	1.0000
K2P	12	-1166.9360	2358.5958	0.0000	1.0000
JC69G	12	-1170.5802	2365.8843	0.0000	1.0000
SYM	16	-1166.5022	2366.2784	0.0000	1.0000
JC69I	12	-1171.0910	2366.9059	0.0000	1.0000
JC69IG	13	-1170.3494	2367.5453	0.0000	1.0000
JC69	11	-1183.7358	2390.0827	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1123.4677	2278.9354	0.4343	0.4343
HKYIG	17	-1123.0865	2280.1729	0.2339	0.6683
HKYI	16	-1124.1644	2280.3287	0.2164	0.8847
GTRG	20	-1121.4909	2282.9817	0.0574	0.9421
GTRIG	21	-1121.1270	2284.2540	0.0304	0.9725
GTRI	20	-1122.2274	2284.4549	0.0275	1.0000
HKY	15	-1138.6546	2307.3092	0.0000	1.0000
GTR	19	-1136.4752	2310.9504	0.0000	1.0000
F81G	15	-1145.5759	2321.1519	0.0000	1.0000
F81I	15	-1146.1919	2322.3839	0.0000	1.0000
F81IG	16	-1145.2848	2322.5695	0.0000	1.0000
K2PG	13	-1152.5323	2331.0647	0.0000	1.0000
K2PI	13	-1153.0759	2332.1518	0.0000	1.0000
K2PIG	14	-1152.2326	2332.4652	0.0000	1.0000
SYMG	17	-1151.9316	2337.8632	0.0000	1.0000
SYMI	17	-1152.4411	2338.8823	0.0000	1.0000
SYMIG	18	-1151.6328	2339.2656	0.0000	1.0000
F81	14	-1158.8413	2345.6825	0.0000	1.0000
K2P	12	-1166.9360	2357.8719	0.0000	1.0000
SYM	16	-1166.5022	2365.0044	0.0000	1.0000
JC69G	12	-1170.5802	2365.1604	0.0000	1.0000
JC69I	12	-1171.0910	2366.1820	0.0000	1.0000
JC69IG	13	-1170.3494	2366.6987	0.0000	1.0000
JC69	11	-1183.7358	2389.4716	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1123.4677	2344.4686	0.6378	0.6378
HKYI	16	-1124.1644	2345.8619	0.3178	0.9556
HKYIG	17	-1123.0865	2349.8020	0.0443	1.0000
GTRG	20	-1121.4909	2364.8982	0.0000	1.0000
GTRI	20	-1122.2274	2366.3714	0.0000	1.0000
HKY	15	-1138.6546	2368.7466	0.0000	1.0000
GTRIG	21	-1121.1270	2370.2663	0.0000	1.0000
F81G	15	-1145.5759	2382.5892	0.0000	1.0000
F81I	15	-1146.1919	2383.8213	0.0000	1.0000
K2PG	13	-1152.5323	2384.3104	0.0000	1.0000
K2PI	13	-1153.0759	2385.3975	0.0000	1.0000
F81IG	16	-1145.2848	2388.1027	0.0000	1.0000
GTR	19	-1136.4752	2388.7711	0.0000	1.0000
K2PIG	14	-1152.2326	2389.8067	0.0000	1.0000
F81	14	-1158.8413	2403.0241	0.0000	1.0000
K2P	12	-1166.9360	2407.0218	0.0000	1.0000
SYMG	17	-1151.9316	2407.4923	0.0000	1.0000
SYMI	17	-1152.4411	2408.5113	0.0000	1.0000
SYMIG	18	-1151.6328	2412.9905	0.0000	1.0000
JC69G	12	-1170.5802	2414.3103	0.0000	1.0000
JC69I	12	-1171.0910	2415.3319	0.0000	1.0000
JC69IG	13	-1170.3494	2419.9445	0.0000	1.0000
SYM	16	-1166.5022	2430.5376	0.0000	1.0000
JC69	11	-1183.7358	2434.5257	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/chr12_5730.phylip = 1 - 444;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5730.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/chr12_5730.phylip = 1 - 444;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5730.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/chr12_5730.phylip = 1 - 444;
 Partition Dummy = 1:591-tomodel-phylip/chr12_5730.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

