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

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

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


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	-813.9064	1658.8363	0.3533	0.3533
F81	14	-815.6435	1660.1806	0.1804	0.5337
HKYI	16	-813.9062	1660.9748	0.1213	0.6550
HKYG	16	-813.9068	1660.9761	0.1212	0.7762
F81I	15	-815.6433	1662.3101	0.0622	0.8384
F81G	15	-815.6440	1662.3115	0.0622	0.9006
HKYIG	17	-813.9092	1663.1290	0.0413	0.9419
F81IG	16	-815.6465	1664.4554	0.0213	0.9631
K2P	12	-820.7580	1666.1770	0.0090	0.9721
GTR	19	-813.5635	1666.7613	0.0067	0.9789
JC69	11	-822.4999	1667.5579	0.0045	0.9834
K2PI	13	-820.7578	1668.2885	0.0031	0.9865
K2PG	13	-820.7585	1668.2898	0.0031	0.9896
GTRI	20	-813.5633	1668.9369	0.0023	0.9919
GTRG	20	-813.5639	1668.9382	0.0023	0.9942
JC69I	12	-822.4997	1669.6604	0.0016	0.9957
JC69G	12	-822.5004	1669.6618	0.0016	0.9973
K2PIG	14	-820.7609	1670.4153	0.0011	0.9984
GTRIG	21	-813.5662	1671.1282	0.0008	0.9991
JC69IG	13	-822.5029	1671.7786	0.0005	0.9997
SYM	16	-820.4578	1674.0780	0.0002	0.9999
SYMI	17	-820.4576	1676.2258	0.0001	0.9999
SYMG	17	-820.4583	1676.2271	0.0001	1.0000
SYMIG	18	-820.4607	1678.3892	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-813.9064	1657.8128	0.3480	0.3480
F81	14	-815.6435	1659.2870	0.1665	0.5145
HKYI	16	-813.9062	1659.8124	0.1280	0.6426
HKYG	16	-813.9068	1659.8137	0.1280	0.7705
F81I	15	-815.6433	1661.2867	0.0613	0.8318
F81G	15	-815.6440	1661.2881	0.0612	0.8930
HKYIG	17	-813.9092	1661.8185	0.0470	0.9400
F81IG	16	-815.6465	1663.2930	0.0225	0.9624
GTR	19	-813.5635	1665.1269	0.0090	0.9714
K2P	12	-820.7580	1665.5160	0.0074	0.9788
JC69	11	-822.4999	1666.9997	0.0035	0.9823
GTRI	20	-813.5633	1667.1266	0.0033	0.9856
GTRG	20	-813.5639	1667.1279	0.0033	0.9889
K2PI	13	-820.7578	1667.5157	0.0027	0.9917
K2PG	13	-820.7585	1667.5170	0.0027	0.9944
JC69I	12	-822.4997	1668.9994	0.0013	0.9957
JC69G	12	-822.5004	1669.0008	0.0013	0.9970
GTRIG	21	-813.5662	1669.1325	0.0012	0.9982
K2PIG	14	-820.7609	1669.5217	0.0010	0.9992
JC69IG	13	-822.5029	1671.0058	0.0005	0.9997
SYM	16	-820.4578	1672.9156	0.0002	0.9998
SYMI	17	-820.4576	1674.9153	0.0001	0.9999
SYMG	17	-820.4583	1674.9166	0.0001	1.0000
SYMIG	18	-820.4607	1676.9214	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
JC69	11	-822.4999	1713.0254	0.6674	0.6674
K2P	12	-820.7580	1715.7258	0.1730	0.8404
F81	14	-815.6435	1717.8651	0.0594	0.8997
JC69I	12	-822.4997	1719.2092	0.0303	0.9300
JC69G	12	-822.5004	1719.2106	0.0303	0.9603
HKY	15	-813.9064	1720.5750	0.0153	0.9756
K2PI	13	-820.7578	1721.9096	0.0079	0.9835
K2PG	13	-820.7585	1721.9109	0.0079	0.9913
F81I	15	-815.6433	1724.0489	0.0027	0.9940
F81G	15	-815.6440	1724.0503	0.0027	0.9967
JC69IG	13	-822.5029	1725.3997	0.0014	0.9981
HKYI	16	-813.9062	1726.7588	0.0007	0.9988
HKYG	16	-813.9068	1726.7601	0.0007	0.9995
K2PIG	14	-820.7609	1728.0998	0.0004	0.9998
F81IG	16	-815.6465	1730.2394	0.0001	1.0000
HKYIG	17	-813.9092	1732.9490	0.0000	1.0000
SYM	16	-820.4578	1739.8620	0.0000	1.0000
GTR	19	-813.5635	1744.6257	0.0000	1.0000
SYMI	17	-820.4576	1746.0458	0.0000	1.0000
SYMG	17	-820.4583	1746.0472	0.0000	1.0000
GTRI	20	-813.5633	1750.8095	0.0000	1.0000
GTRG	20	-813.5639	1750.8109	0.0000	1.0000
SYMIG	18	-820.4607	1752.2361	0.0000	1.0000
GTRIG	21	-813.5662	1756.9996	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/chrZ_8278.phylip = 1 - 485;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_8278.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/chrZ_8278.phylip = 1 - 485;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_8278.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 (JC69)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_8278.phylip = 1 - 485;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_8278.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal);
END;


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

