创建一个函数以在R中生成负二项式随机变量

Question

这是作业，所以rnbinom无济于事。

我创建了两个函数：combr（重复组合）和另一个函数negbinom。

combr <- function(k,r) {
  return(factorial(k+r-1) / (factorial(r) * factorial (k-1)))
  }

negbinom <- function(k,r,p){
  return(combr(k,r)*(p^k)*((1-p)^r))
  }

Negbinom会生成一个号码。如何通过调用n negbinom次来创建函数以生成n个随机数。或者类似地，我如何创建自己的“rnbinom”？