In this paper we consider solutions to the reflection equation related to the higher spin stochastic six vertex model. The corresponding higher spin R-matrix is associated with the affine quantum algebra U q ( sl(2)). The explicit formulas for boundary K-matrices for spins s = 1/2, 1 are well known. We derive difference equations for the generating function of matrix elements of the K-matrix for any spin s and solve them in terms of hypergeometric functions. As a result we derive the explicit formula for matrix elements of the K-matrix for arbitrary spin. In the lower-and upper-triangular cases, the K-matrix simplifies and reduces to simple products of q-Pochhammer symbols.then we obtain K J (y) J J = 1 by using the definition (3.13) of the Φ function. Applying the q-Vandermonde summation formula (A.7) it is easy to check that J j=0 K J (y) l j = 1, (5.5) J j=0 K J (y) l j = J l=0(−vt) l (q −2J ; q 2 ) l (q 2 ; q 2 ) l = (−vtq −2J ; q 2 ) J , (7.19)