Thank you and Daniel Eacker for helping me understand the problem! I now see that the way I bound the array (with the highest dimension last rather than first) was the root of the problem, as it was not compatible with my parametrization in the model and my initial values array. It is clear I need to strengthen my foundational understanding of how to create and manipulate multidimensional arrays so these rookie errors don't reoccur. Thanks again!
This is a follow-up question to my previous post. I am trying to incorporate another sampling level (sites) to this very basic model, and I'm getting the same error as before, indicating that the error could be related to the initial values again. However, I don't understand the problem, as z has three dimensions. z[s,i,f] ~ dbern(psi) # State model And the object zst that I supply as initial values also has three dimensions zst <- apply(y, c(1,2, 4), max) # Avoid data/model/inits #conflict Please...
That was precisely it--thank you! What a silly error. Correct call below: zst <- apply(y, c(1,2), max)
Hello, I'm very new to JAGS/hierarchical modeling. I'm trying to run this basic model, and am getting an error, which I understand to be related to the initial values that I've provided. Since I'm an extreme novice, I would be very grateful for any insights on how to resolve this error. # Bundle data and summarize data bundle > str( win.data <- list(y = y, M1 = dim(y)[1], M2 = dim(y)[2] , J = dim(y)[3] ) ) List of 4 $ y : num [1:256, 1:24, 1:5] 0 0 0 0 0 0 0 0 0 0 ... ..- attr(*, "dimnames")=List...