EDIT: I was considering some ways to get exact equations for the partial sums of some of those 0 radii of convergence functions. One way we can get the partial sums of the geometric series is to notice that $$S_n = \sum_{k=0}^n x^k \implies xS_n - S_n = x^{n+1}-1 \implies S_n = \frac{x^{n+1}-1}{x-1}$$ However, since we have removed the summation, $S_n$ depends on something continuous and so we might expect $S_n$ to agree with the calculated partial sum even at non-integer values.
Since we already have a differential equation representation for many of our sums, it struck me that we can do something similar to obtain exact forms for the partial sum formula. In the case of $S_N = \sum_{n=1}^N (n!)^2 (n+1)x^{-n}$, we have the differential equation $$\frac{d^2}{dx^2}S_N - \frac{1}{x}S_N + \frac{1}{x^2} = (N!)^2 (N+1) x^{-(N+2)}$$
Letting Mathematica solve this gives something in terms of the Meijer-G and Bessel functions. Interestingly though, this does very closely agree with the approximations of the partial sums. However, I think this strategy can't work in general, since finding the differential equation that the function solves is probably not possible for all functions. However, it still might be possible that there is an approximation method that finds approximate differential equations, though I don't have a good idea for how to do this.