The issue isn’t that there is no data. These sorts of issues arise when there is a mismatch between the data you have supplied via pm.set_data() and the rest of the model. Here, it seems like the shape of your observed data observed_eaches and your now-modified indices are incompatible. Indeed, I’m not exactly sure what you are trying to do. Can we see what is in coords?