<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12px"><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span>Olá Manoel,</span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span><br></span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span id="yui_3_16_0_1_1416322549117_21663">obrigada pelo retorno. Não conhecia o GitHub obrigada pela indicação.</span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span><br></span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span id="yui_3_16_0_1_1416322549117_21664">A princípio os valores de m e dp são fixos, não são importantes nesse exemplo, e não tenho nenhum Bug pra solucionar, a função funciona perfeitamente, com a única exceção de que o valor para i=0 não é calculado corretamente.</span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span><br></span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span id="yui_3_16_0_1_1416322549117_21667">É isso mesmo que eu quero, ter os valores de f e p para então realizar a multiplicação que será retornada na matriz. Aqui um exemplo: (usei </span>m=0.1; <span style="font-size: 12px;" id="yui_3_16_0_1_1416322549117_21714">dp=sqrt(0.05))</span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"><span style="font-size: 12px;"><br></span></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style="">> conv(10)</div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style=""> [1] 0.437497149 0.437497149 0.399621300 0.323958199 0.224710500 0.129711990</div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr"></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style=""> [7] 0.061315128 0.023552114 0.007329468 0.001846871</div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style=""><br></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style="">Veja que o primeiro e o segundo valores sãi iguais (o que na verdade corresponde a i=1) e é isso que quero resolver.</div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style=""><br></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style="">Abraços,</div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style=""><br></div><div id="yui_3_16_0_1_1416322549117_21482" dir="ltr" class="" style="">Daniela</div><div></div><div id="yui_3_16_0_1_1416322549117_21483"> </div><div class="signature" id="yui_3_16_0_1_1416322549117_21666"><div id="yui_3_16_0_1_1416322549117_21665"><span style="font-size:16px;" id="yui_3_16_0_1_1416322549117_21696">-- <br><font size="2" id="yui_3_16_0_1_1416322549117_21695">Daniela Rodrigues Recchia</font></span></div><div id="yui_3_16_0_1_1416322549117_21694"><span style="font-size:16px;" id="yui_3_16_0_1_1416322549117_21693"><font size="2" id="yui_3_16_0_1_1416322549117_21692">M.Sc. in Statistics</font><br><var id="yui-ie-cursor"></var></span></div><div style="color:rgb(0, 0, 0);font-size:16px;font-family:arial, helvetica, clean, sans-serif;background-color:transparent;font-style:italic;" id="yui_3_16_0_1_1416322549117_21691"><span style="font-size:16px;"><font size="2"><br></font></span></div><div style="color:rgb(0, 0, 0);font-family:arial, helvetica, clean, sans-serif;font-style:normal;background-color:transparent;" id="yui_3_16_0_1_1416322549117_21690"><span style="font-size:13px;" id="yui_3_16_0_1_1416322549117_21689">Tel. (Mobil): +49 01578 7588382</span></div><div id="yui_3_16_0_1_1416322549117_21688"><span style="font-size:16px;" id="yui_3_16_0_1_1416322549117_21687"><font size="2" id="yui_3_16_0_1_1416322549117_21686">Germany</font></span></div><div style="color:rgb(0, 0, 0);font-size:13px;font-family:arial, helvetica, clean, sans-serif;background-color:transparent;font-style:normal;" id="yui_3_16_0_1_1416322549117_21681"><span style="font-size:16px;"><font size="2"><br></font></span></div><div><span style="font-size:10px;"></span></div><div id="yui_3_16_0_1_1416322549117_21677"><span style="font-size:10px;" id="yui_3_16_0_1_1416322549117_21680"><span style="line-height:18px;font-family:'helvetica neue', arial, sans-serif;font-style:italic;" id="yui_3_16_0_1_1416322549117_21679"><font size="2" id="yui_3_16_0_1_1416322549117_21678">“The standard of success in life is absolutely the amount of joy you feel”</font></span><br></span><span style="font-size:13px;" id="yui_3_16_0_1_1416322549117_21685"><span style="font-size:10px;" id="yui_3_16_0_1_1416322549117_21684"><span style="line-height:18px;font-family:'helvetica neue', arial, sans-serif;font-style:italic;" id="yui_3_16_0_1_1416322549117_21683"><font size="2" id="yui_3_16_0_1_1416322549117_21682">Esther & Jerry Hick</font></span></span><br></span></div><div id="yui_3_16_0_1_1416322549117_21676"><span style="font-size:small;"> </span></div></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"> <font size="2" face="Arial"> Em Terça-feira, 18 de Novembro de 2014 18:01, Manoel Galdino <mcz.fea@gmail.com> escreveu:<br> </font> </div> <br><br> <div class="y_msg_container"><div id="yiv3291562842"><div><div dir="ltr">Seu código tá bem ruim de ler...<div>Recomendo que você utilize o GitHub para postar seus códigos. Além disso, ajudaria se você dissesse o que pretende com a função. Do jeito que está é bem difícil ajudar.</div><div><br clear="none"></div><div>Mas voltando a sua pergunta. Fiz algumas alterações no código (só reorganizei). Veja se está funcionando...</div><div><br clear="none"></div><div><br clear="none"><div><br clear="none"></div><div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;">conv<- function(n) {</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> u <- matrix(nrow=n+1,ncol=1)</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> for (i in 0:n) { </div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> if (i==0) {</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> f <- function(i) {pnorm(i,m,dp)}</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> p<- function(i) {f(i)}</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> } else {</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> f <- function(i) {pnorm(i*1/n,m,dp)}</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> p <- function(i) {f(i)-f(i-1)}<br clear="none"></div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> }</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> I <- function (i) {-p(i)*log2(p(i))}</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> u[i,]<- I(i)}</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> y <-u[0:n,]</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"> y }</div></div></div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"><br clear="none"></div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"><br clear="none"></div><div style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;">Use o browser pra você debugar seu código. <a rel="nofollow" shape="rect" target="_blank" href="http://adv-r.had.co.nz/Exceptions-Debugging.html">No livro do Hadley</a> tem mais explicações sobre como debugar um código.</div><div dir="ltr" style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"><br clear="none"></div><div style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;">abçs</div><div style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;">M</div><div style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;"><br clear="none"></div><div style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;">ps.: Onde que os argumentos 'm' e 'dp' entram na função? Eles são valores globais? Não é a prática mais recomendada...</div><div style="color:rgb(0,0,0);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:12px;">ps.2: a cada iteração do loop, as funções f e p são subscritas, de forma que ao final você terá uma f e uma p, dada por i igual a n. É isso mesmo que você quer?</div></div><div class="yiv3291562842gmail_extra"><br clear="none"><div class="yiv3291562842gmail_quote">On Tue, Nov 18, 2014 at 2:00 PM, Daniela Recchia <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:daniela_recchia@yahoo.com.br" target="_blank" href="mailto:daniela_recchia@yahoo.com.br">daniela_recchia@yahoo.com.br</a>></span> wrote:<br clear="none"><blockquote class="yiv3291562842gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="yiv3291562842yqt2918287992" id="yiv3291562842yqt44159"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12px;"><div dir="ltr"><span>Boa tarde pessoal,</span></div><div dir="ltr"><span><br clear="none"></span></div><div dir="ltr"><span>tenho certeza que meu problema é bem simples mas não estou conseguindo visualizar a solução, criei a função:</span></div><div dir="ltr"><span><br clear="none"></span></div><div dir="ltr">conv<- function(n){u<-matrix(nrow=n+1,ncol=1)</div><div dir="ltr"> for (i in 0:n) { </div><div dir="ltr"> if(i==0) f<-function(i) {pnorm(i,m,dp)}</div><div dir="ltr"> else f<- function(i) {pnorm(i*1/n,m,dp)}</div><div dir="ltr"> if(i==0) p<- function(i) {f(i)}</div><div dir="ltr"> else p<- function(i) {f(i)-f(i-1)}</div><div dir="ltr"> I<- function (i) {-p(i)*log2(p(i))}</div><div dir="ltr"> u[i,]<- I(i)}</div><div dir="ltr"> y<-u[0:n,]; y}</div><div dir="ltr"></div><div dir="ltr"><br clear="none"></div><div dir="ltr">Mas para o caso onde i==0 não está sendo feito corretamente, e preciso idendificar se não tenho problemas na condição <span style="font-size:11.8181819915771px;">f(i)-f(i-1), o que ele retorna no lugar do i=0 é o valor do i=1.</span></div><div dir="ltr"><span style="font-size:11.8181819915771px;"><br clear="none"></span></div><div dir="ltr"><span style="font-size:11.8181819915771px;">Alguém saberia me dar uma dica?</span></div><div dir="ltr"><span style="font-size:11.8181819915771px;"><br clear="none"></span></div><div dir="ltr"><span style="font-size:11.8181819915771px;">Abraços,</span></div><div dir="ltr"><span style="font-size:11.8181819915771px;"><br clear="none"></span></div><div dir="ltr"><span style="font-size:11.8181819915771px;">Daniela</span></div><span class="yiv3291562842HOEnZb"><font color="#888888"></font></span><div></div><div> </div><div><div><span style="font-size:16px;">-- <br clear="none"><font>Daniela Rodrigues Recchia</font></span></div><div><br clear="none"></div><div style="color:rgb(0,0,0);font-size:13px;font-family:arial, helvetica, clean, sans-serif;background-color:transparent;font-style:normal;"><span style="font-size:16px;"><font><br clear="none"></font></span></div><div><span style="font-size:10px;"></span></div><div><span style="font-size:10px;"><span style="line-height:18px;font-family:'helvetica neue', arial, sans-serif;font-style:italic;"><font>“The standard of success in life is absolutely the amount of joy you feel”</font></span><br clear="none"></span><span style="font-size:13px;"><span style="font-size:10px;"><span style="line-height:18px;font-family:'helvetica neue', arial, sans-serif;font-style:italic;"><font>Esther & Jerry Hick</font></span></span><br clear="none"></span></div><div><span style="font-size:small;"> </span></div></div></div></div></div><br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">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.<br clear="none"></blockquote></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div class="yiv3291562842gmail_signature">Manoel Galdino<br clear="none"><a rel="nofollow" shape="rect" target="_blank" href="https://sites.google.com/site/galdinomcz/">https://sites.google.com/site/galdinomcz/</a><br clear="none"></div>
</div></div></div><br><br></div> </div> </div> </div> </div></body></html>