In a related Phys.SE question about supersymmetric Lagrangian $$ \mathcal{L} = - \frac{1}{2} (\partial S)^2 - \frac{1}{2} (\partial P)^2 - \frac{1}{2} \bar{\psi} \partial\!\!\!/ \psi, $$ the fields $S$ and $P$ are said to be Grassmann-even supernumber-valued rather than real- (or complex) valued, so that supersymmetry transformations (with Grassmann-odd $\varepsilon$) $$ \begin{align*} \delta_\varepsilon S &= \bar{\varepsilon} \psi \\ \delta_\varepsilon P &= \bar{\varepsilon} \gamma_5 \psi \\ \delta_\varepsilon \psi &= \partial\!\!\!/ (S + P \gamma_5) \varepsilon \end{align*} $$ can be consistently defined.
My question is: can you still do path integral (or canonical) quantization of Grassmann-even fields ($S$ and $P$) in the same fashion as real (or complex) scalar/pseudoscalar Boson fields? I am asking this because Grassmann-even supernumbers behave differently from real (complex) numbers. For example, the Grassmann-even $ab$ (where $a$ and $b$ are Grassmann-odd) squares to zero (nilpotent) $$ (ab)(ab) = -(aa)(bb) = 0, $$ which is totally different from a real (complex) number.