Articles of numpy

La forma más rápida de generar cadenas delimitadas a partir de una matriz numpy 1d

Tengo un progtwig que necesita convertir muchas matrices de flotadores grandes y unidimensionales de gran tamaño en cadenas delimitadas. Estoy encontrando esta operación bastante lenta en relación con las operaciones matemáticas en mi progtwig y me pregunto si hay una manera de acelerarlo. Por ejemplo, considere el siguiente bucle, que toma 100,000 números aleatorios en […]

Statsmodels: Calcular valores ajustados y R al cuadrado

Estoy ejecutando una regresión de la siguiente manera ( df es un dataframe de pandas ): import statsmodels.api as sm est = sm.OLS(df[‘p’], df[[‘e’, ‘varA’, ‘meanM’, ‘varM’, ‘covAM’]]).fit() est.summary() Lo que me dio, entre otros, un R cuadrado de 0.942 . Entonces quise trazar los y-values y originales y-values ajustados. Para esto, ordené los valores […]

Pandas: restando dos columnas de fecha y el resultado es un entero

Tengo dos columnas en un dataframe de Pandas que son fechas. Busco restar una columna de otra y el resultado es la diferencia en números de días como un entero . Un vistazo a los datos: df_test.head(10) Out[20]: First_Date Second Date 0 2016-02-09 2015-11-19 1 2016-01-06 2015-11-30 2 NaT 2015-12-04 3 2016-01-06 2015-12-08 4 NaT […]

Expresión elegante para el producto de puntos por filas de dos matrices.

Tengo dos matrices numpy 2-d con las mismas dimensiones, A y B, y estoy tratando de calcular el producto puntual de las filas. Yo podría hacer: np.sum(A * B, axis=1) ¿Hay otra manera de hacer esto para que numpy esté haciendo el producto de puntos por filas en un paso en lugar de dos? Tal […]

Python MemoryError al realizar el ajuste con Scikit-learn

Estoy ejecutando Python 2.7 (64 bits) en un sistema Windows 8 de 64 bits con 24 GB de memoria. Al realizar el ajuste de los Sklearn.linear_models.Ridge habituales, el código funciona bien. Problema: Sin embargo, cuando uso Sklearn.linear_models.RidgeCV(alphas=alphas) para el ajuste, me encuentro con el error MemoryError se muestra a continuación en la línea rr.fit(X_train, y_train) […]

El ejemplo de Numpy C-Api da un SegFault

Estoy tratando de entender cómo funciona Python C-Api, y quiero intercambiar arrays numpy entre Python y una extensión C Entonces, comencé este tutorial: http://dsnra.jpl.nasa.gov/software/Python/numpydoc/numpy-13.html Intenté hacer el primer ejemplo allí, un módulo C que calcula la traza de una matriz numpy 2d, fue muy bueno para mí, ya que también quiero hacer operaciones elementales en […]

Usando mkl_set_num_threads con numpy

Estoy tratando de establecer el número de subprocesos para cálculos mkl_set_num_threads con mkl_set_num_threads como este import numpy import ctypes mkl_rt = ctypes.CDLL(‘libmkl_rt.so’) mkl_rt.mkl_set_num_threads(4) pero sigo teniendo una falla de segmentación: Program received signal SIGSEGV, Segmentation fault. 0x00002aaab34d7561 in mkl_set_num_threads__ () from /../libmkl_intel_lp64.so Obtener el número de hilos no es un problema: print mkl_rt.mkl_get_max_threads() ¿Cómo puedo […]

Invirtiendo una matriz booleana numpy usando ~

¿Puedo usar ~A para invertir una matriz numpy de valores booleanos, en lugar de las funciones algo incómodas np.logical_and() y np.invert() ? De hecho, ~ parece funcionar bien, pero no puedo encontrarlo en ningún manual de referencia de nympy, y, lo que es más alarmante, ciertamente no funciona con los escalares (por ejemplo, bool(~True) devuelve […]

¿La configuración de arrays numpy a Ninguna memoria libre?

Tengo cientos de matrices realmente grandes, como las de (600, 800) o (3, 600, 800). Por lo tanto, quiero desasignar la memoria utilizada tan pronto como ya no necesite algo. Pensé: some_matrix = None ¿Debería hacer el trabajo o solo la referencia está establecida en Ninguna, pero en algún lugar de la Memoria el espacio […]

Comparando matrices NumPy para que los NaN se comparen iguales

¿Existe una forma idiomática de comparar dos matrices NumPy que traten a los NaN como iguales entre sí (pero no iguales a nada que no sea ​​un NaN)? Por ejemplo, quiero que las dos matrices siguientes se comparen igual: np.array([1.0, np.NAN, 2.0]) np.array([1.0, np.NAN, 2.0]) y las siguientes dos matrices para comparar desiguales: np.array([1.0, np.NAN, […]