[R-br] Método de Newton Rhapson...

Sexta Fevereiro 8 21:01:54 BRST 2013

O método de Gaus Newton é um método de busca numérica muito utilizado para encontrar os parâmetros de equações não lineares. 
Qual função deu origem as derivadas parciais da figura do site?
Você tem os valores de wi?
Vamos tentar criar uma situação para tentarmos resolver o seu problema.
Como eu só utilizei gaus newton para regressão não linear preciso entender um pouco mais o seu problema.
Att.
Tiago. 

Date: Fri, 8 Feb 2013 14:41:48 -0200
From: andrebvs em bol.com.br
To: r-br em 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 em 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 em bol.com.br
To: r-br em 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 em 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 em 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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130209/b15508df/attachment.html>