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

Input data from file "591-tomodel-phylip/chr2_9158.phylip" (7 taxa, 447 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	16	-1221.3763	2476.0177	0.4343	0.4343
HKYIG	17	-1220.8615	2477.1495	0.2466	0.6808
HKYI	16	-1222.0933	2477.4517	0.2120	0.8928
GTRG	20	-1219.4888	2480.9494	0.0369	0.9297
GTRI	20	-1219.8268	2481.6254	0.0263	0.9560
GTRIG	21	-1218.8907	2481.9555	0.0223	0.9783
K2PG	13	-1228.3559	2483.5525	0.0100	0.9884
K2PIG	14	-1227.6882	2484.3486	0.0067	0.9951
K2PI	13	-1229.1753	2485.1912	0.0044	0.9995
SYMG	17	-1227.9583	2491.3432	0.0002	0.9997
SYMIG	18	-1227.1233	2491.8446	0.0002	0.9999
SYMI	17	-1228.6662	2492.7591	0.0001	1.0000
HKY	15	-1238.8693	2508.8522	0.0000	1.0000
F81IG	16	-1238.3608	2509.9868	0.0000	1.0000
F81G	15	-1240.6150	2512.3436	0.0000	1.0000
JC69IG	13	-1243.5478	2513.9362	0.0000	1.0000
F81I	15	-1241.5995	2514.3127	0.0000	1.0000
GTR	19	-1237.3832	2514.5463	0.0000	1.0000
K2P	12	-1246.0241	2516.7671	0.0000	1.0000
JC69G	12	-1246.6504	2518.0198	0.0000	1.0000
JC69I	12	-1247.5285	2519.7760	0.0000	1.0000
SYM	16	-1245.6011	2524.4674	0.0000	1.0000
F81	14	-1259.5320	2548.0363	0.0000	1.0000
JC69	11	-1265.6031	2553.8131	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1221.3763	2474.7526	0.4109	0.4109
HKYIG	17	-1220.8615	2475.7229	0.2530	0.6639
HKYI	16	-1222.0933	2476.1866	0.2006	0.8645
GTRG	20	-1219.4888	2478.9776	0.0497	0.9142
GTRI	20	-1219.8268	2479.6535	0.0354	0.9497
GTRIG	21	-1218.8907	2479.7814	0.0332	0.9829
K2PG	13	-1228.3559	2482.7119	0.0077	0.9906
K2PIG	14	-1227.6882	2483.3764	0.0055	0.9961
K2PI	13	-1229.1753	2484.3505	0.0034	0.9995
SYMG	17	-1227.9583	2489.9166	0.0002	0.9997
SYMIG	18	-1227.1233	2490.2465	0.0002	0.9999
SYMI	17	-1228.6662	2491.3325	0.0001	1.0000
HKY	15	-1238.8693	2507.7385	0.0000	1.0000
F81IG	16	-1238.3608	2508.7217	0.0000	1.0000
F81G	15	-1240.6150	2511.2299	0.0000	1.0000
GTR	19	-1237.3832	2512.7665	0.0000	1.0000
JC69IG	13	-1243.5478	2513.0955	0.0000	1.0000
F81I	15	-1241.5995	2513.1990	0.0000	1.0000
K2P	12	-1246.0241	2516.0482	0.0000	1.0000
JC69G	12	-1246.6504	2517.3009	0.0000	1.0000
JC69I	12	-1247.5285	2519.0571	0.0000	1.0000
SYM	16	-1245.6011	2523.2023	0.0000	1.0000
F81	14	-1259.5320	2547.0640	0.0000	1.0000
JC69	11	-1265.6031	2553.2062	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1228.3559	2536.0451	0.5844	0.5844
K2PI	13	-1229.1753	2537.6838	0.2576	0.8419
HKYG	16	-1221.3763	2540.3935	0.0664	0.9084
K2PIG	14	-1227.6882	2540.8122	0.0539	0.9623
HKYI	16	-1222.0933	2541.8275	0.0324	0.9947
HKYIG	17	-1220.8615	2545.4664	0.0053	1.0000
SYMG	17	-1227.9583	2559.6601	0.0000	1.0000
GTRG	20	-1219.4888	2561.0288	0.0000	1.0000
SYMI	17	-1228.6662	2561.0760	0.0000	1.0000
GTRI	20	-1219.8268	2561.7047	0.0000	1.0000
SYMIG	18	-1227.1233	2564.0926	0.0000	1.0000
K2P	12	-1246.0241	2565.2789	0.0000	1.0000
GTRIG	21	-1218.8907	2565.9352	0.0000	1.0000
JC69IG	13	-1243.5478	2566.4288	0.0000	1.0000
JC69G	12	-1246.6504	2566.5316	0.0000	1.0000
JC69I	12	-1247.5285	2568.2878	0.0000	1.0000
HKY	15	-1238.8693	2569.2769	0.0000	1.0000
F81G	15	-1240.6150	2572.7683	0.0000	1.0000
F81IG	16	-1238.3608	2574.3626	0.0000	1.0000
F81I	15	-1241.5995	2574.7374	0.0000	1.0000
SYM	16	-1245.6011	2588.8432	0.0000	1.0000
GTR	19	-1237.3832	2590.7151	0.0000	1.0000
JC69	11	-1265.6031	2598.3343	0.0000	1.0000
F81	14	-1259.5320	2604.4999	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_9158.phylip = 1 - 447;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9158.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_9158.phylip = 1 - 447;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9158.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 591-tomodel-phylip/chr2_9158.phylip = 1 - 447;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9158.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

