All variables are labeled (labels are displayed in the dataset spreadsheet) and de-identified in the excel spreadsheet. Missing data were left blank in the excel spreadsheet. Data were analyzed in R software (version 4.1.3) package. Analysis involved descriptive statistics (frequencies and standard deviations) to present sociodemographic characteristics, participants’ experience with COVID-19 and anxiety and depression levels according to the score in GAD-7 and PHQ-9 test. In order to calculate participants’ level of depression each item of the PHQ-9 questionnaire was rated on a 4-point Likert scale and a final total score was calculated. Answers to every item range from 0 (not at all) to 3 (nearly every day). More specifically, “not at all”: 0 points, “Several days”: 1 point, “More than half the days”:2 points and “nearly every day”:3 points. The sum of all 9 items (presented in the study questionnaire file) was used to determine the level of depression for every participant. Total scores ranging from 0 to 4 were rated as “minimal or none”, 5 to 9 as “mild”, 10 to14 as “moderate”, 15 to 19 as “moderately severe”, and 20 to 27 as “severe”. A cut-off point of 10 or above is indicative of major depressive disorder. In order to calculate participants’ level of anxiety, each item of the GAD-7 questionnaire, 7-item questionnaire (presented in the study questionnaire file) was rated on a 4-point Likert scale and answers to every item range from 0 (not at all) to 3 (nearly every day). More specifically, “not at all”: 0 points, “Several days”: 1 point, “More than half the days”:2 points and “nearly every day”:3 points. The sum of all 7 items was used to determine the level of anxiety for every participant. Total scores ranging from 0 to 4 were classified as “not at all”, 5 to 9 as “mildly”, 10 to 14 as “moderately” and ≥15 as “severely”. A cut-off point of 10 or above corresponds to moderate to severe anxiety disorder. After calculating the GAD-7 and PHQ-9 score for each participant, univariate analysis was conducted to evaluate the association between sociodemographic factors and participants’ experience with COVID-19 with the levels of anxiety and depression respectively. Multivariable logistic regression in R was used to determine the association between independent variables (socio-demographic characteristics and participants’ experience with COVID-19) which yielded statistically significant results at the univariate analysis (gender, age, facility location, previous infection, relative/friend hospitalized or deceased, being in a high risk group) with the dichotomous dependent variables determined by the cut-off point of 10 in the GAD-7 and PHQ-9 indicating clinically significant levels of anxiety and depression, respectively. Two-tailed p values of less than 0.05 were deemed statistically significant.