<div dir="ltr">Oi Éder,<div><br></div><div>me desculpe pela demora em responder! </div><div>Muito obrigada pela ajuda! Deu certo!</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 27 de março de 2015 00:16, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Boa noite,</div><div><br></div><div>Acho que é possível simplificar com o uso de scale()...</div><div><br></div><div><span style="font-family:monospace,monospace">### <code r></span></div><div><span style="font-family:monospace,monospace">x <- as.matrix(c(10.373173554, 6.125265914, 3.148949951, 2.722401273, 2.045041946, </span><br></div><div><font face="monospace, monospace">1.498078581, 1.35012278, 1.33545711, 1.190051106, 0.999610112, </font></div><div><font face="monospace, monospace">0.850212256, 0.716720583, 0.600704889, 0.474433269, 0.431826856, </font></div><div><font face="monospace, monospace">0.388576986, 0.266709936, 0.176960029, 0.045935755, 0.02012691, </font></div><div><font face="monospace, monospace">-0.003897135, -0.030136936, -0.074815367, -0.14541063, -0.294553778, </font></div><div><font face="monospace, monospace">-0.453450244, -0.551343912, -0.710794973, -1.222537293, -1.273419528))</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">### Similar ao utilizado no código inicial</font></div><div><font face="monospace, monospace">raiz <- sqrt(sum(x^2))</font></div><div><font face="monospace, monospace">scale(x, center=FALSE, scale=raiz)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">###  Para valores entre 0 e 1</font></div><div><font face="monospace, monospace">scale(x-min(x), center=FALSE, scale=max(x)-min(x))</font></div><div><font face="monospace, monospace">### </code></font></div><div class="gmail_extra"><br clear="all"><div><div><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div></div>
<br><div class="gmail_quote"><div><div class="h5">Em 26 de março de 2015 21:01, Michelle Bau Graczyk <span dir="ltr"><<a href="mailto:mbgraczyk@gmail.com" target="_blank">mbgraczyk@gmail.com</a>></span> escreveu:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Boa noite,<div><br></div><div>Eu estou normalizando autovetores de uma matriz e colocando eles em outra. Para isso eu montei o programa a seguir. A primeira forma, ele só calcula o ultimo elemento da coluna 4 da matriz, e da segunda forma os valores do resultado estão errados.</div><div>Primeira forma:</div><div><div>> a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE)</div><div><br></div><div>> b<-as.matrix(a)</div><div>> matriznorm<-matrix(0,30,4)</div><div>> matriznorm[,1]<-b[,1]# só vou usar este valores do arquivo acima</div><div>> show(matriznorm[,1])</div><div> [1] 10.373173554  6.125265914  3.148949951  2.722401273  2.045041946  1.498078581  1.350122780</div><div> [8]  1.335457110  1.190051106  0.999610112  0.850212256  0.716720583  0.600704889  0.474433269</div><div>[15]  0.431826856  0.388576986  0.266709936  0.176960029  0.045935755  0.020126910 -0.003897135</div><div>[22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973</div><div>[29] -1.222537293 -1.273419528</div><div>> for( i in 1:30){</div><div>+   matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1]</div><div>+ }</div><div>>   </div><div>> somaQuadrados<-sum(matriznorm[,2])</div><div>> raiz<-sqrt(somaQuadrados)</div><div>> matriznorm[,3]<-raiz</div><div>> </div><div>>  for(j in i:30){</div><div>+    matriznorm[j,4]<-matriznorm[j,1]/raiz</div><div>+  }</div><div>> </div><div>> write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt")</div><div>> show(matriznorm)</div><div>              [,1]         [,2]     [,3]        [,4]</div><div> [1,] 10.373173554 1.076027e+02 13.46876  0.00000000</div><div> [2,]  6.125265914 3.751888e+01 13.46876  0.00000000</div><div> [3,]  3.148949951 9.915886e+00 13.46876  0.00000000</div><div> [4,]  2.722401273 7.411469e+00 13.46876  0.00000000</div><div> [5,]  2.045041946 4.182197e+00 13.46876  0.00000000</div><div> [6,]  1.498078581 2.244239e+00 13.46876  0.00000000</div><div> [7,]  1.350122780 1.822832e+00 13.46876  0.00000000</div><div> [8,]  1.335457110 1.783446e+00 13.46876  0.00000000</div><div> [9,]  1.190051106 1.416222e+00 13.46876  0.00000000</div><div>[10,]  0.999610112 9.992204e-01 13.46876  0.00000000</div><div>[11,]  0.850212256 7.228609e-01 13.46876  0.00000000</div><div>[12,]  0.716720583 5.136884e-01 13.46876  0.00000000</div><div>[13,]  0.600704889 3.608464e-01 13.46876  0.00000000</div><div>[14,]  0.474433269 2.250869e-01 13.46876  0.00000000</div><div>[15,]  0.431826856 1.864744e-01 13.46876  0.00000000</div><div>[16,]  0.388576986 1.509921e-01 13.46876  0.00000000</div><div>[17,]  0.266709936 7.113419e-02 13.46876  0.00000000</div><div>[18,]  0.176960029 3.131485e-02 13.46876  0.00000000</div><div>[19,]  0.045935755 2.110094e-03 13.46876  0.00000000</div><div>[20,]  0.020126910 4.050925e-04 13.46876  0.00000000</div><div>[21,] -0.003897135 1.518766e-05 13.46876  0.00000000</div><div>[22,] -0.030136936 9.082349e-04 13.46876  0.00000000</div><div>[23,] -0.074815367 5.597339e-03 13.46876  0.00000000</div><div>[24,] -0.145410630 2.114425e-02 13.46876  0.00000000</div><div>[25,] -0.294553778 8.676193e-02 13.46876  0.00000000</div><div>[26,] -0.453450244 2.056171e-01 13.46876  0.00000000</div><div>[27,] -0.551343912 3.039801e-01 13.46876  0.00000000</div><div>[28,] -0.710794973 5.052295e-01 13.46876  0.00000000</div><div>[29,] -1.222537293 1.494597e+00 13.46876  0.00000000</div><div>[30,] -1.273419528 1.621597e+00 13.46876 -0.09454616</div><div>> </div></div><div>Da segunda forma:</div><div><div>> a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE)</div><div>> b<-as.matrix(a)</div><div>> matriznorm<-matrix(0,30,4)</div><div>> matriznorm[,1]<-b[,1]</div><div>> show(matriznorm[,1])</div><div> [1] 10.373173554  6.125265914  3.148949951  2.722401273  2.045041946  1.498078581  1.350122780</div><div> [8]  1.335457110  1.190051106  0.999610112  0.850212256  0.716720583  0.600704889  0.474433269</div><div>[15]  0.431826856  0.388576986  0.266709936  0.176960029  0.045935755  0.020126910 -0.003897135</div><div>[22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973</div><div>[29] -1.222537293 -1.273419528</div><div>> for( i in 1:30){</div><div>+   matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1]</div><div>+   </div><div>+ somaQuadrados<-sum(matriznorm[,2])</div><div>+ raiz<-sqrt(somaQuadrados)</div><div>+ matriznorm[,3]<-raiz</div><div>+ </div><div>+  for(j in i:30){</div><div>+    matriznorm[j,4]<-matriznorm[j,1]/raiz</div><div>+  }</div><div>+ }</div><div>> write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt")</div><div>> show(matriznorm)</div><div>              [,1]         [,2]     [,3]          [,4]</div><div> [1,] 10.373173554 1.076027e+02 13.46876  1.0000000000</div><div> [2,]  6.125265914 3.751888e+01 13.46876  0.<a href="tel:5084624670" value="+15084624670" target="_blank">5084624670</a></div><div> [3,]  3.148949951 9.915886e+00 13.46876  0.<a href="tel:2528991576" value="+12528991576" target="_blank">2528991576</a></div><div> [4,]  2.722401273 7.411469e+00 13.46876  0.<a href="tel:2135962692" value="+12135962692" target="_blank">2135962692</a></div><div> [5,]  2.045041946 4.182197e+00 13.46876  0.1584251428</div><div> [6,]  1.498078581 2.244239e+00 13.46876  0.1152793173</div><div> [7,]  1.350122780 1.822832e+00 13.46876  0.1033376904</div><div> [8,]  1.335457110 1.783446e+00 13.46876  0.1016853664</div><div> [9,]  1.190051106 1.416222e+00 13.46876  0.0902440218</div><div>[10,]  0.999610112 9.992204e-01 13.46876  0.0755856426</div><div>[11,]  0.850212256 7.228609e-01 13.46876  0.0641564605</div><div>[12,]  0.716720583 5.136884e-01 13.46876  0.0540043425</div><div>[13,]  0.600704889 3.608464e-01 13.46876  0.0452163558</div><div>[14,]  0.474433269 2.250869e-01 13.46876  0.0356888678</div><div>[15,]  0.431826856 1.864744e-01 13.46876  0.0324667092</div><div>[16,]  0.388576986 1.509921e-01 13.46876  0.0292025272</div><div>[17,]  0.266709936 7.113419e-02 13.46876  0.0200398900</div><div>[18,]  0.176960029 3.131485e-02 13.46876  0.0132951406</div><div>[19,]  0.045935755 2.110094e-03 13.46876  0.0034511674</div><div>[20,]  0.020126910 4.050925e-04 13.46876  0.0015121392</div><div>[21,] -0.003897135 1.518766e-05 13.46876 -0.0002927926</div><div>[22,] -0.030136936 9.082349e-04 13.46876 -0.0022641887</div><div>[23,] -0.074815367 5.597339e-03 13.46876 -0.0056207917</div><div>[24,] -0.145410630 2.114425e-02 13.46876 -0.0109238800</div><div>[25,] -0.294553778 8.676193e-02 13.46876 -0.0221227473</div><div>[26,] -0.453450244 2.056171e-01 13.46876 -0.0340370871</div><div>[27,] -0.551343912 3.039801e-01 13.46876 -0.0413498303</div><div>[28,] -0.710794973 5.052295e-01 13.46876 -0.0532327969</div><div>[29,] -1.222537293 1.494597e+00 13.46876 -0.0911767942</div><div>[30,] -1.273419528 1.621597e+00 13.46876 -0.0945461625</div></div><div><br></div><div>Muito obrigada!</div></div>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a 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>
Leia o guia de postagem (<a 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.<br></blockquote></div><br></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" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>