R and SAS code for running multilevel hurdle negative binomial models with cyclical covariates