Описание тега rstanarm
Пакет R, предоставляющий интерфейс для построения и выполнения вывода для байесовских регрессионных моделей.
1
ответ
Подход уникальных перехватов для категориальных переменных в пакете "rstanarm" в R
Справочная информация: McElearth (2016) в своих переосмысленных страницах книги 158-159 использует переменную индекса вместо фиктивного кодирования для переменной 3 категории, называемой "clade", для предсказания "kcal.per.g" (линейная регрессия). В…
21 мар '18 в 20:56
0
ответов
Как обновить объект brmsfit с измененной моделью stan, сгенерированной brms, маргинализуя распределение весов
Я был бы признателен за любую помощь в обновлении моего объекта brmsfit с помощью измененной модели stan, сгенерированной brms, потому что я хочу передать различные столбцы весов вероятности способом, который brms пока не поддерживает, вероятно. Моя…
02 фев '19 в 08:47
1
ответ
Оптимизация гауссовского процесса в Станции
Недавно я столкнулся с гауссовыми моделями процессов и случайно подумал, что они могут стать решением проблемы, над которой я работал в своей лаборатории. У меня есть открытый и связанный вопрос о перекрестной проверке, но я хотел отделить мои вопро…
02 фев '18 в 02:24
1
ответ
Как нам найти разные значения для каждой итерации в rstanarm?
Я написал основной код GLM в rstanarm, Когда я выполняю следующий код, он дает среднее значение переменных-предикторов всех выполненных итераций: summary(TestGLM) print(TestGLM) Я хотел получить значение переменных-предикторов на каждой итерации для…
03 окт '17 в 13:49
3
ответа
Построение апостериорных оценок параметров от нескольких моделей с байесплотом
Я использую отличную библиотеку черчения bayesplot визуализировать апостериорные вероятностные интервалы из моделей, которые я оцениваю с помощью rstanarm, Я хочу графически сравнить рисунки из разных моделей, получив апостериорные интервалы для коэ…
18 окт '18 в 13:54
1
ответ
Как извлечь стандартный код из объекта rstanarm
Есть ли возможность извлечь стандартный код, используемый для выборки MCMC в rstanarm? Я хотел бы сравнить свою собственную параметризацию модели и предыдущих вариантов с той, что использовалась в rstanarm.
26 ноя '18 в 13:01
1
ответ
rstanarm для байесовского иерархического моделирования биномиальных экспериментов
Предположим, есть три биномиальных эксперимента, проведенных в хронологическом порядке. Для каждого эксперимента я знаю количество испытаний и успехов. Чтобы использовать первые два более старых эксперимента, как и в предыдущем третьем эксперименте,…
29 дек '17 в 19:38
1
ответ
Апостериорное прогнозирование на основе группирующей переменной из `stan_glm()` в пакете `rstanarm`?
Мне было интересно, как получить апостериорный прогноз на основе группирующей переменной из stan_glm() в rstanarm пакет? Например, если у меня есть двоичный файл (0, 1) закодированная группирующая переменная называется "vs" по моим данным (база данн…
21 мар '18 в 03:00
1
ответ
Что такое "linear.predictors" как извлекаемый из объекта stan_glm() в пакете "rstanarm"?
Я пишу, чтобы узнать, что "linear.predictors" как возвращено stan_glm() объект. По-видимому, "linear.predictors" не совпадает с предиктором (ами), предоставленными пользователем (документация не помогла). В любом случае, есть ли способ получить знач…
12 мар '18 в 21:31
0
ответов
Регуляризация в stan_lm
Я также пытаюсь научиться использовать stanarm, начиная с простой модели lm: summary(lm(log(rt) ~so, data = rDat)) Основываясь на некоторых подсказках в виньетках rstanarm, я запустил байесовскую версию вышеупомянутого, используя stan_glm, используя…
15 дек '17 в 18:39
1
ответ
Расчет маржинальных эффектов в биномиальном логите с использованием rstanarm
Я пытаюсь получить предельные эффекты, согласно этому посту: http://andrewgelman.com/2016/01/14/rstanarm-and-more/ td <- readRDS("some data") CHAINS <- 1 CORES <- 1 SEED <- 42 ITERATIONS <- 2000 MAX_TREEDEPTH <- 9 md <- td[,.(y,…
11 июл '17 в 14:30
1
ответ
Как поставить разные приоры используя rstanarm
Допустим, у меня есть модель формы y=a_{i} + b_{i,1}*x_{1} + b_{2}*x_{2}, где i=1,2,...,12 и я хотел бы оценить эту модель, используя rstanarm, Можно ли установить разные приоры для каждого перехвата a_{i} (так скажем, первые 4 имеют normal(location…
07 авг '17 в 09:35
1
ответ
Получение стандартизированных коэффициентов из пакета "rstanarm" в R?
Мне было интересно, возможно ли (и, возможно, рекомендуется) получить стандартизированные коэффициенты из stan_glm() в rstanarm пакет? (не нашел ничего конкретного в документации) Могу ли я просто стандартизировать все переменные, как в обычной регр…
17 мар '18 в 02:56
1
ответ
Расчет вероятных интервалов для предельных эффектов в биномиальном логите с использованием rstanarm
В этом методе для расчета предельных эффектов для биномиального логита с использованием rstanarm, /questions/26652385/raschet-marzhinalnyih-effektov-v-binomialnom-logite-s-ispolzovaniem-rstanarm/26652391#26652391 nd <- md nd$x1 <- 0 p0 <- p…
24 янв '18 в 21:28
0
ответов
rstanarm - линейная смешанная модель занимает слишком много времени, чтобы соответствовать
Я подгоняю линейную смешанную модель с stan_lmer с помощью prior_fit_1 <- normal(location = c(2.0, 0.0, 0.0), scale = c(4.0, 9.0, 9.0)) fit_1 <- stan_lmer(response ~ predictor1 + predictor2 + predictor3 (1 + predictor1 + predictor2 | group1) (…
27 июл '18 в 06:45
0
ответов
Определение `newdata` в`posterior_linpred()`из пакета rstanarm R
Мне интересно, почему я получаю ошибку при попытке использовать newdata аргумент posterior_linpred() из пакета rstanarm? Полный код ниже вместе с ошибкой, которую я получаю: PS Я следую аналогичному шаблону, как вариант 2 в этой демонстрации. librar…
11 мар '18 в 19:50
1
ответ
Прогнозирование из полного апостериорного распределения с использованием stan_glmer
Могу я попросить о помощи, пожалуйста? Я использовал биномиальную модель, используя stan_glmer, и выбрал модель, которая, я думаю, лучше всего подходит для данных. Я использовал команду апостериорного прогнозирования, чтобы сравнить мои наблюдаемые …
04 авг '17 в 13:34
1
ответ
Проблемный `rstanarm::stan_lmer` - недопустимое внутреннее подмножество
Я сталкиваюсь со следующей проблемой, которая, кажется, связана с поднабором, сделанным в rstanarm, Смотрите эту связанную проблему Ошибка в xj[i]: неверный тип индекса 'list' Воспроизводимый пример: library(stringi) library(data.table) library(Matr…
21 фев '18 в 19:01
1
ответ
rstanarm Предыдущее местоположение должно быть больше 0
Я пытаюсь вписать линейную модель в Rstanarm, используя иерархическую сжатие до. Однако я получаю сообщение об ошибке, утверждающее, что местоположение предыдущего должно быть больше 0. Error: location > 0 is not TRUE Я немного удивлен, так как h…
26 дек '17 в 16:48
1
ответ
"lines()", действующий как "polygon()" R?
Я пытаюсь нарисовать две черные линии вокруг красной линии регрессии. Но lines() Команда рисует что-то более похожее на polygon() чем простая строка (см. рисунок ниже кода). Мне интересно, есть ли способ просто нарисовать две линии вокруг линии регр…
12 мар '18 в 16:18