Prezados, estou tentando maximizar uma função de verossimilhança para estimar 5 parâmetros. Contudo, meus resultados estão dando NULL.  Sou iniciante na manipulação do R, por isso posso ter cometido erros primários. Favor alguém me ajude.  Meus dados são:  B <-c(1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,1,0,1,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,1,0,1,1,1,1,1,0,0,0,0,0,0,0,0,0)  S <- c(0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0,1,0,1,1,1,0,1,1,1,0,1,1,1,1,1,1,1,0,1,0,0,0,0,0,1,1,1,1,1,1,1,1,1)  A distribuição é de Poisson:  B <- rpois(50, 1)  S <- rpois(50, 1)  A equação que utilizo é:  loglik <- function(param) {  alpha <- param[1]  mu <- param[2]  Eb <- param[3]  Es <- param[4]  delta <- param[5]  ll <- (1-alpha)*exp(-Eb)*((Eb^B)/factorial(B))*exp(-Es)*((Es^S)/factorial(S)) + alpha*delta*exp(-Eb)*((Eb^B)/factorial(B))*exp(-(mu+Es))*((mu+Es)^S/factorial(S)) + alpha*(1-delta)*exp(-(mu+Eb))*((mu+Eb)^B/factorial(B))*exp(-Es)*((Es^S)/factorial(S))  }  res <- maxLik(loglik, grad=NULL, hess=NULL, start=c(0,1), constraints=NULL)  print(res)  coef(res)  Acho que estou errando nos argumentos dessa equação.  Alguém tem alguma sugestão?