Unfortunately, the data file contains the data for a polytomous IRT, and it needs one recode before you can use it. You do not need to do much, just recode the 2's as 1's. This code seems to work, and the presentation explains how to interpret the findings.
GHQ28 <- read.table(file.choose(),
header=TRUE, sep="\t",
na.strings="NA", dec=".",
strip.white=TRUE)
Anxiety<-GHQ28[ ,8:14]
library(ltm)
Anxiety2<-Anxiety
apply(Anxiety2, 2, table)
#recodes 2's to 1's
Anxiety2[Anxiety2==2]<-1
apply(Anxiety2, 2, table)
Result1PL<-rasch(Anxiety2)
summary(Result1PL)
plot(Result1PL)
Result2PL<-ltm(Anxiety2 ~ z1)
summary(Result2PL)
plot(Result2PL)
anova(Result1PL,Result2PL)
#May require more than one run to converge
Result3PL<-tpm(Anxiety2,control=c(iter.qN=10000,
GHk=50,verbose=TRUE),start.val="random")
summary(Result3PL)
plot(Result3PL)
anova(Result2PL,Result3PL)
If you are just starting to learn about item response theory, it might help if you watched an introductory video. There is no R code, but Nick Shryane is slow and deliberate. This video should fill in the gaps for many of you. It's called "What is Item Response Theory?" After 46 minutes, you ought to have a pretty good idea.
As always, I learn a great deal from your comments and emails. Thanks for your help.
No comments:
Post a Comment