Estou tentando ganhar familiaridade com  Teoria de Resposta ao Item no R.

Eu gostaria de gerar um conjunto de itens segundo o graded response model e, a partir deles, um conjunto de respostas (pelo que vi, a função rmvordlogis do pacote ltm faz isso)

Depois, eu gostaria de rodar uma TRI nos dados para ver se ela recupera os parâmetros originais dos itens.

O problema é que eu não consigo gerar os dados.

A função rmvordlogis me dá um alerta deste tipo:

 

 

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