scipy.stats：多元正态拟合？

Question

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

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

fit()

方法似乎没有针对

multivariate_normal

包中的

scipy.stats

分布实现。 我应该使用哪个 python 包来做到这一点？

谢谢您的回答。

Answer 1

今年夏天，

scipy/scipy#18986 在

fit

中添加了

scipy.stats.multivariate_normal

方法。该功能将在 SciPy 1.12.0 中提供，计划于今年年底左右发布。同时，您可以通过从源代码构建或使用夜间轮子来使用该功能。如果您从发布下拉菜单中选择“开发”，您可以在 fit 文档

 中找到 

multivariate_normal 的条目。