Question

我有一个矩阵，我试图通过mirt包的mirt函数运行：

resp.freq <- data.frame(matrix(c(11, 46, 12, 31, 13, 8, 21, 20, 22, 68, 23, 12,
                                 31,  1, 32, 12, 33, 11), nrow = 9, ncol = 2,
                               byrow = T))
dados3 <- matrix(rep(resp.freq$X1, resp.freq$X2), ncol = 1)
dados3 <- data.frame(cbind(as.numeric(substr(dados3, 1, 1)),
                           as.numeric(substr(dados3, 2, 2))))

require(mirt)
m1 <- mirt(dados3, 1, D = 1, SE = TRUE)
coef(m1)

我曾经通过 mirt 0.5.0 运行此操作，并获得以下最终结果：

$X1
        a1    d1     d2
pars 1.913 0.633 -3.152
SE   0.228 0.181  0.322

$X2
        a1    d1     d2
pars 1.659 1.121 -2.488
SE   0.200 0.190  0.261

但是，我的工作站已经更新，并且已经 mirt 1.2.1 ，现在我在运行m1 <- mirt(...时得到以下输出：

Iteration: 1000, Log-Lik: -387.670, Max-Change: 0.00007
EM iterations terminated after 1000 iterations.


Calculating information matrix...
Warning message:
In loadESTIMATEinfo(info = info, ESTIMATE = ESTIMATE, constrain = constrain) :
  Negative SEs set to NaN.

coef(m1)给了我：

$X1
            a1    d1     d2
par     11.515 2.732 -13.98
CI_2.5     NaN 0.685    NaN
CI_97.5    NaN 4.778    NaN

$X2
           a1    d1     d2
par     1.099 0.903 -2.139
CI_2.5  0.712 0.561 -2.611
CI_97.5 1.485 1.245 -1.668

$GroupPars
        MEAN_1 COV_11
par          0      1
CI_2.5      NA     NA
CI_97.5     NA     NA

我已阅读changelog，但无法找到行为更改的原因。我曾试图拼命改变mirt()的参数，但无法达到类似的收敛水平。是什么给了什么？

输出的变化取决于mirt版本

0 个答案: