In this paper, we present a detailed study of the unpolarized nucleon parton distribution function (PDF) employing the approach of parton pseudo-distribution func- tions. 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 ex- trapolation systematic errors of the PDF are considered. In addition to the x dependence of the PDF, we compute their first two moments and compare them with the pertinent phenomenological determinations.