Articles of entrometido

resultado inestable de scipy.cluster.kmeans

El siguiente código proporciona diferentes resultados en cada tiempo de ejecución al agrupar los datos en 3 partes utilizando el método k significa: from numpy import array from scipy.cluster.vq import kmeans,vq data = array([1,1,1,1,1,1,3,3,3,3,3,3,7,7,7,7,7,7]) centroids = kmeans(data,3,100) #with 100 iterations print (centroids) Tres posibles resultados obtenidos fueron: (array([1, 3, 7]), 0.0) (array([3, 7, 1]), 0.0) […]

Construyendo una matriz de filas ‘rodadas’ eficientemente en Numpy

Me gustaría construir una matriz (n,n) partir de una matriz unidimensional, donde cada fila se desplaza (con ajuste) una en relación con la anterior. El siguiente código hace esto: import numpy as np r = np.array([1, 2, 3, 4, 5]) n = len(r) MM = np.zeros((n, n), dtype=r.dtype) for k in range(n): MM[k, :] = […]

Cómo cargar parcialmente una matriz guardada con numpy guardar en python

Tengo una matriz multidimensional con numpy save y solo quiero cargar parcialmente alguna dimensión porque la matriz es muy grande. ¿Cómo puedo hacerlo de forma sencilla? Editar: El contexto es simple y básico: Tienes una matriz de 5 Gb guardada con numpy.save . Pero, solo necesita tener acceso a algunas partes de la matriz A[:,:] […]

Ajustar una curva a una distribución de ley de potencia con curve_fit no funciona

Estoy tratando de encontrar una curva que se ajuste a mis datos que visualmente parecen tener una distribución de la ley de poder. Esperaba utilizar scipy.optimize.curve_fit, pero no importa qué función o normalización de datos intente, obtengo un RuntimeError (parámetros no encontrados o desbordados) o una curva que no se ajusta a mis datos ni […]

¿Kurtosis, sesgo de un gráfico de barras? – Python

¿Cuál es un método eficiente para determinar el sesgo / kurtosis de un gráfico de barras en python? Teniendo en cuenta que los gráficos de barras no están agrupados (a diferencia de los histogtwigs), esta pregunta no tendría mucho sentido, pero lo que estoy tratando de hacer es determinar la simetría de la altura frente […]

Scipy.optimize.root no converge en Python mientras que Matlab fsolve funciona, ¿por qué?

Estoy tratando de encontrar la raíz y de una función llamada f usando Python. Aquí está mi código: def f(y): w,p1,p2,p3,p4,p5,p6 = y[:7] t1 = w – 0.99006633*(p1**0.5) – (-1.010067)*((1-p1)) t2 = w – 22.7235687*(p2**0.5) – (-1.010067)*((1-p2)) t3 = w – 9.71323491*(p3**0.5) – (-1.010067)*((1-p3)) t4 = w – 2.43852877*(p4**0.5) – (-1.010067)*((1-p4)) t5 = w – […]

Cómo leer en una lista de bordes para hacer una matriz dispersa escéptica

Tengo un archivo grande donde cada línea tiene un par de cadenas de 8 caracteres. Algo como: ab1234gh iu9240gh en cada linea Este archivo realmente representa un gráfico y cada cadena es un ID de nodo. Me gustaría leer en el archivo y hacer directamente una matriz de adyacencia escasa dispersa. Luego ejecutaré PCA en […]

Resolviendo la ecuación diferencial de matrices en Python usando Scipy / Numpy-NDSolve equivalente?

Tengo dos matrices numpy: 9×9 y 9×1. Me gustaría resolver la ecuación diferencial en puntos de tiempo discretos, pero tengo problemas para hacer que ODEInt funcione. No estoy seguro de si estoy haciendo lo correcto. Con Mathematica, la ecuación es: Solution = {A[t]} /. NDSolve[{A'[t] == Ab.A[t] && A[0] == A0}, {A[t]}, {t, 0, .5}, […]

Distancia de Hausdorff entre grillas 3D

Tengo varias cuadrículas (matrices numpy [Nk, Ny, Nx]) y me gustaría usar la distancia de Hausdorff como una métrica de similitud de estas cuadrículas. Hay varios módulos en scipy (scipy.spatial.distance.cdist, scipy.spatial.distance.pdist) que permiten calcular la distancia euclidiana entre matrices 2D. Ahora, para comparar las cuadrículas, tengo que elegir una sección transversal (por ejemplo, cuadrícula1 [0 […]

error de tipos no ordenados al importar sklearn

Instalé numpy (1.12.0b1), Scipy (0.18) en Windows. También instalé sci-kit también. Cuando escribí “importar sklearn” en la consola de Python, se produce un error como este: si np_version <(1, 12, 0): TypeError: tipos desordenados: str () <int () ¿Cuál será el problema?