11 Resilience 1

11.1 Intended Learning Outcomes

By the end of this chapter you should be able to:

Explain what resilience is
Convert variables to different data types using as.factor() and as.character()
Begin to problem solve how to reshape data to long-form

11.2 Activity 1: Resilience

For the next set of chapters we're going to use data from a measurement of resilience.

First, take part in this online version of the Brief Resilience Coping Scale by Sinclair and Wallston (2004).
Second, read the Wikipedia entry on resilience. It's quite long so you don't have to read the entire thing (although feel free if you're interested), the main sections to focus on are Overview, Process, and Developing and Sustaining Resilience.
Finally, answer the following questions. Please note that your responses will not save in the browser - if you want to save them, make a note of them somewhere.

Which of these is a strategy suggested to enhance resilience?
Psychological resilience is mainly characterized as a:
Resilient individuals do not experience any increase in stress in response to stressors.
Social support systems are considered external factors influencing a person's level of resilience.
The presence of stress in life is necessary for the practice and development of resilience.

Cognitive-behavioral techniques are suggested for enhancing resilience because they involve actively changing thought patterns and behaviors to foster a more adaptive response to stress.
Psychological resilience is termed a "Process" because it is seen as a dynamic set of interactions between a person and their environment, which can be developed over time.
The statement is "False" because resilient individuals do experience stress; however, they generally manage to return to baseline stress levels more effectively than less resilient individuals.
"True", social support systems are external factors that contribute to a person's resilience by providing emotional and practical resources during times of stress.
"True", the presence of stress provides individuals with opportunities to practice resilience, although the optimal level of stress for this purpose varies among individuals.

11.3 Activity 2: New project

Open RStudio and make a new project for the Resilience chapters (this week and next week).
Create a new Markdown document named "Resilience 1". Delete the default text and then create a new code chunk.
Your environment should be clear but if there are objects in it, remove them by pressing the brush icon.

11.4 Activity 3: Data files

Once you've done all this, it's time to download the files we need and then unzip them and put them in the right place.

First, download the Resilience data zip file to your computer and make sure you know which folder you saved it in.
Then unzip the files and move them into the right folder.

The data we're working with is simulated data for an experiment where people have completed a resilience questionnaire and also take part in a treatment condition. They either are in the control group where they are asked to do a colouring book every night, or they're in an intervention group where they're asked to complete CBT exercises every night. There are 890 participants in total.

The zip file contains three files:

demographic_data contains the participant ID, age, gender, and treatment condition (1 = control, 2 = intervention).
questionnaire_data contains the participant ID and then their responses to each question on the Brief Resilience Scale.
scoring.csv contains the scoring information for the questionnaire as some of the items are reverse-scored (also known as backward scoring).

11.5 Activity 4: Loading the data

In code chunk 1, write and run the code that:

Loads the tidyverse.
Uses read_csv() to load demographic_data.csv into an object named demo_data.
Loads questionnaire_data.csv into an object named q_data.
Loads scoring.csv into an object named scoring

library(package_name)
object_name <- read_csv("file_name.csv")

library(tidyverse)
demo_data <- read_csv("demographic_data.csv")
q_data <- read_csv("questionnaire_data.csv")
scoring <- read_csv("scoring.csv")

Is demographic_data in wide-form or long-form? Is questionnaire_data in wide-form or long-form?

They're both in wide-form at the moment because each dataset has a single row for each participant, with multiple columns for each bit of information.

11.6 Activity 5: Check the data

As always, your first step should be to check the data. Run summary() and str() on both demo_data and q_data and answer the following questions:

What type of data is gender stored as?
What type of data is treatment stored as?
What type of data is participant_ID stored as?
How many missing data points does bounce_back_quickly have (these are coded as NA)?
How many missing data points does age have?

11.7 Activity 6: Convert data types

A useful next step with data processing is to fix any data types as it can prevent problems further down the line. We've done this before but we'll extend it a little now and show you other functions to convert.

First, let's do a conversion you've already done - we'll use factor() to convert treatment to a factor and update the labels. We'll store this in a new object named demo_cleaned but we will just overwrite the existing variable treatment rather than creating a new one.

demo_cleaned <- demo_data %>%
  mutate(treatment = factor(treatment, # the variable you're converting
                            levels = c(1,2), # original values
                            labels = c("control", "intervention"))) # new labels

If you now run summary(demo_cleaned), you'll see that treatment is now coded as a factor and it helpfully tells us how many participants are in each group.

gender is also a categorical variable that we might want to include in our analysis but it requires a slightly different, and simpler, approach because the values are already stored as text, R just doesn't yet know that these represent discrete categories. Instead of using factor() and updating the levels and labels, we can just use as.factor(), which just overwrites whatever is in that column as a factor.

demo_cleaned <- demo_data %>%
  mutate(treatment = factor(treatment, 
                            levels = c(1,2), 
                            labels = c("control", "intervention")),
         gender = as.factor(gender)) # simpler version when no need to update labels

There are also similar functions for other data types. Although it won't really cause us any problems, we can convert participant_ID. It's currently stored as a numeric variable (because the IDs are just numbers) but really it's a character variable (it's just the same as a name).

demo_cleaned <- demo_data %>%
  mutate(treatment = factor(treatment, 
                            levels = c(1,2), 
                            labels = c("control", "intervention")),
         gender = as.factor(gender),
         participant_ID = as.character(participant_ID)) # convert to character

Now create an object named q_cleaned that cleans the q_data object by converting participant_ID to a character variable.

q_cleaned <- q_data %>%
  mutate(participant_ID = as.character(participant_ID))

11.8 Activity 7: Join and reshape the data

We can now do our first join and join together the demographic and questionnaire data so we've got a single wide-form dataset. In this code we create an object named dat_wide that joins together demo_cleaned and q_cleaned.

dat_wide <- inner_join(x = demo_cleaned, y = q_cleaned, by = "participant_ID")

Once we've done that, we need to reshape it into long-form using pivot_longer(). See if you can remember how to do this using just the hints. Reshaping data is probably the hardest conceptual leap you've taken since first learning to code and it only gets easier if you try and think through the problem yourself.

Name the object you're going to create dat_long
Start with the dataset dat_wide
Use pivot_longer()
The first column you want to pivot is named bounce_back_quickly and the last column is named long_time_over_setbacks.
You want to send the names of the columns to a variable called item
You want to send the values to a variable called response

object_name <- starting_data %>%
  pivot_longer(cols = first_var:last_var, 
               names_to = "name", 
               values_to = "value")

dat_long <- dat_wide %>%
  pivot_longer(cols = bounce_back_quickly:long_time_over_setbacks, 
               names_to = "item", 
               values_to = "response")

Finally, now that we've got it all in long-form, we can join on scoring.

Write the code to create an object named dat that joins dat_long to scoring by their common columns.

dat <- inner_join(x = object1, y = object2, by = c("commoncol1", "commoncol2"))

dat <- inner_join(x = dat_long, y = scoring, by = c("item", "response"))

11.9 Activity 8: Summarise and visualise

Now we've got everything together we can calculate the mean scores for each participant and plot some initial findings. First, we'll create an object that calculates the mean score for each participant.

Because we also want this dataset to include the demographic variables, we have to list these in group_by().
There's some missing data (we'll come back to this in the next chapter), so you need to add on na.rm = TRUE to get it to calculate the means.
When you use group_by(), it's always safest to include ungroup() at the end at it prevents future operations you do on that dataset accidentally being grouped by variables you didn't intend.

dat_scores <- dat %>%
  group_by(participant_ID, age, gender, treatment) %>%
  summarise(resilience_score = mean(score, na.rm = TRUE)) %>%
  ungroup()

## `summarise()` has grouped output by 'participant_ID', 'age', 'gender'. You can
## override using the `.groups` argument.

Now that we've got the mean scores, we can look at group differences.For example, scores for different genders:

dat_scores %>%
  group_by(gender) %>%
  summarise(group_score = mean(resilience_score, na.rm = TRUE)) %>%
  ungroup()

gender	group_score
man	3.000816
non-binary	3.098810
woman	3.017898

Or for the treatment conditions:

dat_scores %>%
  group_by(treatment) %>%
  summarise(group_score = mean(resilience_score, na.rm = TRUE))%>%
  ungroup()

treatment	group_score
control	2.810133
intervention	3.197737

Finally, make plot the data however you want. Here's a few options you could try.

A violin-boxplot.

A faceted histogram.

A boxplot with raw data points, and the mean with error bars representing standard error.

I haven't come close to teaching you how to make this one but if you want to stretch yourself, see if you can figure it out and send me (Emily) a DM with your code. I would start with the boxplot, then try and flip it, then add in the data raw data points and the axis labels. The colours and the mean + SE will cause you the most trouble. Happy Googling :)

11.10 Finished

Finally, try knitting the file to HTML and remember to make a note of any mistakes you made and how you fixed them or any other useful information you learned. Then save your Markdown, and quit RStudio completely.

10 Big five personality 2

12 Resilience 2