# What is wrong with this definite integral (expected value of a function of a truncated Normal random variable)?

I am solving a problem involving the following integral: $$\int_{\ln\left(\frac{(b-c)d-a}{1+b-c}\right)}^\infty \ln(a+\exp(x)-b\left|\exp(x)-d\right|-c(\exp(x)-d)) \frac{1}{\sqrt{2\pi}\sigma}\exp\left({-\frac{(x-\mu)^2}{2\sigma^2}}\right) \ dx$$ with \begin{aligned} a\geq0, \\ 0\leq c0. \end{aligned} The lower bound ensures that the argument of the logarithm is positive, so the logarithm is well defined. The integral can be interpreted as the expected value of the following function of a random variable $$X$$, $$\ln(a+\exp(X)-b\left|\exp(X)-d\right|-c(\exp(X)-d)),$$ where $$X\sim \text{truncated}\ N(\mu,\sigma^2)$$ where the left tail of the Normal distribution is truncated at $$\ln\left(\frac{(b-c)d-a}{1+b-c}\right)$$.

I have mostly forgotten the little calculus that I once knew, so I resorted to some online solvers such as www.integral-calculator.com and www.Desmos.com. They said Antiderivative or integral could not be found and Undefined, respectively. I then turned to Sage Math but did not manage to make it work either (this is most likely my own fault; I am just a beginner at Sage Math).

However, I think the integral should be doable and its value should be finite.

• As $$x\downarrow-\infty$$, the logarithm asymptotes along $$\ln(a)$$.
As $$x\uparrow\infty$$, the logarithm asymptotes along $$(1-B-C)x$$.
Thus the function behaves at most as a linear function of $$x$$.
• The (truncated) normal distribution has finite first moment (and many more moments than that).
• Therefore, the integral should be finite.

What surprised me even more was that the solvers failed to calculate a finite value of the integral even when concrete values of $$a,\ b,\ c,\ d,\ \mu,\ \sigma$$ were supplied, at least for some sets of such values (e.g. when $$d$$ is relatively large).

I would appreciate any help getting this integral done. I have plotted all the functions involved in the integral using Desmos here. You can check it out and adjust with the parameters very easily to see how the functions behave. Perhaps that will be helpful in helping me. A screenshot is presented below:

The dashed blue line is the argument of $$\ln(\cdot)$$.
The violet line is $$\ln(a+\exp(x)-b\left|\exp(x)-d\right|-c(\exp(x)-d))$$.
The green line is the density of the random variable w.r.t. which I am integrating.
The red line is the integrand (the product of the logarithm and the density).
The black line is the lower bound of the definite integral.

Update: As suggested in the comments, finding a closed-form expression for the integral may be impossible. However, I think I found an explanation for why the integral failed numerically for concrete values of $$a,\ b,\ c,\ d$$. Though the integrand is asymptotically linear for $$x\rightarrow+\infty$$, its calculation requires exponentiation (followed by a logarithm). When executed naively, exponent of a large number quickly exceeds the limits of standard software, so numerical integration fails. A (probably) naive workaround is to introduce an upper bound $$ub$$ such that conditions A and B are satisfied.

• Condition A: the density $$\frac{1}{\sqrt{2\pi}\sigma}\exp\left({-\frac{(x-\mu)^2}{2\sigma^2}}\right)$$ is numerically equal to zero.
• Condition B: the exponent $$\exp(x)$$ is small enough so that software can handle it without returning infinity or error or the like.

This seems be working OK for me in R software using the function integrate.

• Just out of curiosity : where did you catch this monster ? Could you give a bit of context ? Commented Aug 21, 2020 at 9:53
• Thanks for the explanation ! I am very, very skeptical if you look for a closed form. Cheers :-) Commented Aug 21, 2020 at 9:59
• @ClaudeLeibovici, but do you agree that the integral should be finite? I am not getting a finite value even when I insert concrete numbers for $a$, $b$, $c$ and $d$ and try numerical integration. This is puzzling me... Commented Aug 21, 2020 at 10:02
• Could you give me a set of parameters $(a,b,c,d,\mu,\sigma)$ ? I should try. Commented Aug 21, 2020 at 10:04
• @ClaudeLeibovici, for example, $a=0,\ b=0.1,\ c=0.05,\ d=1,\ \mu=0,\ \sigma=0.05$. I said $a>0$ in the post, but actually it is $a\geq 0$; I will edit that. Commented Aug 21, 2020 at 10:14

Let's denote your integral with $$I$$, so $$I = \int_{\ln\left(\frac{(b-c)d-a}{1+b-c}\right)}^{\infty} \ln(a+\exp(x)-b|\exp(x)-d|-c(\exp(x)-d)) \frac{1}{\sqrt{2\pi}\sigma}\exp\left({-\frac{(x-\mu)^2}{2\sigma^2}}\right) \ dx$$

We have: \begin{align*}&\quad\;a+\exp(X)−b|\exp(X)−d|−c(\exp(X)−d) \\ &= a+d+\exp(X)-d −b|\exp(X)−d|−c(\exp(X)−d) \\ &=a+d+(\exp(X) - d)(1-b\cdot \text{sgn}(X-d) - c) \end{align*}

By the assumptions we have $$0 \le (1-b\cdot \text{sgn}(X-d) - c) \le 1$$

and so we get: $$-d \le (\exp(X) - d)(1-b\cdot \text{sgn}(X-d) - c) \le \exp(X)$$ hence $$a \le a+\exp(X)−b|\exp(X)−d|−c(\exp(X)−d) \le a+d+\exp(X)$$

And we get: $$\ln(a)\left(1-\Phi\left(\ln\left(\frac{(b-c)d-a}{1+b-c}\right)\right)\right) \le I \le E[\ln(a+d+\exp(X))] < +\infty$$ The bounds are not pretty sharp but enough to show your integral is finite hence exists.

I leave it to yourself to show that $$E[\ln(r+\exp(X))]$$ exists for $$r > 0$$ if $$X \sim \mathcal{N}(\mu,\sigma^2)$$ and you will get the result by setting $$r=a+d$$.

• Thank you, this is helpful. Commented Aug 25, 2020 at 9:28