rstanarm - 使用二进制数据更新先验,无任何事件

时间:2016-12-28 02:53:39

标签: r bayesian rstan

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)

0 个答案:

没有答案