Playing with permutations

• Slasher films are notoriously gory and are said to contain recurring biases.
  • One such bias, is that women in slasher films are more likely to survive than men.
• This lab will focus on the statistical question: Are women in slasher films more likely to survive until the end of the film than men?
• To answer this question, we’ll learn how to use permuted data to gauge how likely an event occurs by chance.
• To begin, use the data function to load the slasher data file.
  • The data contains information about 485 characters from a random sample of 50 slasher horror films.

Initial thoughts…

• To familiarize yourself with the data, answer the following:
  • How many variables and observations are contained in the data and what are the possible values of the variables?
  • Which gender had more survivors? Write down a few sentences as to how you came to your conclusion. Be sure to look at both the counts and percentages of survivors in each group before deciding.
  • Calculate the difference between the percentage of females who survived and the percentage of males who survived. Is the difference large enough to conclude that women tend to survive more often than men?

Tally whoa … !

• Something you might have noticed is that these two lines of code aren’t equivalent:

tally(gender ~ survival, data = slasher, margin = TRUE)

tally(survival ~ gender, data = slasher, margin = TRUE)

• The first line of code takes the group of survivors and tells us how many of them were Male or Female.
• The other takes the group of females / males and tells us how many of them Dies or Survives.
• The last question on the previous slide can be answered using the line of code below. Why?
  • • Pro-tip: Include the option format = "percent" to obtain a two-way table with percentages.

tally(survival ~ gender, format = "percent", data = slasher, margin = TRUE)

##           gender
## survival      Female      Male
##   Dies      77.47748  86.69202
##   Survives  22.52252  13.30798
##   Total    100.00000 100.00000

Examining differences

• When we’re comparing the difference between two quantities, such as survival rates of slasher films, it can be difficult to decide how different two values need to be before we can conclude that the difference didn’t just happen by chance.
  • To help us decide when a difference is not due to chance, we’ll use repeated random shuffling.
• By using repeated random shuffling, we’ll estimate how often our actual difference occurs by chance.

Do the shuffle!

• When we shuffle data, we use our original data set as a starting point.
  • Run the following and write down the resulting table on a piece of paper.

tally(survival ~ gender, data = slasher)

• Now run the following to randomly reassign each survival status to each observation. Compare the resulting table to the one you wrote down.

tally(shuffle(survival) ~ gender, 
      data = slasher)

Let’s compare …

• How many people, in total, survived the slasher film before shuffling? How many people survived after shuffling?
• How has shuffling our data changed the percentage of women who survived compared to men who survived?
  • Is the difference in percentages from your shuffled data larger or smaller than the difference from the original data? Interpret what this means.
• Explain why shuffling our data one time is not enough to decide if the difference seen in our actual data occurs by chance or not.

Detecting differences

• To help us decide if the difference in percentages in our actual data occurs by chance or not, we can use the do() function to shuffle our data many times and see how often our actual difference occurred by chance.

Run the following lines of code:

set.seed(7)

shuffled_outcomes <- do(10) * tally(shuffle(survival) ~ gender, format = "percent", 
                                    data = slasher)

View(shuffled_outcomes)

• In how many simulations did a higher percentage of males survive than females?

• What is the largest difference in percentages of survival between males and females?

• What patterns are emerging from these simulations?

• Ten simulations is not enough. Use do, tally and shuffle functions to shuffle the survival variable and tally the percentage of women who survived 500 times. Assign your 500 shuffles the name shuffled_survivors. Use set.seed(1)

Now what?

• The next step to find out how often our actual difference occur by chance is to compare it to the differences in our shuffled data.
• To compute the differences for each shuffle we can use the mutate function.
  • Fill in the blanks to add a new column that contains the difference between Survives.Female and Survives.Male to our shuffled_survivors data.

shuffled_survivors <- mutate(shuffled_survivors, 
            diff = ____ - ____)

Time to decide

• Create a histogram of the differences in our shuffled_survivors data. Based on your plot, answer the following:

  • What was the typical difference in percentages between men and women survivors?

• Include a vertical line in your histogram of the actual difference by running the code below:

add_line(vline = 22.52252 - 13.30798)

• Does the actual difference occur very often by chance alone?
• Does gender play a role in whether or not a character will survive in a horror film? Explain your reasoning.
• If you wanted to survive in a horror film, would you want to play a female character or a male character?

Summary

• By shuffling the survival label, we made it so that the proportion of males and females who survived the slasher film was random.
  • The males and females survived by chance alone.
• If surviving the film occurred purely by chance, then most of the time the difference in survival proportions was close to zero.
  • Notice how most values in the histogram occur close to zero.
• When we look to see how often our actual difference occurs in our shuffled data, if the actual difference doesn’t occur very often then perhaps there is something more going on than just chance alone …

On your own

• Carry out another 500 simulations but this time shuffle the gender variable instead of the survival variable.
  • Include the code set.seed(1) before your 500 simulations to make your answer reproducible.
• Does shuffling the gender variable instead of the survival variable change your answer to the question Does gender play a role in whether or not a character will survive in a horror film?
  • Why or why not?