warning in MCMC: logProb of data node Counts[1]: logProb is -Inf.

8 views
Skip to first unread message

Agnieszka Tomczyk

unread,
May 28, 2026, 3:35:22 PM (5 days ago) May 28
to nimble-users
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!
Reply all
Reply to author
Forward
0 new messages