diff --git a/Project.toml b/Project.toml
index 8283a4f..737874f 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "ApproxManifoldProducts"
 uuid = "9bbbb610-88a1-53cd-9763-118ce10c1f89"
 keywords = ["MM-iSAMv2", "SLAM", "inference", "sum-product", "belief-propagation", "nonparametric", "manifolds", "functional"]
-version = "0.4.8"
+version = "0.4.9"
 
 [deps]
 Base64 = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"
diff --git a/src/Legacy.jl b/src/Legacy.jl
index b00367c..6810bc9 100644
--- a/src/Legacy.jl
+++ b/src/Legacy.jl
@@ -111,7 +111,6 @@ Ndim(x::ManifoldKernelDensity, w...;kw...) = Ndim(x.belief,w...;kw...)
 Npts(x::ManifoldKernelDensity, w...;kw...) = Npts(x.belief,w...;kw...)
 
 getWeights(x::ManifoldKernelDensity, w...;kw...) = getWeights(x.belief, w...;kw...)
-Random.rand(x::ManifoldKernelDensity, d::Integer=1) = rand(x.belief, d)
 
 getKDERange(x::ManifoldKernelDensity, w...;kw...) = getKDERange(x.belief, w...;kw...)
 getKDEMax(x::ManifoldKernelDensity, w...;kw...) = getKDEMax(x.belief, w...;kw...)
diff --git a/src/services/ManifoldKernelDensity.jl b/src/services/ManifoldKernelDensity.jl
index 9e201dd..d61273b 100644
--- a/src/services/ManifoldKernelDensity.jl
+++ b/src/services/ManifoldKernelDensity.jl
@@ -162,6 +162,9 @@ function sample(x::ManifoldKernelDensity{M,B,L,P}, N::Int) where {M,B,L,P}
   vecP, lbls
 end
 
+Random.rand(mkd::ManifoldKernelDensity, N::Integer=1) = sample(mkd, N)[1]
+
+
 
 function resample(x::ManifoldKernelDensity, N::Int)
   pts, = sample(x, N)