Sparse logistic regression for RR Lyrae vs binaries classification

RR Lyrae (RRL) are old, low-mass radially pulsating variable stars in their core helium burning phase. They are popular stellar tracers and primary distance indicators, since they obey to well defined period-luminosity relations in the near-infrared regime. Their photometric identification is not trivial, indeed, RRL samples can be contaminated by eclipsing binaries, especially in large datasets produced by fully automatic pipelines. Interpretable machine-learning approaches for separating eclipsing binaries from RRL are thus needed. Ideally, they should be able to achieve high precision in identifying RRL while generalizing to new data from different instruments. In this paper, we train a simple logistic regression classifier on Catalina Sky Survey (CSS) light curves. It achieves a precision of 87% at 78% recall for the RRL class on unseen CSS light curves. It generalizes on out-of-sample data (ASAS/ASAS-SN light curves) with a precision of 85% at 96% recall. We also considered a L1-regularized version of our classifier, which reaches 90% sparsity in the light-curve features with a limited trade-off in accuracy on our CSS validation set and -- remarkably -- also on the ASAS/ASAS-SN light curve test set. Logistic regression is natively interpretable, and regularization allows us to point out the parts of the light curves that matter the most in classification. We thus achieved both good generalization and full interpretability.