This paper explores application of the Finite Element Method (FEM) and Shear Strength Reduction (SSR) analysis to compute probabilities of failure for slopes. It does so using two probabilistic approaches: the Point Estimate Method (PEM) and limited numbers of Monte Carlo simulations. The paper explains why probabilistic analysis with numerical methods such as the FEM is challenging, and how the PEM and Monte Carlo simulations can be used to calculate the statistical moments of output variables and to estimate slope probability of failure. One of the paper's two examples describes application of probabilistic FEM analysis to determine slope probability of failure due to the random distribution of joints (joint networks).