In this note we prove a monotonicity result with respect to the parameter $\nu$ of the cumulative distribution function for the McKay $I_\nu$ Bessel distribution and uniform upper bounds for a bilinear expression containing modified Bessel function of the first kind $I_\nu$. Certain implications, among others with the Horn function $\Phi_2$ and for the Gaussian hypergeometric function close the exposition.