Tag: scipy

Стандартная ошибка нелинейной регрессии

Я занимаюсь физическими симуляторами в Монте-Карло с Python, и я не могу определить стандартную ошибку для коэффициентов нелинейного наименьшего квадрата. Первоначально я использовал scipy.stats.linregress SciPy для своей модели, так как я думал, что это будет линейная модель, но заметил, что это на самом деле какая-то функция власти. Затем я использовал polyfit NumPy со степенями свободы, […]

«Покраска» одного массива на другой с использованием python / numpy

Я пишу библиотеку для обработки слежения за глазами в Python, и я довольно новичок во всем мире numpy / scipy. По сути, я хочу взять массив значений (x, y) во времени и «нарисовать» некоторую фигуру на холсте в этих координатах. Например, форма может быть размытым кругом. Операция, которую я имею в виду, более или менее […]

Нелинейная регрессия e ^ (- x) с использованием scipy, python, numpy

Код ниже дает мне плоскую линию для линии наилучшего соответствия, а не хорошую кривую вдоль модели e ^ (- x), которая соответствовала бы данным. Может ли кто-нибудь показать мне, как исправить код ниже, чтобы он соответствовал моим данным? import numpy as np import matplotlib.pyplot as plt import scipy.optimize def _eNegX_(p,x): x0,y0,c,k=py = (c * np.exp(-k*(x-x0))) […]

numpy: эффективное чтение большого массива

У меня есть двоичный файл, содержащий плотную матрицу n*m из 32-битных поплавков. Каков наиболее эффективный способ прочитать его в массив numpy упорядоченный по Fortran? Файл имеет размер в несколько гигабайт. Я получаю управление форматом, но он должен быть компактным (т. Е. Длиной около 4*n*m байтов) и должен быть легко создан из кода, отличного от Python. […]

2D Сортировка с NumPy – сортировка 1 строка, а другая – сортировка

Скажем, у меня есть массив NumPy: [[4 9 2] [5 1 3]] Я хочу отсортировать нижнюю строку этого массива, но верхняя строка будет следовать за сортировкой, так что я получаю: [[9 2 4] [1 3 5]] Я знаю, что вы можете сортировать, как это, используя функцию sorted (), но для этого требуется ввод и вывод […]

Геометрическая медиана многомерных точек

У меня есть массив трехмерных точек: a = np.array([[2., 3., 8.], [10., 4., 3.], [58., 3., 4.], [34., 2., 43.]]) Как я могу вычислить геометрическую медиану этих точек?

Пуассоновский доверительный интервал с numpy

Я пытаюсь поместить пуассоновские строки ошибок на гистограмме, которую я создаю с помощью matplotlib, но я не могу найти функцию numpy, которая даст мне 95% -ный доверительный интервал, предполагающий данные о пуассонии. В идеале решение не зависит от scipy, но все будет работать. Существует ли такая функция? Я много узнал о загрузке, но в моем […]

Как создать белое изображение в Python?

После выполнения моей домашней работы я наткнулся на проблему, связанную с Python и манипуляциями с изображениями. Должен сказать, использование образа lib не является вариантом. Итак, здесь from scipy.misc import imread,imsave from numpy import zeros imga = zeros([100,100,3]) h = len(imga) w = len(imga[0]) for y in range(h): for x in range(w): imga[y,x] = [255,255,255] imsave("Result.jpg",imga) […]

Сравнение цепочек цепей с двумя предикатами

В Numpy я могу генерировать логический массив следующим образом: >>> arr = np.array([1, 2, 1, 2, 3, 6, 9]) >>> arr > 2 array([False, False, False, False, True, True, True], dtype=bool) Можно ли объединить сравнения? Например: >>> 6 > arr > 2 array([False, False, False, False, True, False, False], dtype=bool)

эквивалент python qnorm, qf и qchi2 из R

Мне нужен квантиль некоторых распределений в python. В r можно вычислить эти значения с помощью функций qf, qnorm и qchi2. Существует ли какой-либо эквивалент Python этих R-функций? Я смотрел на scipy, но я ничего не нашел.

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