We introduced least absolute shrinkage and selection operator (lasso) in obtaining periodic signals in unevenly spaced time-series data. A very simple formulation with a combination of a large set of sine and cosine functions has been shown to yield a very robust estimate, and the peaks in the resultant power spectra were very sharp. We studied the response of lasso to low signal-to-noise data, asymmetric signals and very closely separated multiple signals. When the length of the observation is sufficiently long, all of them were not serious obstacles to lasso. We analyzed the 100-year visual observations of delta Cep, and obtained a very accurate period of 5.366326(16) d. The error in period estimation was several times smaller than in Phase Dispersion Minimization. We also modeled the historical data of R Sct, and obtained a reasonable fit to the data. The model, however, lost its predictive ability after the end of the interval used for modeling, which is probably a result of chaotic nature of the pulsations of this star. We also provide a sample R code for making this analysis.