Tag: scipy

Как получить доступ к разреженным матричным элементам?

type(A) <class 'scipy.sparse.csc.csc_matrix'> A.shape (8529, 60877) print A[0,:] (0, 25) 1.0 (0, 7422) 1.0 (0, 26062) 1.0 (0, 31804) 1.0 (0, 41602) 1.0 (0, 43791) 1.0 print A[1,:] (0, 7044) 1.0 (0, 31418) 1.0 (0, 42341) 1.0 (0, 47125) 1.0 (0, 54376) 1.0 print A[:,0] #nothing returned Теперь я не понимаю, когда я набираю A […]

Любой способ решить систему связанных дифференциальных уравнений в питоне?

Я работал с sympy и scipy, но не могу найти или выяснить, как решить систему связанных дифференциальных уравнений (нелинейных, первого порядка). Так есть ли способ решить связанные дифференциальные уравнения? Уравнения имеют вид: V11'(s) = -12*v12(s)**2 v22'(s) = 12*v12(s)**2 v12'(s) = 6*v11(s)*v12(s) – 6*v12(s)*v22(s) – 36*v12(s) с начальными условиями для v11 (s), v22 (s), v12 (s).

Внекорневая обработка редких массивов CSR

Как можно применить некоторую функцию параллельно на кусках разреженного массива CSR, сохраненного на диске с помощью Python? Последовательно это можно сделать, например, сохраняя массив CSR с помощью joblib.dump открывая его с помощью joblib.load(.., mmap_mode="r") и обрабатывая куски строк один за другим. Это можно сделать более эффективно с помощью dask ? В частности, если предположить, что […]

Установить pyopencv внутри virtualenv с помощью -no-site-packages

Я пытаюсь установить pyopencv в virtualenv, созданный с помощью опции -no-site-packages: pip install pyopencv Но я получаю следующую ошибку на Ubuntu 10.04.3: CMake Error at CMakeLists.txt:186 (find_package): Could not find a configuration file for package OpenCV. Set OpenCV_DIR to the directory containing a CMake configuration file for OpenCV. The file will have one of the […]

нелинейный метод наименьших квадратов python

Я немного из глубины своей математики, связанной с моей проблемой, поэтому прошу прощения за любую неправильную номенклатуру. Я смотрел на использование scipy функции lesssq, но не уверен, что это правильная функция. У меня есть следующее уравнение: eq = lambda PLP,p0,l0,kd : 0.5*(-1-((p0+l0)/kd) + np.sqrt(4*(l0/kd)+(((l0-p0)/kd)-1)**2)) У меня есть данные (8 наборов) для всех членов, кроме kd […]

Обратный сортировка и argsort в python

Я пытаюсь написать функцию в Python (все еще noob!), Которая возвращает индексы и множество документов, упорядоченных внутренними продуктами их оценок tfidf. Процедура такова: Вычислить вектор внутренних продуктов между doc idx и всеми другими документами Сортировать в порядке убывания Верните «баллы» и индексы от второго до конца (т. Е. Не самого себя) Код, который у меня […]

Сплайн с ограничениями на границе

Я измерил данные по трехмерной сетке, например f(x, y, t) . Я хочу интерполировать и сгладить эти данные в направлении t с помощью сплайнов. В настоящее время я делаю это с помощью scipy.interpolate.UnivariateSpline : import numpy as np from scipy.interpolate import UnivariateSpline # data is my measured data # data.shape is (len(y), len(x), len(t)) data […]

Существует ли функция python (scipy) для определения параметров, необходимых для получения целевой мощности?

В R есть очень полезная функция, которая помогает с определением параметров для двухстороннего t-теста для получения целевой статистической мощности. Эта функция называется power.prop.test . http://stat.ethz.ch/R-manual/R-patched/library/stats/html/power.prop.test.html Вы можете вызвать его, используя: power.prop.test(p1 = .50, p2 = .75, power = .90) И он скажет вам размер выборки, необходимый для получения этой мощности. Это чрезвычайно полезно для сдерживания […]

Игнорирование значений -Inf в массивах с использованием numpy / scipy в Python

У меня есть массив NxM в numpy, который я хотел бы записать в журнал, и игнорировать записи, которые были отрицательными до начала регистрации. Когда я беру журнал отрицательных записей, он возвращает -Inf, поэтому в результате у меня будет матрица с некоторыми значениями -Inf. Затем я хочу суммировать столбцы этой матрицы, но игнорируя значения -Inf – […]

ndimage отсутствует из scipy

Я пытаюсь использовать библиотеку ndimage из scipy, но, по-видимому, отсутствует. Я провел тесты как с numpy, так и с scipy, и результаты были в порядке. Я использую numpy 1.6.1 и scipy 0.10.0, установленные из официальных пакетов на sourceforge. Бег import numpy import scipy import pprint print(scipy.version.version) print(numpy.version.version) img = scipy.ndimage.imread("") дает 0.10.0 1.6.1 Traceback (most […]

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