Hi Guys,
I am getting a below error in my model:
warning: logProb of data node Counts[1]: logProb is -Inf.
What may be the reason?
Model details:
This is a spacial analysis. I am modelling the counts of events (var = Counts), with 2 explanatory variables (var1, var2) and one offset - population (var = Pop). No spacial random effects are added.
My code:
model <- nimbleCode({
beta0 ~ dnorm(0, sd=100)
beta1 ~ dnorm(0, sd=100)
beta2 ~ dnorm(0, sd=100)
for(i in 1:N){
log(mu[i]) <- log(Pop[i]) + beta0 + (beta1*Var1[i]) + (beta2*Var2[i])
Counts[i] ~ dpois(mu[i])
}
})
Data_Nimble <- select(data, Counts, Pop, Var1, Var2) %>%
st_drop_geometry() %>%
as.list()
Constants1 <- list(N = nrow(data))
Inits1 <- list(beta0 = 3,
beta1 = 3,
beta2 = 3)
model.nimble <- nimbleMCMC(model, data=Data_Nimble,
constants=Constants1,
niter = 1000,
inits = Inits1,
samplesAsCodaMCMC = TRUE, summary=TRUE,
monitors=c("beta0","beta1","beta2"),
setSeed = 1994)
Thanks for your help!