In R you can just plug in a vector of probabilities directly to get the likelihood. Ie: p = seq(.01, .99, by = .01) y1 = dbinom(7, 47, p) y2 = dnbinom(

In R you can just plug in a vector of probabilities directly to get the likelihood. Ie: p = seq(.01, .99, by = .01) y1 = dbinom(7, 47, p) y2 = dnbinom(40, 7, p) It doesn't matter much for this use case but that will also be much faster. It does simplify the code quite a bit too though.