The relay feedback system often shows an asymmetric response because of initial transient states, disturbances, and process nonlinearity. To restore a symmetric response of the relay feedback system, an iterative adjustment of input or output bias is required. Instead of trying to obtain a symmetric response, the asymmetric response can be analyzed to estimate the ultimate data or a parametric model of process. However, the asymmetric response causes additional errors in estimating the ultimate properties of the process. An area method is proposed to reduce these errors. Because integrals (areas) of the relay responses reduce the effects of the high-order harmonic terms significantly, the proposed method shows better accuracy in obtaining frequency response data and parametric models compared with previous approaches.