<html><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:10pt"><div><span class="Apple-style-span" style="font-family: Verdana, Geneva, Helvetica, Arial, sans-serif; ">Prezados, estou tentando maximizar uma função de verossimilhança para estimar 5 parâmetros. Contudo, meus resultados estão dando NULL. <br><br>Sou iniciante na manipulação do R, por isso posso ter cometido erros primários. Favor alguém me ajude. <br><br>Meus dados são: <br><br>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) <br><br>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) <br><br>A distribuição é de Poisson: <br><br>B <- rpois(50, 1) <br>S <- rpois(50, 1) <br><br>A equação que utilizo é: <br><br>loglik <- function(param) { <br>alpha
<- param[1] <br>mu <- param[2] <br>Eb <- param[3] <br>Es <- param[4] <br>delta <- param[5] <br>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)) <br>} <br>res <- maxLik(loglik, grad=NULL, hess=NULL, start=c(0,1), constraints=NULL) <br>print(res) <br>coef(res) <br><br>Acho que estou errando nos argumentos dessa equação. <br><br>Alguém tem alguma sugestão?</span><br></div><div style="text-align: left; font-family: arial, helvetica, sans-serif; " align="right"><font color="#000080" face="Geneva, Arial, Sans-serif" size="1"><i></i></font></div><div><font style="font-family: arial, helvetica, sans-serif; " size="1"><font style="color:rgb(0, 0, 0);" size="2"><span
style="font-weight:bold;"></span></font></font></div></div></body></html>