Using data that include the direction and the sense of motion on individual fault surfaces determined by slickenside lineations, it is possible to reconstruct reduced stress tensors that correspond to the orientation of stress axes and to the ratio ϕ =(&sgr;2-&sgr;3)/(&sgr;1-&sgr;3) between principal stress values (&sgr;1≥&sgr;2≥&sgr;3, compression being positive). No assumption is made concerning the orientation of fault planes relative to stress axes, so that reactived faults are taken into account as well as newly created ones. Qualitative and quantitative methods for analysis of fault slip data were developed during the last 10 years. The practical limitation of the methods and the necessity for critical field observations are emphasized. These methods can be applied to focal mechanisms of earthquakes. A more complex analysis of heterogeneous data sets, involving an iterative separation of different stress systems, is also discussed. This analysis enables one to distinguish successive faulting events. Careful qualitative study in the field is in all cases essential.