A recursive formulation is used to retrieve sinusoids from nonuniformly sampled data. The approach combines EDC-type estimators to select the number of sinusoids. The algorithm enables to retrieve the frequencies above the Nyquist frequency. The procedure can efficiently model a time series with linear trend in data. The approach achieves the Cramer–Rao lower bound above the threshold SNR.