We explain the origins of yield stress, shear thinning, and normal stress differences in rigid fiber suspensions. We investigate the interplay between the hydrodynamic, noncontact attractive and repulsive, and interfiber contact interactions. The shear-thinning viscosity and finite yield stress obtained from the Immerse Boundary Method simulations are in quantitative agreement with experiential results from the literature. In this study, we show that attractive interactions result in yield stress and shear thinning rheology in the suspensions of rigid fibers. This is an important finding, given the ongoing discussion regarding the origin of the yield stress for suspensions of fibers. The ability of the proposed model to quantitatively predict the rheology is not limited to only shear thinning and yield stress but also extends to normal stresses.