Statistical Analysis
Descriptive statistics with mean, median, standard deviation, maximum, and minimum value were used to characterize participants. Where variables were nominal or ordinal (e.g., gender), relative and absolute frequencies were reported.For the primary outcome, a generalized linear mixed-effects model of the Gaussian family with log-link was fitted to the data to estimate the change in neck pain-related work productivity loss [The statistical analyses were performed using Stata® Version 15.1 (StataCorp, College Station, Texas, USA) and R® (Boston, USA) statistical software. Significance level was set at alpha = 0.05. We reported all model estimates with corresponding 95% confidence intervals (95% CI). Data analysts were blinded to the identity and group allocation of the participants. The data were analysed on an intention-to-treat basis. The study was registered on ClinicalTrials.gov (NCT04169646, 