Seaborn jointplot定制带宽

2024-10-16 17:20:03 发布

您现在位置:Python中文网/ 问答频道 /正文

我尝试使用Seaborn的jointplot函数来创建一个二维核密度估计。我以前用过R,用'ks'包计算插件带宽,计算KDE如下(我不想使用“经验法则”方法):

library(ks)
data(unicef)
x_vals <- unicef[sample(nrow(unicef), 10), 1]
y_vals <- unicef[sample(nrow(unicef), 10), 1]
mydf <- data.frame(x_vals, y_vals)
H <- Hpi(mydf)
fhat <- kde(x=mydf, H=H)
plot(fhat)

上面的代码段返回带宽H的2x2矩阵

我想在Seaborn的jointplot函数中使用相同的带宽计算。两个问题:i)是否有一种以Python为中心的方法来计算带宽,就像上面的R方法一样?;ii)一旦我有了带宽,如何将其传递给joinpoint,因为它接受的是一个值列表,而不是一个矩阵?

谢谢


Tags: sample方法函数插件data矩阵seabornks