Analysis of the psychological effects of Covid-19 in Europe during the first weeks of the pandemic.The survey on which this study is based was conducted by international researchers around the world, trying to map out all the factors, that might affects people’s psychological wellbeing and their ability to make good decisions during the COVID-19 (or “Coronavirus”) outbreak. Researchers from many countries are collaborating on this project, to help scientists and decision makers to help and communicate.
The project was chosen as it deals with analyzing all the factors that in this emergency situation can have threatening consequences in people’s lives: the isolation, the economic situation and the stress due to the stringent measures to stop the pandemic are part of the side effects of covid-19. The survey on which the next analyzes are based was designed by social scientists and psychologists from more than 50 different universities. The study in question presents more than 150,000 individual respondents from more than 50 different countries. The following analyzes describe a series of variables selected in the survey, taking as a sample the responses provided by citizens located in Europe who answered the survey between March 30, 2020 and May 30, 2020.
For more informations check this Nature article.
As said in the ‘Preamble’ this project is focuses on the responses provided by citizens located in Europe’s countries.
#filtering by EU Countries
<- c("Austria","Belgium","Bulgaria","Croatia","Cyprus","Czech Republic",
target "Denmark","Estonia","Finland","France","Germany","Greece","Hungary",
"Ireland","Italy","Latvia","Lithuania","Luxembourg","Malta",
"Netherlands","Poland","Portugal","Romania","Slovakia","Slovenia",
"Spain","Sweden")
<- filter(cleaned_d, Country %in% target)
EU
#fixing date visualization
$RecordedDate <- as.Date(EU$RecordedDate , format= "%Y-%m-%d") EU
An important point of the project is based on the analysis of perceived individual stress, obtained through the Perceived Stress Scale, based on the answers to 10 questions, invented by psychologists Cohen, Karmak, and Mermelstrein (1983), using indicators of stress responses, for instance, perceived lack of control over events, pressure from mounting difficulties and feeling upset about unexpected changes. Scores are considered moderate above 2.4, and high above 3.7. More informations about PS Scale here
PSS scores are obtained by reversing responses (e.g., 0 = 4, 1 = 3, 2 = 2, 3 = 1 & 4 = 0) to the four positively stated items (items 4, 5, 7, & 8) and then summing across all scale items. A short 4 item scale can be made from questions 2, 4, 5 and 10 of the PSS 10 item scale.
#note that the na values are omitted
$Scale_PSS10_UCLA_4[EU$Scale_PSS10_UCLA_4 == 1] <- 5
EU$Scale_PSS10_UCLA_4[EU$Scale_PSS10_UCLA_4 == 2] <- 4
EU$Scale_PSS10_UCLA_4[EU$Scale_PSS10_UCLA_4 == 4] <- 2
EU$Scale_PSS10_UCLA_4[EU$Scale_PSS10_UCLA_4 == 5] <- 1
EU$Scale_PSS10_UCLA_5[EU$Scale_PSS10_UCLA_5 == 1] <- 5
EU$Scale_PSS10_UCLA_5[EU$Scale_PSS10_UCLA_5 == 2] <- 4
EU$Scale_PSS10_UCLA_5[EU$Scale_PSS10_UCLA_5 == 4] <- 2
EU$Scale_PSS10_UCLA_5[EU$Scale_PSS10_UCLA_5 == 5] <- 1
EU$Scale_PSS10_UCLA_7[EU$Scale_PSS10_UCLA_7 == 1] <- 5
EU$Scale_PSS10_UCLA_7[EU$Scale_PSS10_UCLA_7 == 2] <- 4
EU$Scale_PSS10_UCLA_7[EU$Scale_PSS10_UCLA_7 == 4] <- 2
EU$Scale_PSS10_UCLA_7[EU$Scale_PSS10_UCLA_7 == 5] <- 1
EU$Scale_PSS10_UCLA_8[EU$Scale_PSS10_UCLA_8 == 1] <- 5
EU$Scale_PSS10_UCLA_8[EU$Scale_PSS10_UCLA_8 == 2] <- 4
EU$Scale_PSS10_UCLA_8[EU$Scale_PSS10_UCLA_8 == 4] <- 2
EU$Scale_PSS10_UCLA_8[EU$Scale_PSS10_UCLA_8 == 5] <- 1
EU
<- select(EU,Country, Scale_PSS10_UCLA_1:Scale_PSS10_UCLA_10) %>%
eu_stress ddply( .(Country), summarize,
Rate_PSS10_UCLA1=mean(Scale_PSS10_UCLA_1,na.rm=TRUE),
Rate_PSS10_UCLA2=mean(Scale_PSS10_UCLA_2,na.rm=TRUE),
Rate_PSS10_UCLA3=mean(Scale_PSS10_UCLA_3,na.rm=TRUE),
Rate_PSS10_UCLA4=mean(Scale_PSS10_UCLA_4,na.rm=TRUE),
Rate_PSS10_UCLA5=mean(Scale_PSS10_UCLA_5,na.rm=TRUE),
Rate_PSS10_UCLA6=mean(Scale_PSS10_UCLA_6,na.rm=TRUE),
Rate_PSS10_UCLA7=mean(Scale_PSS10_UCLA_7,na.rm=TRUE),
Rate_PSS10_UCLA8=mean(Scale_PSS10_UCLA_8,na.rm=TRUE),
Rate_PSS10_UCLA9=mean(Scale_PSS10_UCLA_9,na.rm=TRUE),
Rate_PSS10_UCLA10=mean(Scale_PSS10_UCLA_10,na.rm=TRUE))%>%
mutate(total_stress = rowMeans(select(., -Country)))
Levels of stress were moderate or lower in many countries. Poland and Croatia reported the highest levels of stress in Europe, Denmark and Netherlands the lowest.
<- select(filter(EU,Dem_gender=="Male" | Dem_gender=="Female"), Dem_gender,
eu_stress_gender :Scale_PSS10_UCLA_10) %>%
Scale_PSS10_UCLA_1ddply( .(Dem_gender), summarize,
Rate_PSS10_UCLA1=mean(Scale_PSS10_UCLA_1,na.rm=TRUE),
Rate_PSS10_UCLA2=mean(Scale_PSS10_UCLA_2,na.rm=TRUE),
Rate_PSS10_UCLA3=mean(Scale_PSS10_UCLA_3,na.rm=TRUE),
Rate_PSS10_UCLA4=mean(Scale_PSS10_UCLA_4,na.rm=TRUE),
Rate_PSS10_UCLA5=mean(Scale_PSS10_UCLA_5,na.rm=TRUE),
Rate_PSS10_UCLA6=mean(Scale_PSS10_UCLA_6,na.rm=TRUE),
Rate_PSS10_UCLA7=mean(Scale_PSS10_UCLA_7,na.rm=TRUE),
Rate_PSS10_UCLA8=mean(Scale_PSS10_UCLA_8,na.rm=TRUE),
Rate_PSS10_UCLA9=mean(Scale_PSS10_UCLA_9,na.rm=TRUE),
Rate_PSS10_UCLA10=mean(Scale_PSS10_UCLA_10,na.rm=TRUE))%>%
mutate(total_stress = rowMeans(select(., -Dem_gender)))
From what emerges women result more stressed then men. In line with this Psychological article.
<- select(EU, Dem_employment,Scale_PSS10_UCLA_1:Scale_PSS10_UCLA_10) %>%
eu_stress_empl ddply( .(Dem_employment), summarize,
Rate_PSS10_UCLA1=mean(Scale_PSS10_UCLA_1,na.rm=TRUE),
Rate_PSS10_UCLA2=mean(Scale_PSS10_UCLA_2,na.rm=TRUE),
Rate_PSS10_UCLA3=mean(Scale_PSS10_UCLA_3,na.rm=TRUE),
Rate_PSS10_UCLA4=mean(Scale_PSS10_UCLA_4,na.rm=TRUE),
Rate_PSS10_UCLA5=mean(Scale_PSS10_UCLA_5,na.rm=TRUE),
Rate_PSS10_UCLA6=mean(Scale_PSS10_UCLA_6,na.rm=TRUE),
Rate_PSS10_UCLA7=mean(Scale_PSS10_UCLA_7,na.rm=TRUE),
Rate_PSS10_UCLA8=mean(Scale_PSS10_UCLA_8,na.rm=TRUE),
Rate_PSS10_UCLA9=mean(Scale_PSS10_UCLA_9,na.rm=TRUE),
Rate_PSS10_UCLA10=mean(Scale_PSS10_UCLA_10,na.rm=TRUE))%>%
mutate(total_stress = rowMeans(select(., -Dem_employment)))
The most stressed category is the one of the students, as expected, due to isolation and social distancing. This article stated that Polish University’s students recorded an higher level of stress during the Covid-19 pandemic. Notably, Poland is also the most stressed Country.
<- select(filter(EU,RecordedDate=="2020-03-30"|
eu_stress_2months =="2020-04-06"|
RecordedDate=="2020-04-12"|
RecordedDate=="2020-04-18"|
RecordedDate=="2020-04-24"|
RecordedDate=="2020-04-30"|
RecordedDate=="2020-05-05"|
RecordedDate=="2020-05-11"|
RecordedDate=="2020-05-17"|
RecordedDate=="2020-05-23"|
RecordedDate=="2020-05-30"),
RecordedDate:Scale_PSS10_UCLA_10) %>%
RecordedDate,Scale_PSS10_UCLA_1ddply( .(RecordedDate), summarize,
Rate_PSS10_UCLA1=mean(Scale_PSS10_UCLA_1,na.rm=TRUE),
Rate_PSS10_UCLA2=mean(Scale_PSS10_UCLA_2,na.rm=TRUE),
Rate_PSS10_UCLA3=mean(Scale_PSS10_UCLA_3,na.rm=TRUE),
Rate_PSS10_UCLA4=mean(Scale_PSS10_UCLA_4,na.rm=TRUE),
Rate_PSS10_UCLA5=mean(Scale_PSS10_UCLA_5,na.rm=TRUE),
Rate_PSS10_UCLA6=mean(Scale_PSS10_UCLA_6,na.rm=TRUE),
Rate_PSS10_UCLA7=mean(Scale_PSS10_UCLA_7,na.rm=TRUE),
Rate_PSS10_UCLA8=mean(Scale_PSS10_UCLA_8,na.rm=TRUE),
Rate_PSS10_UCLA9=mean(Scale_PSS10_UCLA_9,na.rm=TRUE),
Rate_PSS10_UCLA10=mean(Scale_PSS10_UCLA_10,na.rm=TRUE))%>%
mutate(total_stress = rowMeans(select(., -RecordedDate)))
Level of stress remained fairly stable over the middle of April, with a negligible decrease in the second half of the month. Notably, stress level increased at the end of May and that is odd because with the beginning of summer, people should be less stressed. We do not have any data that will confirm this trend.
The analysis will focus on the cause of stress of the respondents, obtaining the data through asking them questions regarding various factors that represent a source of stress.
<- select(EU, Expl_Distress_1:Expl_Distress_24)%>%
eu_stress_source na.omit()
= data.frame(apply(eu_stress_source,2,function(x)mean(x[x<99]))) %>%
eu_stress_source ::rownames_to_column(var = "col1") %>%
tibble`colnames<-`(c("stress_source", "mean"))
<-(c("No religious activities", "Not knowing how to stop covid-19",
listlab"Feeling ashamed for acting differently","Adapt work to Digital Platforms",
"Access to necessities(food etc..)","Behavior of adults I live with",
"Behavior of childrens I live with","Adapt to social life on digital platforms",
"No trovels outside my Country","Loneliness","Time i spend in proximity to others",
"Children's education","Civil services(Police,sanitations..)","Income",
"Not knowing about developments with Covid-19","Time I spend inside","Work",
"Job prospect","No social activities","Worry over friends and relatives who live far away",
"Not knowing how long the measures will last","Risk of catching covid-19",
"Risk of being hospedalized or dying","National economy"))
The results indicate that most respondents are concerned about the state of the national economy. Economic considerations were followed closely by health-related risks, such as the risks of being hospitalized and of contracting the new disease.
Partecipants were asked how much they trusted six major institutions in relation to the Covid-19 emergency (on a scale from 1 = not at all to 10 = completely), in particular partecipants were asked about their trust towards World Health Organization (WHO), the national governments’ efforts to tackle the COVID19, the Police, the civil service, the national government and the healthcare systems.
Here there are OECD guidelines on measuring institutions trust.
<- select(EU, OECD_insititutions_1:OECD_insititutions_6)%>%
eu_trust na.omit()
= data.frame(apply(eu_trust,2,function(x)mean(x))) %>%
eu_trust ::rownames_to_column(var = "col1") %>%
tibble`colnames<-`(c("trust", "mean"))
<-(c("Country's Governments","Country's civil Service","Country's Police",
lablist"Governments effort against covid-19","WHO","Country's Healthcare system"))
Europeans reported only medium level of trust at all, with the highest levels of confidence for their countries’ healthcare system and the WHO. Notably, Countries governments reported the lowest level of trust.
<- select(EU,Country,OECD_insititutions_1) %>%
eu_mean_oecd1 ddply( .(Country), summarize,
Rate_OECD_insititutions_1=mean(OECD_insititutions_1,na.rm=TRUE))%>%
arrange(desc( Rate_OECD_insititutions_1))
Notably, Poland is the most stressed Nation but it is also one of the country in which the population trust least in their own government. Finland, on the other hand, is the nation with the least stress related to Covid, and it is also one of the countries in which the population places an high level of compliance for their own government.
<- select(EU,Country,OECD_insititutions_4) %>%
eu_mean_oecd4 ddply( .(Country), summarize,
Rate_OECD_insititutions_4=mean(OECD_insititutions_4,na.rm=TRUE))%>%
arrange(desc( Rate_OECD_insititutions_4))
Similar patterns as the trust level in the governments, with Bulgaria that reported the lowest level of trust for the national effort to handle the pandemic.
It was also asked participants to judge the appropriateness of the countries’ measures in response to the COVID-19 on a scale from 0 (too little), to 5 (appropriate, the black dashed line in the graph), to 10 (too much).
<- select(EU,Country,Trust_countrymeasure) %>%
eu_trustinc group_by(Country) %>%
::summarize(Mean = mean(Trust_countrymeasure, na.rm=TRUE)) dplyr
Slovakia and Slovenia were those countries where citizens considered the measures stronger than necessary; in the plot, the vertical mid-line indicates that the measures are appropriate. France and Hungary citizens judged their governments’ response less than appropriate.
<- select(EU,Country,OECD_insititutions_6,Compliance_1:Compliance_6) %>%
eu_conf_corr ddply( .(Country), summarize,
Rate_OECD=mean(OECD_insititutions_6,na.rm=TRUE),
Rate_Compliance_1=mean(Compliance_1,na.rm=TRUE),
Rate_Compliance_2=mean(Compliance_2,na.rm=TRUE),
Rate_Compliance_3=mean(Compliance_3,na.rm=TRUE),
Rate_Compliance_4=mean(Compliance_4,na.rm=TRUE),
Rate_Compliance_5=mean(Compliance_5,na.rm=TRUE))%>%
mutate(total_Compl = rowMeans(select(., -Country,-Rate_OECD)))
= lm( total_Compl ~Rate_OECD, data = eu_conf_corr)
mod summary(mod)
Call:
lm(formula = total_Compl ~ Rate_OECD, data = eu_conf_corr)
Residuals:
Min 1Q Median 3Q Max
-0.27076 -0.09898 0.03277 0.09216 0.26395
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.20673 0.15954 26.367 <2e-16 ***
Rate_OECD 0.03408 0.02464 1.383 0.179
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 0.1386 on 25 degrees of freedom
Multiple R-squared: 0.07106, Adjusted R-squared: 0.0339
F-statistic: 1.912 on 1 and 25 DF, p-value: 0.1789
$coefficients mod
(Intercept) Rate_OECD
4.20672802 0.03407768
From what emerges Poland not having confidence in the measures taken by governments and does not even try to reduce transmission, while Denmark is in the opposite situation. In countries with higher trust towards governments ’efforts, citizens were also more likely to report higher levels of compliance with directives aimed at controlling the spread of the virus. Notably, Portugal reported higher levels of compliance than the model would predict based on the acceptance of government efforts. An outlier in the other direction, Latvia reported lower levels of compliance than those the model predicted, given levels of trust.
r = cor(p$total_stress, p$Rate_OECD_insititutions_1)) (
[1] -0.7951361
^2 r
[1] 0.6322414
In conclusions, Student is the most stressed category, as expected, due to isolation and social distancing; Poland is the most stressed Nation in Europe and as reported in this article Polish University’s Students recorded an highly level of stress during the Covid-19 pandemic. Women result more stressed than men, as stated in this Psychological article. Remarkably, around the Europe the Country’s healthcare system is trusted followed by the WHO. There is a strong negative correlation between an high level of stress and a lack of trust in governments, indeed citizens of northern Europe’s countries are less stressed and they have also an high level of confidence in their governments.
Click Here to start.