The paper investigates the probability of failure of slopes using both traditional and more advanced probabilistic analysis tools. The advanced method, called the random finite-element method, uses elastoplasticity in a finite-element model combined with random field theory in a Monte-Carlo framework. The traditional method, called the first-order reliability method, computes a reliability index which is the shortest distance (in units of directional equivalent standard deviations) from the equivalent mean-value point to the limit state surface and estimates the probability of failure from the reliability index. Numerical results show that simplified probabilistic analyses in which spatial variability of soil properties is not properly accounted for, can lead to unconservative estimates of the probability of failure if the coefficient of variation of the shear strength parameters exceeds a critical value. The influences of slope inclination, factor of safety (based on mean strength values), and cross corr...