<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span>Obrigada pela ajuda. Tem que instalar algum pacote diferente? Vou tentar fazer. att, Fátima</span></div><div><br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div dir="ltr"> <hr size="1"> <font size="2" face="Arial"> <b><span style="font-weight:bold;">De:</span></b> "andrebvs@bol.com.br" <andrebvs@bol.com.br><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quinta-feira, 10 de Outubro de 2013 19:23<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Função-escalonar<br> </font> </div> <div class="y_msg_container"><br><div id="yiv1484456136"><div><div>Com essa função
chega-se ao resultado final da matriz escalonada, porém, não apresenta as matrizes intermediárias até chegar ao resultado final.<br clear="none"><br clear="none">A = matrix(c(-2,4,6,0,0,0,-1,1,1,2,0,2),nr=3)<br clear="none"><br clear="none">reducedRowEchelonForm <- function(A){<br clear="none">n<-nrow(A)<br clear="none">m<-ncol(A)<br clear="none">i<-j<-1<br clear="none">while(i<=n&&j<=m){<br clear="none">while(j<=m){<br clear="none">currentColumn <- A[,j]<br clear="none">currentColumn[1:n < i] <- 0<br clear="none"># find maximum pivot in current column at or below current row<br clear="none">which <- which.max(abs(currentColumn))<br clear="none">pivot <- currentColumn[which]<br clear="none">if(pivot==0){#checkfor0pivot<br clear="none">j<-j+1<br clear="none">next<br clear="none">5<br clear="none">}<br clear="none">if (which > i) A[c(i, which),] <- A[c(which, i),] # exchange rows<br
clear="none">A[i,] <- A[i,]/pivot # pivot<br clear="none">row<-A[i,]<br clear="none">A<-A-outer(A[,j],row) #sweep<br clear="none">A[i,] <- row # restore current row<br clear="none">j<-j+1<br clear="none">break<br clear="none">}<br clear="none">i<-i+1<br clear="none">}<br clear="none">#0rowstobottom<br clear="none">zeros <- which(ap
ply(A[,1:m], 1, function(x) all(x == 0)))<br clear="none">if (length(zeros) > 0){<br clear="none">zeroRows <- A[zeros,]<br clear="none">A <- A[-zeros,]<br clear="none">A <- rbind(A, zeroRows)<br clear="none">rownames(A) <- NULL<br clear="none">}<br clear="none">A<br clear="none">}<br clear="none"><br clear="none">reducedRowEchelonForm(A)<br clear="none"><br clear="none">Att.<br clear="none">André</div>
<div> </div>
<div><br clear="none"><br clear="none"></div>
<hr style="border-top:1px solid #ccc;">
<div class="yiv1484456136yqt5220233419" id="yiv1484456136yqt89743"><div><br clear="none"><b>De:</b> Fernando A. Souza < nandodesouza@gmail.com ><br clear="none"><b>Enviada:</b> Quinta-feira, 10 de Outubro de 2013 12:32<br clear="none"><b>Para:</b> r-br@listas.c3sl.ufpr.br<br clear="none"><b>Assunto:</b> Re: [R-br] Função-escalonar<br clear="none"><br clear="none"></div>
<div class="yiv1484456136moz-cite-prefix">Em 10-10-2013 12:26, <a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-abbreviated" href="">andrebvs@bol.com.br</a> escreveu:</div>
<blockquote>
<div>O escalonameno da função matinv() mostra as operações realizadas até chegar no resulado final?<br clear="none"> <br clear="none"> <i>Att.</i><br clear="none"> <i>André</i></div>
<div> </div>
<hr style="border-top:1px solid #ccc;">
<div><br clear="none"> <b>De:</b> Benilton Carvalho < <a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-abbreviated" href="">beniltoncarvalho@gmail.com</a> ><br clear="none"> <b>Enviada:</b> Quinta-feira, 10 de Outubro de 2013 10:09<br clear="none"> <b>Para:</b> r-br < <a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-abbreviated" href="">r-br@listas.c3sl.ufpr.br</a> ><br clear="none"> <b>Assunto:</b> Re: [R-br] Função-escalonar<br clear="none"> </div>
<div dir="ltr">
<div>
<div>scale() nao e' o q vc procura, Fatima.<br clear="none"> </div>
Ha' algum tempo, qdo o pacote Design ainda era suportado, havia a funcao matinv() que fazia escalonamento por Gauss-Jordan... mas, atualmente, o pacote nao e' suportado... Se vc tiver uma versao antiga do R, vc ainda pode tentar instala'-lo.</div>
b</div>
<div class="yiv1484456136gmail_extra"><br clear="none"> <br clear="none">
<div class="yiv1484456136gmail_quote">Em 10 de outubro de 2013 07:54, Fernando A. Souza <span dir="ltr"><<a rel="nofollow" shape="rect" href="">nandodesouza@gmail.com</a>></span> escreveu:<br clear="none">
<blockquote class="yiv1484456136gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>
<div>Em 10-10-2013 07:12, Maria de Fatima Ferreira Almeida escreveu:</div>
<blockquote>
<div class="yiv1484456136im">
<div style="font-size: 12pt; font-family: arial, helvetica, sans-serif;">
<div>Gostaria de saber se existe alguma função do R que faz escalonamento de matrizes. Alguém poderia me ajudar?</div>
</div>
<br clear="none"><fieldset></fieldset></div>
<pre>_______________________________________________
R-br mailing list
<a rel="nofollow" shape="rect" href="">R-br@listas.c3sl.ufpr.br</a>
<a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
Veja a função scale() do R</div>
<br clear="none"> _______________________________________________<br clear="none"> R-br mailing list<br clear="none"> <a rel="nofollow" shape="rect" href="">R-br@listas.c3sl.ufpr.br</a><br clear="none"> <a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none"> Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</blockquote>
</div>
</div>
<div> </div>
<br clear="none"><fieldset class="yiv1484456136mimeAttachmentHeader"></fieldset><br clear="none">
<pre>_______________________________________________
R-br mailing list
<a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-abbreviated" href="">R-br@listas.c3sl.ufpr.br</a>
<a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-freetext" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-freetext" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<div>Ah sim ! Benilton, você tem razão,.Scale não faz o que ela deseja não! No endereço :<a rel="nofollow" shape="rect" class="yiv1484456136moz-txt-link-freetext" target="_blank" href="http://socserv.socsci.mcmaster.ca/jfox/Courses/R-course/exercises-programming.pdf">http://socserv.socsci.mcmaster.ca/jfox/Courses/R-course/exercises-programming.pdf</a>, há uma função que faz o escalonamento de uma matriz pelo escalonamento gaussiano desenvolvida pelo Jonh Fox. Ela deve resolver seu problema<br clear="none"> </div></div>
</div></div><br><div class="yqt5220233419" id="yqt44194">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div><br><br></div> </div> </div> </div></body></html>