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çoEm 21 de março de 2014 11:16, Wagner Bonat <wbonat@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@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@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 ajudaJoelmir
R-br mailing list
R-br@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@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@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@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.