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) […]
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, :] = […]
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[:,:] […]
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 […]
¿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 […]
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 – […]
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 […]
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}, […]
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 […]
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?