diff --git a/docs/api_reference.rst b/docs/api_reference.rst
index 9fd29c38..1c175ac2 100644
--- a/docs/api_reference.rst
+++ b/docs/api_reference.rst
@@ -15,12 +15,25 @@ This reference provides detailed documentation for user functions in the current
    :members:
 
 
+:mod:`preliz.multidimensional`
+======================================
+
+.. automodule:: preliz.multidimensional
+   :members:
+
+
 :mod:`preliz.predictive`
 ======================================
 
 .. automodule:: preliz.predictive
    :members:
 
+:mod:`preliz.ppls`
+======================================
+
+.. automodule:: preliz.ppls
+   :members:
+
 :mod:`preliz.distributions.continuous`
 ======================================
 
diff --git a/preliz/__init__.py b/preliz/__init__.py
index fd0a6931..dc017ca9 100644
--- a/preliz/__init__.py
+++ b/preliz/__init__.py
@@ -13,10 +13,9 @@
 from .predictive import *
 from .ppls import *
 from .unidimensional import *
+from .multidimensional import *
 
 
-__all__ = ["maxent", "mle", "ppa", "roulette", "quartile"]
-
 __version__ = "0.8.1"
 
 _log = logging.getLogger("preliz")
@@ -36,3 +35,6 @@
 _preliz_style_path = os_path.join(os_path.dirname(__file__), "styles")
 style.core.USER_LIBRARY_PATHS.append(_preliz_style_path)
 style.core.reload_library()
+
+# clean namespace
+del handler, logging, os_path, rcParams, _preliz_style_path, _log
diff --git a/preliz/multidimensional/__init__.py b/preliz/multidimensional/__init__.py
new file mode 100644
index 00000000..c1ba8c28
--- /dev/null
+++ b/preliz/multidimensional/__init__.py
@@ -0,0 +1,3 @@
+from .dirichlet_mode import dirichlet_mode
+
+__all__ = ["dirichlet_mode"]
diff --git a/preliz/unidimensional/dirichlet_mode.py b/preliz/multidimensional/dirichlet_mode.py
similarity index 100%
rename from preliz/unidimensional/dirichlet_mode.py
rename to preliz/multidimensional/dirichlet_mode.py
diff --git a/preliz/unidimensional/__init__.py b/preliz/unidimensional/__init__.py
index 33a1f5f7..58a6f596 100644
--- a/preliz/unidimensional/__init__.py
+++ b/preliz/unidimensional/__init__.py
@@ -4,6 +4,5 @@
 from .quartile import quartile
 from .quartile_int import quartile_int
 from .roulette import roulette
-from .dirichlet_mode import dirichlet_mode
 
-__all__ = ["beta_mode", "maxent", "mle", "roulette", "quartile", "quartile_int", "dirichlet_mode"]
+__all__ = ["beta_mode", "maxent", "mle", "roulette", "quartile", "quartile_int"]