rstanarm - 使用二进制数据更新priors而不会发生任何事件
嗨,我开始探索rstanarm包,同时学习第一道概率课程。我的问题不是作业,而是出于好奇心,因为我读了https://cran.r-project.org/web/packages/rstanarm/vignettes/binomial.html
想象一下,在这个插图中提供的示例中(参见下面的代码),switch变量只有前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)