[R-br] Dificuldade em maximização/minimização de funções! (Repassando help Adailton Alcantara)

[Ivan Bezerra Allaman]

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 em yahoo.com.br/ivanalaman em gmail.com
@
\end{signature}
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120525/b1270b51/attachment.html>