<div dir="ltr">Olá amigos,<div>  Estou criando uma rotina (ainda em teste) para usar nos meus dados originais...Será que poderiam me ajudar pq desse erro:</div><div><br></div><div><div>teste<-c(0,2,0,1,0,1,3,0,0,0,4)</div><div>x <-c(0,2,1,5,3,1,0,1,5,7,2)</div><div><br></div><div>info1<-c()</div><div>info2<-c()</div><div>dzip1= function(y,x,p,beta){</div><div>      for (i in 1:length(y))</div><div>{<span class="" style="white-space:pre">               </span>beta0 = beta[1] </div><div>            beta1 = beta[2] </div><div>            lambda=exp(beta0+beta1*x[i])</div><div>            if (y[i]==0) {info1[i]<-p+((1-p)*exp(-lambda))}</div><div>            if (y[i]>0)  {info2[i]<-((1-p)*exp(-lambda)*(lambda^y[i]))/factorial(y[i])}</div><div>}</div><div>return(na.exclude(c(info1,info2)))</div><div>}</div><div><br></div><div>dzip1(teste,x,0.3,c(0.5,2.5))</div><div><br></div><div><br></div><div>MV1 <- function(p,beta) </div><div>{    D = dzip1(teste,x,p,beta)</div><div>     -sum(log(D))</div><div>}</div><div><br></div><div>library(stats4) </div><div>mle(MV1, start = list(p=0.5, beta=c(0.5,1)),method="L-BFGS-B")</div><div><br></div><div><br></div><div><div><font color="#ff0000">> mle(MV1, start=list(p=0.5, beta=c(0.5,1)),method="L-BFGS-B")</font></div><div><font color="#0000ff">Erro em optim(start, f, method = method, hessian = TRUE, ...) : </font></div><div><font color="#0000ff">  objeto (list) não é coercível para tipo 'double'</font></div></div><div><br></div><div>O mesmo programa funcionou mto bem, sem usar a função de ligação.</div><div>O problema começou qdo coloquei os betas.</div><div>Como não conheço mto a função mle, talvez vcs possam me ajudar....</div><div><br></div><div><br></div><div>Abçs,</div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><br></div><div><br></div><div>------------------------------------------------------------------------------------------</div><font size="2">Vanessa Felix do Nascimento Sergio<br>Graduada em Estatística - ICE/UFJF<br>Mestranda em Saúde Coletiva - IMS/UERJ</font><div><font size="2">Assistente de pesquisa I - Instituto de Pesquisa Econômica Aplicada – IPEA<br>(32) 9142-8428<br>(21) 97557-3956<br>-------------------------------------------------------------------------------------------<br></font><br></div></div></div></div></div>
</div></div>