The z-transform is introduced to analyze a full discretization method for a partial integro-differential equation (PIDE) with a weakly singular kernel. In this method, spectral collocation is used for the spatial discretization, and, for the time stepping, the finite difference method combined with the convolution quadrature rule is considered. The global stability and convergence properties of complete discretization are derived and numerical experiments are reported.