PyMC3 posterior prediction example; how to incorporate a matrix of betas?

You have to make sure that all the shapes work together. Things that come to mind and would be issues:

  • You are creating pm.Data but not using it later in the model, so updating the data won’t change anything in the model, only what ArviZ stored as constant_data
  • You are only updating only treatment and not covariates, those are added, so their shapes need to match for the model to work.
1 Like