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

Input data from file "591-tomodel-phylip/chr1_32246.phylip" (7 taxa, 643 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-1198.7645	2438.8794	0.4259	0.4259
GTRG	20	-1198.7666	2438.8837	0.4250	0.8509
GTRIG	21	-1198.7646	2441.0171	0.1463	0.9972
GTR	19	-1205.2907	2449.8013	0.0018	0.9990
HKYI	16	-1209.8850	2452.6390	0.0004	0.9995
HKYG	16	-1210.0124	2452.8937	0.0004	0.9998
HKYIG	17	-1209.8850	2454.7492	0.0002	1.0000
HKY	15	-1217.7543	2466.2742	0.0000	1.0000
F81I	15	-1228.4997	2487.7650	0.0000	1.0000
F81G	15	-1228.6484	2488.0623	0.0000	1.0000
F81IG	16	-1228.4995	2489.8681	0.0000	1.0000
SYMI	17	-1227.9677	2490.9147	0.0000	1.0000
SYMG	17	-1227.9682	2490.9156	0.0000	1.0000
SYMIG	18	-1227.9677	2493.0315	0.0000	1.0000
F81	14	-1235.6887	2500.0461	0.0000	1.0000
SYM	16	-1234.3771	2501.6233	0.0000	1.0000
K2PI	13	-1241.4662	2509.5111	0.0000	1.0000
K2PG	13	-1241.6544	2509.8874	0.0000	1.0000
K2PIG	14	-1241.4666	2511.6020	0.0000	1.0000
K2P	12	-1249.5204	2523.5360	0.0000	1.0000
JC69I	12	-1258.5997	2541.6947	0.0000	1.0000
JC69G	12	-1258.7339	2541.9631	0.0000	1.0000
JC69IG	13	-1258.5999	2543.7786	0.0000	1.0000
JC69	11	-1265.6929	2553.8042	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1198.7645	2437.5289	0.4217	0.4217
GTRG	20	-1198.7666	2437.5332	0.4208	0.8424
GTRIG	21	-1198.7646	2439.5292	0.1551	0.9976
GTR	19	-1205.2907	2448.5814	0.0017	0.9992
HKYI	16	-1209.8850	2451.7700	0.0003	0.9996
HKYG	16	-1210.0124	2452.0247	0.0003	0.9999
HKYIG	17	-1209.8850	2453.7700	0.0001	1.0000
HKY	15	-1217.7543	2465.5086	0.0000	1.0000
F81I	15	-1228.4997	2486.9995	0.0000	1.0000
F81G	15	-1228.6484	2487.2968	0.0000	1.0000
F81IG	16	-1228.4995	2488.9991	0.0000	1.0000
SYMI	17	-1227.9677	2489.9355	0.0000	1.0000
SYMG	17	-1227.9682	2489.9364	0.0000	1.0000
SYMIG	18	-1227.9677	2491.9354	0.0000	1.0000
F81	14	-1235.6887	2499.3773	0.0000	1.0000
SYM	16	-1234.3771	2500.7543	0.0000	1.0000
K2PI	13	-1241.4662	2508.9324	0.0000	1.0000
K2PG	13	-1241.6544	2509.3087	0.0000	1.0000
K2PIG	14	-1241.4666	2510.9332	0.0000	1.0000
K2P	12	-1249.5204	2523.0407	0.0000	1.0000
JC69I	12	-1258.5997	2541.1994	0.0000	1.0000
JC69G	12	-1258.7339	2541.4678	0.0000	1.0000
JC69IG	13	-1258.5999	2543.1999	0.0000	1.0000
JC69	11	-1265.6929	2553.3859	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1209.8850	2523.2284	0.4408	0.4408
HKYG	16	-1210.0124	2523.4830	0.3881	0.8289
GTRI	20	-1198.7645	2526.8518	0.0720	0.9010
GTRG	20	-1198.7666	2526.8561	0.0719	0.9728
HKYIG	17	-1209.8850	2529.6945	0.0174	0.9902
HKY	15	-1217.7543	2532.5008	0.0043	0.9945
GTRIG	21	-1198.7646	2533.3182	0.0028	0.9973
GTR	19	-1205.2907	2533.4382	0.0027	1.0000
F81I	15	-1228.4997	2553.9917	0.0000	1.0000
F81G	15	-1228.6484	2554.2889	0.0000	1.0000
F81IG	16	-1228.4995	2560.4574	0.0000	1.0000
F81	14	-1235.6887	2561.9033	0.0000	1.0000
SYMI	17	-1227.9677	2565.8599	0.0000	1.0000
SYMG	17	-1227.9682	2565.8609	0.0000	1.0000
K2PI	13	-1241.4662	2566.9923	0.0000	1.0000
K2PG	13	-1241.6544	2567.3686	0.0000	1.0000
SYM	16	-1234.3771	2572.2126	0.0000	1.0000
SYMIG	18	-1227.9677	2572.3260	0.0000	1.0000
K2PIG	14	-1241.4666	2573.4592	0.0000	1.0000
K2P	12	-1249.5204	2576.6345	0.0000	1.0000
JC69I	12	-1258.5997	2594.7932	0.0000	1.0000
JC69G	12	-1258.7339	2595.0616	0.0000	1.0000
JC69IG	13	-1258.5999	2601.2597	0.0000	1.0000
JC69	11	-1265.6929	2602.5135	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32246.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32246.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32246.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32246.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32246.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32246.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

