The aim of this paper is to study the following problem: { ( - Δ ) p , β s u = f ( x , u ) in Ω , u = 0 in ℝ N ∖ Ω , \left\{\begin{aligned} \displaystyle(-\Delta)^{s}_{p,\beta}u&\displaystyle=f(x% ,u)&&\displaystyle\phantom{}\text{in }\Omega,\\ \displaystyle u&\displaystyle=0&&\displaystyle\phantom{}\text{in }\mathbb{R}^{% N}\setminus\Omega,\end{aligned}\right. where Ω is a smooth bounded domain of ℝ N {\mathbb{R}^{N}} containing the origin, ( - Δ ) p , β s u ( x ) := PV ∫ ℝ N | u ( x ) - u ( y ) | p - 2 ( u ( x ) - u ( y ) ) | x - y | N + p s d y | x | β | y | β (-\Delta)^{s}_{p,\beta}u(x):=\mathrm{PV}\int_{\mathbb{R}^{N}}\frac{\lvert u(x)% -u(y)\rvert^{p-2}(u(x)-u(y))}{\lvert x-y\rvert^{N+ps}}\frac{dy}{\lvert x\rvert% ^{\beta}\lvert y\rvert^{\beta}} with 0 ≤ β < N - p s 2 {0\leq\beta<\frac{N-ps}{2}} , 1 < p < N {1<p<N} , s ∈ ( 0 , 1 ) {s\in(0,1)} , and p s < N {ps<N} . The main purpose of this work is to prove the existence of a weak solution under some hypotheses on f . In particular, we will consider two cases: (i) f ( x , σ ) = f ( x ) {f(x,\sigma)=f(x)} ; in this case we prove the existence of a weak solution, that is, in a suitable weighted fractional Sobolev space for all f ∈ L 1 ( Ω ) {f\in L^{1}(\Omega)} . In addition, if f ⪈ 0 {f\gneq 0} , we show that the problem above has a unique entropy positive solution. (ii) f ( x , σ ) = λ σ q + g ( x ) {f(x,\sigma)=\lambda\sigma^{q}+g(x)} , σ ≥ 0 {\sigma\geq 0} ; in this case, according to the values of λ and q , we get the largest class of data g for which the problem above has a positive solution.