[R-br] Pacote ltm

[Felipe Buchbinder]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130827/5c98431f/attachment.html>