<div dir="ltr"><div class="gmail_quote">Prezados, mais uma vez peço desculpas... Não consigo responder diretamente ao tópico dando o reply na mensagem que chega até a minha caixa de e-mail. Será que estou fazendo algum procedimento errado??? </div><div class="gmail_quote"><br></div><div class="gmail_quote">Como preciso muito solucionar meu caso, estou reenviando uma nova mensagem com o mesmo título do tópico anteriormente criado. Se existir a possibilidade de mesclar esse tópico com os outros dois criados eu ficaria agradecido. Se alguém puder me orientar tbm em como responder uma mensagem no R-BR eu agradeço, pois já participei de alguns tópicos aqui e acredito que sempre respondia ao e-mail que chegava na minha conta...  Neste tópico vou colocar as mensagens postadas anteriormente para pelo menos seguir uma linha de raciocínio, e novamente peço para me orientarem como responder nesse tópico atual. Obrigado.</div><div class="gmail_quote"><br></div><div class="gmail_quote">###########################################</div><div class="gmail_quote"><br></div><div class="gmail_quote">A seguir peço que leiam a situação atual do problema:</div><div class="gmail_quote"><br></div><div class="gmail_quote">###########################################</div><div class="gmail_quote"><br></div><div class="gmail_quote"><div dir="ltr"><div>Caros colegas,</div><div><br></div><div>estou com um problema que não consigo resolver... Na verdade acredito que já identifiquei o problemas, mas não sei como contornar esse tipo de situação no R. É o seguinte estou maximizando a função de log-verossimilhança de uma dist. de probabilidade, com o intuito de estimar os parâmetros da dist. Entretanto a otimização não está convergindo para o parâmetro conhecido. Essa dist. possui dois parâmetros, onde um deles o "z" assume valores menores que 2, e o outro parâmetro "eta" assume valores maiores que 0. Consegui simular situações para os casos que tenho 1<z<2 e para 0<z<1, entretanto para valores negativos de z, não estou conseguindo convergência para o parâmetro. Estou usando a função optim() do R, selecionando nessa função o método de Nelder-Mead (já usei outros métodos e nenhum convergiu). O problema que ocorre é que para z negativo, durante as iterações ocorrem o cálculo de logaritmos cujos argumentos estão com valores negativos, tenho quase certeza que o problema da não convergência está aí, aí eu gostaria de saber se tem como colocar alguma restrição durante o processo de iteração para que os argumentos das funções logarítmicas nunca sejam valores negativos. A seguir coloco o código que estou usando:</div><div><br></div><div>##############################################</div><div><br></div><div>##### Fornecendo os parâmetros</div><div>##############################</div><div>n<-50</div><div>z<--1.95</div><div>eta<-0.1</div><div>z.est=vector()</div><div>eta.est=vector()</div><div><br></div><div>p=1</div><div>while (p<1001) {</div><div><br></div><div>###### Gerando valores da varíavel cujos parâmetros são z e eta</div><div><br></div><div>u<-runif(n)<span style="white-space:pre-wrap">      </span></div><div>y = eta*((1-(u)^((1-z)/(2-z)))/(1-z))</div><div><br></div><div>#### Obtendo as estimativas para z e eta</div><div>########################################</div><div><br></div><div>vero <- function(par,x){</div><div>z = par[1]</div><div>eta = par[2]</div><div>saida<-((sum(log(1 - ((1 - z)*x*(1/eta)))))/(1 - z)) + ((n)*log(2 - z)) + (n*log(1/eta))</div><div>return(-saida)</div><div>}</div><div><br></div><div>saida1<-optim(par=c(-1,4),fn=vero,</div><div>      method="Nelder-Mead",x=y</div><div>      )</div><div>z.est[p]<-saida1[1]$par[1]</div><div>eta.est[p]<-saida1[1]$par[2]</div><div>p=p+1</div><div>}</div><div>mean(z.est)</div><div>mean(eta.est)</div><div><br></div><div>##############################################</div><div><br></div><div>Como estou replicando o procedimento mil vezes, irei utilizar a média das 1000 estimativas encontradas para z e para eta... Então fazendo isso tenho a seguinte saída:</div><div><br></div><div>> mean(z.est)</div><div>[1] -30177882</div><div><br></div><div>> mean(eta.est)</div><div>[1] 964132.7</div><div><br></div><div>Conforme observamos as estimativas estão muito distantes dos parâmetros reais definidos no início do código (z = -1.95 e eta = 0.1).</div><div><br></div><div><br></div><div>Dessa forma, desconfio que o problema seja decorrente dos vários erros de logaritmos, que se repetem durante a interação e são da seguinte forma:</div><div><br></div><div>1: In log(1 - ((1 - z) * x * (1/eta))) : NaNs produzidos</div><div><br></div><div><br></div><div>Então tem alguma forma de controlar esse problema que ocorre durante a interação??? Tenho que fazer de alguma forma com que os argumentos fiquem sempre positivos, ou seja: (1-((1-z)*x*(1/eta))) tem que ser sempre positivo...</div><div><br></div><div>Por favor, caso vcs possam me ajudar, eu agradeço imensamente.</div><div><br></div><div>Att,</div><div><br></div><div>Romero.</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>###########################################</div><div>###########################################</div><div>Mensagem anterior enviada por <span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">Wagner bonat:</span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">**************************************</span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div><div style="font-family:arial,sans-serif;font-size:12.7272720336914px"><div>Uma outra opcao é vc restringir o espaco de busca da optim() usando o metodo L-BFGS-B. Vc precisa identificar para quais valores do seu parametro o erro acontece e se certificar que o algoritmo nao caminhe nesta direcao.<br><br></div>Uma segunda opcao é se for possivel calcule pelo menos a primeira derivada da funcao de log-verossimilhanca e a forneca ao algoritmo BFGS ou L-BFGS-B isso ajuda o algoritmo ficar sempre no espaco parametrico, porque a derivada nao deixa ele tentar valores meio que absurdos.<br><br></div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px">Para uma implementacao definitiva eu aconselharia usar um Fisher-score muito eficiente e normalmente partindo de valores iniciais razoaveis vc sempre converge. 'E claro que para isso vc precisa calcular a matriz de derivadas segundas e tomar a esperanca, dependendo da situacao pode ser razoavel.</span><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div><div><div>###########################################</div><div>###########################################</div><div>Mensagem anterior enviada por <span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">Romero (rlmsf)</span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">**************************************</span></div></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div><span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">Olá Walmes e Vinícius,</span><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">obrigado pelas dicas, pelo menos já estou buscando encontrar alguma forma de resolver, pois confesso, que já não estava vendo mais como sair do dilema.  </div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">Vinicius, já tinha tentado fazer a otimização da função de verossimilhança (sem o log), entretanto, não sei o que acontece, mas as estimativas ficam todas iguais aos pontos iniciais, ou seja, para as mil replicações que estou fazendo, todos os resultados são iguais aos pontos iniciais... Com certeza, fazendo dessa forma, não deve estar chegando nem na segunda iteração (não sei o motivo).</div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">Walmes, a sua dica é boa, alguém já tinha falado em algo do tipo, mas eu não tinha entendido bem... acho que agora com sua explicação ficou mais claro, resumindo, por favor veja se é isso: preciso encontrar funções cujos possíveis resultado sejam pertencentes ao domínio dos meus dois parâmetros de interesse, com a exigência de que o argumento dessas funções sejam sempre valores pertencentes aos reais, certo??</div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">Pensando assim, como o parâmetro z é sempre menor que 2 e o parâmetro eta é sempre um número positivo, então teoricamente posso usar a seguinte transformação: </div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">z= 2-exp(x) (argumento x real e resultado da função sempre menor que 2)   e    eta=exp (k) (argumento k real e resultado da função sempre maior que zero). É isso mesmo, vcs concordam???  </div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">Obrigado pela atenção,</div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px">Romero.</div></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><div style="color:rgb(34,34,34);font-family:arial;font-size:small">###########################################</div><div style="color:rgb(34,34,34);font-family:arial;font-size:small">###########################################</div><div style="color:rgb(34,34,34);font-family:arial;font-size:small">Mensagem anterior enviada por <span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">Vinicius Brito</span></div><div style="color:rgb(34,34,34);font-family:arial;font-size:small"><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">**************************************</span></div></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><br></div><div style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:12.7272720336914px"><span style="color:rgb(34,34,34);font-size:12.7272720336914px">Romero, </span><div style="color:rgb(34,34,34);font-size:12.7272720336914px"><br></div><div style="color:rgb(34,34,34);font-size:12.7272720336914px">faz tempo que não tenho esse tipo de problema na minha vida.</div><div style="color:rgb(34,34,34);font-size:12.7272720336914px"><br></div><div style="color:rgb(34,34,34);font-size:12.7272720336914px">A sugestão do Walmes me parece a mais correta.</div><div style="color:rgb(34,34,34);font-size:12.7272720336914px"><br></div><div style="color:rgb(34,34,34);font-size:12.7272720336914px">Outra abordagem que me parece ser razoável e tentar maximizar sem a função de log, pois deveria obter o mesmo resultado.</div><div style="color:rgb(34,34,34);font-size:12.7272720336914px"><br></div><div style="color:rgb(34,34,34);font-size:12.7272720336914px">Ou também usar uma abordagem bayesiana. Determinar uma distribuição a posteriori e então determinar as distribuições condicionais completas para os seus parâmetros. Sem dúvida esse tipo de approach vai te dar muito mais trabalho.</div><div style="color:rgb(34,34,34);font-size:12.7272720336914px"><br></div><div style="color:rgb(34,34,34);font-size:12.7272720336914px">Abs</div><div style="color:rgb(34,34,34);font-size:12.7272720336914px"><br></div><div style="color:rgb(34,34,34);font-size:12.7272720336914px">Vinicius</div></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div><div>###########################################</div><div>###########################################</div><div>Mensagem anterior enviada por <span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">Walmes</span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">**************************************</span></div></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap"><br></span></div><div></div>
</div>
</div><span style="color:rgb(0,0,0);font-family:'trebuchet ms',sans-serif;font-size:13px">Em problemas de otimização o primeiro conselho é ter uma função objetivo com argumentos no domínio dos reais. Tente colocar reparametrizações na sua função objetivo para garantir isso. Por exemplo, no caso da distribuição normal a função de log-verrossimilhança tem que ser escrita com mu=mu e sigma=exp(theta), não se estima diretamente sigma pois sigma>0 mas exp(theta)>0 para todos theta nos reais. A mesma idéia serve para binomial, onde p está no (0,1), não se otimiza em p mas em theta sendo que p=1/(1+exp(theta)) que vem da função de ligação log(p/(1-p))=theta.</span><br style="color:rgb(0,0,0);font-family:'trebuchet ms',sans-serif;font-size:13px"><br style="color:rgb(0,0,0);font-family:'trebuchet ms',sans-serif;font-size:13px"><span style="color:rgb(0,0,0);font-family:'trebuchet ms',sans-serif;font-size:13px">À disposição.</span><br style="color:rgb(0,0,0);font-family:'trebuchet ms',sans-serif;font-size:13px"><span style="color:rgb(0,0,0);font-family:'trebuchet ms',sans-serif;font-size:13px">Walmes.</span><div><font color="#000000" face="trebuchet ms, sans-serif"><br></font></div><div><font color="#000000" face="trebuchet ms, sans-serif"><br></font></div><div><div>###################################################</div><div>###################################################</div><div>Mensagem anterior enviada por <span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">Romero (rlmsf) (1º Tópico)</span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px;white-space:nowrap">**************************************</span><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:12.7272720336914px;white-space:nowrap">*******</span></div><div><span style="color:rgb(102,102,102);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:12.7272720336914px;white-space:nowrap"><br></span></div><div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Caros colegas,</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">estou com um problema que não consigo resolver... Na verdade acredito que já identifiquei o problemas, mas não sei como contornar esse tipo de situação no R. É o seguinte estou maximizando a função de log-verossimilhança de uma dist. de probabilidade, com o intuito de estimar os parâmetros da dist. Entretanto a otimização não está convergindo para o parâmetro conhecido. Essa dist. possui dois parâmetros, onde um deles o "z" assume valores menores que 2, e o outro parâmetro "eta" assume valores maiores que 0. Consegui simular situações para os casos que tenho 1<z<2 e para 0<z<1, entretanto para valores negativos de z, não estou conseguindo convergência para o parâmetro. Estou usando a função optim() do R, selecionando nessa função o método de Nelder-Mead (já usei outros métodos e nenhum convergiu). O problema que ocorre é que para z negativo, durante as iterações ocorrem o cálculo de logaritmos cujos argumentos estão com valores negativos, tenho quase certeza que o problema da não convergência está aí, aí eu gostaria de saber se tem como colocar alguma restrição durante o processo de iteração para que os argumentos das funções logarítmicas nunca sejam valores negativos. A seguir coloco o código que estou usando:</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">##############################################</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">##### Fornecendo os parâmetros</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">##############################</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">n<-50</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">z<--1.95</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">eta<-0.1</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">z.est=vector()</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">eta.est=vector()</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">p=1</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">while (p<1001) {</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">###### Gerando valores da varíavel cujos parâmetros são z e eta</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">u<-runif(n)<span class="" style="white-space:pre">        </span></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">y = eta*((1-(u)^((1-z)/(2-z)))/(1-z))</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">#### Obtendo as estimativas para z e eta</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">########################################</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">vero <- function(par,x){</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">z = par[1]</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">eta = par[2]</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">saida<-((sum(log(1 - ((1 - z)*x*(1/eta)))))/(1 - z)) + ((n)*log(2 - z)) + (n*log(1/eta))</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">return(-saida)</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">}</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">saida1<-optim(par=c(-1,4),fn=vero,</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">      method="Nelder-Mead",x=y</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">      )</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">z.est[p]<-saida1[1]$par[1]</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">eta.est[p]<-saida1[1]$par[2]</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">p=p+1</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">}</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">mean(z.est)</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">mean(eta.est)</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">##############################################</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Como estou replicando o procedimento mil vezes, irei utilizar a média das 1000 estimativas encontradas para z e para eta... Então fazendo isso tenho a seguinte saída:</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">> mean(z.est)</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><span id="skype_c2c_container" class="" dir="ltr" tabindex="-1"><span class="" dir="ltr"><span class="" id="non_free_num_ui"><img width="0" height="0" class="" src="chrome-extension://lifbcibllhkdhoafpjfnlhfpfgnpldfl/call_skype_logo.png"><span class="">[1] -30177882</span><span class=""></span></span></span></span></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">> mean(eta.est)</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">[1] 964132.7</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Conforme observamos as estimativas estão muito distantes dos parâmetros reais definidos no início do código (z = -1.95 e eta = 0.1).</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Dessa forma, desconfio que o problema seja decorrente dos vários erros de logaritmos, que se repetem durante a interação e são da seguinte forma:</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">1: In log(1 - ((1 - z) * x * (1/eta))) : NaNs produzidos</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Então tem alguma forma de controlar esse problema que ocorre durante a interação??? Tenho que fazer de alguma forma com que os argumentos fiquem sempre positivos, ou seja: (1-((1-z)*x*(1/eta))) tem que ser sempre positivo...</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Por favor, caso vcs possam me ajudar, eu agradeço imensamente.</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Att,</div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px"><br></div><div style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13px">Romero.</div></div><div class="gmail_signature"><br></div>
</div></div>