<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span>Adailton,</span></div><div><br><span></span></div><div><span>Procure seguir o guia de postagem e evite anexar arquivos! <br></span></div><div><br><span></span></div><div><span>Repassando então a dúvido de Adailton Alcantara para os nobres amigos do R. </span></div><div><br><span></span></div><div><span>OBS: Segue o link no final da postagem com as funções que o Adailton criou!<br></span></div><div><span><br></span></div><div><span>Prezados,<br><br>Nome: Adeilton Pedro de Alcântara<br><br>Venho pedir ajuda.<br><br>Definição do problema:<br><br>Eu tenho uma superfície (mistura de cinco distribui��es normal bivariada) lambda_23(x,y) e uma reta lambda_1(t). Considerando a reta e a superf�cie como conhcidas ( distribui��es verdadeiras) o meu objetivo � recuperar essas distribui��es maximuzando a
equa��o de log verossimilhan�a ( no programa que enviei denominada matricaalg1). <br><br>Ou seja o meu objetivo � recuperar a estrutura do polin�mio que gerou um determinado conjunto de dados utilizando uma equa��o de log-verossimilhan�a sem a necessidade de escrever um modelo param�trico.<br><br>A minha dificuldade est� no comando nlminb() do R que encontra-se na fun��o metricaalg1. <br><br>Eu n�o consigo lidar com estes comandos de maximiza��o/minimiza��o de fun��es. Quando altero o chute inicial, embora tenha converg�ncia, os valores encontrados s�o completamente diferentes do valor que eu quero encontrar. Ou seja dependendo do chute inicial a fun��o me retorna coeficientes completamente diferentes.<br><br>Eu n�o sei se o problema est� na op��o control.list() ou no chute inicial. Provavelmente n�o.<br><br>O programa que gera o PPNH( Processo de Poisson n�o homog�neo) est�
funcionando. <br><br>A dificuldade est� no uso do comando nlminb().<br>J� utilizei nlm(), optim(), e DEoptim() para maximizar a equa��o de log-verossimilhan�a por�m n�o consigo enxergar o que est� errado .<br><br>O problema deve ser de m�ltiplos m�ximos na verossimilhan�a, pois tenho muitos coeficientes.<br><br><span style="font-weight: bold; background-color: rgb(255, 0, 0);">Já estou a quatro anos e não consigo sair disso.</span><br><br>Atenciosamente,<br><br>Adeilton Pedro de Alcântara.<br><br>Estou enviando um c�digo na linguagem R. O objetivO � recuperar uma reta (lambda_1(t) )e uma mistura de 5 distribui��es normal bivariada (lambda_23(x,y) ).<br><br>A minha dificuldade est� no chute inicial, ou seja os resultados dependem muito do chute inicial. Al�m disso, cada vez que eu modifici o chute inicial tenho um resultado completamente diferente.<br><br>O senhor verifica a fun��o onde eu gero o chute
inicial (CHUTEINICIALalg1) e a fun��o metricaalg1 (fun��o onde tenho a verossimilhan�a). Nesta fun��o estou fazendo a maximiza��o num�rica da equa��o de log-verossimilhan�a utilizaddo o comando "nlminb() do R.<br><br>Onde tem CC<-(as.single(CC)>0) # Coloquei positivo para n�o ficar muito tempo rodando. O senhor pode deixar CC<-(as.single(CC)).<br><br>Para rodar a fun��o basta seguir os comandos a partir de <br><br>#################<br>#Rodando a fun��o#<br>#################<br><br>O arquivo com as rotinas do R est� no anexo deste e-mail<br><br><br>Me ajude, por favor.<br><br>Atenciosamente,<br><br>Adeilton Pedro de Alcântara</span></div><div><br><span></span></div><div><span>http://dl.dropbox.com/u/33619290/codigoAdailton.txt</span></div><div> </div><div style="background-color:transparent;" align="center"><font style="background-color:transparent;" face="comic sans ms" size="2"><font
size="1"><b><br></b></font></font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\begin{signature}</font></font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2"><<>>=</font></font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Prof. Dr. Ivan Bezerra Allaman</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Universidade Estadual de Santa Cruz</font></div><div style="text-align:left;background-color:transparent;" align="center"><font
class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Departamento de Ciências Exatas e Tecnológicas</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Ilhéus/BA - Brasil</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Fone: +55 73 3680-5596</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com</font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">@</font></font></div><div
style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\end{signature}</font></font></div></div></body></html>