<div dir="ltr">Talvez isso ajude !<br><div><br>ll <- function(mu,x){-sum(dnorm(x,mean=mu, sd=1,log=TRUE))}<br>ll(mu = 0, x = rnorm(10))<br>dados <- matrix(rnorm(1000), ncol=100, nrow=10)<br>saida <- list()<br>saida[1][[1]] <- optim(par= 0, f = ll, x = dados[,1], method="Brent", upper = 10, lower = -10)<br>
saida[1][[1]]$par<br>for(i in 2:100){<br>    print(saida[c(i-1)][[1]]$par)<br>    saida[i][[1]] <- optim(par= saida[c(i-1)][[1]]$par, f = ll, x = dados[,i], method="Brent", upper = 10, lower = -10)<br>    }<br>
<br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 21 de março de 2014 15:23, jborssoi . <span dir="ltr"><<a href="mailto:joelmirb@gmail.com" target="_blank">joelmirb@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Olá Wagner,</div><div>é 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).</div>

<div>Imagino algo do tipo:</div><div>for(i in 1:J){<br>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)</div><div>}</div>

<div> </div><div>mas não sei como fazer para que o optim atualize "par".</div><div> </div><div>abraço</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 21 de março de 2014 11:16, Wagner Bonat <span dir="ltr"><<a href="mailto:wbonat@gmail.com" target="_blank">wbonat@gmail.com</a>></span> escreveu:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">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 ??<br>

</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">Em 21 de março de 2014 14:55, Diogo Ferrari <span dir="ltr"><<a href="mailto:diogoferrari@gmail.com" target="_blank">diogoferrari@gmail.com</a>></span> escreveu:<div>

<div><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><p dir="ltr">Nao sei se compreendi bem sua questao, mas acredito q optim faz exatamente isso automaticamente.</p>


<div class="gmail_quote"><div><div>On Mar 21, 2014 9:26 AM, "jborssoi ." <<a href="mailto:joelmirb@gmail.com" target="_blank">joelmirb@gmail.com</a>> wrote:<br type="attribution"></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">


<div><div>
<div dir="ltr"><div>Caros,</div><div>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?</div>




<div> </div><div>Grato pela ajuda</div><div>Joelmir</div></div>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div></div></div><span><font color="#888888"><br>

<br clear="all"><br>-- <br>Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná
</font></span></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div></div></div><br></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><br>-- <br>Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná
</div>