Professores,

 

Estou tentando ganhar familiaridade com a função rmvordlogis do pacote ltm, de Teoria de Resposta ao Item.

 

Para isto, tentei gerar a função tal qual descrita no arquivo de ajuda do pacote. A única diferença é que eu inclui o parâmetro model = "grm''

 

thetas <- lapply(1:5, function(u) c(seq(-1, 1, len = 2), 1.2))
x <- rmvordlogis(10,thetas,model = "grm")

 

A resposta que eu obtive foi a seguinte:

 

Erro em sample(ncatg[j], 1, prob = probs[[j]][i, ]) : 
  probabilidade não positiva

O que pode estar ocasionando este problema e como eu poderia gerar resposas para o caso GRM?

 

Desde já obrigado!

 

Felipe Buchbinder