For $0\leq r\leq 1$ and $0 < p < 1$, define the Kullback-Leibler divergence between the Bernoulli(r) and Bernoulli(p) distributions as: $D(r||p) := r\log\frac{r}{p} + \bar{r}\log\frac{\bar{r}}{\bar{p}}$ where $\bar{r}:=1-r, \bar{p}:=1-p.$

Fix $0<\theta_0,\theta_1<1$ and $p.$

Consider the function

$$\eta(r) := \frac{D(r\theta_1+\bar{r}\theta_0||p\theta_1+\bar{p}\theta_0)}{D(r||p)}.$$

This function can be made continuous at $r=p$ if we define $\eta(p)$ suitably.

Plots for different values of $p,\theta_0,\theta_1$ show that $\eta(r)$ is concave in $r$. The second derivative of $\eta(\cdot)$ seems quite involved though. Perhaps I am ignorant of some sophisticated techniques that could be useful here. I would be grateful for a proof that $\eta(\cdot)$ is concave. Or any ideas that could help in proving such a thing.