The solution I fond is to use taylor expansions for the first moment of functions of random variables. The details can be found here: vignette_taylor_series