[R-br] Iteração com optim

[Wagner Bonat]

Talvez isso ajude !

ll <- function(mu,x){-sum(dnorm(x,mean=mu, sd=1,log=TRUE))}
ll(mu = 0, x = rnorm(10))
dados <- matrix(rnorm(1000), ncol=100, nrow=10)
saida <- list()
saida[1][[1]] <- optim(par= 0, f = ll, x = dados[,1], method="Brent", upper
= 10, lower = -10)
saida[1][[1]]$par
for(i in 2:100){
    print(saida[c(i-1)][[1]]$par)
    saida[i][[1]] <- optim(par= saida[c(i-1)][[1]]$par, f = ll, x =
dados[,i], method="Brent", upper = 10, lower = -10)
    }



Em 21 de março de 2014 15:23, jborssoi . <joelmirb em gmail.com> escreveu:

> Olá Wagner,
> é isso... tenho 100 conjuntos de dados. Eu preciso maximizar o vetor de
> parâmetros para a primeira amostra. Já para a segunda (e assim por diante)
> o optim deve utilizar o vetor obtido para a primeira amostra (e assim por
> diante).
> Imagino algo do tipo:
> for(i in 1:J){
> maxL.rep <- optim(par=???, fn, gr, method="BFGS",
> control=list(fnscale=-1,maxit=1000),y=W1$y[,i], u=W1$u[,i], X=d[,c(1,2)],
> Z=Zi, sigma.e=1)
> }
>
> mas não sei como fazer para que o optim atualize "par".
>
> abraço
>
>
> Em 21 de março de 2014 11:16, Wagner Bonat <wbonat em gmail.com> escreveu:
>
> O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados
>> diferentes e quer estimar para cada um determinado modelo ... Nao consegui
>> entender ... Que metodo vc quer usar para inferencia ??
>>
>>
>> Em 21 de março de 2014 14:55, Diogo Ferrari <diogoferrari em gmail.com>escreveu:
>>
>> Nao sei se compreendi bem sua questao, mas acredito q optim faz
>>> exatamente isso automaticamente.
>>> On Mar 21, 2014 9:26 AM, "jborssoi ." <joelmirb em gmail.com> wrote:
>>>
>>>>  Caros,
>>>> tenho 100 amostras e preciso fazer a estimação de parâmetros usando
>>>> a função optim. Como faço para que, a cada nova estimação, o vetor com o
>>>> chute inicial ("par") seja atualizado com as estimativas obtidas no passo
>>>> anterior?
>>>>
>>>> Grato pela ajuda
>>>> Joelmir
>>>>
>>>> _______________________________________________
>>>> R-br mailing list
>>>> R-br em listas.c3sl.ufpr.br
>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>>> código mínimo reproduzível.
>>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>
>>
>>
>> --
>> Wagner Hugo Bonat
>> LEG - Laboratório de Estatística e Geoinformação
>> UFPR - Universidade Federal do Paraná
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>



-- 
Wagner Hugo Bonat
LEG - Laboratório de Estatística e Geoinformação
UFPR - Universidade Federal do Paraná
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/e0e422aa/attachment.html>