[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