This is the most natural approach. Since you seem to be seeking different methods I will mention a generalization (Wilson Reflection Formula) from a slightly different viewpoint.
An equivalent way to state Wilson's theorem is that any complete system of representatives of nonzero remainders mod $\,p\,$ has product $\equiv -1.\,$ In particular this is true for any sequence of $\,p\,$ consecutive integers, after removing its $\rm\color{#c00}{multiple}$ of $\,p.\,$ Thus the standard remainder sequence when left-shifted by $\,k\,$ i.e. $\,{-}k,\ldots,-1,\require{cancel}\cancel{\color{#c00}0,} 1,2,\ldots,\, p\!-\!1\!-\!k\,$ has product $\,(-1)^k k!\, (p\!-\!1\!-\!k)!\equiv -1\pmod{\!p},\, $ so
$$\qquad\qquad (p\!-\!1\!-\!k)!\, \equiv\, \dfrac{(-1)^k}{k!}\! \pmod{\!p}\qquad [\text{Wilson Reflection Formula}]$$
The OP is the special case $\, k = 1\,$ (leftshift by $1)$
Remark $\ $ The essence of Wilson's theorem is group-theoretical, so if you know a little group theory I highly recommend that you look at some prior posts on the group-theoretic viewpoint, which more clearly highlight the innate involution symmetry (negation/inversion "reflection")