Page 199 #15

Solution

# function to print out the value of n, p, j, and the possion value
bin_poiss <- function(n,p,j){
  #binomial function
  bin <- dbinom(j, n, p) 
  lamda <- p*n
  poiss_value <- dpois(j, lamda)
  print(n)
  print(p)
  print(j)
  print(poiss_value)
}

#Result

# n = 100, p = .03, j = 0 to 20 where 
for(j in seq(0,20)){
  bin_poiss(100, .03, j)
}
## [1] 100
## [1] 0.03
## [1] 0
## [1] 0.04978707
## [1] 100
## [1] 0.03
## [1] 1
## [1] 0.1493612
## [1] 100
## [1] 0.03
## [1] 2
## [1] 0.2240418
## [1] 100
## [1] 0.03
## [1] 3
## [1] 0.2240418
## [1] 100
## [1] 0.03
## [1] 4
## [1] 0.1680314
## [1] 100
## [1] 0.03
## [1] 5
## [1] 0.1008188
## [1] 100
## [1] 0.03
## [1] 6
## [1] 0.05040941
## [1] 100
## [1] 0.03
## [1] 7
## [1] 0.02160403
## [1] 100
## [1] 0.03
## [1] 8
## [1] 0.008101512
## [1] 100
## [1] 0.03
## [1] 9
## [1] 0.002700504
## [1] 100
## [1] 0.03
## [1] 10
## [1] 0.0008101512
## [1] 100
## [1] 0.03
## [1] 11
## [1] 0.0002209503
## [1] 100
## [1] 0.03
## [1] 12
## [1] 5.523758e-05
## [1] 100
## [1] 0.03
## [1] 13
## [1] 1.274713e-05
## [1] 100
## [1] 0.03
## [1] 14
## [1] 2.731529e-06
## [1] 100
## [1] 0.03
## [1] 15
## [1] 5.463057e-07
## [1] 100
## [1] 0.03
## [1] 16
## [1] 1.024323e-07
## [1] 100
## [1] 0.03
## [1] 17
## [1] 1.807629e-08
## [1] 100
## [1] 0.03
## [1] 18
## [1] 3.012715e-09
## [1] 100
## [1] 0.03
## [1] 19
## [1] 4.756919e-10
## [1] 100
## [1] 0.03
## [1] 20
## [1] 7.135379e-11