Bei der Verwendung dieses Codes sind einige seltsame Probleme aufgetreten:

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

Jede Iteration dieser for-Schleife sollte die Binomialwahrscheinlichkeit berechnen, dassi Häufigkeit des Auftretens des Ereignisses. Jede Iteration fasst auch das Ergebnis zusammen. Dies sollte zu dem Ergebnis führenprob Variable nie größer als 1, aber nach 7 oder so für Schleifeniterationen geht alles zur Hölle undprob übersteigt 1.

Ich dachte, es könnte eine Frage der Genauigkeit von Ziffern sein, also versuchte ich es mitRmpfr aber ohne Erfolg blieb das gleiche Problem bestehen.

Ich habe mich gefragt, ob es irgendwelche Tipps oder Pakete gibt, um diese Situation zu überwinden, oder ob ich damit nicht weiterkomme.