Inclusive jet production data are important for constraining the gluon distribution in the global QCD analysis of parton distribution functions. With the addition of recent CDF and D0 Run II jet data, we study a number of issues that play a role in determining the up-to-date gluon distribution and its uncertainty, and produce a new set of parton distributions that make use of that data. We present in detail the general procedures used to study the compatibility between new data sets and the previous body of data used in a global fit. We introduce a new method in which the Hessian matrix for uncertainties is ``rediagonalized to obtain eigenvector sets that conveniently characterize the uncertainty of a particular observable.