[R-br] Dúvida função no R

[Felipe Buchbinder]

Fábio,

Muito obrigado pela ajuda. Mas o *if *já não dá conta do caso em que p.vec
não possui um elemento maior do que seed? E como p.vec é construido dentro
de um loop que adiciona um elemento durante *k* iterações, como ele pode
ter mais do que *k* elementos? Eu achava que já tinha controlado para
evitar esses casos... não?

mais uma vez obrigado,

Felipe

2012/9/13 Fabio Mathias Corrêa <fabio.ufla em yahoo.com.br>

> O problema está nesta parte aqui!
>
> Nem sempre o vetor  p.vec possui um elemento maior que o seed, e quando
> isto ocorre retorna interger(0)
>
> 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.
>
> Precisa corrigir ai!
>
>  temp = which(p.vec>seed) #we use this temporary vector in order to avoid
> warnings should its length be zero
> print(temp)
>   if(length(temp)>0){resp=max(temp)}else{resp=1}
>
> Valeu!
>
>         Fábio Mathias Corrêa
>
>    Universidade Estadual de Santa Cruz
> Departamento de Ciências Exatas e da Terra - DCET
>
>
> Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna
> CEP 45662-900. Ilhéus-Bahia
>
>
> Tel.: 73-3680-5076
>   ------------------------------
> *De:* Felipe Buchbinder <felbuch em gmail.com>
> *Para:* r-br em listas.c3sl.ufpr.br
> *Enviadas:* Quinta-feira, 13 de Setembro de 2012 15:24
> *Assunto:* Re: [R-br] Dúvida função no R
>
> Obrigado pela referência,
>
> Preparei um CMR que envio anexo.
>
> Mais uma vez, obrigado
> Felipe
>
> 2012/9/13 FHRB Toledo <fernandohtoledo em gmail.com>
>
> Felipe,
>
> Para CMR, consulte: https://gist.github.com/1088208
>
> att,
> FH
>
>
> 2012/9/13 Felipe Buchbinder <felbuch em gmail.com>
>
> O k foi definido no início do código, fora da função, como k = 5.
>
> Desculpe a ignorância, mas o que é um CMR? Como eu posso fazer um para
> enviar à lista?
>
> Obrigado pela disposição,
>
> Felipe
>
>
> 2012/9/13 Leonard de Assis <assis.leonard em gmail.com>
>
>  Onde foi definido o k?
>
> []s
> Leonard de Assis
> assis <dot> leonard <at> gmail <dot> com
>
> Em 13/09/2012 12:56, Felipe Buchbinder escreveu:
>
> Prezados professores,
>
>  Estou com um problema em uma função no R.
> A função é a seguinte:
>
>   resp <- function(subject,item,trait)
> {
>
>   p.vec = NULL
>
>   for(category in 1:k)  # k = 5
>   {
>     p = ogive ( theta[subject,trait] , a[item] , b[item,category] , D )
>  # aqui é uma logit (D=1) ou probit (D=1.7)
>     p.vec = c(p.vec,p)
>   }
>
>   #Generates a random number
>
>   seed<-runif(1,0.2,0.8)
>
>   #Selects an answer
>
>   temp = which ( p.vec>seed )
>   if( length (temp) > 0 ) { resp = max (temp) } else { resp = 1 }
>   resp
> }
>
>
>
>  Esta função está inserida em um loop e deveria me retornar um número
> inteiro entre 1 e 5.
> Da primeira vez, isto acontece.
> Mas a partir da segunda, ela começa a me retornar valores como 20, 16,...
>
>  Alguma idéia de porque isso acontece???
>
>  Desde já, obrigado!
>
>  Felipe Buchbinder
>
>
>
>
> _______________________________________________
> R-br mailing listR-br em listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120914/8e0bd9d8/attachment.html>