【问题标题】:How to plot error ellipses from covariance matrix如何从协方差矩阵中绘制误差椭圆
【发布时间】:2017-12-01 16:53:41
【问题描述】:

我用一些参数拟合了我的数据,现在我得到了它们的最优值 popt 和它们的协方差矩阵 pcov。我假设这些不确定性是高斯的,我想绘制如下图所示的误差椭圆图。有没有一些简单的方法可以做到(无需手动完成每个子图)?我正在查看pandas.tools.plotting.scatter_matrix,它似乎有点接近,但我不知道如何将其应用于我的多元高斯。

【问题讨论】:

  • 请提供可重现的小代码
  • 这个情节不是我的(我在某个地方找到的)——这是我想要实现的,我正在寻找这样的代码
  • 这可能是一种误解。 SO 不是代码编写服务。假设您已经尝试过某些事情,您应该以链接到其他没有帮助的问题的形式来说明它是什么,并提供失败尝试的minimal reproducible example。见How to AskSO checklist 也可能有所帮助。
  • 我不是在找人为我编写代码,而是为我指明正确的方向。我花了一些时间寻找一些东西,壁橱里的东西是我提到的pandas 包。但是因为我从不使用它,所以我不确定我想要实现的目标是否已经写在包中。我假设更有经验的用户可以告诉我这是否容易——这是我的问题。
  • 没有神奇的plot_my_errorellipses_from_data_to_grid函数。那么这应该以“过于宽泛”还是“离题(因为要求软件推荐)”而结束?

标签: python pandas matplotlib


【解决方案1】:

我找到了这个corner.py 包,文档here,它可以满足我的需要。

简单示例:

import corner
import numpy as np

ndim, nsamples = 2, 10000
np.random.seed(42)
samples = np.random.randn(ndim * nsamples).reshape([nsamples, ndim])
figure = corner.corner(samples)

【讨论】:

    【解决方案2】:

    如果您使用Seaborn,您可以尝试创建PairGrid 或使用pairplot() 函数。 PairGrid 更灵活,而 pairplot() 是一种使用默认散点矩阵查看某些数据的快速方法。

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 2019-02-15
      • 2021-05-21
      • 1970-01-01
      • 1970-01-01
      • 2016-05-03
      • 2012-12-09
      • 1970-01-01
      • 2014-03-06
      相关资源
      最近更新 更多