Вопросы по теме 'mcmc'
Алгоритм Метрополиса-Гастингса с использованием распределения предложений, отличного от гауссового, в Matlab
В настоящее время я работаю над своим выпускным проектом для получения степени по математике, который основан на обзоре алгоритма Метрополиса-Гастингса и некоторых числовых примерах. До сих пор я получил отличные результаты, используя мое...
1474 просмотров
schedule
24.10.2023
Ошибка в векторе (список, n.chains): недопустимый аргумент «длина»
Я использую R2jags и CODA для получения некоторой диагностической статистики для своих цепей MCMC, но у меня возникают проблемы. Я хочу запустить MCMC следующим образом:
modelfit <- jags(data=jags.data, inits=jags.inits, model.params, n.iter =...
5320 просмотров
schedule
14.06.2022
Вызов R из C из R в mcmc
Поэтому мне приходится запускать большие mcmc, которые невозможны в R. Поэтому я пишу циклы в C и вызываю их из R, используя .C.
Однако внутри этого цикла у меня иногда есть этапы выборки или этапы моделирования, для которых в R есть стандартные...
283 просмотров
schedule
03.03.2023
Обработка нескольких изображений с помощью PyMC3
У меня возникла проблема с обработкой изображений, которую я решил использовать, чтобы поэкспериментировать и узнать больше о PyMC3. Я провел немало времени, возясь с нелинейными решателями и методами грубой силы, и пока ничто меня не радует.
Моя...
894 просмотров
schedule
16.01.2023
Трудности на pymc3 по сравнению с pymc2, когда задействованы дискретные переменные
Я обновляю некоторые расчеты, в которых я использовал pymc2 для pymc3, и у меня возникают некоторые проблемы с поведением сэмплеров, когда у меня есть некоторые дискретные случайные переменные в моей модели. В качестве примера рассмотрим следующую...
1429 просмотров
schedule
08.04.2023
настройка MCMC с логарифмической вероятностью и логарифмической нормой до использования PyMC
Я новичок в pyMC, и я все еще не могу построить структуру своего MCMC с помощью pyMC. Я хотел бы создать цепочку, и я не понимаю, как вместе определить мои параметры и функцию логарифма правдоподобия. Моя функция хи-квадрат определяется следующим...
1039 просмотров
schedule
15.03.2024
Metropolis Hastings с вероятностью Custom Log в Pymc
Я хочу использовать pymc, чтобы использовать цепочку MH для выборки из пользовательской вероятности журнала. Но я не могу заставить его работать и не могу найти достойный пример в Интернете.
def getPymcVariable(data):
def logp(value):
......
618 просмотров
schedule
27.09.2022
MAP в PyMC выдает непредвиденную ошибку
Я не понимаю, почему MAP выдает ошибку, тогда как MCMC работает нормально в том же сценарии? Я пишу ниже соответствующую часть кода.
tau = Uniform('tau', lower=0.01, upper=5, doc='tau')
rv = [ Multinomial("rv"+str(i), count[i],...
132 просмотров
schedule
31.05.2022
Графики гистограмм в pymc, что означают разные аспекты?
Я определил стохастическую случайную величину (и многие другие, но ради этого вопроса достаточно одной)
tau = pm.DiscreteUniform("tau", lower = 0, upper = 74)
После выборки с использованием MCMC, когда я строю трассировку tau , я получаю...
244 просмотров
schedule
05.10.2022
Что подразумевается под настройкой пошаговых методов в pymc
Я пытаюсь поэкспериментировать с разными значениями аргументов MCMC.sample в pymc.
Я просмотрел справочные страницы для MCMC.sample и обнаружил:
tune_interval: int
Шаговые методы будут настраиваться с интервалом в такое количество итераций,...
186 просмотров
schedule
17.06.2023
Как выполнить анализ панельных данных в байесовской модели с помощью pymc
каждый. У меня есть вопрос о том, как проводить анализ панельных данных в байесовской модели с помощью pymc. Данные такие:
..........................................................
User Time x1 x2 x3 Y
1...
1106 просмотров
schedule
08.05.2022
Стандартная версия модели JAGS, которая включает сумму дискретных значений - возможно ли это?
Я пытался запустить эту модель в Стэне. У меня есть работающая версия JAGS (которая возвращает сильно автокоррелированные параметры), и я знаю, как сформулировать ее как CDF двойной экспоненты (с двумя скоростями), которая, вероятно, будет работать...
950 просмотров
schedule
22.04.2023
NUTS pymc3 не может хорошо работать с моей иерархической моделью для байесовских нейронных сетей?
У меня есть иерархическая модель для изучения байесовских сетей с единственным скрытым слоем. Сетевые параметры разделены на 4 группы весов и смещений «вход-скрытый» и «скрытый-выход». Гауссовский априор определяется для каждой группы параметров....
841 просмотров
schedule
20.03.2023
Вызов подпрограмм Stan из программы на C++
Я прочитал здесь , что можно (и я интерпретировал прямо) вызвать Stan подпрограммы из программы на C++.
У меня есть несколько сложных функций логарифмического правдоподобия, которые я закодировал на C++, и я действительно понятия не имею, как я...
1327 просмотров
schedule
31.08.2022
Почему эта иерархическая модель Пуассона не соответствует истинным параметрам сгенерированных данных?
Я пытаюсь подогнать иерархическую регрессию Пуассона для оценки time_delay для каждой группы и в глобальном масштабе. Я не понимаю, применяет ли pymc автоматически функцию ссылки на журнал к mu или мне нужно сделать это явно:
with pm.Model() as...
556 просмотров
schedule
25.05.2022
Управление объектом mcmc.list в R
Я использовал JAGS, вызываемый через rjags, для создания объекта mcmc.list foldD_samples, который содержит мониторы трассировки для большого количества стохастических узлов (>800 узлов).
Теперь я хотел бы использовать R для вычисления довольно...
3487 просмотров
schedule
15.08.2023
MCMCglmm для непрерывных данных
Я имею в виду пакет R MCMCglmm (Обобщенные линейные модели смешанных эффектов с цепью Маркова Монте-Карло), см. Cran.r-project.org/web/packages/MCMCglmm/MCMCglmm.pdf
Хотя MCMCglmm задается как обобщенная модель смешанного эффекта (и поэтому не...
138 просмотров
schedule
26.11.2022
простой Gamma GLM в STAN
Пробую простой Gamma GLM в STAN и R, но он сразу вылетает
генерировать данные:
set.seed(1)
library(rstan)
N<-500 #sample size
dat<-data.frame(x1=runif(N,-1,1),x2=runif(N,-1,1))
#the model
X<-model.matrix(~.,dat)
K<-dim(X)[2] #number...
1167 просмотров
schedule
27.07.2023
Две проблемы при написании скрипта для вычисления марковского совместного распределения (на питоне)
Я новичок в Python, недавно я работаю над каким-то проектом, чтобы выполнить вычисление совместного распределения марковского процесса.
Пример стохастического ядра используется в недавнем исследовании Гамильтона (2005), который исследует...
282 просмотров
schedule
10.09.2023
Гауссова модель смешивания PyMC3
Я следил за примером смешанной модели Гаусса для PyMC3 здесь: https://github.com/pymc-devs/pymc3/blob/master/pymc3/examples/gaussian_mixture_model.ipynb и прекрасно работает с искусственным набором данных.
Я пробовал это с реальным набором...
2527 просмотров
schedule
08.02.2023