Prezados professores,
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
}