In this paper we perform an extensive study of the existence, uniqueness (or multiplicity) and stability of nonnegative solutions to the semilinear elliptic equation −∆u = λu − u p in Ω, with the nonlinear boundary condition ∂u/∂ν = u r on ∂Ω. Here Ω is a smooth bounded domain of IR d with outward unit normal ν, λ is a real parameter and p, r > 0. We also give the precise behavior of solutions for large |λ| in the cases where they exist. The proofs are mainly based on bifurcation techniques, sub-supersolutions and variational methods.