Python scipy.stats：多元正态拟合？

我使用MCMC方法对数据进行了分析。这是仅针对两个参数的结果示例（注意：我的模型有4个参数）：

我想做一个多元拟合，但是与单变量分布不同，

scipy.stats

包中的

multivariable\u normal

分布似乎没有实现

fit（）

方法。 我应该使用哪个python包来实现这一点

---

谢谢你的回答。

这能回答你的问题吗@马克斯皮耶里尼。它部分地回答了这个问题。它适用于正态分布，但我正在寻找一个更通用的多元分布包。

fit

方法（）仅用于单变量分布。