Sample_prior_predictive doesn't work with Negative Binomial likelihood

Questions
2

Does it also fail if you explicitly pass the shape to the negative binomial? Something like shape=bike_data["temp_std"].to_numpy().shape?