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@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 ajuda Joelmir
_______________________________________________ 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.