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

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

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
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
HKYI	16	-582.4122	1198.9746	0.3111	0.3111
GTRI	20	-578.0998	1199.5731	0.2307	0.5418
GTRG	20	-578.6213	1200.6162	0.1369	0.6787
HKYG	16	-583.3819	1200.9140	0.1180	0.7967
HKYIG	17	-582.4123	1201.2531	0.0996	0.8963
GTRIG	21	-578.0933	1201.9124	0.0716	0.9679
F81I	15	-586.4602	1204.8101	0.0168	0.9847
F81G	15	-587.0150	1205.9197	0.0097	0.9944
F81IG	16	-586.4602	1207.0706	0.0054	0.9998
K2PI	13	-594.2182	1215.8583	0.0001	0.9999
K2PG	13	-594.9858	1217.3934	0.0000	0.9999
K2PIG	14	-594.2182	1218.0834	0.0000	0.9999
SYMI	17	-591.1291	1218.6868	0.0000	1.0000
GTR	19	-588.9450	1218.9299	0.0000	1.0000
SYMG	17	-592.0646	1220.5578	0.0000	1.0000
HKY	15	-594.4618	1220.8134	0.0000	1.0000
SYMIG	18	-591.1292	1220.9835	0.0000	1.0000
JC69I	12	-598.0678	1221.3497	0.0000	1.0000
JC69G	12	-598.7588	1222.7315	0.0000	1.0000
JC69IG	13	-598.0680	1223.5580	0.0000	1.0000
F81	14	-597.3797	1224.4065	0.0000	1.0000
K2P	12	-605.8169	1236.8477	0.0000	1.0000
SYM	16	-602.8364	1239.8230	0.0000	1.0000
JC69	11	-609.4240	1241.8713	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-578.0998	1196.1996	0.3021	0.3021
HKYI	16	-582.4122	1196.8244	0.2211	0.5232
GTRG	20	-578.6213	1197.2427	0.1793	0.7025
GTRIG	21	-578.0933	1198.1866	0.1119	0.8144
HKYG	16	-583.3819	1198.7638	0.0838	0.8982
HKYIG	17	-582.4123	1198.8245	0.0813	0.9795
F81I	15	-586.4602	1202.9204	0.0105	0.9900
F81G	15	-587.0150	1204.0300	0.0060	0.9960
F81IG	16	-586.4602	1204.9204	0.0039	0.9999
K2PI	13	-594.2182	1214.4364	0.0000	0.9999
GTR	19	-588.9450	1215.8899	0.0000	0.9999
K2PG	13	-594.9858	1215.9715	0.0000	1.0000
SYMI	17	-591.1291	1216.2582	0.0000	1.0000
K2PIG	14	-594.2182	1216.4364	0.0000	1.0000
SYMG	17	-592.0646	1218.1292	0.0000	1.0000
SYMIG	18	-591.1292	1218.2584	0.0000	1.0000
HKY	15	-594.4618	1218.9236	0.0000	1.0000
JC69I	12	-598.0678	1220.1357	0.0000	1.0000
JC69G	12	-598.7588	1221.5175	0.0000	1.0000
JC69IG	13	-598.0680	1222.1361	0.0000	1.0000
F81	14	-597.3797	1222.7594	0.0000	1.0000
K2P	12	-605.8169	1235.6337	0.0000	1.0000
SYM	16	-602.8364	1237.6728	0.0000	1.0000
JC69	11	-609.4240	1240.8481	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-582.4122	1254.3992	0.5052	0.5052
HKYG	16	-583.3819	1256.3386	0.1916	0.6967
F81I	15	-586.4602	1256.8967	0.1449	0.8416
F81G	15	-587.0150	1258.0063	0.0832	0.9248
HKYIG	17	-582.4123	1259.9977	0.0307	0.9556
K2PI	13	-594.2182	1261.2159	0.0167	0.9723
F81IG	16	-586.4602	1262.4951	0.0088	0.9811
K2PG	13	-594.9858	1262.7510	0.0078	0.9889
JC69I	12	-598.0678	1263.3167	0.0058	0.9947
JC69G	12	-598.7588	1264.6986	0.0029	0.9977
K2PIG	14	-594.2182	1266.8143	0.0010	0.9987
GTRI	20	-578.0998	1268.1680	0.0005	0.9992
JC69IG	13	-598.0680	1268.9156	0.0004	0.9996
GTRG	20	-578.6213	1269.2111	0.0003	0.9999
HKY	15	-594.4618	1272.8999	0.0000	0.9999
F81	14	-597.3797	1273.1373	0.0000	1.0000
GTRIG	21	-578.0933	1273.7534	0.0000	1.0000
SYMI	17	-591.1291	1277.4314	0.0000	1.0000
K2P	12	-605.8169	1278.8148	0.0000	1.0000
SYMG	17	-592.0646	1279.3024	0.0000	1.0000
JC69	11	-609.4240	1280.4307	0.0000	1.0000
SYMIG	18	-591.1292	1283.0300	0.0000	1.0000
GTR	19	-588.9450	1284.2600	0.0000	1.0000
SYM	16	-602.8364	1295.2475	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/chr15_3483.phylip = 1 - 270;
 Partition Dummy = 1:591-tomodel-phylip/chr15_3483.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 (HKYI)]

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


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

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

