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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-579.1801	1190.0684	0.5443	0.5443
HKYG	16	-579.1803	1192.3035	0.1780	0.7224
HKYI	16	-579.1858	1192.3146	0.1771	0.8994
HKYIG	17	-579.1861	1194.5658	0.0574	0.9569
GTR	19	-577.7553	1196.2543	0.0247	0.9816
GTRG	20	-577.7524	1198.5483	0.0078	0.9894
GTRI	20	-577.7568	1198.5570	0.0078	0.9972
GTRIG	21	-577.7561	1200.8721	0.0025	0.9997
F81	14	-588.6143	1206.7180	0.0001	0.9998
K2P	12	-591.9172	1208.9329	0.0000	0.9999
F81G	15	-588.6150	1208.9382	0.0000	0.9999
F81I	15	-588.6169	1208.9419	0.0000	0.9999
K2PG	13	-591.9179	1211.1219	0.0000	1.0000
K2PI	13	-591.9195	1211.1252	0.0000	1.0000
F81IG	16	-588.6176	1211.1780	0.0000	1.0000
K2PIG	14	-591.9202	1213.3297	0.0000	1.0000
SYM	16	-591.1434	1216.2297	0.0000	1.0000
SYMG	17	-591.1441	1218.4818	0.0000	1.0000
SYMI	17	-591.1458	1218.4851	0.0000	1.0000
SYMIG	18	-591.1465	1220.7535	0.0000	1.0000
JC69	11	-600.2988	1223.5240	0.0000	1.0000
JC69I	12	-600.2986	1225.6959	0.0000	1.0000
JC69G	12	-600.2997	1225.6979	0.0000	1.0000
JC69IG	13	-600.2998	1227.8859	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-579.1801	1188.3603	0.4972	0.4972
HKYG	16	-579.1803	1190.3606	0.1829	0.6801
HKYI	16	-579.1858	1190.3717	0.1819	0.8620
HKYIG	17	-579.1861	1192.3722	0.0669	0.9289
GTR	19	-577.7553	1193.5106	0.0379	0.9667
GTRG	20	-577.7524	1195.5048	0.0140	0.9807
GTRI	20	-577.7568	1195.5135	0.0139	0.9946
GTRIG	21	-577.7561	1197.5121	0.0051	0.9997
F81	14	-588.6143	1205.2287	0.0001	0.9998
F81G	15	-588.6150	1207.2300	0.0000	0.9999
F81I	15	-588.6169	1207.2337	0.0000	0.9999
K2P	12	-591.9172	1207.8343	0.0000	1.0000
F81IG	16	-588.6176	1209.2351	0.0000	1.0000
K2PG	13	-591.9179	1209.8357	0.0000	1.0000
K2PI	13	-591.9195	1209.8389	0.0000	1.0000
K2PIG	14	-591.9202	1211.8403	0.0000	1.0000
SYM	16	-591.1434	1214.2868	0.0000	1.0000
SYMG	17	-591.1441	1216.2883	0.0000	1.0000
SYMI	17	-591.1458	1216.2915	0.0000	1.0000
SYMIG	18	-591.1465	1218.2930	0.0000	1.0000
JC69	11	-600.2988	1222.5977	0.0000	1.0000
JC69I	12	-600.2986	1224.5973	0.0000	1.0000
JC69G	12	-600.2997	1224.5993	0.0000	1.0000
JC69IG	13	-600.2998	1226.5997	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-579.1801	1243.7662	0.8791	0.8791
HKYG	16	-579.1803	1249.4604	0.0510	0.9301
HKYI	16	-579.1858	1249.4714	0.0507	0.9808
K2P	12	-591.9172	1252.1591	0.0132	0.9940
HKYIG	17	-579.1861	1255.1657	0.0029	0.9970
F81	14	-588.6143	1256.9409	0.0012	0.9982
K2PG	13	-591.9179	1257.8542	0.0008	0.9989
K2PI	13	-591.9195	1257.8575	0.0008	0.9997
F81G	15	-588.6150	1262.6360	0.0001	0.9998
F81I	15	-588.6169	1262.6397	0.0001	0.9998
JC69	11	-600.2988	1263.2287	0.0001	0.9999
K2PIG	14	-591.9202	1263.5526	0.0000	0.9999
GTR	19	-577.7553	1263.6915	0.0000	1.0000
F81IG	16	-588.6176	1268.3348	0.0000	1.0000
JC69I	12	-600.2986	1268.9221	0.0000	1.0000
JC69G	12	-600.2997	1268.9241	0.0000	1.0000
GTRG	20	-577.7524	1269.3794	0.0000	1.0000
GTRI	20	-577.7568	1269.3881	0.0000	1.0000
SYM	16	-591.1434	1273.3865	0.0000	1.0000
JC69IG	13	-600.2998	1274.6182	0.0000	1.0000
GTRIG	21	-577.7561	1275.0805	0.0000	1.0000
SYMG	17	-591.1441	1279.0817	0.0000	1.0000
SYMI	17	-591.1458	1279.0850	0.0000	1.0000
SYMIG	18	-591.1465	1284.7802	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32208.phylip = 1 - 297;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32208.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32208.phylip = 1 - 297;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32208.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32208.phylip = 1 - 297;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32208.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

