In Section 3.1, scores of the COVID-19 anxiety variable were categorized into quintiles, and the quintiles were dummy-coded with the lowest one used as the reference category [20] Regression models were used to estimate the relationship between the COVID-19 anxiety quintiles and the PHQ-15 subscales (pain, gastrointestinal, cardiopulmonary, and fatigue). Model 1 included the four dummy-coded COVID-19 anxiety variables as predictors of the four PHQ-15 subscales. The regression coefficient for each dummy-coded variable is interpreted as the mean difference between each quintile and the lowest one. Model 1 was also run separately with the total PHQ-15 summed scale score replacing the subscales. In Model 2 the covariates (age, gender, income, pre-existing health problems, and GAD) were included as predictors, with the addition of Italian region mourning for COVID-19 losses factors.