diff --git a/preliz/__init__.py b/preliz/__init__.py
index f7085769..b62485a2 100644
--- a/preliz/__init__.py
+++ b/preliz/__init__.py
@@ -20,7 +20,7 @@
 from preliz.internal.citations import citations
 
 
-__version__ = "0.27.0"
+__version__ = "0.27.1"
 
 
 # Allow legend outside plot in maxent to be included in the saved figure
diff --git a/preliz/ppls/agnostic.py b/preliz/ppls/agnostic.py
index 3f9e6f38..d18eaba1 100644
--- a/preliz/ppls/agnostic.py
+++ b/preliz/ppls/agnostic.py
@@ -262,11 +262,11 @@ def looper(*args, **kwargs):
             model = func(*args, **kwargs)
             model.build()
             if var_to_plot is None:
-                var_to_plot = model.observed_RVs[0].name
+                var_to_plot = model.backend.model.observed_RVs[0].name
 
             with disable_pymc_sampling_logs():
                 idata = model.prior_predictive(iterations)
-            results = idata[group].stack(sample=("chain", "draw"))[var_to_plot].values.T
+            results = idata[group].dataset.stack(sample=("chain", "draw"))[var_to_plot].values.T
             if group == "prior":
                 results = np.atleast_2d(results)
 
@@ -275,8 +275,8 @@ def looper(*args, **kwargs):
                 if var_to_plot is None:
                     var_to_plot = model.observed_RVs[0].name
                 with disable_pymc_sampling_logs():
-                    idata = sample_prior_predictive(samples=iterations)
-                results = idata[group].stack(sample=("chain", "draw"))[var_to_plot].values.T
+                    idata = sample_prior_predictive(draws=iterations)
+                results = idata[group].dataset.stack(sample=("chain", "draw"))[var_to_plot].values.T
             if group == "prior":
                 results = np.atleast_2d(results)