We perform numerical fits of Grand Unified Models based on SO(10), using various combinations of 10-, 120- and 126-dimensional Higgs representations. Both the supersymmetric and non-supersymmetr