<div dir="ltr"><div><div>Professores,<br><br></div>Obrigado pela ajuda. Modifiquei o meu CMR para o seguinte:<br><br>library(ltm)<br>thetas <- lapply(1:5, function(u) c(seq(-1, 1, len = 2), 1.2))<br>x <- rmvordlogis(10,thetas,model = "grm",IRT=F)<br>
grm(x)<br><br></div>Porém, o grm não recupera os parâmetros -1, 1 e 1.2 para os itens. Ao invés disso, ele encontrou:<br><br><span class="" style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229)"><pre tabindex="0" class="" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-width:initial;outline-color:initial;border-style:none;border-width:initial;border-color:initial;white-space:pre-wrap!important;margin:0px;line-height:1.2">
Call:
grm(data = x)
Coefficients:
Extrmt1 Extrmt2 Dscrmn
Item 1 -0.457 0.525 4.973
Item 2 -0.250 0.363 2.051
Item 3 -0.500 0.314 2.589
Item 4 -0.346 0.600 0.961
Item 5 -0.161 -8.819 -0.251
Log.Lik: -46.384</pre></span><div><br>E, de vez em quando, ele retorna a seguinte mensagem:<br><br><span class="" style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229)"><pre tabindex="0" class="" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-width:initial;outline-color:initial;border-style:none;border-width:initial;border-color:initial;white-space:pre-wrap!important;margin:0px;line-height:1.2">
<span class="" style="color:rgb(197,6,11)">Warning messages:
</span><span class="" style="color:rgb(197,6,11)">1: glm.fit: algorithm did not converge
</span><span class="" style="color:rgb(197,6,11)">2: glm.fit: fitted probabilities numerically 0 or 1 occurred
</span><span class="" style="color:rgb(197,6,11)">3: glm.fit: algorithm did not converge
</span><span class="" style="color:rgb(197,6,11)">4: glm.fit: fitted probabilities numerically 0 or 1 occurred
</span><span class="" style="color:rgb(197,6,11)">5: glm.fit: algorithm did not converge
</span><span class="" style="color:rgb(197,6,11)">6: glm.fit: fitted probabilities numerically 0 or 1 occurred<br><br></span><br></pre></span>Estes resultados são de se esperar? Por que o grm parece não recuperar os parâmetros originais?<br>
<br></div><div>Mais uma vez, obrigado,<br><br></div><div>Felipe<br> </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/9/8 Alexandre Serpa <span dir="ltr"><<a href="mailto:serpa.alexandre@gmail.com" target="_blank">serpa.alexandre@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Felipe,<div><br></div><div>Na linha de comando <span style="font-family:arial,sans-serif;font-size:13px">"x <- rmvordlogis(10,thetas,model = "grm")" você deve adicionar IRT=F.</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Você pode eliminar o modelo na geração dos dados randômicos "(model="grm")" e especificar a parametrização ao chamar o comando "grm".</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Desse modo, vai funcionar "x<-rmvordlogis(10,thetas)".</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Abs</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div class="gmail_extra"><br clear="all"><div>Alexandre Serpa<br><br>Psicólogo<br>Especialista em Métodos Computacionais Estatísticos<br>Mestre em Educação<div>Doutorando em Psicologia<br>email: <a href="mailto:serpa.alexandre@gmail.com" target="_blank">serpa.alexandre@gmail.com</a></div>
</div>
<br><br><div class="gmail_quote">Em 8 de setembro de 2013 15:38, Ivan Bezerra Allaman <span dir="ltr"><<a href="mailto:ivanalaman@yahoo.com.br" target="_blank">ivanalaman@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div class="h5">
<div dir="ltr"><div>A mensagem de erro é extremamente óbvia!!! e o seu exemplo não é reproduzível pois faltou colocar "library(ltm)".<br><br></div><div>O problema está na especificação do argumento "model"!!! Ao abrir a função tem-se:<br>
<br>probs <- if (model == "grm") {<br> gammas <- lapply(thetas, function(x) {<br> nx <- length(x)<br> if (IRT) <br> cbind(plogis(x[nx] * (z - matrix(x[-nx], n, <br>
nx - 1, TRUE))), 1)<br> else cbind(plogis(matrix(x[-nx], n, nx - 1, TRUE) - <br> x[nx] * z), 1)<br> })<br> lapply(gammas, function(x) {<br> nc <- ncol(x)<br>
cbind(x[, 1], x[, 2:nc] - x[, 1:(nc - 1)])<br> })<br> }<br><br></div><div>Aqui, este objeto é calculado de modo que contém valores negativos em suas listas. Logo, e claramente, a função "sample" irá acusar erro, uma vez que,<br>
<br>X <- matrix(0, n, p)<br> for (j in 1:p) {<br> for (i in 1:n) X[i, j] <- sample(ncatg[j], 1, prob = probs[[j]][i, <br> ])<br> }<br><br></div><div>Este objeto "probs" é justamente usado no argumento "prob" da função "sample" e como sabemos, não existe probabilidade negativa segundo os axiomas de Kolmogorov.<br>
<br></div><div>Detectei o erro, porém não posso lhe dar a solução pois TRI não é algo que entendo!!<br><br></div><div>(s,f,p)<span><font color="#888888"><br></font></span></div><span><font color="#888888"><div>
Allaman<br></div><div><br></div><div><br></div><div><br></div></font></span></div>
<br></div></div><div class="im">_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div></blockquote></div><br></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>