We analyse the recent precision measurements of the lepton-hadron deep inelastic scattering at CERN and SLAC to extract model independent constraints among the nucleon matrix elements of the twist-4 operators. We also study a parameterization of these matrix elements and point out the possibility that the matrix elements of the quark-gluon mixed operator has a negative value of the order of \(- (400 \pm 100\ {\rm MeV})^2\) at 5 GeV\(^2\) renormalization scale.