Articles of numpy scipy

Scipy expit: comportamiento inesperado. NaNs

Noté que algunos nanos aparecían inesperadamente, en mis datos. (y expandiendo y ocultando todo lo que tocaron) Hicieron una investigación cuidadosa y produjeron un ejemplo de trabajo mínimo: >>> import numpy >>> from scipy.special import expit >>> expit(709) 1.0 >>> expit(710) nan Expit es el logit inverso. Documentación Scipy aquí . Lo que nos dice: […]

Elemento sabio poder de la matriz scipy.sparse

¿Cómo puedo elevar una matriz scipy.sparse a un poder, de forma elemental? numpy.power debería, según su manual , hacer esto, pero falla en matrices dispersas: >>> X <1353×32100 sparse matrix of type '’ with 144875 stored elements in Compressed Sparse Row format> >>> np.power(X, 2) Traceback (most recent call last): File “”, line 1, in […]

punto scipy.sparse extremadamente lento en Python

El siguiente código ni siquiera terminará en mi sistema: import numpy as np from scipy import sparse p = 100 n = 50 X = np.random.randn(p,n) L = sparse.eye(p,p, format=’csc’) XTdot(L).dot(X) ¿Hay alguna explicación de por qué esta multiplicación de matrices está colgando?

¿Cómo usar scipy.integrate para obtener el volumen de una esfera truncada?

Estoy luchando con el uso de scipy.integrate, usé tplquad, pero ¿cómo puedo usar la integrate para obtener el volumen de la esfera (truncada)? Gracias import scipy from scipy.integrate import quad, dblquad, tplquad from math import* from numpy import * R = 0.025235 #radius theta0 = acos(0.023895) #the angle from the edge of truncated plane to […]

Sustituya la difusión numpy usando scipy.sparse.csc_matrix

Tengo en mi código la siguiente expresión: a = (b / x[:, np.newaxis]).sum(axis=1) donde b es un ndarray de forma (M, N) , y x es un ndarray de forma (M,) . Ahora, b es en realidad escasa, por lo que para la eficiencia de la memoria me gustaría sustituir en scipy.sparse.csc_matrix o csr_matrix . […]

minimizando una función multivariable y diferenciable usando scipy.optimize

Estoy intentando minimizar la siguiente función con scipy.optimize : cuyo gradiente es este (Para aquellos que estén interesados, esta es la función de probabilidad de un modelo de Bradley-Terry-Luce para comparaciones por pares. Muy relacionado con la regresión logística). Es bastante claro que agregar una constante a todos los parámetros no cambia el valor de […]

La implementación sin problemas de MATLAB (promedio móvil de n puntos) en NumPy / Python

La función smooth de Matlab, por defecto, suaviza los datos utilizando una media móvil de 5 puntos. ¿Cuál sería la mejor manera de hacer lo mismo en python? Por ejemplo, si esta es mi información. 0 0.823529411764706 0.852941176470588 0.705882352941177 0.705882352941177 0.676470588235294 0.676470588235294 0.500000000000000 0.558823529411765 0.647058823529412 0.705882352941177 0.705882352941177 0.617647058823529 0.705882352941177 0.735294117647059 0.735294117647059 0.588235294117647 0.588235294117647 1 0.647058823529412 […]

¿Cómo instalar numpy y scipy en OS X?

Hola chicos, soy nuevo en Mac, así que tengan paciencia conmigo. Estoy usando el leopardo de nieve 10.6.4 en este momento. Quiero instalar numpy y scipy, así que descargué los archivos python2.6, numpy y scipy dmg de su sitio oficial. Sin embargo, estoy teniendo problemas para importar números: Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/numpy/core/multiarray.so: no matching architecture in universal wrapper […]

Iterando sobre la matriz escasa dispersa por columna

Estoy tratando de averiguar cómo iterar a través de una matriz escasa dispersa por columna. Estoy tratando de calcular la sum de cada columna, luego ponderar los miembros de esa columna por esa sum. Lo que quiero hacer es básicamente: for i=0 to #columns for j=0 to #rows sum=sum+matrix[i,j] for j=0to #rows matrix[i,j]=matrix[i,j]/sum Todos los […]

¿Cómo leer un sistema de ecuaciones diferenciales de un archivo de texto para resolver el sistema con scipy.odeint?

Tengo un sistema grande (> 2000 ecuaciones) de EDO que quiero resolver con el método de Python Scipy. Tengo tres problemas que quiero resolver (¿quizás tendré que hacer 3 preguntas diferentes?). Para simplificar, los explicaré aquí con un modelo de juguete, pero tenga en cuenta que mi sistema es grande. Supongamos que tengo el siguiente […]