  The relation between the Star Formation Rate (SFR) and stellar mass (${\rm
M}_{\star}$) of galaxies represents a fundamental constraint on galaxy
formation and has been studied extensively both in observations and
cosmological simulations. However, the observed amplitude has not been
successfully reproduced in simulations, indicating either that the halo
accretion history and baryonic physics are poorly modeled or that observations
contain biases. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation
of $z\sim1-4 $ galaxies and display the inconsistency between observed
relations that are obtained using different techniques. We employ cosmological
hydrodynamic simulations and compare these with a range of observed SFR$-{\rm
M}_{\star}$ relations. We find that numerical results are consistent with
observations that use Spectral Energy Distribution (SED) techniques to estimate
star formation rates and dust corrections. On the contrary, simulations are not
able to reproduce results that were obtained by combining only UV and IR
luminosities. These imply SFRs at a fixed stellar mass that are larger almost
by a factor of 5 than those of SED measurements for $z \sim1.5-4$. Furthermore,
we find remarkable agreement between the numerical results from various authors
who have employed different cosmological codes and run simulations with
different resolutions. This is interesting for two reasons. A) simulations can
produce realistic populations of galaxies within representative cosmological
volumes even at relatively modest resolutions. B) It is likely that current
numerical codes that rely on similar subgrid multiphase Inter-Stellar Medium
(ISM) models and are tuned to reproduce statistical properties of galaxies,
produce similar results for the SFR$-{\rm M}_{\star}$ relation by construction,
regardless of resolution, box size and, to some extent, the adopted feedback
prescriptions.
