Вопросы по теме 'stan'
Определение новых функций и дистрибутивов на языке BUGS/JAGS/STAN
Я новичок в мире статистического анализа и недавно заинтересовался языком моделирования BUGS/JAGS/STAN. Что меня действительно удивляет, так это то, что я не видел ни одного примера определения новых функций или дистрибутивов, чтобы избежать...
2034 просмотров
schedule
12.06.2023
Стандартная ошибка: требуется объявление неограниченной переменной. нашел симплекс
В Стэне я получаю следующую ошибку:
SYNTAX ERROR, MESSAGE(S) FROM PARSER:
require unconstrained variable declaration. found simplex.
ERROR at line 48
46: for (j in 1:records) {
47: real phenology_predictor;
48:...
1055 просмотров
schedule
01.08.2023
Как получить стандартные ошибки оценок максимального правдоподобия в STAN?
Я использую оптимизацию максимального правдоподобия в Stan, но, к сожалению, функция optimizing() не сообщает о стандартных ошибках:
> MLb4c <- optimizing(get_stanmodel(fitb4c), data = win.data, init = inits)
STAN OPTIMIZATION COMMAND...
3817 просмотров
schedule
01.07.2023
Как преобразовать файл модели STAN в файл Graphviz DOT или другое графическое представление?
У меня есть файл STAN, описывающий иерархическую модель. Я хотел бы визуализировать эту иерархию со всеми параметрами, преобразовав код STAN в файл Graphviz DOT. Подойдет и другое графическое представление.
Рассмотрим следующий небольшой пример:...
433 просмотров
schedule
14.04.2022
Код Rstan для простой многомерной линейной модели
Я пытаюсь использовать Rstan, чтобы соответствовать модели из примеров Байесовских идей и анализа данных Кристенсена, Джонсона, Бранскума и Хэнсона: Введение для ученых и статистиков. Авторы используют WinBUGS, поэтому необходима некоторая адаптация....
2988 просмотров
schedule
23.06.2022
Модели смесей в Stan - векторизация
Сейчас я изучаю Стэна и хотел реализовать простую модель смеси.
В справочнике (stan-reference-2.14.0) уже есть решение:
data {
int<lower=1> K; // number of mixture components
int<lower=1> N; // number of data points
real y[N];...
1252 просмотров
schedule
05.05.2022
Указание целочисленной скрытой переменной в stan
Я изучаю байесовский анализ данных. Я пытаюсь воспроизвести учебные пособия от Trond Reitan от stan, которые изначально были созданы WinBugs.
В частности, у меня есть следующие данные и модель
weta.windata<-list(numdet=c(0, 0, 1, 0, 0, 0,...
1626 просмотров
schedule
19.01.2023
Полиномиальный логит с использованием rstanarm
Я знаю, что вы можете подобрать биномиальную логит-модель с пакетом rstanarm , используя stan_glm и установив family = binomial(link="logit") , но дает ли rstanarm возможность подобрать полиномиальную логит-модель?
Я нашел это...
1239 просмотров
schedule
04.10.2022
Гиперприоры для иерархических моделей со Стэном
Я ищу подходящую модель для оценки множественных вероятностей биномиальных данных с помощью Стэна. Я использовал бета-априоры для каждой вероятности, но я читал об использовании гиперприоров для объединения информации и поощрения сокращения оценок....
866 просмотров
schedule
09.09.2022
Подавлять запущенные сообщения от RStan в выводе Rmarkdown HTML
Ниже приведен мой код r.
```{r message=FALSE, warning=FALSE, cache=0,eval=TRUE, error=FALSE}
stan_m1 <- rethinking::map2stan(
alist(
y ~ dbinom(n, p),
logit(p) <- alpha + bP*P + bA*A + bV*V,
alpha ~ dnorm(0, 10),
bP ~...
2038 просмотров
schedule
12.11.2022
Как использовать stan в rmarkdown
Я хотел бы получить оценочные коэффициенты модели, используя rstan в rnotebook
У меня есть следующий блок stan :
```{stan output.var="rats"}
data {
int<lower=0> N;
int<lower=0> T;
real x[T];
real y[N,T];
real xbar;
}...
2744 просмотров
schedule
22.05.2022
Проблемный `rstanarm::stan_lmer` - недопустимое внутреннее подмножество
Я столкнулся со следующей проблемой, которая, по-видимому, связана с подмножеством, выполненным в rstanarm . См. эту связанную проблему
Ошибка в xj[i]: неверный тип индекса 'список'
воспроизводимый пример:
library(stringi)...
86 просмотров
schedule
20.08.2022
Значения гистограммы искажены при использовании `ppc_stat()`
Я использую следующие данные:
speed <- c(28, 26, 33, 24, 34, -44, 27, 16, 40, -2, 29, 22, 24, 21, 25, 30, 23, 29, 31, 19, 24, 20, 36, 32, 36, 28, 25, 21, 28, 29, 37, 25, 28, 26, 30, 32, 36, 26, 30, 22, 36, 23, 27, 27, 28, 27, 31, 27, 26, 33,...
66 просмотров
schedule
14.01.2024
Реализация экспоненциальной общей линейной модели в Stan / RStan
У меня такая модель:
Я хочу узнать, как реализовать эту модель.
Исследование и попытка
В качестве примера я видел следующий пуассоновский GLM:
```{stan output.var="PoissonGLMQR"}
data{
int<lower=1> n;...
618 просмотров
schedule
08.05.2022
Как rstan хранит апостериорные образцы для отдельных цепочек?
Я хотел бы понять, как вывод extract в rstan упорядочивает апостериорные образцы. Я понимаю, что могу просмотреть апостериорные образцы из каждой цепочки, используя as.array ,
stanfit <- sampling(
model,
data = stan.data)
fitarray...
281 просмотров
schedule
05.02.2024
Заказанная оценка пробита в стандарте
Я пытаюсь воспроизвести упорядоченную пробит-модель JAGS из книги Джона Крушке «Проведение байесовского анализа» (стр. 676) в Стэне:
Модель JAGS:
model {
for ( i in 1:Ntotal ) {
y[i] ~ dcat( pr[i,1:nYlevels] )
pr[i,1] <-...
156 просмотров
schedule
26.12.2022
Задание априорного распределения на матрице в rstan
У меня возникли проблемы с получением байесовской модели смешанных эффектов для получения стационарных и хорошо перемешанных цепочек. Я создал свои собственные данные, поэтому знаю, какие параметры должна извлекать модель. К сожалению, из-за того,...
283 просмотров
schedule
25.02.2022
Преобразование модели Jags в модель стана
Я пытаюсь преобразовать модель Jags в модель Stan.
Зубцы:
model{
for (i in 1:n){
theta[i] ~ dbeta(u*s, s-u*s)
y[i] ~ dbin(theta[i],N[i])
}
u ~ dunif(0,1)
s ~ dlnorm(4,0.25)
}
Стэн:
data {...
255 просмотров
schedule
07.06.2022
Есть ли функция R, которая может помочь мне изменить размер точек диаграммы рассеяния на парижском графике?
Я использовал pairs(fit, pars = c("mu", "sigma", "lambda"), include = TRUE) для создания гистограмм и диаграмм рассеяния между каждой парой переменных моей модели стана. Сейчас точки на диаграмме рассеяния очень большие, например:
Я хочу...
103 просмотров
schedule
14.10.2022
альтернатива заполнению данных в stan
Мне нужно передать данные стана, где векторы одной и той же переменной имеют разную длину (из-за отсутствия данных). Я видел здесь , что нет возможности создать массив с векторами разной длины, но поскольку с момента публикации прошло некоторое...
11 просмотров
schedule
17.06.2023