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
_______________________________________________
R-br mailing list
[email protected]
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.

Responder a