Tag: scipy

неустойчивый результат от scipy.cluster.kmeans

Следующий код дает разные результаты в каждой рабочей среде, при кластеризации данных на 3 части с использованием метода k-меток: 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) Были получены три возможных результата: (array([1, 3, 7]), 0.0) (array([3, 7, 1]), 0.0) (array([7, 3, 1]), 0.0) […]

Ограничьте мышь Маяви на вращающуюся Землю вокруг своей оси

Используя iPython Notebook, я смог вывести земной шар Земли с кодом вроде: from mayavi import mlab from mayavi.sources.builtin_surface import BuiltinSurface ocean_blue = (0.4, 0.5, 1.0) r = 6371 # km sphere = mlab.points3d(0, 0, 0, name='Globe', scale_mode='none', scale_factor=r * 2.0, color=ocean_blue, resolution=50) sphere.actor.property.specular = 0.20 sphere.actor.property.specular_power = 10 continents_src = BuiltinSurface(source='earth', name='Continents') continents_src.data_source.on_ratio = 1 […]

Как получить имя распространения из spicy.stats. замороженное распределение?

Доступ к имени замороженного дистрибутива При создании замороженного дистрибутива из пакета scipy.stats , как можно получить доступ к имени дистрибутива после того, как экземпляр распространения будет заморожен? При попытке получить доступ к атрибуту .name возникает ошибка, так как она больше не является атрибутом переменной rv . import scipy.stats as stats # Get the name of […]

Получение квантилей SciPy для соответствия функции Stata xtile

Я унаследовал некоторый старый код Stata (Stata11), который использует функцию xtile для классификации наблюдений в векторе по их квантилям (в данном случае только стандартные 5 квинтилей, 20%, 40%, 60%, 80%, 100% ). Я пытаюсь реплицировать часть кода в Python, и я использую функцию mquantiles() для вычисления. Насколько я могу судить по документации Stata и поиску […]

передача аргументов функции для подгонки

Я пытаюсь установить функцию, которая принимает в качестве входных 2 независимых переменных x, y и 3 параметров, которые будут найдены a, b, c. Это мой тестовый код: import numpy as np from scipy.optimize import curve_fit def func(x,y, a, b, c): return a*np.exp(-b*(x+y)) + cy= x = np.linspace(0,4,50) z = func(x,y, 2.5, 1.3, 0.5) #works ok […]

Функция качания python pandas с двумя аргументами

Любовь моего новичка к python проходит тяжелый процесс … Мне нужно вычислить функцию в скользящем окне фиксированной длины (скажем: 5). Функция требует двух параметров. Я хорошо знаю ответ здесь, который почти идентичен, но я продолжаю получать ошибки. Мой код прост: import numpy as np import pandas as pd import scipy as sp import scipy.stats df […]

Альтернативы сплавленного типа в цитоне

Я работаю над переписыванием модуля python, первоначально написанного на C, используя python-C api для Cython. Модуль также использует NumPy. Главная задача проекта – поддерживать текущую скорость модуля, а также работать для всех типов данных Numpy. Я думаю использовать плавный тип данных, чтобы сделать его общим, но я беспокоюсь из-за его узкого места воздействия на производительность. […]

Использование scipy для минимизации функции, которая также принимает невариантные параметры

Я хочу использовать модуль scipy.optimize для минимизации функции. Скажем, моя функция f(x,a) : def f(x,a): return a*x**2 При фиксированном a я хочу минимизировать f(x,a) по x . С scipy я могу импортировать, например, функцию fmin (у меня есть старый scipy: v.0.9.0), дать начальное значение x0 и затем оптимизировать ( документацию ): from scipy.optimize import fmin […]

Нет модуля с именем scipy.stats – Почему, несмотря на то, что scipy установлен

Как использовать python и scipy для получения случайной переменной poissio? Ничего себе .. Я установил scipy, и для документов, которые я получил, нет модуля с именем scipy.stats? Я на ubuntu 12.04. Итак, иди сюда http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.poisson.html ubuntu@ubuntu:~/Downloads$ sudo apt-get install python-scipy Reading package lists… Done Building dependency tree Reading state information… Done python-scipy is already the […]

Как установить распределение weibull на данные с помощью python?

Я ищу, чтобы найти наилучшие подходящие параметры Weibull для набора данных, используя Python 3.4. import scipy.stats as ss list1 = [] list2 = [] for x in range(0, 10): list1.append(ss.exponweib.pdf(x, a=1, c=2.09, scale=10.895, loc=0)) list2.append(ss.weibull_min.pdf(x, c=2.09, loc=0, scale=10.895)) if list1[x]-list2[x] < .000000001: list1[x]=list2[x] if list1 == list2: print("true") print(ss.distributions.weibull_min.fit(list1, floc=0)) print(ss.distributions.weibull_min.fit(list1, loc=0)) print(ss.distributions.weibull_min.fit(list1, floc=0)) print(ss.distributions.exponweib.fit(list1, […]

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