Edit: This function is only convex on its positive values, as pointed out in the comments. No such function exists - see my other answer.
A solution obtained via some modifications to RJTK's answer:
We set up our function so that it passes through $(0,0), (1/2,1), (3/4,e), (7/8,e^e),(15/16,e^{e^e}),$ and so on. Thus, after $n$ applications of log it will have a domain of $(1-2^{-n},1)$.
To interpolate between these points, we let $a_n=1-2^{-n}$ and define the function from $a_n$ to $a_{n+1}$ to be $\exp(\exp(\ldots\exp(g(\frac{x-a_n}{a_{n+1}-a_n})))$, where there are $n$ applications of the exponential and $g$ is a strictly convex increasing function on $[0,1]$ passing through $(0,0)$ and $(1,1)$.
Since $\log(f(x)) = f(2x-1)$ for all $x$ such that this is defined, we need only check convexity of $f$ and the convexity of each application of log will follow.
Clearly, $f$ is continuous and strictly convex on $(a_n,a_{n+1})$. To see that it is convex at each $a_n$, we want to show that the slope at the start of the interval $(a_n,a_{n+1})$ exceeds that at the end of the interval $(a_{n-1},a_n)$. But $f$ has the property that $f'(x)=2f(x)\cdot f'(2x-1)$. So, by induction and the fact that $f$ is monotonically increasing, it suffices to check convexity at $x=\frac12$.
Here is where $g$ must be chosen; we need that $g'(1) \le 2g'(0)$, by the previous identity. Taking $g(x)=\frac{3x+x^2}{4}$ suffices for this.
(Originally, I tried defining $f(x)$ piecewise between integers, but this forced $g'(1) \le g'(0)$, i.e., that $g$ not be strictly convex. Allowing for the slope to double at each successive piece gets around this, albeit in a somewhat hacky way.)