Tag: scipy

Среднее взвешенное скользящее среднее в python

У меня есть данные, отобранные по существу случайными интервалами. Я хотел бы вычислить взвешенное скользящее среднее с использованием numpy (или другого пакета python). У меня есть грубая реализация скользящей средней, но у меня возникают проблемы с поиском хорошего способа сделать взвешенную скользящую среднюю, так что значения к центру бункера взвешиваются больше, чем значения по краям. […]

Определение доверительных интервалов для оценки максимального правдоподобия

Я пытаюсь написать код для создания доверительных интервалов для количества разных книг в библиотеке (а также создать информативный сюжет). Мой двоюродный брат учится в начальной школе, и каждую неделю ему дают книгу своего учителя. Затем он читает его и возвращает его вовремя, чтобы получить еще одну на следующей неделе. Через некоторое время мы начали замечать, […]

Внедрение «фильтра Kurtosis» с использованием scipys generic_filter

У меня есть массив 5000*5000 numpy, на котором я хочу рассчитать Kurtosis для окон размером 25. Я попытался использовать функцию generic_filter в generic_filter найденную в ndimage.filters например: import numpy as np from scipy.stats import kurtosis from scipy.ndimage.filters import generic_filter mat = np.random.random_sample((5000, 5000)) kurtosis_filter = generic_filter(mat, kurtosis, size=25, mode='reflect') Это никогда не заканчивается, и я […]

Классифицировать массив Python ближайшим «семенным» регионом?

У меня есть растр экологических местообитаний, который я преобразовал в двумерный массив numpy Python (example_array ниже). У меня также есть массив, содержащий «семенные» регионы с уникальными значениями (seed_array ниже), которые я хотел бы использовать для классификации областей среды обитания. Я бы хотел «вырастить» мои семенные регионы в «области моей среды обитания», чтобы местам обитания присваивался […]

Научные вычисления в Python для программистов MATLAB

Мне было интересно, знает ли кто-нибудь хороший учебник или вводный текст по научным вычислениям на Python / SciPy для программистов MATLAB. Я думал о чем-то в духе текста Дэвида Хибелера на R для программистов MATLAB, который обеспечивает бок о бок ссылку между обоими языками. Чтобы уточнить, моя текущая среда выбора – это Python (x, y) […]

Надежный алгоритм определения ширины пиков

Я спросил, как программно судить диапазоны спектра, и @detly предложил использовать FWHM (полная ширина с половиной максимума), чтобы определить ширину пиков. Я обыскал вокруг и обнаружил, что FWHM можно использовать для подгонки моделей (я на самом деле это непрофессионал!), Особенно модели Guassian. В частности, 2.354 * sigma – это ширина для модели Гуаса. Я смотрю […]

Самый быстрый способ вычислить k наибольших собственных значений и соответствующих собственных векторов с numpy

Я имею большую NxN плотную симметричную матрицу и хочу, чтобы собственные векторы соответствовали k наибольшим собственным значениям. Каков наилучший способ их поиска (желательно, используя numpy, но, возможно, вообще используя blas / atlas / lapack, если это единственный способ пойти)? В общем случае N намного больше k (скажем N> 5000, k <10). Кажется, что у Numpy […]

Как обрабатывать зависимость от scipy в setup.py

Я пытаюсь создать setup.py для проекта, который зависит от SciPy. Следующий setup.py воспроизводит это: setup( name='test', version='0.1', install_requires=['scipy'] ) При установке этого с использованием python setup.py develop он генерирует следующую ошибку: ImportError: No module named numpy.distutils.core Однако, когда я устанавливаю scipy, используя pip , он установил его с колеса, и он работает отлично. Итак, мои […]

Ограниченная линейная регрессия в Python

У меня есть классическая проблема линейной регрессии формы: y = X b где y – вектор отклика X – матрица входных переменных, а b – вектор подходящих параметров, которые я ищу. Python предоставляет b = numpy.linalg.lstsq( X , y ) для решения задач этой формы. Однако, когда я использую это, я имею тенденцию получать либо […]

Как заставить нулевой перехват в линейной регрессии?

Я немного извиняюсь, так что извиняюсь, если на этот вопрос уже был дан ответ, я посмотрел и не смог найти конкретно то, что искал. У меня есть некоторые более или менее линейные данные формы x = [0.1, 0.2, 0.4, 0.6, 0.8, 1.0, 2.0, 4.0, 6.0, 8.0, 10.0, 20.0, 40.0, 60.0, 80.0] y = [0.50505332505407008, 1.1207373784533172, […]

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