Estimating covariate functions associated to multivariate risks: a level set approach