Вопросы по теме 'mcmc'
Пример кода ядра со стороны cuRand Mersenne twister __device__
Я работаю над реализацией алгоритма Монте-Карло цепи Маркова (MCMC) на NVIDIA CUDA GPU. Алгоритм CPU MCMC использует высококачественный генератор случайных чисел Mersenne Twister, и я хотел бы использовать его в написанных мною ядрах графического...
3525 просмотров
schedule
05.07.2022
Изучение дискретных параметров HMM в PyMC
Я пытаюсь узнать параметры простого дискретного HMM, используя PyMC. Я моделирую дождливо-солнечную модель со страницы Wiki на HMM. Модель выглядит следующим образом:
Я использую следующие приоры.
theta_start_state ~ beta(20,10)...
1450 просмотров
schedule
20.05.2024
R : функция для создания распределения смеси
Мне нужно создать образцы из смешанного дистрибутива
40% выборок исходят из гауссова (среднее значение = 2, стандартное отклонение = 8)
20% образцов поступают из Коши (местоположение = 25, масштаб = 2)
40% выборок исходят от гауссова...
9155 просмотров
schedule
07.07.2022
Стохастические переменные в pymc
Я встречал такие термины, как
runiform, rbinomial etc
во многих местах.
Я нигде не мог найти о них. Я вижу только их использование.
Что они обозначают и чем отличаются от uniform, binomial соответственно
209 просмотров
schedule
29.05.2022
PyMC: установка ограничений при подгонке моделей
Я пытаюсь установить ограничения при подгонке переменных с помощью подхода MCMC с PyMC. Например, я определил следующие стохастические модели в PyMC.
import pymc as pm
a=pm.Uniform('a',lower=0.,upper=1.,value=0.2)...
405 просмотров
schedule
19.03.2023
Составная модель в PyMC
Я пытаюсь использовать PyMC 2.3 для получения оценки параметра составной модели . Под «составным» я подразумеваю случайную величину, которая следует за распределением, параметром которого является другая случайная величина. («вложенные» или...
775 просмотров
schedule
25.10.2022
Как определить общую детерминированную функцию в PyMC
В моей модели мне нужно получить значение моей детерминированной переменной из набора родительских переменных, используя сложную функцию Python.
Возможно ли это сделать?
Ниже приведен код pyMC3, который показывает, что я пытаюсь сделать в...
1350 просмотров
schedule
13.04.2023
Как рассчитать смоделированные значения при построении графика несоответствия для соответствия?
Я пытаюсь построить график несоответствия для проверки соответствия после получения значений наилучшего соответствия с помощью MCMC с использованием pymc. Мой код выглядит так:
import pymc
import numpy as np
import matplotlib.pyplot as plt,...
153 просмотров
schedule
23.12.2023
Metropolis Hastings для модели линейной регрессии
Я пытаюсь реализовать алгоритм Метрополиса-Гастингса для простой линейной регрессии на C (без использования других библиотек (boost, Eigen и т.д.) и без двумерных массивов)*. Для лучшего тестирования кода/оценки графиков трассировки я переписал код...
1289 просмотров
schedule
27.05.2024
использовать ведущего для выборки функций распределения вероятностей?
У меня есть два массива (pdf_#), каждый из которых содержит вероятность параметра (val_#). Как я могу выбрать из этих дистрибутивов для создания совместного апостериорного распределения (например, как в угловой участок здесь )? Для ведущего похоже,...
367 просмотров
schedule
26.05.2024
LaplacesDemon: когда я должен суммировать априорную плотность?
Я перехожу с JAGS на LaplacesDemon и пытаюсь переписать некоторые из своих кодов. Я прочитал Учебное пособие по LaplacesDemon и Примеры LaplacesDemon виньетки, и я немного смущен некоторыми примерами в виньетки.
В простом примере из...
73 просмотров
schedule
03.12.2022
Как интерпретировать апостериорную вероятность регрессионной переменной?
Я использую следующую модель в WINBUGS для запуска иерархической байесовской регрессии, где beta являются моими ковариатами:
Если я изменю эту модель, добавив следующий код:
# posterior probabilities of Positive beta's
p.beta0 <-...
370 просмотров
schedule
03.12.2023
Как исправить «Узел несовместим с родителями» в R2jags::jags
Я работаю с R-пакетом R2jags. После запуска кода, который я прикрепляю ниже, R выдал сообщение об ошибке: «Узел несовместим с родителями».
Я пытался решить это. Однако сообщение об ошибке сохраняется. Я использую следующие переменные:
i)...
205 просмотров
schedule
12.02.2023
Как пакетировать преобразованное (масштабированное и квантованное) бета-распределение в вероятность тензорного потока
Я пытаюсь подогнать бета-распределение к результатам опроса с дискретными оценками (1, 2, 3, 4, 5) . Для этого мне нужен рабочий log_prob вероятности бета-версии TensorFlow. Однако существует проблема с обработкой пакетной обработки в бета-версии....
140 просмотров
schedule
02.04.2023
Сообщение об ошибке в JAGS при определении ограничений для априорных точек Попытка переопределить узел
Я пытаюсь взять образцы из предыдущих дистрибутивов и ограничить некоторые суммы и различия, чтобы они были больше 0, однако, когда я запускаю этот код, который обычно работает в OpenBugs / WinBugs, я получаю следующую ошибку «Попытка переопределить...
40 просмотров
schedule
25.02.2022
Guided-LDA без свернутой выборки Гиббса
Я использую пакет Python Guided-LDA со свернутой выборкой Гиббса [ подробнее ]. Но мне нужен управляемый LDA без свернутой выборки Гиббса.
Как я могу отделить управляемую LDA-выборку и свернутую выборку Гиббса или как получить результаты без...
111 просмотров
schedule
28.01.2024
MCMCglmm: извлечение задней поверхности смоделированных фиксированных эффектов
Мне нужно извлечь апостериорные оценки SE для каждого фиксированного эффекта из моей модели.
В иллюстративных целях набор данных, аналогичный тому, который я использую, будет набором данных ChickWeight в базе R .
Я извлекаю апостериорные оценки...
24 просмотров
schedule
03.05.2022
Прогнозировать значения из сложной модели Rjags
Я впервые работаю с R2Jag, цепочками MCM и байесовскими моделями, и у меня возникают проблемы с вычислением прогнозируемых значений для моей модели. Модель основана на исследовании Hallmann et al. 2017 г. (...
37 просмотров
schedule
27.06.2022
Статистика Гельмана-Рубина для модели MCMCglmm в R
У меня есть многомерная модель с этой (приблизительной) формой:
library(MCMCglmm)
mod.1 <- MCMCglmm(
cbind(OFT1, MIS1, PC1, PC2) ~
trait-1 +
trait:sex +
trait:date,
random = ~us(trait):squirrel_id + us(trait):year,
rcov...
73 просмотров
schedule
25.06.2022