Soit $(X,Y)$ un couple aléatoire à valeurs dans $\Lambda\times J$, où $\Lambda\subset\mathbbR^d$ et $J\subset\mathbbR$ sont supposés bornés. Nous construisons un estimateur plug-in des ensembles de niveaux de la fonction de régression $r$ de $Y$ sur $X$, à partir d'un estimateur à noyaux de $r$. Nous obtenons une vitesse de convergence du même ordre que celle obtenue par Cadre dans le cas de la densité. Nous discutons ensuite les résultats obtenus sur des données simulées.