  We study the global Q^2 dependence of large x, F_2 nucleon structure function
data, with the aim of providing a perturbative-QCD based, quantitative analysis
of parton-hadron duality. As opposed to previous analyses at fixed x, we use a
framework in fixed W^2. We uncover a breakdown of the twist-4 approximation
with a renormalon type improvement at O(1/Q^4) which, by affecting the initial
evolution of parton distributions, will have consequences for pQCD analyses
also at large x and very large Q^2.
