Obrigado a todos pela ajuda. Abs
Em 29 de junho de 2015 10:01, Rafael Garcia Cunha <[email protected]> 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 <[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. > -- * 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.
