We address the issue of nonlinear modes in a two-dimensional waveguide array, spatially distributed in the Lieb lattice geometry, and modeled by a saturable nonlinear Schr\"odinger equation. In particular, we analyze the existence and stability of vortex-type solutions finding localized patterns with symmetric and asymmetric profiles, ranging from topological charge S=1 to S=3. By taking into account the presence of anisotropy, which is inherent to experimental realization of waveguide arrays, we identify different stability behaviors according to their topological charge. Our findings might give insight on experimental feasibility to observe these kind of vortex states.