I would like to estimate the function value of the sigmoid over an expectation, that is: \begin{equation} \sigma(\mathbb{E}_{p(x)}[f(x)]), \end{equation} where $\sigma(x) = \frac{1}{1 + e^{-x}}$, and $p(x)$ is the one we only access its samples but cannot evaluate its density.
To estimate $\sigma(\mathbb{E}_{p(x)}[f(x)])$, we could use $\sigma(\frac{1}{n}\sum_{i=1}^n f(x_i))$, but it is biased. My question is, how could we define an unbiased estimator.