Matplotlib - 根据 log10 值计算的箱线图，但以对数刻度显示

Question

我认为这是一个简单的问题，但我似乎仍然想不出一个简单的解决方案。我有一组分子丰度数据，其值范围为多个数量级。我想用boxplots (box-and-whiskers plots) 来表示这些丰度，并且由于值的范围很广，我希望在对数刻度上计算这些框。 我知道我可以只计算数据的 log10 并将其发送到 matplotlib 的boxplot，但这不会在以后的图中保留对数刻度。

所以我的问题基本上是这样的： 当我根据我的值的 log10 计算出箱线图时，我如何转换之后的绘图以对数刻度显示，而不是与 log10 值呈线性关系？ 我可以更改刻度标签以部分解决此问题，但我不知道如何将对数刻度带回绘图。

或者有另一种更直接的方法来绘制这个。可能已经包含此选项的其他软件包？

非常感谢您的帮助。

Answer 1

我建议不要在原始值上绘制箱线图并将 y 轴设置为对数，因为箱线图函数并非设计为跨数量级工作，并且您可能会得到太多异常值（取决于您的数据，当然）。

相反，您可以绘制数据的对数并手动调整 y 标签。

这是一个非常粗略的例子：

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.ticker import MultipleLocator, FormatStrFormatter

np.random.seed(42)

values = 10 ** np.random.uniform(-3, 3, size=100)

fig = plt.figure(figsize=(9, 3))


ax = plt.subplot(1, 3, 1)

ax.boxplot(np.log10(values))
ax.set_yticks(np.arange(-3, 4))
ax.set_yticklabels(10.0**np.arange(-3, 4))
ax.set_title('log')

ax = plt.subplot(1, 3, 2)

ax.boxplot(values)
ax.set_yscale('log')
ax.set_title('raw')

ax = plt.subplot(1, 3, 3)

ax.boxplot(values, whis=[5, 95])
ax.set_yscale('log')
ax.set_title('5%')

plt.show()

右图显示了原始值的箱线图。这会导致许多异常值，因为最大晶须长度计算为四分位间距（盒子高度）的倍数（默认值：1.5），不会跨越数量级。

或者，您可以指定在给定的百分位范围内绘制胡须： ax.boxplot(values, whis=[5, 95]) 在这种情况下，您会得到固定数量的异常值 (5%) 高于和低于。

Answer 2

你可以使用 plt.yscale:

plt.boxplot(data); plt.yscale('log')