I want to study the impact of heat waves (0 for non-heat wave days, 1 for heat wave days) on mortality during the warm period over a 20-year span. The code I use to create the crossbasis function using the dlnm package in R is: