2.5 Pacote irtoys

O pacote irtoys possui uma função chamada est() que possui características similares ao software BILOG quanto as opções de estimação. Inclusive, é possível utilizar as ferramentas de estimação do BILOG caso esse esteja instalado em seu computador.

Nesse material, não será abordado o BILOG. Será utilizado o pacote ltm como ferramenta de estimação dos parâmetros dos modelos de interesse.

Para mais informações sobre a função est() veja o arquivo de ajuda do pacote (help(est)).

O ajuste de modelos com três parâmetros pode ser obtido da seguinte maneira:

O objeto manha.par contém os valores dos parâmetros a, b e c do modelo com 3 parâmetros. Para visualizas os primeiros parâmetros digite:

Para obter o θ ou a proficiência ou habilidade de cada indivíduo, após a estimação dos parâmetros do modelo, pode-se utilizar a função eap(). É necessário fornecer o arquivo de respostas e o objeto com os parêmetros estimados do modelo:

O objeto manha.sco contém o vetor θ de cada indivíduo.

Um objeto final, com os scores, identificação e posição dos respondentes pode ser criado com os seguintes comandos:

Pode-se visualizar o resultado ordenado pelo número de acertos,

Se for de interesse, pode-se alterar a escala do score com a função score.transform() do pacote CTT. Por exemplo, pode-se mudar a escala N(0,1) para uma escala N(500,10):

Para obtenção da curva característica pode-se obter um plot do objeto que contém o resultado da estimação dos parâmetros.

Figura 3: Curva característica dos itens.

  > plot(manha.tpm,item=1:30,sub='SARESP-manhã',legend=F)

Figura 4: Curva característica dos itens.

  > par(mfrow=c(2,2))
  > plot(manha.tpm,items=1:5)
  > plot(manha.tpm,type="IIC",items=1:5)
  > plot(manha.tpm,type="IIC",items=0)
  > par(mfrow=c(1,1))

A estimativa da proficiência de cada indivíduo pode ser obtida com a função factor.scores():

Figura 5: Proficiência para o turno manhã.

  > plot(manha.prof)