drunk <- function(n) {
  t <- 0:n
  p <- dbinom(t, n, 1/4)
  x <- 2*t-n
  plot(x, p, xlim=c(-60,60), ylim=c(0,1))
  points(x[p==max(p)], p[p==max(p)], pch="X")
  Sys.sleep(0.2)
  x[p==max(p)][1]
}
mle <- sapply(0:60, drunk)
plot(0:60, mle)