In this paper, we obtain gradient estimates of the positive solutions to weighted p-Laplacian type equations with a gradient-dependent nonlinearity of the form
0.1
$${\rm div }( \vert x \vert ^\sigma \vert \nabla u \vert ^{p-2}\nabla u) = \vert x \vert ^{-\tau }u^q \vert \nabla u \vert ^m\quad {\rm in}\;\Omega^*: = \Omega {\rm \setminus }\{ 0\} .$$ Here,
$\Omega \subseteq {\open R}^N$ denotes a domain containing the origin with
$N\ges 2$, whereas
$m,q\in [0,\infty )$,
$1<p\les N+\sigma $ and
$q>\max \{p-m-1,\sigma +\tau -1\}$. The main difficulty arises from the dependence of the right-hand side of (0.1) on x, u and
$ \vert \nabla u \vert $, without any upper bound restriction on the power m of
$ \vert \nabla u \vert $. Our proof of the gradient estimates is based on a two-step process relying on a modified version of the Bernstein's method. As a by-product, we extend the range of applicability of the Liouville-type results known for (0.1).