confidence intervals around expected total scores from multi-group analysis

[Cathy Labrish]

Hi Phil, 

Hope this finds you well.  Am trying to get a plot of expected total scores with 95% confidence from a multigroup analysis using the mirt plot method.  Confidence intervals are not appearing.  Can you suggest what I am going incorrectly.  I am conducting a DIF/DTF analysis on a measure of 23 items, and setting anchors as the items listed in invar.final.set1.  

Here is the code:

group <- pooledat$group

itcol <- colnames(pooledat)[-24]

invar.final.set1 <- itcol[c(2,8,13,14,16,18,20,22,23)]
dif.final.set1 <- c(1,3:7,9:12,15,17,19,21)

remove(modMG_final.set1)
modMG_final.set1 <- multipleGroup(pooledat[itcol], 1, group, SE=T,
                 invariance = c(invar.final.set1,
                                           'free_mean',
                                           'free_var'),
                            technical = list(NCYCLES = 2000))

plot(modMG_final.set1,MI=2000)

data is attached.