Different approaches (Mean Field theory, Monte Carlo simulations and a MF-MC hybrid model) to predict order parameters of model-solutes in nematic phases have been tested and compared among them and with experimental data (when available). Peculiarities of the methods are discussed in order to explore their potentialities in tackling the orientational problem. (C) 2003 Elsevier B.V.. All rights reserved.