Я прочитал файл Matlab, содержащий большое количество массивов в качестве набора данных, в Python, сохраняющий словарь Matlab под макетом имен переменных, используя команду: mat = loadmat('Sample Matlab Extract.mat') Есть ли способ, которым я могу использовать функцию Python для записи в функции csv, чтобы сохранить эту переменную словаря Matlab, которую я читал в Python, в виде […]
У меня есть кривая в 3D-пространстве. Я хочу использовать кусочно-кубическую интерполяцию, сохраняющую фигуру, на ней аналогично pchip в matlab. Я исследовал функции, предоставляемые в scipy.interpolate, например interp2d, но функции работают для некоторых структур кривой, а не для точек данных, которые у меня есть. Любые идеи о том, как это сделать? Вот данные: x,y,z 0,0,0 0,0,98.43 […]
У меня есть нелинейная модель, которая выглядит так: Темная сплошная линия подходит для модели, а серая часть – это сырые данные. Краткая версия вопроса: как я могу получить вероятность соответствия этой модели, поэтому я могу выполнить тест отношения правдоподобия? Предположим, что остаток нормально распределен. Я относительно новичок в статистике, и мои текущие мысли таковы: Получите […]
Я работаю над Windows-машиной и версией python версии 2.7. Я уже установил numpy, и теперь я пытаюсь установить scipy. Я попытался установить его, используя команду easy install, и загрузив zip-файл, а затем запустив файл setup.py. Когда я пытаюсь установить его с помощью easy_install c:\Python27\Scripts>easy_install scipy Searching for scipy Reading https://pypi.python.org/simple/scipy/ Best match: scipy 0.13.3 Downloading […]
Учитывая следующую марковскую матрицу: import numpy, scipy.linalg A = numpy.array([[0.9, 0.1],[0.15, 0.85]]) Стационарная вероятность существует и равна [.6, .4] . Это легко проверить, если взять большую мощность матрицы: B = A.copy() for _ in xrange(10): B = numpy.dot(B,B) Здесь B[0] = [0.6, 0.4] . Все идет нормально. Согласно википедии : Стационарный вектор вероятности определяется как […]
Я не знаю, с чего начать с этой проблемы, поскольку у меня не было большого опыта в этом, но требуется решить эту часть проекта с помощью компьютера. У меня есть ODE 2-го порядка: m = 1220 k = 35600 g = 17.5 a = 450000 и b составляет от 1000 до 10000 с шагом 500. […]
У меня есть трехмерное изображение мозга (давайте назовем его вспышкой), и на данный момент это 263 x 256 x 185. Я хочу изменить его размер другого изображения (назовите его whole_brain_bravo); 256 x 256 x 176 и (надеюсь) использовать интерполяцию lanczos для повторной выборки (Image.ANTIALIAS). Моя (неудачная) попытка: from scipy import ndimage as nd import nibabel […]
Я хочу поднять 2-мерный array numpy, назовем его A , мощностью некоторого числа n , но до сих пор мне не удалось найти функцию или оператор для этого. Я знаю, что я мог бы применить его к типу matrix и использовать тот факт, что тогда (подобно тому, что было бы поведением в Matlab), A**n делает […]
У меня две разные машины с scipy 0.12 и PIL. На одной машине, когда я пытаюсь прочитать файл .png, он возвращает массив целых чисел с размером (wxhx 3): In[2]: from scipy.ndimage.io import imread In[3]: out = imread(png_file) In[4]: out.shape Out[4]: (750, 1000, 4) На другой машине, используя тот же файл изображения, возвращается объект PIL.PngImagePlugin.PngImageFile завернутый […]
У меня есть коллекция O (N) NxN scipy.sparse.csr_matrix , и каждая разреженная матрица имеет порядок N элементов. Я хочу добавить все эти матрицы вместе, чтобы получить обычный массив NxN numpy. (N порядка 1000). Расположение ненулевых элементов внутри матриц таково, что полученная сумма, конечно, не является разреженной (фактически фактически нет нулевых элементов). На данный момент я […]