A came across following log-likelihood function with variables $x$ and $y$ with domain $x+y\leq0$ and some constant $c > 1$: \begin{equation} \log\mathcal{L}(x,y) = \log\left(1 + \log_c\left(1-\frac{1-c^{-1}}{e^{-x}+1}\right)+\log_c\left(1-\frac{1-c^{-1}}{e^{-y}+1}\right)\right) + \ldots \end{equation} I was able to show that all terms of this log-likelihood function are concave except for the first shown term. Hence, if this term was also concave, the whole log-likelihood function would be concave. Plotting this first term for various $c$ suggests that it is concave over the domain $x+y\leq 0$. See for example for $c=2$:
Does anyone know how to prove that it is really concave?