如何在RStudio中通过此操作计算大数？

Question

我正在执行此操作，并且当过程给我NaN时，

    Choques<-c(1:10)
    print (Choques)
    pr<-0
    n<-3818
    p<-0.040633627

    for(i in Choques) {
    pr[i]<-(factorial(n)/(factorial(Choques[i])*factorial(n-Choques[i])))*p^Choques[i]*(1-p)^(n-Choques[i]) 
    print (pr[i])
    }

但是将变量n的值更改为较小的数字，例如20，如果它显示数字，我想知道是否有任何方法可以通过数字来更改NaN，我认为数字太大。

Answer 1

替换

factorial(n)/(factorial(Choques[i])*factorial(n-Choques[i]))

with

choose(n, Choques[i])