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 |
Partager