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

Регрессия PyMC многих регрессий?
Я не так давно использую PyMC, но был доволен тем, как быстро мне удалось получить линейную регрессию с нуля (этот код должен работать без изменений в IPython): import pandas as pd from numpy import * import pymc data=pd.DataFrame(rand(40))...
922 просмотров
schedule 07.05.2022

Решение обратных задач с помощью PyMC
Предположим, у нас есть априорное значение для X (например, X ~ Gaussian) и прямой оператор y = f(x) . Предположим далее, что мы наблюдали y с помощью эксперимента, и этот эксперимент можно повторять бесконечно. Предполагается, что выходной...
1728 просмотров
schedule 18.08.2022

Как смоделировать смесь 3 нормалей в PyMC?
Существует вопрос о перекрестной проверке о том, как использовать PyMC для подгонки двух нормальных распределений к данным. Ответ Cam.Davidson.Pilon состоял в том, чтобы использовать распределение Бернулли для присвоения данных одному из две...
2198 просмотров
schedule 03.09.2022

pymc 3.0 Прогнозирование апостериорного распределения
Я конвертирую очень простой пример из pymc 2.3 в pymc 3.0 и не могу понять, как выбрать (или получить MAP) из прогнозирующего апостериорного распределения. Следуя предложению в документах (7.3 Проверка модели и диагностика: качество подгонки), я могу...
1015 просмотров
schedule 01.06.2023

Регрессия с «однонаправленным» шумом
Я хотел бы оценить параметры простой линейной функции и член шума с гамма-распределением по данным. (Примечание. Это дополнительный вопрос https://stats.stackexchange.com/questions/88676/regression-with-unidirectional-noise , но упрощенный и более...
126 просмотров
schedule 03.06.2022

Сопоставьте нелинейную функцию с данными/наблюдениями с помощью pyMCMC/pyMC
Я пытаюсь сопоставить некоторые данные с функцией Гаусса (и более сложной). Я создал небольшой пример ниже. Мой первый вопрос: правильно ли я делаю? Мой второй вопрос: как добавить ошибку в направлении x, т. е. в положении x...
8907 просмотров

PyMC, детерминированные узлы в циклах
Я немного новичок в Python и PyMC и быстро прогрессирую. Но меня просто смущает использование установки детерминированных значений 2D-матрицы. У меня есть модель ниже, которую я не могу правильно разобрать. Проблема связана с установкой значения...
902 просмотров
schedule 26.10.2022

Процесс Дирихле в PyMC 3
Я хотел бы реализовать пример процесса Дирихле, указанный в Реализация процессов Дирихле для байесовских полупараметрических моделей (источник: здесь ) в PyMC 3. В примере вероятности поломки палки вычисляются с помощью декоратора...
897 просмотров
schedule 12.08.2022

Свернутая выборка/отдельные шаги Метрополис-Гастингс
Моя модель имеет три параметра, скажем, theta_1 , theta_2 и nu . Я хочу отобрать theta_1 , theta_2 из апостериорной с маргинализацией nu (что можно сделать аналитически), т. е. из p( theta_1, theta_2 | D) вместо p(theta_1, theta_2, nu...
449 просмотров
schedule 28.03.2023

Ошибка при запуске PyMC — значение стохастика выходит за рамки его поддержки или запрещает текущие значения его родителей.
Я пытаюсь использовать PyMC, чтобы взять в качестве входных данных четыре различных предиктора состояния здоровья и объединить их, чтобы получить общую апостериорную вероятность того, что у пациента есть состояние, учитывая подмножество предикторов,...
1281 просмотров
schedule 24.03.2024

Количество параметров в MCMC
Я хочу взять образец из моего апостериорного распределения, используя пакет pymc . Мне интересно, есть ли ограничение на количество измерений, которые может обрабатывать такой алгоритм. Моя логарифмическая вероятность - это сумма 3 гауссианов и 1...
229 просмотров
schedule 20.06.2022

качество подгонки в pymc и построение расхождений
Я использую PYMC 2.3.4. Я нашел потрясающий. Теперь я хотел бы немного улучшить подгонку и расхождения графиков, как показано в разделе 7.3 документации ( https://pymc-devs.github.io/pymc/modelchecking.html ) В документации сказано, что для графика...
456 просмотров
schedule 07.12.2022

Как получить скорость принятия из базы данных pymc pickle
Я использую приведенный ниже метод, чтобы получить уровень принятия после запуска MCMC с pymc (вдохновленный этот ответ ): MDL.step_method_dict[next(iter(MDL.stochastics))][0].ratio (или есть более простой способ?) Но если я сохраню MDL в...
236 просмотров
schedule 23.07.2022

Понимание кода с декораторами в pymc, Python
Я пытался получить некоторую информацию об иерархических моделях и нашел этот замечательный пост: https://sl8r000.github.io/ab_testing_statistics/use_a_hierarchical_model/ В нем, закройте середину поста, автор делится каким-то кодом. Есть часть,...
121 просмотров

PyMC DiscreteMetropolis в дискретной плавающей сетке
В настоящее время я пытаюсь решить проблему из астрофизики, которую можно упростить следующим образом: Я хотел подогнать линейную модель (скажем, y = a + b*x ) к наблюдаемым данным, и я хочу использовать PyMC для характеристики апостериорных...
147 просмотров
schedule 27.02.2023

Ошибка при построении модели PyMC с использованием pymc.Matplot
Я получаю сообщение об ошибке при построении модели PyMC с использованием pymc.Matplot после успешной подгонки #!/usr/bin/env python2 import pylab import pymc as pm import numpy as np import seir as mod reload(mod) mc = pm.MCMC(mod)...
162 просмотров
schedule 29.04.2023

Прогнозирование апостериорности новых данных в байесовской линейной регрессии с использованием PyMC3
Использование PyMC3 для выполнения байесовской линейной регрессии. Я построил свою модель и хочу предсказать апостериорное значение для новых значений X, используя ту же модель. Я пытаюсь следовать инструкциям на веб-сайте документации:...
1618 просмотров

Удобный способ преобразований в pymc3
В некоторых случаях pymc3 автоматически выполняет преобразования переменных (см. Каковы интервальные преобразования в pymc3 для равномерного распределения? , что означает вывод find_MAP в pymc3? ). В ответе...
88 просмотров
schedule 30.12.2022

Правильная спецификация вложенной модели логита в pymc3
Я работаю над моделью вложенной логистической регрессии с 3 результатами, представляющими варианты A, B или C. Первый уровень представляет выбор между A и B или C, а второй уровень представляет выбор между B и C. некоторые выдуманные данные приведены...
247 просмотров
schedule 14.02.2023

PyMC3 - различия в способах передачи наблюдений в модель - ›разница в результатах?
Я пытаюсь понять, есть ли какая-либо значимая разница в способах передачи данных в модель - агрегированных или отдельных испытаний (обратите внимание, что это будет разумный вопрос только для определенных распределений, например, биномиальных)....
262 просмотров
schedule 09.11.2022