
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr2_24827.phylip" (7 taxa, 304 characters)

Minimum AIC  model: SYMI
Minimum AICc model: SYMI
Minimum BIC  model: K2P


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
SYMI	17	-657.1336	1350.4071	0.2053	0.2053
SYMG	17	-657.2607	1350.6612	0.1808	0.3860
SYM	16	-658.6125	1351.1204	0.1437	0.5297
K2PI	13	-662.5147	1352.2846	0.0803	0.6100
K2PG	13	-662.5993	1352.4538	0.0738	0.6838
SYMIG	18	-657.2736	1352.9473	0.0576	0.7414
K2P	12	-664.0524	1353.1769	0.0514	0.7928
GTRI	20	-655.3824	1353.7331	0.0389	0.8317
GTRG	20	-655.5133	1353.9947	0.0341	0.8658
GTR	19	-656.8292	1354.3344	0.0288	0.8947
K2PIG	14	-662.6101	1354.6736	0.0243	0.9190
HKYI	16	-660.3949	1354.6853	0.0242	0.9431
HKYG	16	-660.4691	1354.8337	0.0224	0.9656
HKY	15	-662.0184	1355.7034	0.0145	0.9801
GTRIG	21	-655.5257	1356.3280	0.0106	0.9907
HKYIG	17	-660.4808	1357.1015	0.0072	0.9980
JC69I	12	-668.5997	1362.2715	0.0005	0.9985
JC69G	12	-668.6785	1362.4291	0.0005	0.9990
JC69	11	-670.0963	1363.0967	0.0004	0.9994
JC69IG	13	-668.6879	1364.6310	0.0002	0.9995
F81I	15	-666.5399	1364.7464	0.0002	0.9997
F81G	15	-666.6190	1364.9047	0.0001	0.9998
F81	14	-668.0355	1365.5243	0.0001	1.0000
F81IG	16	-666.6285	1367.1525	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMI	17	-657.1336	1348.2673	0.2153	0.2153
SYMG	17	-657.2607	1348.5213	0.1896	0.4049
SYM	16	-658.6125	1349.2249	0.1334	0.5383
SYMIG	18	-657.2736	1350.5473	0.0689	0.6071
GTRI	20	-655.3824	1350.7649	0.0618	0.6689
GTRG	20	-655.5133	1351.0265	0.0542	0.7231
K2PI	13	-662.5147	1351.0294	0.0541	0.7772
K2PG	13	-662.5993	1351.1986	0.0497	0.8269
GTR	19	-656.8292	1351.6584	0.0395	0.8664
K2P	12	-664.0524	1352.1047	0.0316	0.8980
HKYI	16	-660.3949	1352.7898	0.0224	0.9204
HKYG	16	-660.4691	1352.9382	0.0208	0.9413
GTRIG	21	-655.5257	1353.0514	0.0197	0.9610
K2PIG	14	-662.6101	1353.2203	0.0181	0.9791
HKY	15	-662.0184	1354.0367	0.0120	0.9911
HKYIG	17	-660.4808	1354.9616	0.0076	0.9987
JC69I	12	-668.5997	1361.1993	0.0003	0.9990
JC69G	12	-668.6785	1361.3570	0.0003	0.9993
JC69	11	-670.0963	1362.1926	0.0002	0.9995
F81I	15	-666.5399	1363.0797	0.0001	0.9996
F81G	15	-666.6190	1363.2381	0.0001	0.9998
JC69IG	13	-668.6879	1363.3758	0.0001	0.9999
F81	14	-668.0355	1364.0711	0.0001	1.0000
F81IG	16	-666.6285	1365.2570	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-664.0524	1396.7091	0.6272	0.6272
K2PI	13	-662.5147	1399.3508	0.1674	0.7946
K2PG	13	-662.5993	1399.5200	0.1538	0.9484
JC69	11	-670.0963	1403.0799	0.0259	0.9743
K2PIG	14	-662.6101	1405.2587	0.0087	0.9831
JC69I	12	-668.5997	1405.8037	0.0066	0.9897
JC69G	12	-668.6785	1405.9613	0.0061	0.9958
SYM	16	-658.6125	1408.6974	0.0016	0.9974
HKY	15	-662.0184	1409.7921	0.0009	0.9983
SYMI	17	-657.1336	1411.4568	0.0004	0.9987
JC69IG	13	-668.6879	1411.6972	0.0003	0.9991
SYMG	17	-657.2607	1411.7108	0.0003	0.9994
HKYI	16	-660.3949	1412.2622	0.0003	0.9997
HKYG	16	-660.4691	1412.4107	0.0002	0.9999
F81	14	-668.0355	1416.1094	0.0000	0.9999
SYMIG	18	-657.2736	1417.4538	0.0000	1.0000
HKYIG	17	-660.4808	1418.1511	0.0000	1.0000
F81I	15	-666.5399	1418.8351	0.0000	1.0000
F81G	15	-666.6190	1418.9935	0.0000	1.0000
GTR	19	-656.8292	1422.2819	0.0000	1.0000
F81IG	16	-666.6285	1424.7294	0.0000	1.0000
GTRI	20	-655.3824	1425.1054	0.0000	1.0000
GTRG	20	-655.5133	1425.3671	0.0000	1.0000
GTRIG	21	-655.5257	1431.1090	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24827.phylip = 1 - 304;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24827.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24827.phylip = 1 - 304;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24827.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24827.phylip = 1 - 304;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24827.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

