Вопросы по теме 'pymc3'
Как сделать усеченный нормальный априор: преобразование pymc2 в pymc3
Как в pymc3 настроить усеченный нормальный априор? В pymc2 это довольно просто (см. Ниже), но в pymc3, похоже, больше нет доступного усеченного нормального распределения.
Pymc2:
TruncatedNormal('gamma_own_%i_' % i, mu=go, tau=v_gamma_inv,...
1498 просмотров
schedule
09.11.2021
Как использовать распределение StudentT в pymc3?
Я не уверен, считается ли это вопросом или сообщением об ошибке. Я разместил суть GitHub здесь: https://gist.github.com/jbwhit/a9012e04b0f48e582c22
Я нашел этот вопрос ( pymc3: иерархическая модель с несколькими наблюдаемыми переменными ) быть...
1061 просмотров
schedule
10.11.2021
Многомерная нормальная модель с отдельными средними, ошибка несоответствия размеров
Я хочу смоделировать многомерную (в примере двумерную) нормаль с отдельными средними и ковариационной матрицей, почему-то я просто не могу добиться соответствия размеров. Модель относительно проста, вот процесс генерации данных:
import numpy as...
633 просмотров
schedule
24.09.2021
Как использовать массив распределений в качестве аргумента для нового распределения в pymc3
Я пытаюсь смоделировать следующую проблему:
Перед некоторым моментом переключения (выбранным из DiscreteUniform distribution) значения выбираются из Uniform distribution с границами early_low , early_high . После этого переключения во времени...
163 просмотров
schedule
14.09.2021
Отследить отдельные размеры с помощью графика трассировки PYMC3?
Я использовал PYMC3 для выполнения вывода по модели байесовской логистической регрессии. Я хочу найти апостериор по весам $ \ beta \ in \ mathbb {R} ^ K $ с учетом гауссовского априорного значения $ \ mathcal {N} \ sim (0, I) $.
При этом...
1128 просмотров
schedule
22.10.2021
PyMC3 намного медленнее, чем PyMC в выборке Metropolis
Я пытаюсь сравнить скорость выборки между PyMC и PyMC3.
PyMC:
p1 = pymc.Normal('p1', 10, 0.5)
p2 = pymc.Gamma('p2', 11, 5)
p3 = pymc.Normal('p3', p1, p2)
model = pymc.Model([p1, p2, p3])
mcmc = pymc.MCMC(model)...
965 просмотров
schedule
14.09.2021
Как согласовать TFP с результатами PyMC3 MCMC?
При тестировании TFP я попытался взять образец из апостериорного распределения модели сопряженных нормальных норм (известная дисперсия), то есть
x | mu ~ нормальный (mu, 1.)
mu ~ Нормальный (4., 2.)
Пробоотборник tf.mcmc.RandomWalkMetropolis...
226 просмотров
schedule
10.09.2021
Прогнозирование базовой байесовской линейной регрессии с помощью PyMC3
Я хотел бы использовать мою модель PyMC3 LR, чтобы получить диапазон HPD 80% для значения прогнозируемой переменной y по мере появления новых данных. Таким образом, экстраполируйте достоверное распределение значений для y на новое значение x ,...
1356 просмотров
schedule
06.03.2022
Создание прогнозов на основе предполагаемых параметров в pymc3
Я сталкиваюсь с распространенной проблемой, мне интересно, может ли кто-нибудь помочь. Я часто хотел бы использовать pymc3 в двух режимах: обучение (т.е. фактически выполнение вывода по параметрам) и оценка (т.е. использование предполагаемых...
7853 просмотров
schedule
01.04.2022
Биномиальное распределение с параметром loc в pymc3
Я хотел бы использовать биномиальное распределение со сдвигом на параметр loc (как в scipy ) в модели pymc3. Например.:
with pm.Model() as m1:
prob = pm.Beta('prob',alpha=2,beta=2)
x = pm.Binomial('x',n=20,p=prob,loc=5)
Но...
149 просмотров
schedule
13.05.2022
PyMC3 и Theano — код Theano, который работает, перестает работать после импорта pymc3
Некоторый простой теано-код, который отлично работает, перестает работать, когда я импортирую pymc3
Вот несколько фрагментов, чтобы воспроизвести ошибку:
#Initial Theano Code (this works)
import theano.tensor as tsr
x = tsr.dscalar('x')
y...
1743 просмотров
schedule
01.07.2022
Старая трассировка группировки в стиле PyMC3, построенная с помощью Arviz
У меня есть старый пост в блоге, где я тренирую модель PyMC3. Вы можете найти сообщение в блоге здесь , но суть модели показана ниже. .
with pm.Model() as model:
mu_intercept = pm.Normal('mu_intercept', mu=40, sd=5)
mu_slope =...
100 просмотров
schedule
05.07.2022
Почему я получаю несоответствие размеров в моей иерархической модели PyMC3?
По сути, это пример «Несколько монет от нескольких монетных дворов / бейсболистов» из Doing Bayesian Data Analysis, Second Edition (DBDA2). Я считаю, что у меня есть код PyMC3, который функционально эквивалентен, но один работает, а другой нет....
454 просмотров
schedule
08.07.2022
Построение многомерного нормального распределения с вероятностными параметрами в PyMC3
Я хочу построить многомерную нормальную модель в PyMC3, в которой матрица среднего значения и точности включает вероятностные переменные. h предназначен для использования в качестве скрытой переменной в более крупном проекте, к которому принадлежит...
448 просмотров
schedule
20.07.2022
Инкрементальные байесовские обновления с многомерными параметрами
Я пытаюсь использовать PYMC3 для байесовской модели, где я хотел бы неоднократно обучать свою модель новым невидимым данным. Я думаю, мне нужно будет обновлять априорные значения апостериорной частью ранее обученной модели каждый раз, когда я вижу...
79 просмотров
schedule
03.08.2022
Простая байесовская сеть с проблемой подбрасывания монеты
Я пытаюсь реализовать байесовскую сеть и решить проблему регрессии с помощью PYMC3. В моей модели у меня есть честная монета в качестве родительского узла. Если родительским узлом является H, дочерний узел выбирает нормальное распределение N(5,0.2);...
248 просмотров
schedule
27.08.2022
Распараллеливание PyMC
Может ли кто-нибудь дать общие инструкции о том, как можно распараллелить код PyMC MCMC . Я пытаюсь запустить LASSO регрессию, следуя приведенному примеру здесь . Я где-то читал, что параллельная выборка выполняется по умолчанию, но мне все еще...
4574 просмотров
schedule
08.10.2022
Определение стохастических и детерминированных переменных с помощью pymc3
Я пытаюсь использовать собственные переменные stochastic и deterministic с pymc3 , но старый опубликованный рецепт для pymc2.3 объяснял, как мы можем параметризовать наши переменные, больше не работает. Например, я пытался использовать этот...
2466 просмотров
schedule
08.04.2023
доступ к значению случайной величины в PyMC3
В PyMC2 есть методы random() и value() для генерации случайного значения и получения текущего значения случайных переменных. Есть ли способ сделать то же самое в PyMC3?
p = pm.Dirichlet('p', theta=np.array([1., 1., 1.]))
p.random()
p.value
945 просмотров
schedule
12.01.2023
Можно ли использовать операторы if в детерминированной функции PyMC?
После прочтения книги Кэма Дэвидсона-Пилона Вероятностное программирование и байесовские методы для хакеров , я решил попробовать свои силы в решении проблемы обучения со скрытой марковской моделью (HMM) с помощью PyMC. Пока что код не...
933 просмотров
schedule
20.07.2023