Welcome @Sebastian_Botero! Here’s one solution (include a binary indicator for whether there is missing data for a sampling occasion): Ragged data and mixture (occupancy) model
You can additionally provide a fill value for your observation matrix y
like -1, to make sure you get an error if you accidentally pass it to a bernoulli or binomial pmf.