小编Dav*_*rsh的帖子

我完全不熟悉seaborn，因此很抱歉，如果这是一个简单的问题，但是我在文档的任何地方都找不到关于如何在kdeplot中控制由n_levels绘制的级别的描述。这是一个例子：

import seaborn as sns
import numpy as np
import matplotlib.pyplot as plt

x,y=np.random.randn(2,10000)

fig,ax=plt.subplots()
sns.kdeplot(x,y, shade=True,shade_lowest=False, ax=ax,n_levels=3,cmap="Reds")
plt.show()

这是结果图：

我希望能够知道显示了哪些置信度，这样我就可以将图标记为“阴影区域显示（a，b，c）百分比置信区间”。我会天真地假设n_levels以某种方式与高斯中的等效“ sigmas”相关，但是从示例看来，情况并非如此。

理想情况下，我希望能够通过将元组传递给kdeplot来指定显示的间隔，例如：

levels=[68,95,99]

并绘制这些置信区域。

编辑：感谢@Goyo和@tom，我想我可以澄清我的问题，并半路寻找我要的答案。如前所述，n_levels传递给plt.cotourfas levels，因此可以传递列表。但是sns.kdeplot绘制PDF时，PDF中的值不对应于我要寻找的置信区间（因为这些对应于PDF的积分）。

我需要做的是通过sns.kdeplot在x,y集成（和标准化）PDF的值，然后我就可以进入如n_levels=[0.68,0.95,0.99,1]。

编辑2：我现在已经解决了这个问题。见下文。我使用二维规范化直方图定义置信区间，然后将其作为级别传递给规范化kde图。抱歉，我可以制作一个返回级别的函数，但我将所有内容明确键入。

import numpy as np
import scipy.optimize
import matplotlib.pyplot as plt
import seaborn as sns

# Generate some random data
x,y=np.random.randn(2,100000)

# Make a 2d normed histogram
H,xedges,yedges=np.histogram2d(x,y,bins=40,normed=True)

norm=H.sum() # Find the norm of the sum
# …