Given $x \sim N(0, I)$ and $f : \mathbb{R}^{n} \longrightarrow \mathbb{R}$ which is $L$-Lipschitz, we have that $f - \mathbb{E}f$ is a sub-Gaussian random variable, specifically that $$ \| f(x) - \mathbb{E} f(x) \|_{\psi_2} \leq C L \:. $$
I want to show that if in addition $f \geq 0$, then for all $p > 0$, the more general statement $$ \| f(x) - (\mathbb{E} f(x)^p)^{1/p} \|_{\psi_2} \leq C_p L $$ holds, where $C_p$ is a constant which depends on $p$. This is Exercise 5.2.5 from Vershynin's book: http://www-personal.umich.edu/~romanv/teaching/2015-16/626/HDP-book.pdf (and the notation above is using his notation).
Using the fact that if $X$ is sub-gaussian then $$ \| X\|_p \leq C \|X\|_{\psi_2} \sqrt{p} \:, \:\: p \geq 1 \:, $$ (Eq 2.16 in Vershynin's HDP book), I was able to show that (just by triangle inequality) $$ \| X - \|X\|_p \|_{\psi_2} \leq (C_1 + C_2 \sqrt{p}) \| X \|_{\psi_2} \:. $$
I wanted then to apply this fact to the sub-Gaussian r.v. $Z := f(x) - \mathbb{E} f(x)$, but this doesn't work for $p \leq 1$ since the constant does not come out of the $p$-th moment of $Z$, so I don't think this approach is correct. Furthermore, it only applies when $p \geq 1$.
Any other hints/suggestions?