Bonjour tout le monde,
j'ai un problème dans l'exécution de ce script et j'ai besoin de votre aide :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
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 : Sélectionner tout - Visualiser dans une fenêtre à part
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.