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: ivanala...@yahoo.com.br/ivanala...@gmail.com
@
\end{signature}_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
