Вопросы по теме 'pymc'
Как pymc представляет функцию априорного распределения и правдоподобия?
Если pymc реализует алгоритм Метрополиса-Гастингса для получения выборок из апостериорной плотности по интересующим параметрам, то для того, чтобы решить, переходить ли к следующему состоянию в цепи Маркова, он должен быть в состоянии оценить что-то,...
2264 просмотров
schedule
18.10.2022
Как реализовать в PyMC функцию случайной величины, которую может взять MCMC Metropolis?
Если у вас есть случайная величина $X$ и функция $f$, вы можете определить $y=f(X)$ как новую случайную величину с функцией плотности вероятности следующим образом:
$p(y)=(f^{-1})'(y)p(x)$. Подробнее см. здесь .
Теперь я определил случайную...
1152 просмотров
schedule
31.03.2023
Как построить модель, требующую внешнего пакета в PyMC3?
Я не уверен, что это вопрос PyMC3 или вопрос Theano. Я долгое время использовал PyMC2, чтобы подогнать космологию к данным о сверхновых. Для этого требуются некоторые запутанные интегралы (см., например, http://arxiv.org/abs/astroph/9905116 )...
449 просмотров
schedule
18.06.2023
Проблема с библиотекой pymc
Я пытаюсь запустить следующий код:
import pymc as pm
alpha = 1.0/count_data.mean() #count_data is the variable that holds txtc
lambda_1 = pm.Exponential("lambda_1", alpha)
lambda_2 = pm.Exponential("lambda_2", alpha)
tau =...
627 просмотров
schedule
02.02.2023
Общие модели смесей в pymc
У меня есть дистрибутив с несколькими горбами. Я хотел бы попробовать подобрать несколько различных типов распределений для каждого горба, гауссовского, экспоненциального, распределения Вейбуйля и т. д. Однако в нынешнем виде кажется, что мне нужно...
146 просмотров
schedule
12.02.2023
Проблемы с запуском линейной регрессии Pymc (масштабирование входных параметров?)
Следуя этому примеру , выполните довольно простую байесовскую линейную регрессию с использованием PYMC3 ( я надеюсь, что учусь) я получаю начальный пример для запуска, но затем пытаюсь использовать свои собственные данные и получаю:
ValueError:...
464 просмотров
schedule
22.08.2022
PyMC: выборка шаг за шагом?
Я хотел бы знать, почему сэмплер невероятно медленный при пошаговой выборке. Например, если я запускаю:
mcmc = MCMC(model)
mcmc.sample(1000)
выборка быстрая. Однако, если я запускаю:
mcmc = MCMC(model)
for i in arange(1000):...
812 просмотров
schedule
24.09.2022
Как я могу создать модель Маркова с дискретным состоянием с помощью pymc?
Я пытаюсь понять, как правильно создать модель цепи Маркова с дискретным состоянием с помощью pymc .
В качестве примера (см. в nbviewer ) давайте сделаем цепь длиной T = 10, в которой марковское состояние бинарное, начальное распределение...
2840 просмотров
schedule
05.07.2023
Импортировать модуль с параметром в Python
Можно ли импортировать модуль с каким-либо параметром в Python?
Под параметром я подразумеваю то, что в модуле существует переменная, которая не инициализирована в этом модуле, но я использую эту переменную в этом модуле. Короче говоря, я хочу,...
27588 просмотров
schedule
05.10.2022
Предупреждение pymc MAP: значение Stochastic tau не является ни числовым, ни массивом с dtype с плавающей запятой. Рекомендовать метод подгонки fmin (по умолчанию)
Я посмотрел на аналогичный вопрос здесь
предупреждение pymc: значение не является ни числовым, ни массив с плавающей запятой dtype
но ответов нет, может ли кто-нибудь сказать мне, должен ли я игнорировать это предупреждение или что делать в...
179 просмотров
schedule
07.02.2024
TypeError: hasattr(): имя атрибута должно быть строкой в pymc
Я просмотрел следующие ссылки, но ни одна из них не дает решения, которое я ищу.
https://github.com/pymc-devs/pymc/issues/125
Ошибка PyMC: hasattr(): имя атрибута должно быть строкой
Мне нужно написать функцию, которая с учетом априорных...
673 просмотров
schedule
08.05.2022
Понимание модели катастрофы в PyMC
I start learning PyMC and strungle to understand the very first tutorial´s example.
disasters_array = \
np.array([ 4, 5, 4, 0, 1, 4, 3, 4, 0, 6, 3, 3, 4, 0, 2, 6,
3, 3, 5, 4, 5, 3, 1, 4, 4, 1, 5, 5, 3, 4, 2, 5,...
158 просмотров
schedule
21.09.2022
Использование pymc3 для соответствия распределению Стьюдента
Не уверен, что я делаю что-то глупое или у pymc3 есть ошибка, но пытаясь подогнать T-распределение к нормальному, я получаю количество степеней свободы (от 0,18 до 0,25, я ожидаю чего-то высокого, 4-5 по крайней мере). Конечно, я получаю ту же...
1130 просмотров
schedule
15.01.2023
байесовская PCA с использованием PyMC
Я пытаюсь реализовать байесовский PCA, используя библиотеку PyMC для python. Но я застрял там, где определяю координаты более низкого измерения...
Модель
x = Wz + e
где x — вектор наблюдения, W — матрица преобразования, а z — вектор...
1072 просмотров
schedule
13.08.2022
Интерпретация апостериорных распределений прогона MCMC с помощью pymc
pymc отлично! Он действительно открыл мой мир для MCMC, так что спасибо за его кодирование.
В настоящее время я использую pymc для оценки некоторых параметров и доверительных интервалов путем подгонки функции к наблюдениям. Для большинства...
742 просмотров
schedule
23.07.2023
AIC и BIC смешанной модели PyMC
Я использую PyMC, чтобы подогнать некоторые данные к прямой линии. В данных есть выбросы, поэтому я адаптировал некоторый код (третий пример по ссылке) написано Джейком Вандерпласом для его учебника. Метод использует векторную переменную qi для...
1422 просмотров
schedule
13.05.2023
Как использовать pymc для параметризации вероятностной графической модели?
Как можно использовать pymc для параметризации вероятностной графической модели?
Предположим, у меня есть PGM с двумя узлами X и Y . Допустим, X->Y — это график.
Причем X принимает два значения {0,1} , а Y также принимает два...
549 просмотров
schedule
25.01.2023
pymc3: как моделировать коррелированные точки пересечения и наклона в многоуровневой линейной регрессии
В примере Pymc3 для многоуровневой линейной регрессии (пример здесь , с набором данных radon из работы Гельмана и др. (2007)), точки пересечения (для разных округов) и наклоны (для квартир с подвалом и без) имеют нормальный априор. Как я могу...
884 просмотров
schedule
12.08.2022
Надежная байесовская корреляция с PyMC3
Это дополнительный вопрос по анализу байесовской корреляции, как описано в этот пример для PyMC2 .
Я успешно перенес ненадежный подход, в котором используется многомерное нормальное распределение, на PyMC3, но я борюсь с надежной версией, в...
515 просмотров
schedule
08.07.2023
что такое итерация в pymc
У меня есть модель в pymc, упрощенная
import pymc as pm
a = pm.Normal('a', 0,5)
b = pm.Normal('b', 0,5)
cnt = 0
@pm.deterministic
def mu(a=a, b=b):
global cnt
cnt += 1
return a + b*x # x and y were loaded from csv
y_hat =...
111 просмотров
schedule
10.04.2023