[R-br] Erro EMV

[Vanessa Felix]

Olá amigos,
  Estou criando uma rotina (ainda em teste) para usar nos meus dados
originais...Será que poderiam me ajudar pq desse erro:

teste<-c(0,2,0,1,0,1,3,0,0,0,4)
x <-c(0,2,1,5,3,1,0,1,5,7,2)

info1<-c()
info2<-c()
dzip1= function(y,x,p,beta){
      for (i in 1:length(y))
{ beta0 = beta[1]
            beta1 = beta[2]
            lambda=exp(beta0+beta1*x[i])
            if (y[i]==0) {info1[i]<-p+((1-p)*exp(-lambda))}
            if (y[i]>0)
 {info2[i]<-((1-p)*exp(-lambda)*(lambda^y[i]))/factorial(y[i])}
}
return(na.exclude(c(info1,info2)))
}

dzip1(teste,x,0.3,c(0.5,2.5))


MV1 <- function(p,beta)
{    D = dzip1(teste,x,p,beta)
     -sum(log(D))
}

library(stats4)
mle(MV1, start = list(p=0.5, beta=c(0.5,1)),method="L-BFGS-B")


> mle(MV1, start=list(p=0.5, beta=c(0.5,1)),method="L-BFGS-B")
Erro em optim(start, f, method = method, hessian = TRUE, ...) :
  objeto (list) não é coercível para tipo 'double'

O mesmo programa funcionou mto bem, sem usar a função de ligação.
O problema começou qdo coloquei os betas.
Como não conheço mto a função mle, talvez vcs possam me ajudar....


Abçs,
-- 


------------------------------------------------------------------------------------------
Vanessa Felix do Nascimento Sergio
Graduada em Estatística - ICE/UFJF
Mestranda em Saúde Coletiva - IMS/UERJ
Assistente de pesquisa I - Instituto de Pesquisa Econômica Aplicada – IPEA
(32) 9142-8428
(21) 97557-3956
-------------------------------------------------------------------------------------------
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150721/5fed946a/attachment.html>