You do, in fact, have the following statement:
Lemma Let $H$ be a open, simply connected domain bounded by a Jordan curve $\partial H$. If $f$ is an analytic function on $H$ that extends continuously to $\partial H$ and is injective on $\partial H$, then $f$ is univalent in $H$ and $f(H)$ is the inner domain of the Jordan curve $f(\partial H)$.
(See C. Pommerenke, Univalent Functions; Lemma 1.1 and Corollary 9.5.)
This implies that step (1) in the proof you wrote down carries over directly to the case of arbitrary analytic functions.
Step (2) can be formulated thus:
By step one it suffices to consider the zero sets of $P(z e^{i\theta}) - P(z e^{-i\theta})$, when $z \neq 0$ and $\theta \in (0,\pi/2]$. This is equivalent, in the region concerned, to looking at the zero sets of the function
$ Q(\theta;z) = \frac{P(z e^{i\theta}) - P(z e^{-i\theta})}{z (e^{i\theta} - e^{-i\theta})} $
using that the denominator is non-zero in the region concerned. Since the numerator, a difference of two holomorphic functions (in $z$), vanishes at $z = 0$, we have that $Q(\theta;z)$ is also holomorphic in $z$. Using that \lim_{\theta\to 0}Q(\theta;z) = P'(z) and Q(\theta;0) = Q(\theta';0) (so locally univalent implies that $Q(\theta;0) \neq 0$), you get that $P$ being univalent is equivalent to $Q$ having no zeroes.
Indeed, if you look at the Zbl review of Dieudonné's 1931 paper (I don't have the paper handy to check, unfortunately), the "infinite sum" version of what you want is described. (I assume that is how the theorem is stated in the paper, but I am not sure. Most citations to the criterion, however, cites this other slightly earlier paper which only deals with the polynomials.)
Thanks to the French national library, I found a digital copy of the paper in question ("Sur les fonctions univalentes"), from which I quote
La condition nécessaire et suffisante pour que la fonction (1) soit univalente dans le cercle unité, est que l'équation en $x$ $ > \phi(x,\theta) \equiv 1 + a_2 x \frac{\sin 2\theta}{\sin \theta} + > \ldots + a_n x^{n-1} \frac{\sin n\theta}{\sin \theta} + \ldots = 0 $ n'ait aucune racine dans le cercle unité, quelle que soit la valeur de $\theta$ comprise entre $0$ et $2\pi$.
where he defined the function in expression (1) to be
$ f(z) = z + a_2 z^2 + \ldots + a_n z^n + \ldots $