Regularity estimates for elliptic equations involving fractional Neumann boundary conditions