Let 1 < p < +∞ and let Ω ⊂ ℝN be either a ball or an annulus. We continue the analysis started in [Boscaggin, Colasuonno, Noris, ESAIM Control Optim. Calc. Var. (2017)], concerning quasilinear Neumann problems of the type $-\Delta _pu = f(u),\quad u > 0\,{\rm in }\,\Omega ,\quad \partial _\nu u = 0\,{\rm on }\,\partial \Omega .$
$-\Delta _pu = f(u),\quad u > 0\,{\rm in }\,\Omega ,\quad \partial _\nu u = 0\,{\rm on }\,\partial \Omega .$