[R-br] Gerar e analizar dados por TRI

[Alexandre Serpa]

Felipe,

Na linha de comando  "x <- rmvordlogis(10,thetas,model = "grm")" você deve
adicionar IRT=F.

Você pode eliminar o modelo na geração dos dados randômicos "(model="grm")"
e especificar a parametrização ao chamar o comando "grm".

Desse modo, vai funcionar "x<-rmvordlogis(10,thetas)".

Abs


Alexandre Serpa

Psicólogo
Especialista em Métodos Computacionais Estatísticos
Mestre em Educação
Doutorando em Psicologia
email:   serpa.alexandre em gmail.com


Em 8 de setembro de 2013 15:38, Ivan Bezerra Allaman <
ivanalaman em yahoo.com.br> escreveu:

> A mensagem de erro é extremamente óbvia!!! e o seu exemplo não é
> reproduzível pois faltou colocar "library(ltm)".
>
> O problema está na especificação do argumento "model"!!! Ao abrir a função
> tem-se:
>
> probs <- if (model == "grm") {
>             gammas <- lapply(thetas, function(x) {
>                 nx <- length(x)
>                 if (IRT)
>                   cbind(plogis(x[nx] * (z - matrix(x[-nx], n,
>                     nx - 1, TRUE))), 1)
>                 else cbind(plogis(matrix(x[-nx], n, nx - 1, TRUE) -
>                   x[nx] * z), 1)
>             })
>             lapply(gammas, function(x) {
>                 nc <- ncol(x)
>                 cbind(x[, 1], x[, 2:nc] - x[, 1:(nc - 1)])
>             })
>         }
>
> 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,
>
> X <- matrix(0, n, p)
>         for (j in 1:p) {
>             for (i in 1:n) X[i, j] <- sample(ncatg[j], 1, prob =
> probs[[j]][i,
>                 ])
>         }
>
> 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.
>
> Detectei o erro, porém não posso lhe dar a solução pois TRI não é algo que
> entendo!!
>
> (s,f,p)
> Allaman
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130908/0928d1f2/attachment.html>