Lab 2E  The Horror Movie Shuffle
Lab 2E  The Horror Movie Shuffle
Directions: Follow along with the slides, completing the questions in blue on your computer, and answering the questions in red in your journal.
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 theslasher
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
orFemale
. 
The other takes the group of females / males and tells us how many of them
Dies
orSurvives
. 
The last question on the previous slide can be answered using the line of code below. Why?
 Protip: Include the option
format = "percent"
to obtain a twoway 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
 Protip: Include the option
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 the code above and perform 500 shuffles.
Assign
your 500 shuffles the same nameshuffled_outcomes
. Useset.seed(1)
.
Now what?

The next step to find out how often our actual difference occurs 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
andSurvives.Male
to ourshuffled_outcomes
data.shuffled_outcomes < mutate(shuffled_outcomes, diff = ____  ____)
Time to decide

Create a
histogram
of thedifference
s in ourshuffled_outcomes
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 thesurvival
variable.– Include the code
set.seed(1)
before your 500 simulations to make your answer reproducible. 
Does shuffling the
gender
variable instead of thesurvival
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?