We perform the NLO QCD fit to the combined deep inelastic scattering (DIS)
data at high x from the SLAC, BCDMS, and NMC collaborations. The model-inde
pendent x shape of the high-twist contribution to structure function F-L is
extracted. The twist-4 contribution to F-L is found to be in a qualitative
agreement with the predictions of the infrared renormalon model. The twist
-6 contribution exhibits a weak trend to negative values, although on the w
hole, it is compatible with zero within the errors.