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?