Eu tenho encontrado alguns problemas estranhos ao usar este código:

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)
}

Cada iteração deste loop for deve calcular a probabilidade binomial que,i O número de ocorrências do evento ocorre de acordo com a frequência. Cada iteração também resume o resultado. Isso deve resultar naprob variável nunca superior a 1, mas depois de 7 ou mais para iterações de loop, tudo vai para o inferno eprob Excede 1.

Eu pensei que poderia ser uma questão de dígitos de precisão, então eu tentei usarRmpfr mas sem sucesso, o mesmo problema persistiu.

Eu queria saber se há alguma dica ou pacote para superar essa situação, ou se estou preso a isso.