# File "Rdimag": a simple Monte Carlo program to estimate probabilities
#   related to Joe DiMaggio's hitting streaks, for the book "Knock on Wood"
#   by Jeffrey S. Rosenthal.

games = 1736
prob = 0.7919133
nrep = 10^5
longsum = longest = num56 = 0

for (count in 1:nrep) {
  thelist = rbinom(games, 1, prob)
  cur = long = 0
  for (i in 1:games) {
    if (thelist[i]==1) {
      cur = cur+1
      if (cur > long)
        long = cur
    } else {
      cur = 0
    }
  }
  longsum = longsum + long
  if (long > longest)
    longest = long
  if (long >= 56)
    num56 = num56 + 1
}
print(longsum/nrep)
print(longest)
print(num56/nrep)