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

我正在执行此操作,当该过程给我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,我认为数字太大。

评论
  • Leaco
    Leaco 回复

    更换

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

    choose(n, Choques[i])