latent regression model in mirt

242 views
Skip to first unread message

Balal Ezanloo

unread,
Apr 14, 2018, 10:01:57 AM4/14/18
to mirt-package
Hi Phil

I know "mirt" can handle below model. i can run it by below syntax, but i am not sure about it. also where is the regression coefficient of "ksi" on Z in output?any help?

covdata <- # a binary matrix for 1000 persone in 5 variables(v1:v5)
z <- rowMeans(covdata)
data <- # another binary matrix for the same 1000 persone in 5 variables

LRM.2pl.out <- mirt(data , 1, '2PL', covdata=covdata, formula = ~ z)


Phil Chalmers

unread,
Apr 14, 2018, 11:55:46 AM4/14/18
to Balal Ezanloo, mirt-package
Xi is set to 1 for the 2PL model for identification. Otherwise, you could freely estimate the latent variable variance term if you set one of the item slope parameters  to 1 (this is a common strategy in SEM). 

Phil

--
You received this message because you are subscribed to the Google Groups "mirt-package" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mirt-package+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Balal Ezanloo

unread,
Apr 14, 2018, 1:23:34 PM4/14/18
to mirt-p...@googlegroups.com
Thanks

Did I introduce the model correctly to the "mirt"? i am not sure about it.

Also,  in the output did i find correctly the regression coefficient of "ksi" on Z, as below?
 
                F1
  (Intercept)  0.000
r.p.comleted 1.872

Phil

To unsubscribe from this group and stop receiving emails from it, send an email to mirt-package...@googlegroups.com.

Balal Ezanloo

unread,
Apr 16, 2018, 3:25:45 AM4/16/18
to mirt-package
Hi Phil

How can i add an intercept of Z in "formula" part

LRM.2pl.out <- mirt(data , 1, '2PL', covdata=covdata, formula = ~ z) 

Phil Chalmers

unread,
Apr 16, 2018, 9:13:26 PM4/16/18
to Balal Ezanloo, mirt-package
Only by modifying the starting values explicitly via the pars input. But, why would you want to? This is fixed at 0 for identification reasons.

Phil

To unsubscribe from this group and stop receiving emails from it, send an email to mirt-package+unsubscribe@googlegroups.com.

balal izanloo

unread,
Apr 17, 2018, 1:22:48 AM4/17/18
to Phil Chalmers, mirt-package
Thanks

I have a model that presented hear. other peoples run it by "mplus" but i want to run it by "mirt". "mirt" results is not the same as "mplus" one. i want to add intercepts to my syntax in mirt and compare the results.




Reply all
Reply to author
Forward
0 new messages