Вопросы по теме '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 просмотров

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