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
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.
Faça com um laço for por exemplo chute <- ... for(...){ op <- optim(chute, ...) chute <- op$par } À disposição. Walmes.
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á
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.
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@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@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.
_______________________________________________ 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á
participantes (4)
-
Diogo Ferrari -
jborssoi . -
Wagner Bonat -
walmes .