Adailton, Procure seguir o guia de postagem e evite anexar arquivos! Repassando então a dúvido de Adailton Alcantara para os nobres amigos do R. OBS: Segue o link no final da postagem com as funções que o Adailton criou! Prezados, Nome: Adeilton Pedro de Alcântara Venho pedir ajuda. Definição do problema: 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). 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. A minha dificuldade est� no comando nlminb() do R que encontra-se na fun��o metricaalg1. 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. Eu n�o sei se o problema est� na op��o control.list() ou no chute inicial. Provavelmente n�o. O programa que gera o PPNH( Processo de Poisson n�o homog�neo) est� funcionando. A dificuldade est� no uso do comando nlminb(). 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 . O problema deve ser de m�ltiplos m�ximos na verossimilhan�a, pois tenho muitos coeficientes. Já estou a quatro anos e não consigo sair disso. Atenciosamente, Adeilton Pedro de Alcântara. 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) ). 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. 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. Onde tem CC<-(as.single(CC)>0) # Coloquei positivo para n�o ficar muito tempo rodando. O senhor pode deixar CC<-(as.single(CC)). Para rodar a fun��o basta seguir os comandos a partir de ################# #Rodando a fun��o# ################# O arquivo com as rotinas do R est� no anexo deste e-mail Me ajude, por favor. Atenciosamente, Adeilton Pedro de Alcântara http://dl.dropbox.com/u/33619290/codigoAdailton.txt \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}