quap错误（alist（紫色~dbinom（总计，p），总计<- Total[ID], p <- a[ID], : non-finite value supplied by optim

Question

library(rethinking)

G <- c("D","ND")
Purple <- c(670, 601)
Pink <- c(315,604)
Total <- c(985, 1205)
GID <- c(1,2)
Data <- data.frame(G, Purple, Pink, Total, GID)

model_1a <- quap(
  alist(
    Purple ~ dbinom(Total,p) , 
    Total <- Total[ID] ,
    p <- a[ID] ,
    a[ID] ~ dbeta(10,10)    
  ) , data = Data)
precis(model_1a, depth=2)

当我运行这个问题时，它会抛出错误：

# Error in quap(alist(Purple ~ dbinom(Total, p), Total <- Total[ID], p <- a[ID],  : non-finite value supplied by optim".

尽管我的教授在他的一次讲座中做了同样的事情，但他的模型可以运行，而我的没有。

我只是想按照教授发布的帮助视频完成我的实验室作业，我遇到了这个问题。

Answer 1

组装所有依赖项需要很长时间，但在运行代码并查看

quap

调用中的参数后，我注意到没有

ID

的定义，并想知道它是否可能是三重错别字（

ID 

->

GID

):

model_1a <- quap(
    alist(
        Purple ~ dbinom(Total,p) , 
        Total <- Total[GID] ,
        p <- a[GID] ,
        a[GID] ~ dbeta(10,10)    
    ) , data = Data)
precis(model_1a, depth=2)
#-----------------
     mean   sd 5.5% 94.5%
a[1] 0.68 0.01 0.65  0.70
a[2] 0.50 0.01 0.48  0.52

Seems to have been so.