Improve example on expanded sample_posterior_predictive

Sample posterior_predictive relies on other groups to infer volatility, so it's important to keep the whole InferenceData.

Also there's a shape error if `sample_dims` is not specified

Author: ricardoV94

pymc/sampling/forward.py

@@ -711,9 +711,14 @@ def sample_posterior_predictive(
 
     .. code:: python
 
-        expanded_data = idata.posterior.expand_dims(pred_id=5)
+        expanded_idata = idata.copy()
+        expanded_idata.posterior = idata.posterior.expand_dims(pred_id=5)
         with model:
-            idata.extend(pm.sample_posterior_predictive(expanded_data))
+            pm.sample_posterior_predictive(
+                expanded_idata,
+                sample_dims=["chain", "draw", "pred_id"],
+                extend_inferencedata=True,
+            )
 
 
     """