Re: [R-br] Método de Jacobi
Da uma olhada neste site e verifique se os códigos te atendem. http://morotalab.org/Mrode2005/mme/mme.html#section0002 Att. Tiago. Wecsley Prates <woprates@gmail.com> escreveu: Olá a todos... Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro. Essa programação me fornece um erro do tipo Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero. Como resolvo??? Agradeço a atenção a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis for (k in 1:length(x)-1){ for (i in k+1:length(x)){ m <- a[i,k]/a[k,k] a[i,k] <- 0 for(j in k+1:length(x)){ a[i,j] <- a[i,j]-m*a[k,j] b[i] <- b[i]-m*b[k] } } } x[length(x)] <- b[length(x)]/a[length(x),length(x)] for (k in length(x)-1:1){ s <- 0 for (j in k+1:length(x)){ s <- s+a[k,j]*x[j] } x[k] <- (b[k]-s)/a[k,k] } -- * Wecsley O. Prates* *Doutorando em Estatística - Universidade Federal de Minas Gerais* _______________________________________________ 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.
Confira se - for (k in 1:length(x)-1) - for (k in length(x)-1:1) estão fazendo o que é necessário. Por exemplo, se* x <- 1:10*; o primeiro "for" retorna o vetor *0 1 2 3 4 5 6 7 8 9*, e o segundo retorna o valor *9*. Talvez o que você esteja interessado seja - for (k in 1:(length(x)-1)) - for (k in (length(x)-1):1) Em 29 de junho de 2015 09:47, Tiago Souza Marçal < tiagosouzamarcal@hotmail.com> escreveu:
Da uma olhada neste site e verifique se os códigos te atendem.
http://morotalab.org/Mrode2005/mme/mme.html#section0002
Att.
Tiago.
Wecsley Prates <woprates@gmail.com> escreveu:
Olá a todos...
Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro.
Essa programação me fornece um erro do tipo
Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero.
Como resolvo???
Agradeço a atenção
a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis
for (k in 1:length(x)-1){
for (i in k+1:length(x)){
m <- a[i,k]/a[k,k]
a[i,k] <- 0 for(j in k+1:length(x)){
a[i,j] <- a[i,j]-m*a[k,j]
b[i] <- b[i]-m*b[k] } } }
x[length(x)] <- b[length(x)]/a[length(x),length(x)]
for (k in length(x)-1:1){
s <- 0 for (j in k+1:length(x)){
s <- s+a[k,j]*x[j] }
x[k] <- (b[k]-s)/a[k,k] }
-- * Wecsley O. Prates*
*Doutorando em Estatística - Universidade Federal de Minas Gerais *
_______________________________________________ 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.
-- Rafael Garcia Cunha YNWA
Corrigindo: Por exemplo, se* x <- 1:10*; *1:lenght(x)-1* retorna o vetor *0 1 2 3 4 5 6 7 8 9*, e *length(x)-1:1* retorna o valor *9* Em 29 de junho de 2015 09:58, Rafael Garcia Cunha <rafareds@gmail.com> escreveu:
Confira se
- for (k in 1:length(x)-1) - for (k in length(x)-1:1)
estão fazendo o que é necessário. Por exemplo, se* x <- 1:10*; o primeiro "for" retorna o vetor *0 1 2 3 4 5 6 7 8 9*, e o segundo retorna o valor *9*. Talvez o que você esteja interessado seja
- for (k in 1:(length(x)-1)) - for (k in (length(x)-1):1)
Em 29 de junho de 2015 09:47, Tiago Souza Marçal < tiagosouzamarcal@hotmail.com> escreveu:
Da uma olhada neste site e verifique se os códigos te atendem.
http://morotalab.org/Mrode2005/mme/mme.html#section0002
Att.
Tiago.
Wecsley Prates <woprates@gmail.com> escreveu:
Olá a todos...
Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro.
Essa programação me fornece um erro do tipo
Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero.
Como resolvo???
Agradeço a atenção
a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis
for (k in 1:length(x)-1){
for (i in k+1:length(x)){
m <- a[i,k]/a[k,k]
a[i,k] <- 0 for(j in k+1:length(x)){
a[i,j] <- a[i,j]-m*a[k,j]
b[i] <- b[i]-m*b[k] } } }
x[length(x)] <- b[length(x)]/a[length(x),length(x)]
for (k in length(x)-1:1){
s <- 0 for (j in k+1:length(x)){
s <- s+a[k,j]*x[j] }
x[k] <- (b[k]-s)/a[k,k] }
-- * Wecsley O. Prates*
*Doutorando em Estatística - Universidade Federal de Minas Gerais *
_______________________________________________ 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.
-- Rafael Garcia Cunha
YNWA
-- Rafael Garcia Cunha YNWA
Obrigado a todos pela ajuda. Abs Em 29 de junho de 2015 10:01, Rafael Garcia Cunha <rafareds@gmail.com> escreveu:
Corrigindo: Por exemplo, se* x <- 1:10*; *1:lenght(x)-1* retorna o vetor *0 1 2 3 4 5 6 7 8 9*, e *length(x)-1:1* retorna o valor *9*
Em 29 de junho de 2015 09:58, Rafael Garcia Cunha <rafareds@gmail.com> escreveu:
Confira se
- for (k in 1:length(x)-1) - for (k in length(x)-1:1)
estão fazendo o que é necessário. Por exemplo, se* x <- 1:10*; o primeiro "for" retorna o vetor *0 1 2 3 4 5 6 7 8 9*, e o segundo retorna o valor *9*. Talvez o que você esteja interessado seja
- for (k in 1:(length(x)-1)) - for (k in (length(x)-1):1)
Em 29 de junho de 2015 09:47, Tiago Souza Marçal < tiagosouzamarcal@hotmail.com> escreveu:
Da uma olhada neste site e verifique se os códigos te atendem.
http://morotalab.org/Mrode2005/mme/mme.html#section0002
Att.
Tiago.
Wecsley Prates <woprates@gmail.com> escreveu:
Olá a todos...
Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro.
Essa programação me fornece um erro do tipo
Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero.
Como resolvo???
Agradeço a atenção
a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis
for (k in 1:length(x)-1){
for (i in k+1:length(x)){
m <- a[i,k]/a[k,k]
a[i,k] <- 0 for(j in k+1:length(x)){
a[i,j] <- a[i,j]-m*a[k,j]
b[i] <- b[i]-m*b[k] } } }
x[length(x)] <- b[length(x)]/a[length(x),length(x)]
for (k in length(x)-1:1){
s <- 0 for (j in k+1:length(x)){
s <- s+a[k,j]*x[j] }
x[k] <- (b[k]-s)/a[k,k] }
-- * Wecsley O. Prates*
*Doutorando em Estatística - Universidade Federal de Minas Gerais *
_______________________________________________ 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.
-- Rafael Garcia Cunha
YNWA
-- Rafael Garcia Cunha
YNWA
_______________________________________________ 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.
-- * Wecsley O. Prates* *Doutorando em Estatística - Universidade Federal de Minas Gerais*
participantes (3)
-
Rafael Garcia Cunha -
Tiago Souza Marçal -
Wecsley Prates