As in this previous question, let a $k-$almost prime be a positive integer having exactly $k$ prime factors, not necessarily distinct. Let $\mathbb{P}_k$ be the set of the $k$-almost primes and let $$ \rho_k(n):=\sum\limits_{\substack{q\in \mathbb{P_k}\\q\le n}}\frac1q. $$ The answer to that question states that (I presume for fixed $k$) the asymptotic estimate for $\rho_k(n)$ has leading term $$ \rho_k(n) \asymp \frac{1}{k!}(\log \log n)^k, $$ and error term of $O((\log \log n)^{k-1}).$
How would a rigorous proof showing the error term proceed?
Is there a reference to such a result in the literature?