binomial.py

Created by numworks

Created on April 11, 2018

397 Bytes

Ce script contient deux fonctions pour le calcul des probabilités binomiales : binomialpdf(k,N,p) renvoie p(X=k) et binomialcdf(k,N,P) renvoie p(X≤k) pour N répétitions et une probabilité de succès p.

from math import *
def factorial(n):
  n=int(n)
  if n<0:
    return None
  res=1
  for i in range(1,n+1):
    res=res*i
  return res
def n_choose_k(n,k):
  return int(factorial(n)/(factorial(k)*factorial(n-k)))
def binomialpdf(k,N,p):
  return n_choose_k(N,k)*(p**k)*((1-p)**(N-k))
def binomialcdf(k,N,p):
  res=0
  for i in range(k+1):
    res=res+binomialpdf(i,N,p)
  return res