We compared the method of least squares (LS), Pope’s iterative data snooping (IDS) and Huber’s M-estimator (HU) in realistic leveling networks, for which the heights or the vertical displacements of points are known. The study was conducted using the Monte Carlo simulation, in which one repeatedly generates sets of observations related to the measurement data, then calculates values of the estimators and, finally, assesses it with respect to the real coordinates. To simulate outliers we used popular mixture models with two or more normal distributions. It is shown that for small, strong networks robust methods IDS and HU are more accurate than LS, but for large, weak networks occurring in practice there is no significant difference between the considered methods in the accuracy of the solution.