I don't know a reference, but let me discuss a strategy for a proof of the comparison that I think should work.
By Artin's comparison, if we take a constant sheaf $\mathbb Z/\ell^n$ and then pull back to the cover where we invert $f$ and adjoin an $m$th root of $f$, pushforward to $X$, and then pull back to the zero locus of $f$, the resulting complex of sheaves will not depend on if we perform these operations in the worlds of analytic or étale sheaves.
We want to use this as a building block for the comparison. The way to do this to observe there is a natural map from the above complex for $m_1$ to the above complex for $m_2$ if $m_1$ divides $m_2$, and so we can take a colimit as $m$ goes to $\infty$. After this, we can take the limit as $n$ goes to $\infty$. This gives nearby cycles as a $\mathbb Z_\ell$-sheaf, and then we can tensor with $\mathbb Q_\ell$. If we do this in the étale world, we recover the usual construction of nearby cycles.
Indeed, the only difference between what I said and the usual construction is that the usual construction involves pulling back to the algebraic closure of some local field, i.e. (in characteristic zero) to the inverse limit of all these coverings (together with étale neighborhoods of zero in $\mathbb A^1$, which are irrelevant after taking $i^*$). But the pushforward from the inverse limit of coverings is the colimit of the pushforward from coverings by a general result (applied one open set at a time, since it's stated in terms of global sections rather than sheaves) with hypotheses that are easily checked in this case.
So it suffices to check that this construction agrees with the usual construction in the analytic world, where we adjoin a logarithm of $f$ instead of an $m$'th root, and where we work with $\mathbb Q$ coefficients. I would attempt this in two steps:
First, like in the étale world, check that the cohomology of the colimits of the finite covers agrees with the cohomology of the infinite cover. I would do this by expressing the finite cover as the quotient of the infinite cover by $\log f \mapsto \log f + 2\pi i m$ so the cohomology of the finite cover can be computed from the cohomology of the infinite cover and the action of the operator $\log f \mapsto \log f + 2 \pi i m$, i.e. the $m$th power of a generator of monodromy. Thinking about the action of this operator should give the statement.
Next, check that the $\mathbb Q_\ell$-cohomology constructed this way isjust the usual $\mathbb Q$-cohomology tensored with $\mathbb Q$. This should just require using the universal coefficient theorem to calculate the $\mathbb Z/\ell^n$ cohomology in terms of the $\mathbb Z$ cohomology and use that the $\mathbb Z$-cohomology is a finitely generated abelian group (surely that is known in this setting even with the logarithm map) to see that the inverse limit over $n$ agrees with the tensor product with $\mathbb Z_\ell$, then tensor both sides with $\mathbb Q$.
Everything I say should work with stacks but now finding references to individual steps will be more difficult.
In the infinite-dimensional case, it depends what you mean, but the étale side would usually be a limit of finite-dimensional things so the question is if the analytic side is as well.