  We provide sufficient conditions on the profile $\varphi$, on the sequence of
random variables $\varepsilon_j>0$ and on the sequence of random vectors
$y_j\in\mathbb{R}^n$ such that
$\mathscr{E}\left(\frac{1}{\varepsilon_j^n(\omega)}\int_{z\in\mathbb{R}^n}\varphi\left(\frac{|x-z-y_j(\omega)|}{\varepsilon_j(\omega)}\right)f(z)
dz\right)\longrightarrow f(x)$ when $j\to\infty$ for almost every
$x\in\mathbb{R}^n$, $f\in L^p(\mathbb{R}^n)$, $1\leq p\leq\infty$, where
$\mathscr{E}$ denotes the expectation, $\varepsilon_j$ tends to
$0\in\mathbb{R}$ in law and $y_j$ tends to $\mathbf{0}\in\mathbb{R}^n$ in law.
