<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)<br></div><div>Allaman<br></div><div><br></div><div><br></div><div><br></div></div>