We prove a generalized resolvent estimate of Stokes equations with nonhomogeneous Robin boundary condition and divergence condition in the $L_q$ framework $(1 < q < \infty)$ in a domain of $\bm{R}^n$ ( $n \geq 2$ ) that is a bounded domain or the exterior of a bounded domain. The Robin condition consists of two conditions: $\nu\cdot u=0$ and $\alpha u + \beta(T(u,p)\nu - \nu) = h$ on the boundary of the domain with $\alpha, \beta \geq 0$ and $\alpha + \beta = 1$ , where $u$ denotes a velocity vector, $p$ a pressure, $T(u, p)$ the stress tensor for the Stokes flow, and $\nu$ the unit outer normal to the boundary of the domain. It presents the slip condition when $\beta =1$ and the non-slip one when $\alpha = 1$ , respectively.