[R-sig-ME] Multi-level Rasch Model Per Douglas Bates' paper
Simon Harmel
@|m@h@rme| @end|ng |rom gm@||@com
Wed May 13 17:07:38 CEST 2020
Hi All!
I'm following this paper <https://www.jstatsoft.org/article/view/v020i02> (
https://www.jstatsoft.org/article/view/v020i02) by Prof. Bates where after
fitting the model (*pp. 14-15*), they obtain what they call *item
easiness* *"from
the estimates of the fixed effects and the conditional modes of the random
effects."*
In short, I wonder how to obtain item easiness estimates for each of my
models (m1 & m2) below? *Thank you, Simon*
library(glmmTMB)
dat <- read.csv('https://raw.githubusercontent.com/ilzl/i/master/d.csv')
form11 <- y ~ item_type + (1 | item_id) + (1 | person_id)
form22 <- y ~ item_type + gender + (1 | item_id) + (1 | person_id)
m1 <- glmmTMB(form11, data = subset(dat, person_id <= 40),
family = beta_family())
m2 <- glmmTMB(form22, data = subset(dat, person_id <= 40),
family = beta_family())
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list