fix mixture ACOV

philchalmers · Aug 11, 2018 · 1f39651 · 1f39651
1 parent 7102031
commit 1f39651
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git a/R/EMstep.group.R b/R/EMstep.group.R
@@ -451,6 +451,7 @@ EM.group <- function(pars, constrain, Ls, Data, PrepList, list, Theta, DERIV, so
                     any.prior=FALSE, itemclass=10L,
                     dat=matrix(sapply(1L:length(pars),
                                       function(g) sum(pars[[g]][[J+1L]]@rr)), 1L))
+            mixtype@item.Q <- matrix(1, nrow = mixtype@ncat, ncol = mixtype@nfact)
             deriv <- Deriv(mixtype, Theta = matrix(1, 1L, length(pars)), estHess=TRUE)
             h[mixtype@parnum, mixtype@parnum] <- deriv$hess
         } else mixtype <- NULL

diff --git a/tests/tests/test-04-multipleGroup.R b/tests/tests/test-04-multipleGroup.R
@@ -63,7 +63,7 @@ test_that('one factor', {
     set.seed(1)
     mod_mixture <- suppressWarnings(multipleGroup(dat, 1, itemtype = 'Rasch', GenRandomPars = TRUE,
                                  verbose = FALSE, dentype = 'mixture-2', SE=TRUE))
-    expect_equal(extract.mirt(mod_mixture, 'condnum'), 115.5747, tolerance=1e-4)
+    expect_equal(extract.mirt(mod_mixture, 'condnum'), 112.2736, tolerance=1e-4)
     so <- summary(mod_mixture, verbose=FALSE)
     expect_equal(so[[1]]$class_proportion, .512121, tolerance=1e-4)