Rocking response of rigid bodies with rectangular footprint, freely standing on horizontal rigid plane is studied analytically. Bodies are subjected to simulated single component of horizontal earthquakes. The effect of baseline correction, applied to simulated excitations, on the rocking response is first examined. The sensitiveness of rocking motion to the details of earthquakes and geometric properties of rigid bodies is investigated. Due to the demonstrated sensitivity of rocking response to these factors, prediction of rocking stability must be made in the framework of probability theory. Therefore, using a large number of simulated earthquakes, the effects of duration and shape of intensity function of simulated earthquakes on overturning probability of rigid bodies are studied. In the case when a rigid body is placed on any floor of a building, the corresponding probability is compared to that of a body placed on the ground. For this purpose, several shear frames are employed. Finally, the viability of the energy balance equation, which was introduced by Housner in 1963 and widely used by nuclear power industry to estimate the rocking stability of bodies, is evaluated. It is found that the equation is robust. Examples are also given to show how this equation can be used.