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?