Вопросы по теме 'pymc3'

Перенос кода PyMC2 в PyMC3 — иерархическая модель для спортивной аналитики
Я попробовал следующий код, но столкнулся с проблемами. Я думаю, что проблема в .values, но как мне закодировать это как объект Theano? Ниже приведен мой источник данных home_team,away_team,home_score,away_score Wales,Italy,23,15...
1001 просмотров
schedule 29.03.2024

сохранить подграфики трассировки PyMC3 в файл изображения
Я пытаюсь очень просто построить подграфики, сгенерированные функцией трассировки PyMC3 (см. здесь ) в файл. Функция генерирует numpy.ndarray (2d) подграфиков. Мне нужно переместить или скопировать эти подграфики в matplotlib.figure, чтобы...
1564 просмотров
schedule 12.11.2022

Какой степпер Монте-Карло pymc3 я могу использовать для пользовательского категориального распределения?
Я работаю над реализацией скрытых цепей Маркова в pymc3. Я продвинулся довольно далеко в реализации скрытых состояний. Ниже я показываю простую марковскую цепь с двумя состояниями: import numpy as np import pymc3 as pm import theano.tensor as...
377 просмотров

PyMC3 - Пуассон работает с точкой переключения, экспоненциальный - нет
Следующая модель (взятая из Байесовских методов для хакеров) работает с Пуассоном. count_data = np.loadtxt("data/txtdata.csv") n_count_data = len(count_data) with pm.Model() as model: alpha = 1.0/count_data.mean() # Recall count_data is the...
939 просмотров
schedule 17.03.2024

Сезонные переменные PYMC3
Я относительно новичок в PYMC3 и пытаюсь реализовать временной ряд байесовской структуры (BSTS) без регрессоров, например, модель подходит здесь в R. Модель выглядит следующим образом: Я могу реализовать локальный линейный тренд с помощью...
969 просмотров
schedule 06.04.2024

Получение апостериорного распределения разницы между двумя переменными с использованием PYMC3
Теперь предположим, что мы смотрим на дневные цены двух акций, A и B. Предыстория проста: все цены имеют нормальное распределение, причем mu_A и mu_B равномерно распределены на [10 100], а sigma_A и sigma_B также равномерно распределены на [1,100]....
624 просмотров
schedule 27.12.2023

как рассчитать журнал после GP по трассе в pymc3
Вариант использования Предположим, у меня есть наблюдение y_0 в X_0 , которое я хотел бы смоделировать с помощью гауссовского процесса с гиперпараметрами theta . Предположим, что затем я определяю распределение в гиперпараметрах theta путем...
222 просмотров
schedule 19.01.2024

Распределение Вишарта в pymc3
Я делаю байесовский анализ с pymc3 (графическая модель с большим количеством априорных значений), как показано ниже: import pymc3 as pm import theano.tensor as T with pm.Model() as smk: a_u =...
138 просмотров
schedule 03.12.2022

PyMC3 Normal с дисперсией на столбец
Я пытаюсь определить переменную pymc3.Normal со следующим значением: mu: import numpy as np import pymc3 as pm mx = np.array([[0.25 , 0.5 , 0.75 , 1. ], [0.25 , 0.333, 0.25 , 0. ], [0.25 , 0.167, 0. , 0....
235 просмотров
schedule 05.11.2022

Как получить конкретную вероятность определенного значения в распределении?
Я новичок в PyMC3, я начал три недели назад, чтобы ознакомиться с ним, и в настоящее время я занимаюсь обучением, поэтому у меня есть одно сомнение в моей программе. Извините, если это глупый вопрос. У меня есть эта модель: with pm.Model() as...
66 просмотров
schedule 19.02.2024

Суммируйте апостериорную часть одного параметра из массива с помощью arviz
Я оцениваю модель, используя библиотеку pyMC3 в python. В моей реальной модели есть четыре массива параметров, два из которых содержат более 170 000 параметров. Обобщение этого массива параметров требует слишком больших вычислительных ресурсов на...
21 просмотров
schedule 12.10.2022