[R] ability estimate with GRM of IRT
dimitris.rizopoulos at med.kuleuven.be
Fri Jan 19 09:20:34 CET 2007
you need to use the 'resp.patterns' argument of factor.scores(); look
at ?factor.scores for more info, e.g.,
fit <- grm(Environment)
factor.scores(fit, resp.patterns = Environment)
I hope it helps.
ps, it'd be better not to call your fitted GRM model 'grm'; check
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
----- Original Message -----
From: "SHEN,FENG" <fshen at ufl.edu>
To: <r-help at stat.math.ethz.ch>
Sent: Friday, January 19, 2007 5:30 AM
Subject: [R] ability estimate with GRM of IRT
> Hi my friends,
> I have an issue with ability estimates when running GRM of IRT. I
> have responses from 242 subjects but got 183 ability estimates.
> Below is what I did to get the estimates.
> 1) I have a csv file "P1.csv" and I imported it into R and loaded
> the "ltm" package by doing:
> 2) I created a subset that included columns 2 to 9 for the
> analysis by doing:
> 3) I converted the subset into data.frame format by doing:
> 4) I checked the descriptive stats for the s1df by doing:
> And it was confirmed that 242 subjects' responses were imported.
> 5)I ran GRM on the s1df dataset by doing:
> grm<-grm(s1df, Hessian=T)
> 6)Finally, I ran the ability estimated by doing:
> And I got 183 factor-scores for observed response patterns.
> Besides, the 183 estimates are not for the first 183 subjects of
> the 242 because the response patterns of No 1 on the factor-scores
> list do not match those of No 1 in the s1 dataset and the same
> thing holds true for the rest of the response patterns.
> Could you help me find out what the problem is?
> Many thanks in advance!
> R-help at stat.math.ethz.ch mailing list
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help