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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1105.6074	2253.0981	0.3463	0.3463
GTRIG	21	-1104.9737	2254.0239	0.2180	0.5644
GTRI	20	-1106.5278	2254.9390	0.1380	0.7023
HKYG	16	-1110.9204	2255.0496	0.1305	0.8329
HKYIG	17	-1110.2031	2255.7693	0.0911	0.9240
HKYI	16	-1111.5268	2256.2624	0.0712	0.9952
F81G	15	-1116.1067	2263.2777	0.0021	0.9973
F81IG	16	-1115.4288	2264.0664	0.0014	0.9987
F81I	15	-1116.6387	2264.3417	0.0013	1.0000
GTR	19	-1116.6677	2273.0355	0.0000	1.0000
HKY	15	-1123.2400	2277.5443	0.0000	1.0000
K2PG	13	-1129.5342	2285.8720	0.0000	1.0000
F81	14	-1128.6195	2286.1682	0.0000	1.0000
K2PIG	14	-1128.7401	2286.4094	0.0000	1.0000
K2PI	13	-1129.9135	2286.6306	0.0000	1.0000
SYMG	17	-1127.9225	2291.2080	0.0000	1.0000
SYMIG	18	-1127.1798	2291.8864	0.0000	1.0000
SYMI	17	-1128.4378	2292.2385	0.0000	1.0000
JC69G	12	-1134.0687	2292.8246	0.0000	1.0000
JC69IG	13	-1133.3578	2293.5192	0.0000	1.0000
JC69I	12	-1134.4568	2293.6007	0.0000	1.0000
K2P	12	-1143.3297	2311.3466	0.0000	1.0000
SYM	16	-1140.5214	2314.2516	0.0000	1.0000
JC69	11	-1147.5538	2317.6878	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1105.6074	2251.2147	0.3677	0.3677
GTRIG	21	-1104.9737	2251.9475	0.2549	0.6226
GTRI	20	-1106.5278	2253.0556	0.1465	0.7691
HKYG	16	-1110.9204	2253.8407	0.0989	0.8680
HKYIG	17	-1110.2031	2254.4062	0.0746	0.9426
HKYI	16	-1111.5268	2255.0536	0.0539	0.9965
F81G	15	-1116.1067	2262.2134	0.0015	0.9980
F81IG	16	-1115.4288	2262.8575	0.0011	0.9991
F81I	15	-1116.6387	2263.2774	0.0009	1.0000
GTR	19	-1116.6677	2271.3353	0.0000	1.0000
HKY	15	-1123.2400	2276.4800	0.0000	1.0000
K2PG	13	-1129.5342	2285.0684	0.0000	1.0000
F81	14	-1128.6195	2285.2390	0.0000	1.0000
K2PIG	14	-1128.7401	2285.4802	0.0000	1.0000
K2PI	13	-1129.9135	2285.8270	0.0000	1.0000
SYMG	17	-1127.9225	2289.8450	0.0000	1.0000
SYMIG	18	-1127.1798	2290.3596	0.0000	1.0000
SYMI	17	-1128.4378	2290.8755	0.0000	1.0000
JC69G	12	-1134.0687	2292.1374	0.0000	1.0000
JC69IG	13	-1133.3578	2292.7157	0.0000	1.0000
JC69I	12	-1134.4568	2292.9135	0.0000	1.0000
K2P	12	-1143.3297	2310.6594	0.0000	1.0000
SYM	16	-1140.5214	2313.0427	0.0000	1.0000
JC69	11	-1147.5538	2317.1076	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1110.9204	2320.1820	0.5421	0.5421
HKYI	16	-1111.5268	2321.3948	0.2956	0.8377
F81G	15	-1116.1067	2324.4084	0.0655	0.9032
HKYIG	17	-1110.2031	2324.8938	0.0514	0.9546
F81I	15	-1116.6387	2325.4724	0.0385	0.9931
F81IG	16	-1115.4288	2329.1988	0.0060	0.9991
GTRG	20	-1105.6074	2334.1413	0.0005	0.9996
GTRI	20	-1106.5278	2335.9821	0.0002	0.9998
HKY	15	-1123.2400	2338.6749	0.0001	0.9999
K2PG	13	-1129.5342	2338.9707	0.0000	0.9999
GTRIG	21	-1104.9737	2339.0204	0.0000	0.9999
K2PI	13	-1129.9135	2339.7293	0.0000	1.0000
JC69G	12	-1134.0687	2341.8933	0.0000	1.0000
JC69I	12	-1134.4568	2342.6695	0.0000	1.0000
F81	14	-1128.6195	2343.2876	0.0000	1.0000
K2PIG	14	-1128.7401	2343.5288	0.0000	1.0000
JC69IG	13	-1133.3578	2346.6179	0.0000	1.0000
GTR	19	-1116.6677	2350.1156	0.0000	1.0000
SYMG	17	-1127.9225	2360.3326	0.0000	1.0000
K2P	12	-1143.3297	2360.4154	0.0000	1.0000
SYMI	17	-1128.4378	2361.3631	0.0000	1.0000
JC69	11	-1147.5538	2362.7172	0.0000	1.0000
SYMIG	18	-1127.1798	2364.9936	0.0000	1.0000
SYM	16	-1140.5214	2379.3840	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3414.phylip = 1 - 467;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3414.phylip = 1 - 467;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3414.phylip = 1 - 467;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3414.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

