我使用gaussian_kde生成二维密度估计。有没有办法可以计算这个PDF的水平曲线，其具有固定数量的概率质量高于给定高度的属性？
例如，我想确定表面z = z_0下面的PDF体积等于0.5的值z_0。用root或fsolve解决这个问题的简单尝试通常不会收敛。
我已经尝试通过在覆盖大部分空间的网格网格上评估kde来形成一个z值数组，然后使用root和/或fsolve来查找需要a的函数的根值z_0并返回超过z_0的z值之和与期望概率p之间的差值。尽管此功能是平滑且单调的，但root和fsolve并不始终如一地融合在唯一的解决方案上。