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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-780.6547	1594.7999	0.2438	0.2438
HKYG	16	-780.6808	1594.8521	0.2375	0.4813
HKY	15	-781.8852	1595.0820	0.2117	0.6931
HKYIG	17	-780.6827	1597.0468	0.0793	0.7723
GTRG	20	-777.6482	1597.6233	0.0594	0.8318
GTRI	20	-777.6599	1597.6466	0.0587	0.8905
GTR	19	-778.9644	1598.0282	0.0485	0.9390
GTRIG	21	-777.6466	1599.8598	0.0194	0.9584
K2P	12	-788.1796	1601.2047	0.0099	0.9684
K2PI	13	-787.1172	1601.2235	0.0098	0.9782
K2PG	13	-787.1425	1601.2741	0.0096	0.9877
K2PIG	14	-787.1449	1603.4342	0.0033	0.9910
SYM	16	-785.1641	1603.8187	0.0027	0.9937
SYMG	17	-784.0696	1603.8204	0.0027	0.9964
SYMI	17	-784.0814	1603.8441	0.0026	0.9990
SYMIG	18	-784.0687	1606.0216	0.0009	0.9999
F81	14	-791.8834	1612.9111	0.0000	0.9999
F81I	15	-790.8591	1613.0297	0.0000	1.0000
F81G	15	-790.8966	1613.1046	0.0000	1.0000
F81IG	16	-790.8993	1615.2890	0.0000	1.0000
JC69	11	-797.9557	1618.6250	0.0000	1.0000
JC69I	12	-796.9769	1618.7994	0.0000	1.0000
JC69G	12	-797.0170	1618.8796	0.0000	1.0000
JC69IG	13	-797.0193	1621.0278	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-780.6547	1593.3095	0.2267	0.2267
HKYG	16	-780.6808	1593.3617	0.2209	0.4476
HKY	15	-781.8852	1593.7705	0.1800	0.6276
GTRG	20	-777.6482	1595.2964	0.0839	0.7115
GTRI	20	-777.6599	1595.3197	0.0830	0.7945
HKYIG	17	-780.6827	1595.3655	0.0811	0.8756
GTR	19	-778.9644	1595.9287	0.0612	0.9368
GTRIG	21	-777.6466	1597.2931	0.0309	0.9677
K2PI	13	-787.1172	1600.2344	0.0071	0.9748
K2PG	13	-787.1425	1600.2849	0.0069	0.9817
K2P	12	-788.1796	1600.3592	0.0067	0.9884
SYMG	17	-784.0696	1602.1391	0.0027	0.9912
SYMI	17	-784.0814	1602.1627	0.0027	0.9939
K2PIG	14	-787.1449	1602.2898	0.0025	0.9964
SYM	16	-785.1641	1602.3283	0.0025	0.9989
SYMIG	18	-784.0687	1604.1373	0.0010	0.9999
F81I	15	-790.8591	1611.7182	0.0000	0.9999
F81	14	-791.8834	1611.7667	0.0000	1.0000
F81G	15	-790.8966	1611.7931	0.0000	1.0000
F81IG	16	-790.8993	1613.7986	0.0000	1.0000
JC69	11	-797.9557	1617.9115	0.0000	1.0000
JC69I	12	-796.9769	1617.9539	0.0000	1.0000
JC69G	12	-797.0170	1618.0341	0.0000	1.0000
JC69IG	13	-797.0193	1620.0387	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-788.1796	1647.7042	0.7145	0.7145
K2PI	13	-787.1172	1651.5248	0.1058	0.8202
K2PG	13	-787.1425	1651.5754	0.1031	0.9233
HKY	15	-781.8852	1652.9518	0.0518	0.9752
HKYI	16	-780.6547	1656.4362	0.0091	0.9842
HKYG	16	-780.6808	1656.4884	0.0088	0.9931
K2PIG	14	-787.1449	1657.5257	0.0053	0.9983
JC69	11	-797.9557	1661.3111	0.0008	0.9991
HKYIG	17	-780.6827	1662.4376	0.0005	0.9996
JC69I	12	-796.9769	1665.2989	0.0001	0.9997
JC69G	12	-797.0170	1665.3791	0.0001	0.9998
SYM	16	-785.1641	1665.4550	0.0001	0.9999
F81	14	-791.8834	1667.0026	0.0000	0.9999
SYMG	17	-784.0696	1669.2113	0.0000	1.0000
SYMI	17	-784.0814	1669.2349	0.0000	1.0000
GTR	19	-778.9644	1670.8917	0.0000	1.0000
F81I	15	-790.8591	1670.8995	0.0000	1.0000
F81G	15	-790.8966	1670.9744	0.0000	1.0000
JC69IG	13	-797.0193	1671.3291	0.0000	1.0000
GTRG	20	-777.6482	1674.2049	0.0000	1.0000
GTRI	20	-777.6599	1674.2281	0.0000	1.0000
SYMIG	18	-784.0687	1675.1549	0.0000	1.0000
F81IG	16	-790.8993	1676.9253	0.0000	1.0000
GTRIG	21	-777.6466	1680.1470	0.0000	1.0000



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

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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_7982.phylip = 1 - 382;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7982.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 (K2P)]

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

