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

Input data from file "591-tomodel-phylip/chr2_23621.phylip" (7 taxa, 443 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	-1227.0838	2487.4446	0.3738	0.3738
HKYI	16	-1227.3297	2487.9365	0.2923	0.6661
HKYIG	17	-1227.0813	2489.6026	0.1271	0.7931
GTRG	20	-1224.0642	2490.1190	0.0981	0.8913
GTRI	20	-1224.3379	2490.6664	0.0746	0.9659
GTRIG	21	-1224.0642	2492.3231	0.0326	0.9985
HKY	15	-1233.9854	2499.0949	0.0011	0.9996
GTR	19	-1230.8842	2501.5650	0.0003	1.0000
K2PI	13	-1240.2284	2507.3052	0.0000	1.0000
K2PG	13	-1240.5548	2507.9581	0.0000	1.0000
K2PIG	14	-1240.0825	2509.1463	0.0000	1.0000
SYMI	17	-1239.0316	2513.5032	0.0000	1.0000
SYMG	17	-1239.1679	2513.7757	0.0000	1.0000
SYMIG	18	-1238.9239	2515.4609	0.0000	1.0000
K2P	12	-1245.8067	2516.3390	0.0000	1.0000
SYM	16	-1244.3197	2521.9164	0.0000	1.0000
F81I	15	-1258.8779	2548.8799	0.0000	1.0000
F81G	15	-1259.1062	2549.3366	0.0000	1.0000
F81IG	16	-1261.4640	2556.2051	0.0000	1.0000
F81	14	-1263.8826	2556.7464	0.0000	1.0000
JC69I	12	-1268.5512	2561.8280	0.0000	1.0000
JC69G	12	-1268.8450	2562.4156	0.0000	1.0000
JC69IG	13	-1268.5596	2563.9677	0.0000	1.0000
JC69	11	-1273.4052	2569.4230	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1227.0838	2486.1676	0.3386	0.3386
HKYI	16	-1227.3297	2486.6595	0.2648	0.6034
GTRG	20	-1224.0642	2488.1285	0.1270	0.7304
HKYIG	17	-1227.0813	2488.1626	0.1249	0.8553
GTRI	20	-1224.3379	2488.6759	0.0966	0.9519
GTRIG	21	-1224.0642	2490.1283	0.0467	0.9987
HKY	15	-1233.9854	2497.9707	0.0009	0.9996
GTR	19	-1230.8842	2499.7683	0.0004	1.0000
K2PI	13	-1240.2284	2506.4567	0.0000	1.0000
K2PG	13	-1240.5548	2507.1096	0.0000	1.0000
K2PIG	14	-1240.0825	2508.1650	0.0000	1.0000
SYMI	17	-1239.0316	2512.0632	0.0000	1.0000
SYMG	17	-1239.1679	2512.3357	0.0000	1.0000
SYMIG	18	-1238.9239	2513.8477	0.0000	1.0000
K2P	12	-1245.8067	2515.6134	0.0000	1.0000
SYM	16	-1244.3197	2520.6394	0.0000	1.0000
F81I	15	-1258.8779	2547.7558	0.0000	1.0000
F81G	15	-1259.1062	2548.2125	0.0000	1.0000
F81IG	16	-1261.4640	2554.9281	0.0000	1.0000
F81	14	-1263.8826	2555.7651	0.0000	1.0000
JC69I	12	-1268.5512	2561.1024	0.0000	1.0000
JC69G	12	-1268.8450	2561.6900	0.0000	1.0000
JC69IG	13	-1268.5596	2563.1193	0.0000	1.0000
JC69	11	-1273.4052	2568.8105	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1227.0838	2551.6648	0.5305	0.5305
HKYI	16	-1227.3297	2552.1566	0.4149	0.9454
HKYIG	17	-1227.0813	2557.7533	0.0253	0.9707
HKY	15	-1233.9854	2559.3743	0.0112	0.9819
K2PI	13	-1240.2284	2559.6731	0.0097	0.9916
K2PG	13	-1240.5548	2560.3260	0.0070	0.9986
K2P	12	-1245.8067	2564.7362	0.0008	0.9994
K2PIG	14	-1240.0825	2565.4750	0.0005	0.9999
GTRG	20	-1224.0642	2569.9999	0.0001	1.0000
GTRI	20	-1224.3379	2570.5473	0.0000	1.0000
GTRIG	21	-1224.0642	2576.0933	0.0000	1.0000
GTR	19	-1230.8842	2577.5462	0.0000	1.0000
SYMI	17	-1239.0316	2581.6539	0.0000	1.0000
SYMG	17	-1239.1679	2581.9264	0.0000	1.0000
SYM	16	-1244.3197	2586.1366	0.0000	1.0000
SYMIG	18	-1238.9239	2587.5320	0.0000	1.0000
F81I	15	-1258.8779	2609.1593	0.0000	1.0000
F81G	15	-1259.1062	2609.6160	0.0000	1.0000
JC69I	12	-1268.5512	2610.2252	0.0000	1.0000
JC69G	12	-1268.8450	2610.8128	0.0000	1.0000
F81	14	-1263.8826	2613.0751	0.0000	1.0000
JC69	11	-1273.4052	2613.8397	0.0000	1.0000
JC69IG	13	-1268.5596	2616.3357	0.0000	1.0000
F81IG	16	-1261.4640	2620.4252	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_23621.phylip = 1 - 443;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23621.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_23621.phylip = 1 - 443;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23621.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/chr2_23621.phylip = 1 - 443;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23621.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

