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 <[email protected]> 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 < > [email protected]> 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 <[email protected]> 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 >> [email protected] >> 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 [email protected] 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.
