<div dir="ltr">Confira se <div><ul><li><span style="font-size:12.8000001907349px">for (k in 1:length(x)-1)</span><br></li><li><span style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">for (k in length(x)-1:1)</span></span></li></ul><div><span style="font-size:12.8000001907349px">estão fazendo o que é necessário. Por exemplo, se<b> x <- 1:10</b>; o primeiro "for" retorna o vetor <b>0 1 2 3 4 5 6 7 8 9</b>, e o segundo retorna o valor <b>9</b>. Talvez o que você esteja interessado seja</span></div></div><div><div><ul><li><span style="font-size:12.8000001907349px">for (k in 1:(length(x)-1))</span></li><li><span style="font-size:12.8000001907349px">for (k in (length(x)-1):1)</span></li></ul></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 29 de junho de 2015 09:47, Tiago Souza Marçal <span dir="ltr"><<a href="mailto:tiagosouzamarcal@hotmail.com" target="_blank">tiagosouzamarcal@hotmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div>
<div>Da uma olhada neste site e verifique se os códigos te atendem.<br>
<br>
<a href="http://morotalab.org/Mrode2005/mme/mme.html#section0002" target="_blank">http://morotalab.org/Mrode2005/mme/mme.html#section0002</a><br>
<br>
Att.<br>
<br>
Tiago.<br>
<br>
Wecsley Prates <<a href="mailto:woprates@gmail.com" target="_blank">woprates@gmail.com</a>> escreveu:<br>
<br>
</div><div><div class="h5">
<div>
<div dir="ltr">Olá a todos...
<div><br>
</div>
<div>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.</div>
<div><br>
</div>
<div>Essa programação me fornece um erro do tipo</div>
<div><br>
</div>
<div>
<div>Erro em a[i, j] <- a[i, j] - m * a[k, j] : </div>
<div>  substituto tem comprimento zero.</div>
</div>
<div><br>
</div>
<div>Como resolvo???</div>
<div><br>
</div>
<div>Agradeço a atenção</div>
<div><br>
</div>
<div>
<div>a = matriz dos coeficientes</div>
<div>b <- vetor das constantes</div>
<div>x <- chute inicial do vetor das variáveis</div>
<div><br>
</div>
<div>for (k in 1:length(x)-1){</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>for (i in k+1:length(x)){</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>m <- a[i,k]/a[k,k]</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>a[i,k] <- 0</div>
<div><span style="white-space:pre-wrap"></span>for(j in k+1:length(x)){</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>a[i,j] <- a[i,j]-m*a[k,j]</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>b[i] <- b[i]-m*b[k]</div>
<div><span style="white-space:pre-wrap"></span>}</div>
<div><span style="white-space:pre-wrap"></span>}</div>
<div>}</div>
<div><br>
</div>
<div>x[length(x)] <- b[length(x)]/a[length(x),length(x)]</div>
<div><br>
</div>
<div>for (k in length(x)-1:1){</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>s <- 0</div>
<div><span style="white-space:pre-wrap"></span>for (j in k+1:length(x)){</div>
<div><br>
</div>
<div><span style="white-space:pre-wrap"></span>s <- s+a[k,j]*x[j]</div>
<div><span style="white-space:pre-wrap"></span>}</div>
<div><span style="white-space:pre-wrap"><br>
</span></div>
<div><span style="white-space:pre-wrap"></span>x[k] <- (b[k]-s)/a[k,k]</div>
<div>}</div>
<div><br>
</div>
<div><br>
</div>
-- <br>
<div>
<div style="text-align:left"><b><i><font face="verdana, sans-serif"><font color="#000099"> </font><font color="#006600">                      <u><font size="4">Wecsley O. Prates</font></u></font></font></i></b></div>
<div style="text-align:left"><b><i><font face="verdana, sans-serif"><font color="#006600"><u><font size="4"><br>
</font></u></font></font></i></b></div>
<div style="text-align:left"><b><i><font color="#006600" face="verdana, sans-serif"><u>Doutorando em Estatística - Universidade Federal de Minas Gerais<br>
</u></font></i></b><span style="color:rgb(0,102,0);font-family:verdana,sans-serif"> </span></div>
<div><font color="#006600"><b><i></i></b> </font></div>
</div>
</div>
</div>
</div>
</div></div></div>

<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div>Rafael Garcia Cunha</div><div><br></div><div>YNWA</div><div><br></div></div>
</div>