Date: Fri, 8 Feb 2013 14:41:48 -0200
From: andrebvs@bol.com.br
To: r-br@listas.c3sl.ufpr.br
Subject: Re: [R-br] Método de Newton Rhapson...
Obg Tiago!
O método de Gaus Newton é equivalente ao Newton Rhapson?
Para meu caso, eu tenho 3 parâmetros, como adaptar seu exemplo ao meu caso pelo Gaus Newton? Não estou sabendo montar as expressões no R para chegar aos resultados numéricos.
Segue no link abaixo, a função densidade de probabilidade da distribuição de valores extremos (GVE), bem como suas derivadas parcais em relação a cada parâmetro.
https://www.transferbigfiles.com/f1220ccb-8cb3-46c6-8a68-3c95fc24ebc4?rid=CnKlj7VjIyqbgcMXXXqFfg2
desde já agradeço!
Att.
André
Em 08/02/2013 07:01, Tiago Souza Marçal < tiagosouzamarcal@hotmail.com > escreveu:
André se você quiser a resolução especificamente por este método da uma olhado no dicas ridículas que o Walmes postou algo a respeito.
Mas eu estou trabalhando com o método de GAUSS-NEWTON e vou te passar o script para caso você queira utiliza-lo.
No site abaixo tem um material explicando o método.
www.inf.ufsc.br/~ogliari/arquivos/regressao_nao_linear.ppt
No exemplo a seguir objetivo era encontrar os betas de um modelo exponencial.
dia<-c(2,5,7,10,14,19,26,31,34,38,45,52,53,60,65) #variável independente
diag<-c(54,50,45,37,35,25,20,16,18,13,8,11,8,4,6) #variável dependente
d3
bi<-c(55,-0.02) # Chutes iniciais eles devem ser bem feitos para garantir que haja convergência para um minimo global.
sqresi<-crossprod(diag-c(d3(dia,bi[1],bi[2]))) # O erro e a diferença entre o valores observados e a estimativa dos valores observados pelos betas dos chutes iniciais.
i
while(i < 10){
est
fx
d
sqresf
bf
bi
x
sqresi
cat(paste(formatC(c(sqresf, bi), digits=6, format="f"), collapse="\t"), "\n")
i
}
Espero ter contribuído.
Att.
Tiago.
Date: Thu, 7 Feb 2013 21:59:17 -0200
From: andrebvs@bol.com.br
To: r-br@listas.c3sl.ufpr.br
Subject: [R-br] Método de Newton Rhapson...
Olá colegas!
Gostaria de saber, como posso encontrar a solução númerica no R de cada parâmetro (µ, σ, ξ) da expressão da imagem abaixo (no link), através do método Newton Rhapson.
https://www.transferbigfiles.com/e7ec41f6-f49b-4463-a031-85159d12bc83?rid=NjTECXd92TZ3IxdEyePzfw2
desde já agradeço!
Att.
André
_______________________________________________ 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.
_______________________________________________
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.