我正在尝试使用Rethinking R软件包来拟合一个非常简单的模型来估计疾病的患病率，如第6页上的https://www.ecdc.europa.eu/sites/default/files/documents/Methodology-estimating-point-prevalence%20-SARS-CoV-2-infection-pooled-RT-PCR-testing.pdf所述，
这是我的代码:
quap(alist(
p ~ dbeta(.3, .3),
p_test ~ 1 - dbinom(0, s, p), # I tried also p_test <- 1 - dbinom(0, s, p)
k ~ dbinom(w, p_test)
), data = list(s = 10, k = 30, w = 200))
但我收到错误:
Error in pars[[i]] : subscript out of bounds in quap我究竟做错了什么？
最佳答案
重新定义alist中的定义可以解决问题。
quap(
alist(
k ~ dbinom(w, p_test),
p ~ dbeta(.3, .3),
p_test ~ 1 - dbinom(0, s, p)
),
data = list(s = 10, k = 30, w = 200)
)
此代码返回:
Quadratic approximate posterior distribution
Formula:
k ~ dbinom(w, p_test)
p ~ dbeta(0.3, 0.3)
p_test ~ 1 - dbinom(0, s, p)
Posterior means:
p
0.01575975
Log-likelihood: -2.55
关于r - 使用Rethinking R软件包拟合贝叶斯模型时出错，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/62212666/