If you want to compare the asymptotic behaviour of two functions, average values don't turn out to be an efficient tool.
Let $f,g$ be those two functions and let's assume they are positive for the sake of simplicity. The simplest way to compare them would be to compute the ratio $\displaystyle\lim_{x\rightarrow\infty}\frac{f(x)}{g(x)}$. If the result is above 1, then $f$ grows asymptotically faster than $g$; if the result is under 1, then $g$ grows faster; and if it is equal to 1, the two functions are asymptotically equal.
If you want to study the $f,g$ individually, you can determine their leading terms (using the big-O notation) and then compare them. You can also compare their growth rates for large values, i.e. the asymptotic behaviour of their derivatives.
The problem with average values is that you must compute them on finite intervals, otherwise the integrals won't converge, unless the functions are $L^1$, but then it means that $\displaystyle \lim_{x\rightarrow\infty}f(x) = \lim_{x\rightarrow\infty}g(x) = 0$, hence a trivial asymptotic behaviour. One way to circumvent this problem would be a continuous version of moving average, but that seems troublesome, since the integral of the gamma function doesn't possess a closed form anyway.
In conclusion, I would recommend you the first method I mentioned, taking into account the asymptotic behaviour of the factorial given by Stirling's formula, namely $n! \sim \sqrt{2\pi} \,n^{n+1/2}e^{-n}$.