Partial answer
Let $g(x)=xf(x)-4$.
Our functional equation becomes
$$\frac{g\circ g(x)+4}{g(x)}=4+\frac1x$$
If there exists an invertible function $\phi(x):\mathbb R^+\to\mathbb R^+$ such that $g(x)=\phi(\phi^{-1}(x)+1)$, direct substitution gives
$$\frac{\phi(z+2)+4}{\phi(z+1)}=4+\frac1{\phi(z)}$$
or $$\phi(z+2)=\left[4+\frac1{\phi(z)}\right]\phi(z+1)-4$$
by substituting $z=\phi^{-1}(x)$.
Clearly, this recurrence relation extends a $\phi(x)$ defined arbitrarily on $(0,2)$ to the whole $\mathbb R^+$. By computing $\phi^{-1}(x)$ this method generates a large class of solution to the functional equation.
(Of course there are certain restrictions on $\phi(x)$ on $(0,2)$ so that $\phi(x)$ is invertible.)
The solution $f(x)=4+\frac4x$ corresponds to $g(x)=4x$ and $\phi(x)=k\cdot 4^x$ where $k$ is a positive constant.
A feature of this special solution is that $\phi(x)$ is continuous. In contrast, if rather arbitrary values are assigned to $\phi(x)$ on $(0,2)$, a discontinuous solution may yield. Below is the graph of $\phi(x)$ where $\phi(x) = x+1$ on $(0,2)$.
