Вопросы по теме 'rstan'
RStan: определение трехуровневой модели случайных уклонов?
Я работал над трехуровневой моделью RStan, в которой повторяющиеся широкополосные измерения (идентификатор года = yrid) вложены в местные органы власти (идентификатор LA = lay), которые, наконец, вложены в регионы (идентификатор региона = rnid)....
508 просмотров
schedule
18.03.2023
Параметр преобразования стандартной ошибки в int
Я пытаюсь провести анализ точки переключения с помощью STAN. У меня есть вектор данных y , который имеет две разные последовательности гауссовских случайных величин. Цель состоит в том, чтобы найти апостериорное распределение того, когда мог...
1381 просмотров
schedule
02.12.2022
Извлечь BFMI из объекта stanfit в rstan
После установки модели в стан с использованием stan() , как я могу извлечь BFMI для каждой цепочки? Из этого разговора: https://groups.google.com/forum/#!topic/stan-dev/uJhsapVwlk8 , похоже, что BFMI в какой-то момент отображается при использовании...
128 просмотров
schedule
10.04.2022
Ошибка в readRDS(file.rds)
Резюме:
Я установил пакет rstan , теперь я столкнулся с этой ошибкой:
Ошибка в readRDS(file.rds): неизвестный входной формат
Описание:
> traceback()
5: readRDS(file.rds)
4: is(obj <- readRDS(file.rds), "stanmodel")
3:...
5009 просмотров
schedule
14.04.2023
Прогнозирование на основе полного апостериорного распределения с использованием stan_glmer
Могу я попросить о помощи?
Я подобрал биномиальную модель, используя stan_glmer, и выбрал модель, которая, на мой взгляд, лучше всего соответствует данным. Я использовал команду апостериорного предсказания, чтобы сравнить мои наблюдаемые данные с...
634 просмотров
schedule
13.02.2022
Как ставить разные приоры с помощью rstanarm
Допустим, у меня есть модель формы y=a_{i} + b_{i,1}*x_{1} + b_{2}*x_{2} , где i=1,2,...,12 , и я хотел бы оценить эту модель, используя rstanarm .
Можно ли установить разные приоритеты для каждого перехвата a_{i} (так, скажем, первые 4 имеют...
653 просмотров
schedule
30.09.2022
Предупреждения, выдаваемые на этапе компиляции stan
Я получаю два предупреждения на этапе компиляции cpp со всеми программами stan, которые я отправляю.
C:/Larry/R/win-library/3.4/BH/include/boost/config/compiler/gcc.hpp:186:0: предупреждение: "BOOST_NO_CXX11_RVALUE_REFERENCES" переопределено #...
151 просмотров
schedule
12.07.2023
График гауссовского процесса в R с использованием rstan
Я пытаюсь понять, где я ошибаюсь с rstan . Я нашел обходной путь, но кажется, что должен быть лучший вариант для графического отображения розыгрыша сзади, чем то, что я придумал.
Я пытаюсь научиться использовать rstan для моделирования...
298 просмотров
schedule
02.07.2023
Оптимизация гауссовского процесса в Stan / rstan
Недавно я столкнулся с гауссовскими моделями процессов и случайно подумал, что они могут быть решением проблемы, над которой я работал в своей лаборатории. У меня есть открытый и связанный с этим вопрос о перекрестной проверке, но я хотел отделить...
508 просмотров
schedule
23.06.2023
Построение модели переменного пересечения в RStan - сообщение об ошибке
Я пытаюсь построить модель переменного перехвата в RStan , пытаясь предсказать использование женских противозачаточных средств в одном регионе Индии на основе возраста, количества детей и того, живут ли они в городской местности, с перехваты...
548 просмотров
schedule
09.01.2023
Получение стандартизированных коэффициентов из пакета rstanarm в R?
Мне интересно, возможно ли (и, возможно, рекомендуется) получить стандартизированные коэффициенты от stan_glm() в пакете rstanarm ? ( в документации ничего конкретного не нашел )
Могу ли я просто стандартизировать все переменные, как в...
427 просмотров
schedule
16.04.2023
апостериорный прогноз на основе группирующей переменной из `stan_glm()` в пакете `rstanarm`?
Мне было интересно, как получить апостериорный прогноз на основе переменной группировки из stan_glm() в пакете rstanarm ?
Например, если у меня есть бинарная группирующая переменная с кодом (0, 1) , называемая "vs" в моих данных (базовые...
412 просмотров
schedule
29.05.2023
Rstan не работает с многоядерной поддержкой
Я установил R версии 3.5.1 (на настольный ПК с Windows 7 x64), а затем Rtools 3.5 и rstan (через install.packages())
Rstan терпит неудачу с многоядерной поддержкой, например:
library(rstan)
rstan_options(auto_write = TRUE)
options(mc.cores = 3)...
265 просмотров
schedule
04.04.2023
Определение априора как для случайных эффектов, так и для дисперсии случайных эффектов в модели со смешанными эффектами (R brms)
Я хочу соответствовать подсчетам модели Пуассона GLMM. У меня 121 субъект ( subject ), я наблюдаю 8 счетов Пуассона ( count ) на субъект: они соответствуют 2 типам событий ( event ) x 4 периодам ( period ).
'data.frame': 968 obs. of 4...
583 просмотров
schedule
25.05.2022
Обзор модели оценки параметров полиномиальной регрессии Стэна
У меня есть следующая модель полиномиальной регрессии:
Версия изображения
Версия LaTeX
$ Y_i | \ mu_i, \ sigma ^ 2 \ sim \ text {Normal} (\ mu_i, \ sigma ^ 2), i = 1, \ dots, n \ \ text {независимый} $
$ \ mu_i = \ alpha + \ beta_1...
282 просмотров
schedule
06.06.2022
Разработка иерархической версии модели нелинейной кривой роста в Stan
Следующая модель является моделью 1 Приса и Бейнса (1978, Annals of Human Biology ) и используется для описания роста человека.
Мой код Stan для этой модели выглядит следующим образом:
```{stan output.var="test"}
data {...
594 просмотров
schedule
07.05.2023
Неэффективная выборка простой биномиальной GLM с использованием rstan MCMC
Я пытаюсь подогнать биномиальный GLM с помощью пакета rethinking (опирается на rstan MCMC).
Модель подходит, но выборка неэффективна, и Рат указывает на то, что что-то пошло не так. Я не понимаю причину этой проблемы с подгонкой.
Это...
148 просмотров
schedule
28.04.2023
байесовский анализ выживаемости с RST
Я хочу построить код stan (rstan) для анализа выживаемости с использованием распределения Вейбулла. Но мой стандартный код всегда не работает. Если кто-нибудь знает, как справиться с моей проблемой, пожалуйста, научите меня.
Мои данные такие...
517 просмотров
schedule
31.07.2023
Как эффективно распараллелить brms::brm?
Краткое описание проблемы
Я подгоняю модель brms::brm_multiple() к большому набору данных, в котором отсутствующие данные были введены с помощью пакета mice . Размер набора данных делает использование параллельной обработки очень желательным....
1131 просмотров
schedule
20.06.2023
Как указать матрицу предикторов для блока данных stan?
Уважаемое сообщество stackoverflow. Я хочу использовать переменные от w1 до w10 в качестве матрицы предикторов matrix[N, W] weights; в моей модели стана. Я не уверен, как это сделать.
кадр данных
(dat <- data.frame(
id = c(1, 1, 1, 1, 1,...
152 просмотров
schedule
01.07.2023