download.file("http://www.openintro.org/stat/data/kobe.RData", destfile = "kobe.RData")
load("kobe.RData")
head(kobe)
## vs game quarter time
## 1 ORL 1 1 9:47
## 2 ORL 1 1 9:07
## 3 ORL 1 1 8:11
## 4 ORL 1 1 7:41
## 5 ORL 1 1 7:03
## 6 ORL 1 1 6:01
## description basket
## 1 Kobe Bryant makes 4-foot two point shot H
## 2 Kobe Bryant misses jumper M
## 3 Kobe Bryant misses 7-foot jumper M
## 4 Kobe Bryant makes 16-foot jumper (Derek Fisher assists) H
## 5 Kobe Bryant makes driving layup H
## 6 Kobe Bryant misses jumper M
kobe$basket[1:9]
## [1] "H" "M" "M" "H" "H" "M" "M" "M" "M"
# Exercise 1
kobe_streak <- calc_streak(kobe$basket)
barplot(table(kobe_streak))
# Exercise 2
outcomes <- c("heads", "tails")
sample(outcomes, size = 1, replace = TRUE)
## [1] "heads"
sim_fair_coin <- sample(outcomes, size = 100, replace = TRUE)
sim_fair_coin
## [1] "heads" "tails" "heads" "heads" "heads" "tails" "heads" "tails"
## [9] "tails" "tails" "tails" "tails" "tails" "tails" "heads" "tails"
## [17] "tails" "tails" "tails" "tails" "heads" "tails" "heads" "heads"
## [25] "tails" "heads" "tails" "heads" "tails" "heads" "heads" "tails"
## [33] "tails" "tails" "heads" "heads" "tails" "tails" "tails" "tails"
## [41] "heads" "heads" "tails" "heads" "tails" "tails" "tails" "tails"
## [49] "tails" "heads" "tails" "heads" "heads" "heads" "heads" "tails"
## [57] "heads" "tails" "tails" "tails" "heads" "heads" "tails" "heads"
## [65] "tails" "tails" "tails" "tails" "tails" "heads" "tails" "heads"
## [73] "tails" "heads" "heads" "heads" "heads" "tails" "tails" "heads"
## [81] "heads" "heads" "heads" "heads" "tails" "tails" "tails" "tails"
## [89] "heads" "heads" "tails" "tails" "heads" "heads" "heads" "heads"
## [97] "heads" "tails" "heads" "heads"
table(sim_fair_coin)
## sim_fair_coin
## heads tails
## 47 53
sim_unfair_coin <- sample(outcomes, size = 100, replace = TRUE, prob = c(0.2, 0.8))
# Exercise 3
outcomes <- c("H", "M")
sim_basket <- sample(outcomes, size = 1, replace = TRUE)
# Exercise 4
kobe$basket
## [1] "H" "M" "M" "H" "H" "M" "M" "M" "M" "H" "H" "H" "M" "H" "H" "M" "M"
## [18] "H" "H" "H" "M" "M" "H" "M" "H" "H" "H" "M" "M" "M" "M" "M" "M" "H"
## [35] "M" "H" "M" "M" "H" "H" "H" "H" "M" "H" "M" "M" "H" "M" "M" "H" "M"
## [52] "M" "H" "M" "H" "H" "M" "M" "H" "M" "H" "H" "M" "H" "M" "M" "M" "H"
## [69] "M" "M" "M" "M" "H" "M" "H" "M" "M" "H" "M" "M" "H" "H" "M" "M" "M"
## [86] "M" "H" "H" "H" "M" "M" "H" "M" "M" "H" "M" "H" "H" "M" "H" "M" "M"
## [103] "H" "M" "M" "M" "H" "M" "H" "H" "H" "M" "H" "H" "H" "M" "H" "M" "H"
## [120] "M" "M" "M" "M" "M" "M" "H" "M" "H" "M" "M" "M" "M" "H"
sim_basket
## [1] "M"
# On Your Own
# Simulating the independent shooter
outcomes <- c("H", "M")
sim_basket <- sample(outcomes, size = 133, replace = TRUE, prob = c(0.45, 0.55))
table(sim_basket)
## sim_basket
## H M
## 49 84
sim_streak <- calc_streak(sim_basket)
barplot(table(sim_streak))
max(sim_streak)
## [1] 4
# 1. The distribution is skewed right. Typical streak length is 0. Longest streak is 6.
# 2. I would expect the similar distribution, since the probability is the same. The shape of the distribution
#should still be right skewed.
# Kobe's distribution would be very similar because the probability is the same. I would not conclude that the
#hot hand model fits Kobe's shooting pattern. It appears to me that shooting is more independent than dependent.