  We calculate all twist-two, three and four parton distribution functions
involving two quark correlations using the wavefunction of the MIT bag model.
The distributions are evolved up to experimental scales and combined to give
the various nucleon structure functions. Comparisons with recent experimental
data on higher twist structure functions at moderate values of $Q^{2}$ give
good agreement with the calculated structure functions.
