In this paper, we present a detailed study of the unpolarized nucleon parton distribution function (PDF) employing the approach of parton pseudo-distribution functions. We perform a systematic analysis using three lattice ensembles at two volumes, with lattice spacings a = 0.127 fm and a = 0.094 fm, for a pion mass of roughly 400 MeV. With two lattice spacings and two volumes, both continuum limit and infinite volume extrapolation systematic errors of the PDF are estimated. In addition to the x dependence of the PDF, we compute their first two moments and compare them with the pertinent phenomenological determinations.