Проблема с панелью ошибок Python matplotlib

Учитывая эти массивы numpy

x = [0 1 2 3 4 5 6 7 8 9]
y = [[ 0.        ]
     [-0.02083473]
     [ 0.08819923]
     [ 0.9454764 ]
     [ 0.80604627]
     [ 0.82189822]
     [ 0.73613942]
     [ 0.64519742]
     [ 0.56973868]
     [ 0.612912  ]]
c = [[ 0.          0.        ]
     [-0.09127286  0.04960341]
     [-0.00300709  0.17940555]
     [ 0.82319693  1.06775586]
     [ 0.74512774  0.8669648 ]
     [ 0.75177669  0.89201975]
     [ 0.63606087  0.83621797]
     [ 0.57786173  0.7125331 ]
     [ 0.46722312  0.67225423]
     [ 0.54951714  0.67630685]]

Я хочу построить график x,y с планками погрешностей, используя значения в c. Я старался

plt.errorbar(x, y, yerr=c)

Но интерпретатор дает мне эту ошибку:

File "C:\Python\32\lib\site-packages\matplotlib\axes.py", line 3846, in vlines
  for thisx, (thisymin, thisymax) in zip(x,Y)]
File "C:\Python\32\lib\site-packages\matplotlib\axes.py", line 3846, in <listcomp>
  for thisx, (thisymin, thisymax) in zip(x,Y)]
ValueError: too many values to unpack (expected 2)

Значение x в zip равно

[0 1 2 3 4 5 6 7 8 9]

и значение Y в zip равно

[[[ 0.          0.        ]
  [ 0.07043814 -0.11210759]
  [ 0.09120632  0.08519214]
  [ 0.12227947  1.76867333]
  [ 0.06091853  1.55117401]
  [ 0.07012153  1.57367491]
  [ 0.10007855  1.3722003 ]
  [ 0.06733568  1.22305915]
  [ 0.10251555  1.0369618 ]
  [ 0.06339486  1.16242914]]

 [[ 0.          0.        ]
  [-0.07043814  0.02876869]
  [-0.09120632  0.26760478]
  [-0.12227947  2.01323226]
  [-0.06091853  1.67301107]
  [-0.07012153  1.71391797]
  [-0.10007855  1.57235739]
  [-0.06733568  1.35773052]
  [-0.10251555  1.2419929 ]
  [-0.06339486  1.28921885]]]

Я почитал, и похоже, что мой код должен быть правильным (глупое предположение, но я не могу найти доказательства обратного... пока), но похоже, что errorbar это не так. как 2d массив. В документации сказано, что yerr может быть двумерным массивом, при этом первый столбец представляет собой минимальную ошибку, а второй — максимальную.

Что я здесь делаю не так?


person JesseBuesking    schedule 30.01.2013    source источник


Ответы (1)


Были некоторые проблемы с кодом, который я исправил ниже, поэтому он работает без проблем.

import numpy
import pylab

arr = numpy.asarray

x = arr([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])   # put comma between numbers
y = arr([[ 0.        ],                   # make it vector
     [-0.02083473],
     [ 0.08819923],
     [ 0.9454764 ],
     [ 0.80604627],
     [ 0.82189822],
     [ 0.73613942],
     [ 0.64519742],
     [ 0.56973868],
     [ 0.612912  ]]).flatten()
c = arr([[ 0.        ,  0.        ],
     [-0.09127286,  0.04960341],
     [-0.00300709,  0.17940555],
     [ 0.82319693,  1.06775586],
     [ 0.74512774,  0.8669648 ],
     [ 0.75177669,  0.89201975],
     [ 0.63606087,  0.83621797],
     [ 0.57786173,  0.7125331 ],
     [ 0.46722312,  0.67225423],
     [ 0.54951714,  0.67630685]]).T      # transpose
pylab.errorbar(x, y, yerr=c)
pylab.show()

и результат:

введите здесь описание изображения

Удачи.

person Developer    schedule 30.01.2013
comment
Ага! Следует отметить, что передаваемые значения ошибок являются смещениями от фактического значения. Доверительные отношения, которые я предоставил, были фактическими значениями y (отсюда и огромные планки погрешностей). - person JesseBuesking; 30.01.2013