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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-755.4754	1553.2028	0.3189	0.3189
HKYI	16	-760.1144	1553.6718	0.2522	0.5711
HKYG	16	-760.8092	1555.0614	0.1259	0.6970
GTRG	20	-756.5140	1555.2801	0.1129	0.8099
GTRIG	21	-755.4755	1555.4348	0.1045	0.9143
HKYIG	17	-760.1145	1555.8567	0.0846	0.9989
F81I	15	-767.7628	1566.7955	0.0004	0.9993
HKY	15	-768.1411	1567.5520	0.0002	0.9995
F81G	15	-768.2255	1567.7209	0.0002	0.9998
F81IG	16	-767.7627	1568.9684	0.0001	0.9999
GTR	19	-764.5546	1569.1412	0.0001	1.0000
F81	14	-775.1167	1579.3416	0.0000	1.0000
K2PI	13	-782.1167	1591.1913	0.0000	1.0000
K2PG	13	-782.5389	1592.0357	0.0000	1.0000
K2PIG	14	-782.1165	1593.3412	0.0000	1.0000
SYMI	17	-780.3893	1596.4062	0.0000	1.0000
SYMG	17	-780.9280	1597.4836	0.0000	1.0000
SYMIG	18	-780.3892	1598.6024	0.0000	1.0000
JC69I	12	-788.1757	1601.1703	0.0000	1.0000
JC69G	12	-788.5837	1601.9863	0.0000	1.0000
JC69IG	13	-788.1758	1603.3096	0.0000	1.0000
K2P	12	-789.4566	1603.7321	0.0000	1.0000
SYM	16	-787.8235	1609.0900	0.0000	1.0000
JC69	11	-795.3600	1613.4111	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-755.4754	1550.9508	0.3691	0.3691
HKYI	16	-760.1144	1552.2288	0.1948	0.5639
GTRIG	21	-755.4755	1552.9509	0.1358	0.6997
GTRG	20	-756.5140	1553.0281	0.1306	0.8303
HKYG	16	-760.8092	1553.6184	0.0972	0.9276
HKYIG	17	-760.1145	1554.2291	0.0717	0.9992
F81I	15	-767.7628	1565.5257	0.0003	0.9995
HKY	15	-768.1411	1566.2822	0.0002	0.9996
F81G	15	-768.2255	1566.4511	0.0002	0.9998
GTR	19	-764.5546	1567.1091	0.0001	0.9999
F81IG	16	-767.7627	1567.5254	0.0001	1.0000
F81	14	-775.1167	1578.2334	0.0000	1.0000
K2PI	13	-782.1167	1590.2334	0.0000	1.0000
K2PG	13	-782.5389	1591.0778	0.0000	1.0000
K2PIG	14	-782.1165	1592.2330	0.0000	1.0000
SYMI	17	-780.3893	1594.7785	0.0000	1.0000
SYMG	17	-780.9280	1595.8559	0.0000	1.0000
SYMIG	18	-780.3892	1596.7784	0.0000	1.0000
JC69I	12	-788.1757	1600.3514	0.0000	1.0000
JC69G	12	-788.5837	1601.1674	0.0000	1.0000
JC69IG	13	-788.1758	1602.3517	0.0000	1.0000
K2P	12	-789.4566	1602.9132	0.0000	1.0000
SYM	16	-787.8235	1607.6470	0.0000	1.0000
JC69	11	-795.3600	1612.7200	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-760.1144	1615.8504	0.6357	0.6357
HKYG	16	-760.8092	1617.2401	0.3173	0.9531
HKYIG	17	-760.1145	1621.8270	0.0320	0.9851
F81I	15	-767.7628	1625.1709	0.0060	0.9911
HKY	15	-768.1411	1625.9275	0.0041	0.9952
F81G	15	-768.2255	1626.0964	0.0038	0.9990
GTRI	20	-755.4754	1630.4778	0.0004	0.9994
F81IG	16	-767.7627	1631.1470	0.0003	0.9997
GTRG	20	-756.5140	1632.5551	0.0001	0.9999
F81	14	-775.1167	1633.9023	0.0001	1.0000
GTRIG	21	-755.4755	1636.4543	0.0000	1.0000
K2PI	13	-782.1167	1641.9259	0.0000	1.0000
GTR	19	-764.5546	1642.6598	0.0000	1.0000
K2PG	13	-782.5389	1642.7704	0.0000	1.0000
K2PIG	14	-782.1165	1647.9019	0.0000	1.0000
JC69I	12	-788.1757	1648.0677	0.0000	1.0000
JC69G	12	-788.5837	1648.8836	0.0000	1.0000
K2P	12	-789.4566	1650.6294	0.0000	1.0000
JC69IG	13	-788.1758	1654.0443	0.0000	1.0000
JC69	11	-795.3600	1656.4598	0.0000	1.0000
SYMI	17	-780.3893	1662.3765	0.0000	1.0000
SYMG	17	-780.9280	1663.4539	0.0000	1.0000
SYMIG	18	-780.3892	1668.3527	0.0000	1.0000
SYM	16	-787.8235	1671.2686	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/chr8_1617.phylip = 1 - 394;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1617.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 (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_1617.phylip = 1 - 394;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1617.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 BIC model (HKYI)]

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

