Я столкнулся с некоторыми странными проблемами при использовании этого кода:

positions<-c(58256)
occurrencies<-c(30)
frequency<-c(11/5531777)
length<-c(4)

prob<-c(0)
for(i in 0:(occurrencies-1))
{
  pow<-frequency^i
  pow1<-(1-frequency)^(positions-i)
  bin<-choose(positions, i)
  prob<<-prob+(bin*pow*pow1)
}

Каждая итерация цикла for должна вычислять биномиальную вероятность того, чтоi количество случаев возникновения события с учетом частоты. Каждая итерация также подводит итог. Это должно привести кprob переменная никогда не превышает 1, но после 7 или около того для итераций цикла все идет в ад иprob превышает 1.

Я подумал, что это может быть вопрос точности цифр, поэтому я попытался использоватьRmpfr но безрезультатно - та же проблема сохраняется.

Мне было интересно, есть ли какие-либо советы или пакеты для преодоления этой ситуации, или я застрял с этим.