<div dir="ltr"><div>Professores,</div><div> </div><div>Estou tentando ganhar familiaridade com a função rmvordlogis do pacote ltm, de Teoria de Resposta ao Item.</div><div> </div><div>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'' </div>
<div> </div><div>thetas <- lapply(1:5, function(u) c(seq(-1, 1, len = 2), 1.2))<br>x <- rmvordlogis(10,thetas,model = "grm")</div><div> </div><div>A resposta que eu obtive foi a seguinte:</div><div> </div>
<span style="font:13px/15px "Lucida Console";color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:pre-wrap;border-collapse:separate;font-size-adjust:none;font-stretch:normal;background-color:rgb(225,226,229)"><pre style="margin:0px;line-height:1.2;font-family:"Lucida Console";font-size:10pt!important;white-space:pre-wrap!important;outline-style:none" tabindex="0">
<span>Erro em sample(ncatg[j], 1, prob = probs[[j]][i, ]) : 
  probabilidade não positiva</span></pre><div></div></span><div>O que pode estar ocasionando este problema e como eu poderia gerar resposas para o caso GRM?</div><div> </div><div>Desde já obrigado!</div><div> </div><div>
Felipe Buchbinder</div><div> </div><div><br></div></div>