rstanarm - 在没有任何事件的情况下使用二进制数据更新先验
嗨,我在上第一门概率课程时开始探索 rstanarm 包。我的问题不是作业,而是在我阅读https://cran.r-project.org/web/packages/rstanarm/vignettes/binomial.html后出于好奇
想象一下,在这个小插图中提供的示例(参见下面的代码)中,开关变量只有前 10 个观察值,没有事件(即所有值都为零)。尝试在没有任何事件的情况下运行此模型会产生错误,但据我所知,即使没有实例,也应该有一种更新先验的方法。
我的问题:是否可以使用 rstanarm 来模拟一个有 10 个没有事件的观察的案例(所有行都用于零值的开关)?
t_prior <- student_t(df = 7, location = 0, scale = 2.5) fit1 <- stan_glm(switch ~ dist100, data = wells, family = binomial(link = "logit"), prior = t_prior, prior_intercept = t_prior, chains = 4, cores = 2, seed = 1234, iter = 100)