The Interferometric Synthetic Aperture Radar (InSAR) technique has been widely used to obtain the ground surface deformation of geohazards (e.g., mining subsidence and landslides). As one of the inherent errors in the interferometric phase, the digital elevation model (DEM) error is usually estimated with the help of an a priori deformation model. However, it is difficult to determine an a priori deformation model that can fit the deformation time series well, leading to possible bias in the estimation of DEM error and the deformation time series. In this paper, we propose a method that can construct an adaptive deformation model, based on a set of predefined functions and the hypothesis testing theory in the framework of the small baseline subset InSAR (SBAS-InSAR) method. Since it is difficult to fit the deformation time series over a long time span by using only one function, the phase time series is first divided into several groups with overlapping regions. In each group, the hypothesis testing theory is employed to adaptively select the optimal deformation model from the predefined functions. The parameters of adaptive deformation models and the DEM error can be modeled with the phase time series and solved by a least square method. Simulations and real data experiments in the Pingchuan mining area, Gaunsu Province, China, demonstrate that, compared to the state-of-the-art deformation modeling strategy (e.g., the linear deformation model and the function group deformation model), the proposed method can significantly improve the accuracy of DEM error estimation and can benefit the estimation of deformation time series.