3.1. Somatic Symptoms
Correlations for all the variables included in Models 1 and 2 are shown in Table 1. The PHQ-15 subscales and the total score were positively and significantly correlated with COVID-19 anxiety, as well as with GAD-7 total score. These correlations are in line with results from the first wave of the UK study [20]. Figure 1 shows the relationship between COVID-19 anxiety and somatic symptoms.
Table 2 shows the estimates from the regression models. In Model 1, considering COVID-19 first quintile as no-anxiety, we can observe significant differences between the baseline and the other quantiles both in PHQ-15 subscales and in its total score. In the sample, the aggravation of the somatic symptoms between quintiles seems to be characterized by a more linear increase than in the UK data, with a similar wider effect over the 5th quintile.
In Model 2, the control variables were taken into account. The effect of COVID-19 anxiety on somatic symptoms was weaker but similar to Model 1, especially referred to the total scale score. The presence of pre-existing health problems and high scores on GAD-7 significantly worsen specific and general somatic symptoms, confirming UK’s findings. Mourning for COVID-19 losses also has a detrimental effect on psychological health. No effect of household income was found. In terms of regions, Campania showed higher scores of somatic symptoms than the other regions. Campania scores were high in the PHQ-15 pain subscale and in the gastrointestinal subscale, but also overall.