Loi binomiale, méthode du maximum de vraisemblance

Bonjour tout le monde,
j'ai un problème dans l'exécution de ce script et j'ai besoin de votre aide :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
x=c(2236, 1787, 1916, 1761, 1894,
    2135, 2103, 1659, 1468 ,2067,
     2660, 2183, 1798, 1928, 2079,
     2625, 2030, 1971, 1905 ,1957,
     2767, 2169 ,1833, 1863, 2211,
     2451, 1855 ,1577, 1993, 2052,
     2486 ,2062, 1971, 2129 ,1994,
     2514 ,1991 ,2105 ,1859, 2062,
     2339 ,1720 ,1760, 1668, 1734,
    2223, 1852, 2091, 1867, 1950,
    2546, 2080, 1782, 1874 ,1797,
  2329, 1801, 2605, 2264, 1721)
M=max(x)+round(var(x))
log_likelihood <- function(params) { -sum(dbinom(x, params[1], params[2], log=TRUE))}
fit2 <-optim(c(M,.5),log_likelihood)    
fit2

J'ai ce message d'erreur

Code:

1
2
Messages d'avis :
1: In dbinom(x, params[1], params[2], log = TRUE) : production de NaN                                         (fois 50)

Le but de ce script est d'estimer les paramètres de la loi binomiale ajustant l'échantillon par la méthode du maximum de vraisemblance.
Merci d'avance.