Tag: scipy

создание errorbars не обрезается в matplotlib с Python

Я использую matplotlib в Python для построения строки с errorbars следующим образом: plt.errorbar(xvalues, up_densities, yerr=ctl_sds, fmt='-^', lw=1.2, markersize=markersize, markeredgecolor=up_color, color=up_color, label="My label", clip_on=False) plt.xticks(xvalues) Я установил тики по оси х, используя «xticks». Тем не менее, строки ошибок последней точки в значениях x (т.е. xvalues ​​[-1]) обрезаются справа – это означает, что отображается только половина строки […]

Эффективно вычислить 50 Гц содержимого сигнала

Постановка задачи У меня длинный сигнал (454912 выборок), и я хотел бы вычислить оценку 50Гц в нем. Скорость здесь важнее точности. Ожидается, что количество 50 Гц будет меняться со временем. Значение должно быть репрезентативным для всей записи, например, для среднего значения. контекст Сигнал записывается с электрода ЭЭГ. Когда электроды ЭЭГ плохо контактируют с кожей головы, […]

Как вычислить остатки точечного процесса в python

Я пытаюсь воспроизвести работу из http://jheusser.github.io/2013/09/08/hawkes.html в python, за исключением разных данных. Я написал код, чтобы имитировать процесс Пуассона, а также процесс Хоукса, который они описывают. Чтобы сделать модель Hawkes MLE, я определяю функцию логарифмического правдоподобия как def loglikelihood(params, data): (mu, alpha, beta) = params tlist = np.array(data) r = np.zeros(len(tlist)) for i in xrange(1,len(tlist)): […]

Интерпретация значений scipy.stats.entropy

Я пытаюсь использовать scipy.stats.entropy для оценки расхождения Kullback-Leibler (KL) между двумя дистрибутивами. В частности, я хотел бы использовать KL как метрику, чтобы решить, насколько совместимы два дистрибутива. Однако я не могу интерпретировать значения KL. Например: t1 = numpy.random.normal (-2.5,0.1,1000) t2 = numpy.random.normal (-2.5,0.1,1000) scipy.stats.entropy (t1, t2) +0,0015539217193737955 Затем, t1 = numpy.random.normal (-2.5,0.1,1000) t2 = numpy.random.normal […]

Python / Cython: использование SciPy с Cython

Учебник Cython показывает хороший пример использования Numpy с Cython. Тем не менее, у меня есть код, который использует пакет scipy.stats, и при попытке скомпилировать код I ошибки, такие как: dvi.pyx:7:8: 'scipy.stats.pxd' not found Я боюсь, что Scipy не поддерживается Cython (?). Может кто-нибудь прокомментировать использование scipy с Cython или указать мне в сторону некоторых ресурсов […]

Как использовать dorpi5 или dop853 в Python

Я просмотрел scipy.integrate.ode, но я не могу узнать, как использовать эти методы интеграции, dorpi5 и dop853 . Я хотел бы попытаться интегрировать интеграцию ode python и mathematica в мой код python с этими двумя методами, чтобы увидеть, как это влияет на результаты, но не знает, как это сделать.

Нарезка массивных массивов

У меня есть массив 1D numpy и некоторые значения смещения / длины. Я хотел бы извлечь из этого массива все записи, которые попадают в смещение, offset + length, которые затем используются для создания нового «уменьшенного» массива из исходного, который состоит только из тех значений, которые были выбраны парами смещения / длины , Для одной пары […]

Дискретный лапласиан (эквивалент del2) в Python

Мне нужен эквивалент Python / Numpy для Matlab (Octave) дискретного оператора Лапласа (функция) del2 (). Я попробовал пару решений Python, ни один из которых, похоже, не соответствует выходу del2. На Октаве у меня есть image = [3 4 6 7; 8 9 10 11; 12 13 14 15;16 17 18 19] del2(image) это дает результат […]

Индексирование массива python с массивом python с избыточными элементами

У меня возникла проблема с индексированием массива. Предположим, у вас есть массив a и другой массив b, который вы хотите использовать для использования в качестве индекса для a, чтобы присвоить некоторые значения позиции, обозначенной b элементами. a=numpy.zeros(5) print a [ 0. 0. 0. 0. 0.] Теперь я хотел бы увеличить второй элемент дважды b=numpy.array([1,1]) a[b]+=1. […]

Многопанельные временные ряды линий и заполненные контуры с использованием matplotlib?

Если бы я хотел создать комбинированное изображение, подобное показанному ниже ( исходный источник здесь ), не могли бы вы указать мне объекты matplotlib, которые мне нужно собрать? Я пытаюсь работать с объектами AxesImage, и я также загрузил SciKits Timeseries, но мне это нужно, или может ли это так же просто использовать strptime, mktime и strftime […]

Python - лучший язык программирования в мире.