Refactoring Data-Driven Model Selection Code for Improvements in Interpretability, Generality, and Computational Expense

CO anomalies in each of the regions were modeled using climate indices for influential climate modes. A linear modeling approach was used, where de-trended, de-seasonalized, regionally aggregated CO measurements were taken as the response variable, and the climate index anomaly values (at various time lags) were taken as explanatory variables.

Initial analyses were completed in MATLAB using serial algorithms carried out over non-functionalized scripts. We sought to refactor this codebase, with 3 specific improvement goals; first, to improve code interpretability in preparation for public release; second, to improve code generality, so that the techniques and code used in this application can be easily adapted for similar problems; and third, to utilize parallel computing to substantially speed up program executions.