Obese-discrimination-arrests

Obese-discrimination-arrests.R

snf <- read.csv("http://www.jakeporway.com/teaching/data/snf_2.csv", as.is=TRUE)

# First create a subset of the data only consisting of “good” weights and heights.
clean.subset = snf[snf$height > 40 & snf$weight > 90 & snf$weight < 400, ]
# Add a BMI variable to our dataset, where BMI is computed as (weight)*703/(height*height).
bmi = (clean.subset$weight)*703/(clean.subset$height*clean.subset$height)
# What percentage of people with BMI’s greater than or equal to 30 who were stopped were 
# ultimately arrested?
obese.subset = clean.subset[bmi >= 30, ]
obese.arrested.subset = which(obese.subset$arrested == 1)
length(which(obese.subset$arrested == 1)) / nrow(obese.subset)
# Answer: 
[1] 0.0777903