<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><div>O problema está nesta parte aqui!</div><div><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; ">Nem sempre o vetor p.vec possui um elemento maior que o seed, e quando isto ocorre retorna interger(0)</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; ">Outro problema é que o which pega a posição do elemento, então, quando vc pega o max(temp), se o vetor tem tamanho 20, por exemplo, ele vai retornar o 20 e não o
elemento.</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; ">Precisa corrigir ai!</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "> temp = which(p.vec>seed) #we use this temporary vector in order to avoid warnings should its length be zero</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>print(temp)</div><div>
if(length(temp)>0){resp=max(temp)}else{resp=1}</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><br></div></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">Valeu! </div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> Fábio Mathias Corrêa<br><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> Universidade Estadual de Santa Cruz<br>Departamento de Ciências Exatas e da Terra - DCET</div><br><br><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna<br>CEP 45662-900. Ilhéus-Bahia</div><div style="font-family: 'times new roman', 'new
york', times, serif; font-size: 12pt; "><br><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">Tel.: 73-3680-5076<br> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Felipe Buchbinder <felbuch@gmail.com><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, 13 de Setembro de 2012 15:24<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Dúvida função no R<br> </font> </div> <br><div id="yiv717438356">Obrigado pela referência,<div><br></div><div>Preparei um CMR que envio anexo.</div><div><br></div><div>Mais uma vez,
obrigado </div><div>Felipe<br><br><div class="yiv717438356gmail_quote">2012/9/13 FHRB Toledo <span dir="ltr"><<a rel="nofollow" ymailto="mailto:fernandohtoledo@gmail.com" target="_blank" href="mailto:fernandohtoledo@gmail.com">fernandohtoledo@gmail.com</a>></span><br>
<blockquote class="yiv717438356gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div>Felipe,</div><div><br></div>Para CMR, consulte: <a rel="nofollow" target="_blank" href="https://gist.github.com/1088208">https://gist.github.com/1088208</a><div>
<br></div><div>att,</div><div>FH<div><div class="yiv717438356h5"><br><br><div class="yiv717438356gmail_quote">2012/9/13 Felipe Buchbinder <span dir="ltr"><<a rel="nofollow" ymailto="mailto:felbuch@gmail.com" target="_blank" href="mailto:felbuch@gmail.com">felbuch@gmail.com</a>></span><br>
<blockquote class="yiv717438356gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">O k foi definido no início do código, fora da função, como k = 5.<div><br></div><div>Desculpe a ignorância, mas o que é um CMR? Como eu posso fazer um para enviar à lista?</div>
<div><br></div><div>Obrigado pela disposição,</div>
<div><br>Felipe<div><div><br><br><div class="yiv717438356gmail_quote">2012/9/13 Leonard de Assis <span dir="ltr"><<a rel="nofollow" ymailto="mailto:assis.leonard@gmail.com" target="_blank" href="mailto:assis.leonard@gmail.com">assis.leonard@gmail.com</a>></span><br><blockquote class="yiv717438356gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>
<div>Onde foi definido o k?<br>
<pre>[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
</pre>
Em 13/09/2012 12:56, Felipe Buchbinder escreveu:<br>
</div>
<blockquote type="cite"><div><div>Prezados professores,
<div><br>
</div>
<div>Estou com um problema em uma função no R.</div>
<div>A função é a seguinte:</div>
<div><br>
</div>
<blockquote style="margin:0 0 0 40px;border:none;padding:0px;">
<blockquote style="margin:0 0 0 40px;border:none;padding:0px;">
<blockquote class="yiv717438356gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex;">resp
<- function(subject,item,trait)<br>
{<br>
<br>
p.vec = NULL<br>
<br>
for(category in 1:k) # k = 5<br>
{<br>
p = ogive ( theta[subject,trait] , a[item] ,
b[item,category] , D ) # aqui é uma logit (D=1) ou probit
(D=1.7)<br>
p.vec = c(p.vec,p)<br>
}<br>
<br>
#Generates a random number<br>
<br>
seed<-runif(1,0.2,0.8)<br>
<br>
#Selects an answer<br>
<br>
temp = which ( p.vec>seed ) <br>
if( length (temp) > 0 ) { resp = max (temp) } else {
resp = 1 }<br>
resp<br>
}</blockquote>
</blockquote>
</blockquote>
<div><br>
</div>
<div><br>
</div>
<div>Esta função está inserida em um loop e deveria me retornar um
número inteiro entre 1 e 5.</div>
<div>Da primeira vez, isto acontece.</div>
<div>Mas a partir da segunda, ela começa a me retornar valores
como 20, 16,...</div>
<div><br>
</div>
<div>Alguma idéia de porque isso acontece???</div>
<div><br>
</div>
<div>Desde já, obrigado!</div>
<div><br>
</div>
<div>Felipe Buchbinder</div>
<div><br>
</div>
<div><br>
</div>
<br>
<fieldset></fieldset>
<br>
</div></div><div><pre>_______________________________________________
R-br mailing list
<a rel="nofollow" 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>
<a rel="nofollow" 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" 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>
</div></blockquote>
<br>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" 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>
<a rel="nofollow" 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>
Leia o guia de postagem (<a rel="nofollow" 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></blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" 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>
<a rel="nofollow" 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>
Leia o guia de postagem (<a rel="nofollow" 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></blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a rel="nofollow" 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>
<a rel="nofollow" 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>
Leia o guia de postagem (<a rel="nofollow" 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></blockquote></div><br></div>
</div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" 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><br> </div> </div> </div></div></body></html>