download.file("http://www.openintro.org/stat/data/kobe.RData", destfile = "kobe.RData")
load("kobe.RData")
head(kobe)
##    vs game quarter time                                             description
## 1 ORL    1       1 9:47                 Kobe Bryant makes 4-foot two point shot
## 2 ORL    1       1 9:07                               Kobe Bryant misses jumper
## 3 ORL    1       1 8:11                        Kobe Bryant misses 7-foot jumper
## 4 ORL    1       1 7:41 Kobe Bryant makes 16-foot jumper (Derek Fisher assists)
## 5 ORL    1       1 7:03                         Kobe Bryant makes driving layup
## 6 ORL    1       1 6:01                               Kobe Bryant misses jumper
##   basket
## 1      H
## 2      M
## 3      M
## 4      H
## 5      H
## 6      M
kobe$basket[1:9]
## [1] "H" "M" "M" "H" "H" "M" "M" "M" "M"

#Excercise 1) A streak of 1 means that there is a “H” followed by a “M”, while a streak of 0 is has 0 “H”

kobe_streak <- calc_streak(kobe$basket)
barplot(table(kobe_streak))

#Excercise 2) The bar plot is skewed to the left and unimolar. With the majority (mole) of his streaks being 0. His longest streak is 4.

outcomes <- c("heads", "tails")
sample(outcomes, size = 1, replace = TRUE)
## [1] "tails"
sample(outcomes, size = 1, replace = TRUE)
## [1] "heads"
sample(outcomes, size = 1, replace = TRUE)
## [1] "tails"
sample(outcomes, size = 1, replace = TRUE)
## [1] "tails"
sample(outcomes, size = 1, replace = TRUE)
## [1] "heads"
sample(outcomes, size = 1, replace = TRUE)
## [1] "tails"
sim_fair_coin <- sample(outcomes, size = 100, replace = TRUE)
sim_fair_coin
##   [1] "tails" "tails" "tails" "heads" "heads" "heads" "heads" "tails" "tails"
##  [10] "heads" "heads" "tails" "heads" "heads" "tails" "tails" "tails" "tails"
##  [19] "tails" "tails" "heads" "tails" "tails" "heads" "tails" "tails" "heads"
##  [28] "tails" "heads" "heads" "heads" "tails" "tails" "heads" "tails" "tails"
##  [37] "tails" "tails" "tails" "tails" "heads" "heads" "heads" "tails" "tails"
##  [46] "tails" "tails" "heads" "tails" "tails" "tails" "tails" "tails" "tails"
##  [55] "tails" "heads" "heads" "heads" "tails" "heads" "tails" "tails" "tails"
##  [64] "heads" "tails" "heads" "heads" "heads" "heads" "tails" "tails" "heads"
##  [73] "tails" "heads" "tails" "heads" "heads" "heads" "tails" "tails" "tails"
##  [82] "heads" "heads" "heads" "tails" "tails" "heads" "heads" "heads" "heads"
##  [91] "heads" "heads" "tails" "heads" "heads" "heads" "heads" "tails" "heads"
## [100] "tails"
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))
sim_unfair_coin
##   [1] "heads" "tails" "tails" "heads" "tails" "tails" "tails" "tails" "tails"
##  [10] "tails" "tails" "tails" "tails" "tails" "tails" "tails" "heads" "tails"
##  [19] "tails" "tails" "tails" "tails" "tails" "tails" "tails" "heads" "tails"
##  [28] "tails" "tails" "tails" "tails" "tails" "tails" "tails" "tails" "heads"
##  [37] "tails" "tails" "tails" "tails" "heads" "tails" "tails" "tails" "tails"
##  [46] "tails" "tails" "tails" "tails" "tails" "tails" "tails" "tails" "heads"
##  [55] "tails" "tails" "tails" "heads" "heads" "tails" "tails" "heads" "tails"
##  [64] "heads" "heads" "tails" "heads" "heads" "tails" "tails" "tails" "tails"
##  [73] "tails" "tails" "tails" "heads" "tails" "tails" "tails" "tails" "tails"
##  [82] "tails" "tails" "heads" "tails" "heads" "tails" "tails" "heads" "tails"
##  [91] "tails" "tails" "heads" "heads" "tails" "heads" "tails" "tails" "tails"
## [100] "tails"
table(sim_unfair_coin)
## sim_unfair_coin
## heads tails 
##    21    79

#Excercise 3) 22 heads

?sample
## starting httpd help server ... done
outcomes <- c("H", "M")
sim_basket <- sample(outcomes, size = 1, replace = TRUE)

#Excercise 4) sim_basket = sample(outcomes,size=133,replace=TRUE,prob=c(0.45,0.55))

outcomes=c("H","M")
sim_basket = sample(outcomes,size=133,replace=TRUE,prob=c(0.45,0.55))
sim_basket
##   [1] "M" "M" "M" "H" "M" "M" "H" "M" "H" "M" "H" "M" "H" "M" "M" "M" "M" "M"
##  [19] "M" "M" "M" "H" "M" "H" "M" "M" "M" "H" "H" "M" "M" "M" "H" "H" "M" "M"
##  [37] "H" "M" "H" "M" "H" "H" "M" "H" "H" "M" "M" "M" "M" "M" "M" "M" "M" "H"
##  [55] "H" "H" "H" "M" "H" "M" "H" "M" "H" "H" "M" "M" "M" "H" "H" "M" "H" "M"
##  [73] "H" "H" "H" "H" "M" "M" "H" "M" "M" "H" "M" "M" "M" "M" "H" "H" "H" "H"
##  [91] "M" "H" "H" "H" "H" "M" "M" "H" "M" "M" "M" "H" "H" "M" "M" "H" "M" "M"
## [109] "M" "H" "H" "H" "H" "H" "M" "M" "H" "H" "M" "H" "M" "M" "H" "H" "M" "M"
## [127] "H" "H" "M" "H" "H" "M" "H"
table(sim_basket)
## sim_basket
##  H  M 
## 61 72

#On your own #Comparing Kobe to an Independent Shooter

calc_streak(sim_basket)
##  [1] 0 0 0 1 0 1 1 1 1 0 0 0 0 0 0 0 1 1 0 0 2 0 0 2 0 1 1 2 2 0 0 0 0 0 0 0 4 1
## [39] 1 2 0 0 2 1 4 0 1 0 1 0 0 0 4 4 0 1 0 0 2 0 1 0 0 5 0 2 1 0 2 0 2 2 1
Kobe_Sim_Streak=calc_streak(sim_basket)
barplot(table(Kobe_Sim_Streak))

#1) Again, like before the streak length bar plot is unimolar with his longest streak being 4 shots.Unlike the first bar plot, it appears his streaks are a lot long in the simulation compared to the actual game.

#2) If we were to run the sample again the distribution would come out slightly different because this is probablitiy something will happen not an exact science. There is variance involved with these samples, so every sample will be slightly different and never exactly the same.

#3) Both have similar barplots, unimolar and skewed to the left. Both have similar max shooting streaks with 4. A little different but that will come with variance.