Translate mirt parameters into suitable structure for plink package
Source:R/read.mirt.R
read.mirt.RdThis function exports item parameters from the mirt package to the
plink package.
Arguments
- x
a single object (or list of objects) returned from
mirt, bfactor, or a single object returned bymultipleGroup- as.irt.pars
if
TRUE, the parameters will be output as anirt.parsobject- ...
additional arguments to be passed to
coef()
Author
Phil Chalmers rphilip.chalmers@gmail.com
Examples
# \donttest{
## unidimensional
library(plink)
data <- expand.table(LSAT7)
(mod1 <- mirt(data, 1))
#>
Iteration: 1, Log-Lik: -2668.786, Max-Change: 0.18243
Iteration: 2, Log-Lik: -2663.691, Max-Change: 0.13637
Iteration: 3, Log-Lik: -2661.454, Max-Change: 0.10231
Iteration: 4, Log-Lik: -2659.430, Max-Change: 0.04181
Iteration: 5, Log-Lik: -2659.241, Max-Change: 0.03417
Iteration: 6, Log-Lik: -2659.113, Max-Change: 0.02911
Iteration: 7, Log-Lik: -2658.812, Max-Change: 0.00456
Iteration: 8, Log-Lik: -2658.809, Max-Change: 0.00363
Iteration: 9, Log-Lik: -2658.808, Max-Change: 0.00273
Iteration: 10, Log-Lik: -2658.806, Max-Change: 0.00144
Iteration: 11, Log-Lik: -2658.806, Max-Change: 0.00118
Iteration: 12, Log-Lik: -2658.806, Max-Change: 0.00101
Iteration: 13, Log-Lik: -2658.805, Max-Change: 0.00042
Iteration: 14, Log-Lik: -2658.805, Max-Change: 0.00025
Iteration: 15, Log-Lik: -2658.805, Max-Change: 0.00026
Iteration: 16, Log-Lik: -2658.805, Max-Change: 0.00023
Iteration: 17, Log-Lik: -2658.805, Max-Change: 0.00023
Iteration: 18, Log-Lik: -2658.805, Max-Change: 0.00021
Iteration: 19, Log-Lik: -2658.805, Max-Change: 0.00019
Iteration: 20, Log-Lik: -2658.805, Max-Change: 0.00017
Iteration: 21, Log-Lik: -2658.805, Max-Change: 0.00017
Iteration: 22, Log-Lik: -2658.805, Max-Change: 0.00015
Iteration: 23, Log-Lik: -2658.805, Max-Change: 0.00015
Iteration: 24, Log-Lik: -2658.805, Max-Change: 0.00013
Iteration: 25, Log-Lik: -2658.805, Max-Change: 0.00013
Iteration: 26, Log-Lik: -2658.805, Max-Change: 0.00011
Iteration: 27, Log-Lik: -2658.805, Max-Change: 0.00011
Iteration: 28, Log-Lik: -2658.805, Max-Change: 0.00010
#>
#> Call:
#> mirt(data = data, model = 1)
#>
#> Full-information item factor analysis with 1 factor(s).
#> Converged within 1e-04 tolerance after 28 EM iterations.
#> mirt version: 1.44.3
#> M-step optimizer: BFGS
#> EM acceleration: Ramsay
#> Number of rectangular quadrature: 61
#> Latent density type: Gaussian
#>
#> Log-likelihood = -2658.805
#> Estimated parameters: 10
#> AIC = 5337.61
#> BIC = 5386.688; SABIC = 5354.927
#> G2 (21) = 31.7, p = 0.0628
#> RMSEA = 0.023, CFI = NaN, TLI = NaN
plinkpars <- read.mirt(mod1)
plot(plinkpars)
plot(mod1, type = 'trace')
# graded
mod2 <- mirt(Science, 1)
#>
Iteration: 1, Log-Lik: -1629.361, Max-Change: 0.50660
Iteration: 2, Log-Lik: -1617.374, Max-Change: 0.25442
Iteration: 3, Log-Lik: -1612.894, Max-Change: 0.16991
Iteration: 4, Log-Lik: -1610.306, Max-Change: 0.10461
Iteration: 5, Log-Lik: -1609.814, Max-Change: 0.09162
Iteration: 6, Log-Lik: -1609.534, Max-Change: 0.07363
Iteration: 7, Log-Lik: -1609.030, Max-Change: 0.03677
Iteration: 8, Log-Lik: -1608.988, Max-Change: 0.03200
Iteration: 9, Log-Lik: -1608.958, Max-Change: 0.02754
Iteration: 10, Log-Lik: -1608.878, Max-Change: 0.01443
Iteration: 11, Log-Lik: -1608.875, Max-Change: 0.00847
Iteration: 12, Log-Lik: -1608.873, Max-Change: 0.00515
Iteration: 13, Log-Lik: -1608.872, Max-Change: 0.00550
Iteration: 14, Log-Lik: -1608.872, Max-Change: 0.00318
Iteration: 15, Log-Lik: -1608.871, Max-Change: 0.00462
Iteration: 16, Log-Lik: -1608.871, Max-Change: 0.00277
Iteration: 17, Log-Lik: -1608.870, Max-Change: 0.00145
Iteration: 18, Log-Lik: -1608.870, Max-Change: 0.00175
Iteration: 19, Log-Lik: -1608.870, Max-Change: 0.00126
Iteration: 20, Log-Lik: -1608.870, Max-Change: 0.00025
Iteration: 21, Log-Lik: -1608.870, Max-Change: 0.00285
Iteration: 22, Log-Lik: -1608.870, Max-Change: 0.00108
Iteration: 23, Log-Lik: -1608.870, Max-Change: 0.00022
Iteration: 24, Log-Lik: -1608.870, Max-Change: 0.00059
Iteration: 25, Log-Lik: -1608.870, Max-Change: 0.00014
Iteration: 26, Log-Lik: -1608.870, Max-Change: 0.00068
Iteration: 27, Log-Lik: -1608.870, Max-Change: 0.00065
Iteration: 28, Log-Lik: -1608.870, Max-Change: 0.00019
Iteration: 29, Log-Lik: -1608.870, Max-Change: 0.00061
Iteration: 30, Log-Lik: -1608.870, Max-Change: 0.00012
Iteration: 31, Log-Lik: -1608.870, Max-Change: 0.00012
Iteration: 32, Log-Lik: -1608.870, Max-Change: 0.00058
Iteration: 33, Log-Lik: -1608.870, Max-Change: 0.00055
Iteration: 34, Log-Lik: -1608.870, Max-Change: 0.00015
Iteration: 35, Log-Lik: -1608.870, Max-Change: 0.00052
Iteration: 36, Log-Lik: -1608.870, Max-Change: 0.00010
plinkpars <- read.mirt(mod2)
plot(plinkpars)
plot(mod2, type = 'trace')
# gpcm
mod3 <- mirt(Science, 1, itemtype = 'gpcm')
#>
Iteration: 1, Log-Lik: -1689.735, Max-Change: 1.30401
Iteration: 2, Log-Lik: -1618.015, Max-Change: 0.28548
Iteration: 3, Log-Lik: -1615.635, Max-Change: 0.27916
Iteration: 4, Log-Lik: -1613.549, Max-Change: 0.13014
Iteration: 5, Log-Lik: -1613.306, Max-Change: 0.08712
Iteration: 6, Log-Lik: -1613.172, Max-Change: 0.09075
Iteration: 7, Log-Lik: -1612.802, Max-Change: 0.06897
Iteration: 8, Log-Lik: -1612.774, Max-Change: 0.05440
Iteration: 9, Log-Lik: -1612.756, Max-Change: 0.03767
Iteration: 10, Log-Lik: -1612.733, Max-Change: 0.03981
Iteration: 11, Log-Lik: -1612.724, Max-Change: 0.02653
Iteration: 12, Log-Lik: -1612.718, Max-Change: 0.03424
Iteration: 13, Log-Lik: -1612.700, Max-Change: 0.00646
Iteration: 14, Log-Lik: -1612.697, Max-Change: 0.01774
Iteration: 15, Log-Lik: -1612.695, Max-Change: 0.01250
Iteration: 16, Log-Lik: -1612.692, Max-Change: 0.01716
Iteration: 17, Log-Lik: -1612.690, Max-Change: 0.01090
Iteration: 18, Log-Lik: -1612.689, Max-Change: 0.01170
Iteration: 19, Log-Lik: -1612.687, Max-Change: 0.00599
Iteration: 20, Log-Lik: -1612.687, Max-Change: 0.00160
Iteration: 21, Log-Lik: -1612.687, Max-Change: 0.00208
Iteration: 22, Log-Lik: -1612.687, Max-Change: 0.02543
Iteration: 23, Log-Lik: -1612.685, Max-Change: 0.00161
Iteration: 24, Log-Lik: -1612.685, Max-Change: 0.00099
Iteration: 25, Log-Lik: -1612.685, Max-Change: 0.00044
Iteration: 26, Log-Lik: -1612.685, Max-Change: 0.00027
Iteration: 27, Log-Lik: -1612.685, Max-Change: 0.00102
Iteration: 28, Log-Lik: -1612.685, Max-Change: 0.00049
Iteration: 29, Log-Lik: -1612.685, Max-Change: 0.00020
Iteration: 30, Log-Lik: -1612.685, Max-Change: 0.00019
Iteration: 31, Log-Lik: -1612.685, Max-Change: 0.01678
Iteration: 32, Log-Lik: -1612.684, Max-Change: 0.00133
Iteration: 33, Log-Lik: -1612.684, Max-Change: 0.00078
Iteration: 34, Log-Lik: -1612.684, Max-Change: 0.00030
Iteration: 35, Log-Lik: -1612.684, Max-Change: 0.00090
Iteration: 36, Log-Lik: -1612.684, Max-Change: 0.00027
Iteration: 37, Log-Lik: -1612.684, Max-Change: 0.00014
Iteration: 38, Log-Lik: -1612.684, Max-Change: 0.00083
Iteration: 39, Log-Lik: -1612.684, Max-Change: 0.00024
Iteration: 40, Log-Lik: -1612.684, Max-Change: 0.00012
Iteration: 41, Log-Lik: -1612.684, Max-Change: 0.01511
Iteration: 42, Log-Lik: -1612.683, Max-Change: 0.00108
Iteration: 43, Log-Lik: -1612.683, Max-Change: 0.00105
Iteration: 44, Log-Lik: -1612.683, Max-Change: 0.00038
Iteration: 45, Log-Lik: -1612.683, Max-Change: 0.00111
Iteration: 46, Log-Lik: -1612.683, Max-Change: 0.00022
Iteration: 47, Log-Lik: -1612.683, Max-Change: 0.00015
Iteration: 48, Log-Lik: -1612.683, Max-Change: 0.00063
Iteration: 49, Log-Lik: -1612.683, Max-Change: 0.00021
Iteration: 50, Log-Lik: -1612.683, Max-Change: 0.00009
plinkpars <- read.mirt(mod3)
plot(plinkpars)
plot(mod3, type = 'trace')
# nominal
mod4 <- mirt(Science, 1, itemtype = 'nominal')
#>
Iteration: 1, Log-Lik: -2231.749, Max-Change: 3.43133
Iteration: 2, Log-Lik: -1647.755, Max-Change: 0.78640
Iteration: 3, Log-Lik: -1629.925, Max-Change: 0.48800
Iteration: 4, Log-Lik: -1621.546, Max-Change: 0.44557
Iteration: 5, Log-Lik: -1616.776, Max-Change: 0.30268
Iteration: 6, Log-Lik: -1613.789, Max-Change: 0.29869
Iteration: 7, Log-Lik: -1610.323, Max-Change: 0.20615
Iteration: 8, Log-Lik: -1609.697, Max-Change: 0.28886
Iteration: 9, Log-Lik: -1609.348, Max-Change: 0.13388
Iteration: 10, Log-Lik: -1609.105, Max-Change: 0.11615
Iteration: 11, Log-Lik: -1608.975, Max-Change: 0.10205
Iteration: 12, Log-Lik: -1608.881, Max-Change: 0.09345
Iteration: 13, Log-Lik: -1608.596, Max-Change: 0.06321
Iteration: 14, Log-Lik: -1608.565, Max-Change: 0.03301
Iteration: 15, Log-Lik: -1608.548, Max-Change: 0.02973
Iteration: 16, Log-Lik: -1608.507, Max-Change: 0.01790
Iteration: 17, Log-Lik: -1608.500, Max-Change: 0.02676
Iteration: 18, Log-Lik: -1608.494, Max-Change: 0.02509
Iteration: 19, Log-Lik: -1608.472, Max-Change: 0.00996
Iteration: 20, Log-Lik: -1608.468, Max-Change: 0.00382
Iteration: 21, Log-Lik: -1608.467, Max-Change: 0.00338
Iteration: 22, Log-Lik: -1608.464, Max-Change: 0.03036
Iteration: 23, Log-Lik: -1608.462, Max-Change: 0.00312
Iteration: 24, Log-Lik: -1608.462, Max-Change: 0.00244
Iteration: 25, Log-Lik: -1608.461, Max-Change: 0.00209
Iteration: 26, Log-Lik: -1608.460, Max-Change: 0.00238
Iteration: 27, Log-Lik: -1608.460, Max-Change: 0.00206
Iteration: 28, Log-Lik: -1608.460, Max-Change: 0.02332
Iteration: 29, Log-Lik: -1608.458, Max-Change: 0.00177
Iteration: 30, Log-Lik: -1608.458, Max-Change: 0.00144
Iteration: 31, Log-Lik: -1608.458, Max-Change: 0.00192
Iteration: 32, Log-Lik: -1608.458, Max-Change: 0.00042
Iteration: 33, Log-Lik: -1608.458, Max-Change: 0.00168
Iteration: 34, Log-Lik: -1608.458, Max-Change: 0.00024
Iteration: 35, Log-Lik: -1608.458, Max-Change: 0.00416
Iteration: 36, Log-Lik: -1608.458, Max-Change: 0.00063
Iteration: 37, Log-Lik: -1608.457, Max-Change: 0.00059
Iteration: 38, Log-Lik: -1608.457, Max-Change: 0.00021
Iteration: 39, Log-Lik: -1608.457, Max-Change: 0.00124
Iteration: 40, Log-Lik: -1608.457, Max-Change: 0.00057
Iteration: 41, Log-Lik: -1608.457, Max-Change: 0.00036
Iteration: 42, Log-Lik: -1608.457, Max-Change: 0.00023
Iteration: 43, Log-Lik: -1608.457, Max-Change: 0.02043
Iteration: 44, Log-Lik: -1608.456, Max-Change: 0.00064
Iteration: 45, Log-Lik: -1608.456, Max-Change: 0.00021
Iteration: 46, Log-Lik: -1608.456, Max-Change: 0.00104
Iteration: 47, Log-Lik: -1608.456, Max-Change: 0.00061
Iteration: 48, Log-Lik: -1608.456, Max-Change: 0.00021
Iteration: 49, Log-Lik: -1608.456, Max-Change: 0.00017
Iteration: 50, Log-Lik: -1608.456, Max-Change: 0.00168
Iteration: 51, Log-Lik: -1608.456, Max-Change: 0.00018
Iteration: 52, Log-Lik: -1608.456, Max-Change: 0.00017
Iteration: 53, Log-Lik: -1608.456, Max-Change: 0.00250
Iteration: 54, Log-Lik: -1608.456, Max-Change: 0.00019
Iteration: 55, Log-Lik: -1608.456, Max-Change: 0.00018
Iteration: 56, Log-Lik: -1608.456, Max-Change: 0.02338
Iteration: 57, Log-Lik: -1608.455, Max-Change: 0.00081
Iteration: 58, Log-Lik: -1608.455, Max-Change: 0.00080
Iteration: 59, Log-Lik: -1608.455, Max-Change: 0.00022
Iteration: 60, Log-Lik: -1608.455, Max-Change: 0.00084
Iteration: 61, Log-Lik: -1608.455, Max-Change: 0.00040
Iteration: 62, Log-Lik: -1608.455, Max-Change: 0.00033
Iteration: 63, Log-Lik: -1608.455, Max-Change: 0.00026
Iteration: 64, Log-Lik: -1608.455, Max-Change: 0.00126
Iteration: 65, Log-Lik: -1608.455, Max-Change: 0.00015
Iteration: 66, Log-Lik: -1608.455, Max-Change: 0.00083
Iteration: 67, Log-Lik: -1608.455, Max-Change: 0.00022
Iteration: 68, Log-Lik: -1608.455, Max-Change: 0.00018
Iteration: 69, Log-Lik: -1608.455, Max-Change: 0.00014
Iteration: 70, Log-Lik: -1608.455, Max-Change: 0.00151
Iteration: 71, Log-Lik: -1608.455, Max-Change: 0.00007
plinkpars <- read.mirt(mod4)
plot(plinkpars)
plot(mod4, type = 'trace')
## multidimensional
data <- expand.table(LSAT7)
(mod1 <- mirt(data, 2))
#>
Iteration: 1, Log-Lik: -2674.021, Max-Change: 0.20368
Iteration: 2, Log-Lik: -2658.770, Max-Change: 0.12215
Iteration: 3, Log-Lik: -2655.896, Max-Change: 0.07195
Iteration: 4, Log-Lik: -2654.819, Max-Change: 0.03372
Iteration: 5, Log-Lik: -2654.630, Max-Change: 0.02022
Iteration: 6, Log-Lik: -2654.553, Max-Change: 0.01189
Iteration: 7, Log-Lik: -2654.493, Max-Change: 0.00789
Iteration: 8, Log-Lik: -2654.475, Max-Change: 0.00763
Iteration: 9, Log-Lik: -2654.460, Max-Change: 0.00744
Iteration: 10, Log-Lik: -2654.379, Max-Change: 0.00615
Iteration: 11, Log-Lik: -2654.368, Max-Change: 0.00600
Iteration: 12, Log-Lik: -2654.357, Max-Change: 0.00589
Iteration: 13, Log-Lik: -2654.297, Max-Change: 0.00550
Iteration: 14, Log-Lik: -2654.288, Max-Change: 0.00545
Iteration: 15, Log-Lik: -2654.280, Max-Change: 0.00540
Iteration: 16, Log-Lik: -2654.232, Max-Change: 0.00523
Iteration: 17, Log-Lik: -2654.224, Max-Change: 0.00515
Iteration: 18, Log-Lik: -2654.218, Max-Change: 0.00509
Iteration: 19, Log-Lik: -2654.179, Max-Change: 0.00497
Iteration: 20, Log-Lik: -2654.173, Max-Change: 0.00492
Iteration: 21, Log-Lik: -2654.167, Max-Change: 0.00487
Iteration: 22, Log-Lik: -2654.134, Max-Change: 0.00466
Iteration: 23, Log-Lik: -2654.130, Max-Change: 0.00460
Iteration: 24, Log-Lik: -2654.125, Max-Change: 0.00455
Iteration: 25, Log-Lik: -2654.097, Max-Change: 0.00437
Iteration: 26, Log-Lik: -2654.093, Max-Change: 0.00434
Iteration: 27, Log-Lik: -2654.089, Max-Change: 0.00431
Iteration: 28, Log-Lik: -2654.066, Max-Change: 0.00420
Iteration: 29, Log-Lik: -2654.062, Max-Change: 0.00416
Iteration: 30, Log-Lik: -2654.058, Max-Change: 0.00413
Iteration: 31, Log-Lik: -2654.038, Max-Change: 0.00404
Iteration: 32, Log-Lik: -2654.034, Max-Change: 0.00401
Iteration: 33, Log-Lik: -2654.031, Max-Change: 0.00399
Iteration: 34, Log-Lik: -2654.012, Max-Change: 0.00391
Iteration: 35, Log-Lik: -2654.009, Max-Change: 0.00389
Iteration: 36, Log-Lik: -2654.006, Max-Change: 0.00387
Iteration: 37, Log-Lik: -2653.989, Max-Change: 0.00373
Iteration: 38, Log-Lik: -2653.986, Max-Change: 0.00372
Iteration: 39, Log-Lik: -2653.983, Max-Change: 0.00370
Iteration: 40, Log-Lik: -2653.966, Max-Change: 0.00366
Iteration: 41, Log-Lik: -2653.964, Max-Change: 0.00363
Iteration: 42, Log-Lik: -2653.961, Max-Change: 0.00362
Iteration: 43, Log-Lik: -2653.945, Max-Change: 0.00355
Iteration: 44, Log-Lik: -2653.943, Max-Change: 0.00353
Iteration: 45, Log-Lik: -2653.940, Max-Change: 0.00352
Iteration: 46, Log-Lik: -2653.925, Max-Change: 0.00347
Iteration: 47, Log-Lik: -2653.922, Max-Change: 0.00345
Iteration: 48, Log-Lik: -2653.920, Max-Change: 0.00344
Iteration: 49, Log-Lik: -2653.905, Max-Change: 0.00340
Iteration: 50, Log-Lik: -2653.902, Max-Change: 0.00339
Iteration: 51, Log-Lik: -2653.900, Max-Change: 0.00338
Iteration: 52, Log-Lik: -2653.884, Max-Change: 0.00335
Iteration: 53, Log-Lik: -2653.882, Max-Change: 0.00334
Iteration: 54, Log-Lik: -2653.879, Max-Change: 0.00333
Iteration: 55, Log-Lik: -2653.864, Max-Change: 0.00331
Iteration: 56, Log-Lik: -2653.862, Max-Change: 0.00330
Iteration: 57, Log-Lik: -2653.859, Max-Change: 0.00329
Iteration: 58, Log-Lik: -2653.844, Max-Change: 0.00326
Iteration: 59, Log-Lik: -2653.842, Max-Change: 0.00325
Iteration: 60, Log-Lik: -2653.839, Max-Change: 0.00324
Iteration: 61, Log-Lik: -2653.824, Max-Change: 0.00321
Iteration: 62, Log-Lik: -2653.822, Max-Change: 0.00320
Iteration: 63, Log-Lik: -2653.819, Max-Change: 0.00319
Iteration: 64, Log-Lik: -2653.804, Max-Change: 0.00315
Iteration: 65, Log-Lik: -2653.802, Max-Change: 0.00315
Iteration: 66, Log-Lik: -2653.800, Max-Change: 0.00314
Iteration: 67, Log-Lik: -2653.785, Max-Change: 0.00309
Iteration: 68, Log-Lik: -2653.783, Max-Change: 0.00309
Iteration: 69, Log-Lik: -2653.780, Max-Change: 0.00308
Iteration: 70, Log-Lik: -2653.766, Max-Change: 0.00303
Iteration: 71, Log-Lik: -2653.764, Max-Change: 0.00302
Iteration: 72, Log-Lik: -2653.762, Max-Change: 0.00301
Iteration: 73, Log-Lik: -2653.748, Max-Change: 0.00295
Iteration: 74, Log-Lik: -2653.746, Max-Change: 0.00294
Iteration: 75, Log-Lik: -2653.743, Max-Change: 0.00293
Iteration: 76, Log-Lik: -2653.730, Max-Change: 0.00287
Iteration: 77, Log-Lik: -2653.728, Max-Change: 0.00286
Iteration: 78, Log-Lik: -2653.726, Max-Change: 0.00285
Iteration: 79, Log-Lik: -2653.713, Max-Change: 0.00281
Iteration: 80, Log-Lik: -2653.711, Max-Change: 0.00281
Iteration: 81, Log-Lik: -2653.709, Max-Change: 0.00281
Iteration: 82, Log-Lik: -2653.698, Max-Change: 0.00283
Iteration: 83, Log-Lik: -2653.696, Max-Change: 0.00282
Iteration: 84, Log-Lik: -2653.694, Max-Change: 0.00282
Iteration: 85, Log-Lik: -2653.683, Max-Change: 0.00282
Iteration: 86, Log-Lik: -2653.681, Max-Change: 0.00282
Iteration: 87, Log-Lik: -2653.679, Max-Change: 0.00281
Iteration: 88, Log-Lik: -2653.669, Max-Change: 0.00281
Iteration: 89, Log-Lik: -2653.667, Max-Change: 0.00280
Iteration: 90, Log-Lik: -2653.665, Max-Change: 0.00279
Iteration: 91, Log-Lik: -2653.655, Max-Change: 0.00278
Iteration: 92, Log-Lik: -2653.654, Max-Change: 0.00277
Iteration: 93, Log-Lik: -2653.652, Max-Change: 0.00276
Iteration: 94, Log-Lik: -2653.643, Max-Change: 0.00274
Iteration: 95, Log-Lik: -2653.642, Max-Change: 0.00273
Iteration: 96, Log-Lik: -2653.640, Max-Change: 0.00272
Iteration: 97, Log-Lik: -2653.632, Max-Change: 0.00269
Iteration: 98, Log-Lik: -2653.631, Max-Change: 0.00268
Iteration: 99, Log-Lik: -2653.630, Max-Change: 0.00267
Iteration: 100, Log-Lik: -2653.622, Max-Change: 0.00263
Iteration: 101, Log-Lik: -2653.621, Max-Change: 0.00262
Iteration: 102, Log-Lik: -2653.619, Max-Change: 0.00261
Iteration: 103, Log-Lik: -2653.612, Max-Change: 0.00256
Iteration: 104, Log-Lik: -2653.611, Max-Change: 0.00255
Iteration: 105, Log-Lik: -2653.610, Max-Change: 0.00254
Iteration: 106, Log-Lik: -2653.604, Max-Change: 0.00249
Iteration: 107, Log-Lik: -2653.603, Max-Change: 0.00248
Iteration: 108, Log-Lik: -2653.602, Max-Change: 0.00247
Iteration: 109, Log-Lik: -2653.596, Max-Change: 0.00242
Iteration: 110, Log-Lik: -2653.595, Max-Change: 0.00241
Iteration: 111, Log-Lik: -2653.594, Max-Change: 0.00239
Iteration: 112, Log-Lik: -2653.589, Max-Change: 0.00234
Iteration: 113, Log-Lik: -2653.588, Max-Change: 0.00233
Iteration: 114, Log-Lik: -2653.587, Max-Change: 0.00232
Iteration: 115, Log-Lik: -2653.582, Max-Change: 0.00226
Iteration: 116, Log-Lik: -2653.581, Max-Change: 0.00225
Iteration: 117, Log-Lik: -2653.580, Max-Change: 0.00224
Iteration: 118, Log-Lik: -2653.576, Max-Change: 0.00218
Iteration: 119, Log-Lik: -2653.575, Max-Change: 0.00217
Iteration: 120, Log-Lik: -2653.575, Max-Change: 0.00216
Iteration: 121, Log-Lik: -2653.571, Max-Change: 0.00210
Iteration: 122, Log-Lik: -2653.570, Max-Change: 0.00209
Iteration: 123, Log-Lik: -2653.569, Max-Change: 0.00208
Iteration: 124, Log-Lik: -2653.566, Max-Change: 0.00202
Iteration: 125, Log-Lik: -2653.565, Max-Change: 0.00201
Iteration: 126, Log-Lik: -2653.565, Max-Change: 0.00200
Iteration: 127, Log-Lik: -2653.561, Max-Change: 0.00194
Iteration: 128, Log-Lik: -2653.561, Max-Change: 0.00193
Iteration: 129, Log-Lik: -2653.560, Max-Change: 0.00192
Iteration: 130, Log-Lik: -2653.557, Max-Change: 0.00186
Iteration: 131, Log-Lik: -2653.557, Max-Change: 0.00185
Iteration: 132, Log-Lik: -2653.556, Max-Change: 0.00184
Iteration: 133, Log-Lik: -2653.554, Max-Change: 0.00179
Iteration: 134, Log-Lik: -2653.553, Max-Change: 0.00178
Iteration: 135, Log-Lik: -2653.553, Max-Change: 0.00177
Iteration: 136, Log-Lik: -2653.550, Max-Change: 0.00171
Iteration: 137, Log-Lik: -2653.550, Max-Change: 0.00170
Iteration: 138, Log-Lik: -2653.550, Max-Change: 0.00169
Iteration: 139, Log-Lik: -2653.547, Max-Change: 0.00164
Iteration: 140, Log-Lik: -2653.547, Max-Change: 0.00163
Iteration: 141, Log-Lik: -2653.547, Max-Change: 0.00162
Iteration: 142, Log-Lik: -2653.545, Max-Change: 0.00157
Iteration: 143, Log-Lik: -2653.544, Max-Change: 0.00148
Iteration: 144, Log-Lik: -2653.544, Max-Change: 0.00148
Iteration: 145, Log-Lik: -2653.543, Max-Change: 0.00196
Iteration: 146, Log-Lik: -2653.543, Max-Change: 0.00138
Iteration: 147, Log-Lik: -2653.542, Max-Change: 0.00137
Iteration: 148, Log-Lik: -2653.541, Max-Change: 0.00143
Iteration: 149, Log-Lik: -2653.541, Max-Change: 0.00137
Iteration: 150, Log-Lik: -2653.541, Max-Change: 0.00138
Iteration: 151, Log-Lik: -2653.539, Max-Change: 0.00134
Iteration: 152, Log-Lik: -2653.539, Max-Change: 0.00137
Iteration: 153, Log-Lik: -2653.539, Max-Change: 0.00131
Iteration: 154, Log-Lik: -2653.538, Max-Change: 0.00113
Iteration: 155, Log-Lik: -2653.538, Max-Change: 0.00130
Iteration: 156, Log-Lik: -2653.537, Max-Change: 0.00130
Iteration: 157, Log-Lik: -2653.537, Max-Change: 0.00121
Iteration: 158, Log-Lik: -2653.536, Max-Change: 0.00109
Iteration: 159, Log-Lik: -2653.536, Max-Change: 0.00126
Iteration: 160, Log-Lik: -2653.536, Max-Change: 0.00102
Iteration: 161, Log-Lik: -2653.536, Max-Change: 0.00119
Iteration: 162, Log-Lik: -2653.535, Max-Change: 0.00096
Iteration: 163, Log-Lik: -2653.535, Max-Change: 0.00213
Iteration: 164, Log-Lik: -2653.535, Max-Change: 0.00194
Iteration: 165, Log-Lik: -2653.534, Max-Change: 0.00196
Iteration: 166, Log-Lik: -2653.533, Max-Change: 0.00105
Iteration: 167, Log-Lik: -2653.532, Max-Change: 0.00118
Iteration: 168, Log-Lik: -2653.532, Max-Change: 0.00111
Iteration: 169, Log-Lik: -2653.532, Max-Change: 0.00094
Iteration: 170, Log-Lik: -2653.531, Max-Change: 0.00093
Iteration: 171, Log-Lik: -2653.531, Max-Change: 0.00093
Iteration: 172, Log-Lik: -2653.531, Max-Change: 0.00165
Iteration: 173, Log-Lik: -2653.530, Max-Change: 0.00157
Iteration: 174, Log-Lik: -2653.530, Max-Change: 0.00181
Iteration: 175, Log-Lik: -2653.530, Max-Change: 0.00098
Iteration: 176, Log-Lik: -2653.529, Max-Change: 0.00161
Iteration: 177, Log-Lik: -2653.529, Max-Change: 0.00163
Iteration: 178, Log-Lik: -2653.529, Max-Change: 0.00106
Iteration: 179, Log-Lik: -2653.529, Max-Change: 0.00074
Iteration: 180, Log-Lik: -2653.528, Max-Change: 0.00177
Iteration: 181, Log-Lik: -2653.528, Max-Change: 0.00143
Iteration: 182, Log-Lik: -2653.528, Max-Change: 0.00145
Iteration: 183, Log-Lik: -2653.528, Max-Change: 0.00150
Iteration: 184, Log-Lik: -2653.527, Max-Change: 0.00088
Iteration: 185, Log-Lik: -2653.527, Max-Change: 0.00035
Iteration: 186, Log-Lik: -2653.527, Max-Change: 0.00029
Iteration: 187, Log-Lik: -2653.527, Max-Change: 0.00031
Iteration: 188, Log-Lik: -2653.527, Max-Change: 0.00032
Iteration: 189, Log-Lik: -2653.527, Max-Change: 0.00033
Iteration: 190, Log-Lik: -2653.527, Max-Change: 0.00034
Iteration: 191, Log-Lik: -2653.527, Max-Change: 0.00034
Iteration: 192, Log-Lik: -2653.527, Max-Change: 0.00034
Iteration: 193, Log-Lik: -2653.527, Max-Change: 0.00034
Iteration: 194, Log-Lik: -2653.527, Max-Change: 0.00034
Iteration: 195, Log-Lik: -2653.527, Max-Change: 0.00034
Iteration: 196, Log-Lik: -2653.527, Max-Change: 0.00033
Iteration: 197, Log-Lik: -2653.527, Max-Change: 0.00033
Iteration: 198, Log-Lik: -2653.527, Max-Change: 0.00033
Iteration: 199, Log-Lik: -2653.526, Max-Change: 0.00033
Iteration: 200, Log-Lik: -2653.526, Max-Change: 0.00032
Iteration: 201, Log-Lik: -2653.526, Max-Change: 0.00032
Iteration: 202, Log-Lik: -2653.526, Max-Change: 0.00032
Iteration: 203, Log-Lik: -2653.526, Max-Change: 0.00032
Iteration: 204, Log-Lik: -2653.526, Max-Change: 0.00032
Iteration: 205, Log-Lik: -2653.526, Max-Change: 0.00032
Iteration: 206, Log-Lik: -2653.526, Max-Change: 0.00031
Iteration: 207, Log-Lik: -2653.526, Max-Change: 0.00031
Iteration: 208, Log-Lik: -2653.526, Max-Change: 0.00031
Iteration: 209, Log-Lik: -2653.526, Max-Change: 0.00031
Iteration: 210, Log-Lik: -2653.526, Max-Change: 0.00030
Iteration: 211, Log-Lik: -2653.526, Max-Change: 0.00030
Iteration: 212, Log-Lik: -2653.526, Max-Change: 0.00030
Iteration: 213, Log-Lik: -2653.526, Max-Change: 0.00030
Iteration: 214, Log-Lik: -2653.526, Max-Change: 0.00030
Iteration: 215, Log-Lik: -2653.526, Max-Change: 0.00030
Iteration: 216, Log-Lik: -2653.526, Max-Change: 0.00029
Iteration: 217, Log-Lik: -2653.526, Max-Change: 0.00029
Iteration: 218, Log-Lik: -2653.526, Max-Change: 0.00029
Iteration: 219, Log-Lik: -2653.526, Max-Change: 0.00029
Iteration: 220, Log-Lik: -2653.526, Max-Change: 0.00029
Iteration: 221, Log-Lik: -2653.526, Max-Change: 0.00029
Iteration: 222, Log-Lik: -2653.526, Max-Change: 0.00028
Iteration: 223, Log-Lik: -2653.526, Max-Change: 0.00028
Iteration: 224, Log-Lik: -2653.526, Max-Change: 0.00028
Iteration: 225, Log-Lik: -2653.525, Max-Change: 0.00028
Iteration: 226, Log-Lik: -2653.525, Max-Change: 0.00028
Iteration: 227, Log-Lik: -2653.525, Max-Change: 0.00028
Iteration: 228, Log-Lik: -2653.525, Max-Change: 0.00027
Iteration: 229, Log-Lik: -2653.525, Max-Change: 0.00027
Iteration: 230, Log-Lik: -2653.525, Max-Change: 0.00027
Iteration: 231, Log-Lik: -2653.525, Max-Change: 0.00027
Iteration: 232, Log-Lik: -2653.525, Max-Change: 0.00027
Iteration: 233, Log-Lik: -2653.525, Max-Change: 0.00027
Iteration: 234, Log-Lik: -2653.525, Max-Change: 0.00133
Iteration: 235, Log-Lik: -2653.525, Max-Change: 0.00148
Iteration: 236, Log-Lik: -2653.525, Max-Change: 0.00026
Iteration: 237, Log-Lik: -2653.525, Max-Change: 0.00023
Iteration: 238, Log-Lik: -2653.525, Max-Change: 0.00024
Iteration: 239, Log-Lik: -2653.525, Max-Change: 0.00122
Iteration: 240, Log-Lik: -2653.525, Max-Change: 0.00041
Iteration: 241, Log-Lik: -2653.525, Max-Change: 0.00026
Iteration: 242, Log-Lik: -2653.525, Max-Change: 0.00025
Iteration: 243, Log-Lik: -2653.525, Max-Change: 0.00124
Iteration: 244, Log-Lik: -2653.525, Max-Change: 0.00134
Iteration: 245, Log-Lik: -2653.524, Max-Change: 0.00033
Iteration: 246, Log-Lik: -2653.524, Max-Change: 0.00021
Iteration: 247, Log-Lik: -2653.524, Max-Change: 0.00115
Iteration: 248, Log-Lik: -2653.524, Max-Change: 0.00036
Iteration: 249, Log-Lik: -2653.524, Max-Change: 0.00024
Iteration: 250, Log-Lik: -2653.524, Max-Change: 0.00116
Iteration: 251, Log-Lik: -2653.524, Max-Change: 0.00022
Iteration: 252, Log-Lik: -2653.524, Max-Change: 0.00022
Iteration: 253, Log-Lik: -2653.524, Max-Change: 0.00022
Iteration: 254, Log-Lik: -2653.524, Max-Change: 0.00112
Iteration: 255, Log-Lik: -2653.524, Max-Change: 0.00072
Iteration: 256, Log-Lik: -2653.524, Max-Change: 0.00022
Iteration: 257, Log-Lik: -2653.524, Max-Change: 0.00022
Iteration: 258, Log-Lik: -2653.524, Max-Change: 0.00108
Iteration: 259, Log-Lik: -2653.524, Max-Change: 0.00090
Iteration: 260, Log-Lik: -2653.524, Max-Change: 0.00063
Iteration: 261, Log-Lik: -2653.524, Max-Change: 0.00043
Iteration: 262, Log-Lik: -2653.524, Max-Change: 0.00106
Iteration: 263, Log-Lik: -2653.524, Max-Change: 0.00031
Iteration: 264, Log-Lik: -2653.524, Max-Change: 0.00021
Iteration: 265, Log-Lik: -2653.524, Max-Change: 0.00103
Iteration: 266, Log-Lik: -2653.523, Max-Change: 0.00056
Iteration: 267, Log-Lik: -2653.523, Max-Change: 0.00039
Iteration: 268, Log-Lik: -2653.523, Max-Change: 0.00101
Iteration: 269, Log-Lik: -2653.523, Max-Change: 0.00020
Iteration: 270, Log-Lik: -2653.523, Max-Change: 0.00020
Iteration: 271, Log-Lik: -2653.523, Max-Change: 0.00020
Iteration: 272, Log-Lik: -2653.523, Max-Change: 0.00099
Iteration: 273, Log-Lik: -2653.523, Max-Change: 0.00067
Iteration: 274, Log-Lik: -2653.523, Max-Change: 0.00020
Iteration: 275, Log-Lik: -2653.523, Max-Change: 0.00019
Iteration: 276, Log-Lik: -2653.523, Max-Change: 0.00097
Iteration: 277, Log-Lik: -2653.523, Max-Change: 0.00084
Iteration: 278, Log-Lik: -2653.523, Max-Change: 0.00058
Iteration: 279, Log-Lik: -2653.523, Max-Change: 0.00040
Iteration: 280, Log-Lik: -2653.523, Max-Change: 0.00095
Iteration: 281, Log-Lik: -2653.523, Max-Change: 0.00028
Iteration: 282, Log-Lik: -2653.523, Max-Change: 0.00020
Iteration: 283, Log-Lik: -2653.523, Max-Change: 0.00093
Iteration: 284, Log-Lik: -2653.523, Max-Change: 0.00054
Iteration: 285, Log-Lik: -2653.523, Max-Change: 0.00037
Iteration: 286, Log-Lik: -2653.523, Max-Change: 0.00092
Iteration: 287, Log-Lik: -2653.523, Max-Change: 0.00020
Iteration: 288, Log-Lik: -2653.523, Max-Change: 0.00018
Iteration: 289, Log-Lik: -2653.523, Max-Change: 0.00090
Iteration: 290, Log-Lik: -2653.523, Max-Change: 0.00078
Iteration: 291, Log-Lik: -2653.523, Max-Change: 0.00055
Iteration: 292, Log-Lik: -2653.523, Max-Change: 0.00088
Iteration: 293, Log-Lik: -2653.523, Max-Change: 0.00079
Iteration: 294, Log-Lik: -2653.523, Max-Change: 0.00055
Iteration: 295, Log-Lik: -2653.522, Max-Change: 0.00087
Iteration: 296, Log-Lik: -2653.522, Max-Change: 0.00081
Iteration: 297, Log-Lik: -2653.522, Max-Change: 0.00057
Iteration: 298, Log-Lik: -2653.522, Max-Change: 0.00017
Iteration: 299, Log-Lik: -2653.522, Max-Change: 0.00086
Iteration: 300, Log-Lik: -2653.522, Max-Change: 0.00061
Iteration: 301, Log-Lik: -2653.522, Max-Change: 0.00017
Iteration: 302, Log-Lik: -2653.522, Max-Change: 0.00017
Iteration: 303, Log-Lik: -2653.522, Max-Change: 0.00084
Iteration: 304, Log-Lik: -2653.522, Max-Change: 0.00075
Iteration: 305, Log-Lik: -2653.522, Max-Change: 0.00052
Iteration: 306, Log-Lik: -2653.522, Max-Change: 0.00036
Iteration: 307, Log-Lik: -2653.522, Max-Change: 0.00083
Iteration: 308, Log-Lik: -2653.522, Max-Change: 0.00026
Iteration: 309, Log-Lik: -2653.522, Max-Change: 0.00017
Iteration: 310, Log-Lik: -2653.522, Max-Change: 0.00081
Iteration: 311, Log-Lik: -2653.522, Max-Change: 0.00049
Iteration: 312, Log-Lik: -2653.522, Max-Change: 0.00034
Iteration: 313, Log-Lik: -2653.522, Max-Change: 0.00080
Iteration: 314, Log-Lik: -2653.522, Max-Change: 0.00020
Iteration: 315, Log-Lik: -2653.522, Max-Change: 0.00016
Iteration: 316, Log-Lik: -2653.522, Max-Change: 0.00079
Iteration: 317, Log-Lik: -2653.522, Max-Change: 0.00065
Iteration: 318, Log-Lik: -2653.522, Max-Change: 0.00045
Iteration: 319, Log-Lik: -2653.522, Max-Change: 0.00078
Iteration: 320, Log-Lik: -2653.522, Max-Change: 0.00059
Iteration: 321, Log-Lik: -2653.522, Max-Change: 0.00041
Iteration: 322, Log-Lik: -2653.522, Max-Change: 0.00077
Iteration: 323, Log-Lik: -2653.522, Max-Change: 0.00047
Iteration: 324, Log-Lik: -2653.522, Max-Change: 0.00033
Iteration: 325, Log-Lik: -2653.522, Max-Change: 0.00076
Iteration: 326, Log-Lik: -2653.522, Max-Change: 0.00021
Iteration: 327, Log-Lik: -2653.522, Max-Change: 0.00015
Iteration: 328, Log-Lik: -2653.522, Max-Change: 0.00074
Iteration: 329, Log-Lik: -2653.522, Max-Change: 0.00056
Iteration: 330, Log-Lik: -2653.522, Max-Change: 0.00039
Iteration: 331, Log-Lik: -2653.522, Max-Change: 0.00074
Iteration: 332, Log-Lik: -2653.521, Max-Change: 0.00043
Iteration: 333, Log-Lik: -2653.521, Max-Change: 0.00030
Iteration: 334, Log-Lik: -2653.521, Max-Change: 0.00073
Iteration: 335, Log-Lik: -2653.521, Max-Change: 0.00015
Iteration: 336, Log-Lik: -2653.521, Max-Change: 0.00072
Iteration: 337, Log-Lik: -2653.521, Max-Change: 0.00086
Iteration: 338, Log-Lik: -2653.521, Max-Change: 0.00060
Iteration: 339, Log-Lik: -2653.521, Max-Change: 0.00041
Iteration: 340, Log-Lik: -2653.521, Max-Change: 0.00071
Iteration: 341, Log-Lik: -2653.521, Max-Change: 0.00054
Iteration: 342, Log-Lik: -2653.521, Max-Change: 0.00037
Iteration: 343, Log-Lik: -2653.521, Max-Change: 0.00070
Iteration: 344, Log-Lik: -2653.521, Max-Change: 0.00042
Iteration: 345, Log-Lik: -2653.521, Max-Change: 0.00029
Iteration: 346, Log-Lik: -2653.521, Max-Change: 0.00069
Iteration: 347, Log-Lik: -2653.521, Max-Change: 0.00014
Iteration: 348, Log-Lik: -2653.521, Max-Change: 0.00014
Iteration: 349, Log-Lik: -2653.521, Max-Change: 0.00068
Iteration: 350, Log-Lik: -2653.521, Max-Change: 0.00069
Iteration: 351, Log-Lik: -2653.521, Max-Change: 0.00048
Iteration: 352, Log-Lik: -2653.521, Max-Change: 0.00013
Iteration: 353, Log-Lik: -2653.521, Max-Change: 0.00068
Iteration: 354, Log-Lik: -2653.521, Max-Change: 0.00054
Iteration: 355, Log-Lik: -2653.521, Max-Change: 0.00014
Iteration: 356, Log-Lik: -2653.521, Max-Change: 0.00013
Iteration: 357, Log-Lik: -2653.521, Max-Change: 0.00066
Iteration: 358, Log-Lik: -2653.521, Max-Change: 0.00061
Iteration: 359, Log-Lik: -2653.521, Max-Change: 0.00042
Iteration: 360, Log-Lik: -2653.521, Max-Change: 0.00029
Iteration: 361, Log-Lik: -2653.521, Max-Change: 0.00065
Iteration: 362, Log-Lik: -2653.521, Max-Change: 0.00020
Iteration: 363, Log-Lik: -2653.521, Max-Change: 0.00014
Iteration: 364, Log-Lik: -2653.521, Max-Change: 0.00064
Iteration: 365, Log-Lik: -2653.521, Max-Change: 0.00041
Iteration: 366, Log-Lik: -2653.521, Max-Change: 0.00029
Iteration: 367, Log-Lik: -2653.521, Max-Change: 0.00064
Iteration: 368, Log-Lik: -2653.521, Max-Change: 0.00020
Iteration: 369, Log-Lik: -2653.521, Max-Change: 0.00014
Iteration: 370, Log-Lik: -2653.521, Max-Change: 0.00063
Iteration: 371, Log-Lik: -2653.521, Max-Change: 0.00041
Iteration: 372, Log-Lik: -2653.521, Max-Change: 0.00028
Iteration: 373, Log-Lik: -2653.521, Max-Change: 0.00062
Iteration: 374, Log-Lik: -2653.521, Max-Change: 0.00020
Iteration: 375, Log-Lik: -2653.521, Max-Change: 0.00013
Iteration: 376, Log-Lik: -2653.521, Max-Change: 0.00061
Iteration: 377, Log-Lik: -2653.521, Max-Change: 0.00039
Iteration: 378, Log-Lik: -2653.521, Max-Change: 0.00027
Iteration: 379, Log-Lik: -2653.521, Max-Change: 0.00061
Iteration: 380, Log-Lik: -2653.520, Max-Change: 0.00018
Iteration: 381, Log-Lik: -2653.520, Max-Change: 0.00013
Iteration: 382, Log-Lik: -2653.520, Max-Change: 0.00060
Iteration: 383, Log-Lik: -2653.520, Max-Change: 0.00042
Iteration: 384, Log-Lik: -2653.520, Max-Change: 0.00029
Iteration: 385, Log-Lik: -2653.520, Max-Change: 0.00060
Iteration: 386, Log-Lik: -2653.520, Max-Change: 0.00027
Iteration: 387, Log-Lik: -2653.520, Max-Change: 0.00019
Iteration: 388, Log-Lik: -2653.520, Max-Change: 0.00059
Iteration: 389, Log-Lik: -2653.520, Max-Change: 0.00012
Iteration: 390, Log-Lik: -2653.520, Max-Change: 0.00058
Iteration: 391, Log-Lik: -2653.520, Max-Change: 0.00069
Iteration: 392, Log-Lik: -2653.520, Max-Change: 0.00048
Iteration: 393, Log-Lik: -2653.520, Max-Change: 0.00034
Iteration: 394, Log-Lik: -2653.520, Max-Change: 0.00058
Iteration: 395, Log-Lik: -2653.520, Max-Change: 0.00043
Iteration: 396, Log-Lik: -2653.520, Max-Change: 0.00030
Iteration: 397, Log-Lik: -2653.520, Max-Change: 0.00057
Iteration: 398, Log-Lik: -2653.520, Max-Change: 0.00032
Iteration: 399, Log-Lik: -2653.520, Max-Change: 0.00022
Iteration: 400, Log-Lik: -2653.520, Max-Change: 0.00056
Iteration: 401, Log-Lik: -2653.520, Max-Change: 0.00056
Iteration: 402, Log-Lik: -2653.520, Max-Change: 0.00042
Iteration: 403, Log-Lik: -2653.520, Max-Change: 0.00011
Iteration: 404, Log-Lik: -2653.520, Max-Change: 0.00055
Iteration: 405, Log-Lik: -2653.520, Max-Change: 0.00060
Iteration: 406, Log-Lik: -2653.520, Max-Change: 0.00011
Iteration: 407, Log-Lik: -2653.520, Max-Change: 0.00054
Iteration: 408, Log-Lik: -2653.520, Max-Change: 0.00055
Iteration: 409, Log-Lik: -2653.520, Max-Change: 0.00011
Iteration: 410, Log-Lik: -2653.520, Max-Change: 0.00054
Iteration: 411, Log-Lik: -2653.520, Max-Change: 0.00056
Iteration: 412, Log-Lik: -2653.520, Max-Change: 0.00011
Iteration: 413, Log-Lik: -2653.520, Max-Change: 0.00053
Iteration: 414, Log-Lik: -2653.520, Max-Change: 0.00055
Iteration: 415, Log-Lik: -2653.520, Max-Change: 0.00011
Iteration: 416, Log-Lik: -2653.520, Max-Change: 0.00053
Iteration: 417, Log-Lik: -2653.520, Max-Change: 0.00054
Iteration: 418, Log-Lik: -2653.520, Max-Change: 0.00011
Iteration: 419, Log-Lik: -2653.520, Max-Change: 0.00052
Iteration: 420, Log-Lik: -2653.520, Max-Change: 0.00054
Iteration: 421, Log-Lik: -2653.520, Max-Change: 0.00010
Iteration: 422, Log-Lik: -2653.520, Max-Change: 0.00052
Iteration: 423, Log-Lik: -2653.520, Max-Change: 0.00053
Iteration: 424, Log-Lik: -2653.520, Max-Change: 0.00010
Iteration: 425, Log-Lik: -2653.520, Max-Change: 0.00051
Iteration: 426, Log-Lik: -2653.520, Max-Change: 0.00053
Iteration: 427, Log-Lik: -2653.520, Max-Change: 0.00010
Iteration: 428, Log-Lik: -2653.520, Max-Change: 0.00050
Iteration: 429, Log-Lik: -2653.520, Max-Change: 0.00052
Iteration: 430, Log-Lik: -2653.520, Max-Change: 0.00010
Iteration: 431, Log-Lik: -2653.520, Max-Change: 0.00050
Iteration: 432, Log-Lik: -2653.520, Max-Change: 0.00051
Iteration: 433, Log-Lik: -2653.520, Max-Change: 0.00010
Iteration: 434, Log-Lik: -2653.520, Max-Change: 0.00049
Iteration: 435, Log-Lik: -2653.520, Max-Change: 0.00051
Iteration: 436, Log-Lik: -2653.520, Max-Change: 0.00010
#>
#> Call:
#> mirt(data = data, model = 2)
#>
#> Full-information item factor analysis with 2 factor(s).
#> Converged within 1e-04 tolerance after 436 EM iterations.
#> mirt version: 1.44.3
#> M-step optimizer: BFGS
#> EM acceleration: Ramsay
#> Number of rectangular quadrature: 31
#> Latent density type: Gaussian
#>
#> Log-likelihood = -2653.52
#> Estimated parameters: 14
#> AIC = 5335.039
#> BIC = 5403.748; SABIC = 5359.283
#> G2 (17) = 21.13, p = 0.2205
#> RMSEA = 0.016, CFI = NaN, TLI = NaN
plinkpars <- read.mirt(mod1)
plinkpars
#> An object of class "irt.pars"
#> Slot "pars":
#> a1 a2
#> [1,] -2.0071991 0.8703411 2.6479286 0
#> [2,] -0.8488431 -0.5221152 0.7876390 0
#> [3,] -2.1529666 -1.8364974 2.4829603 0
#> [4,] -0.7559009 -0.0280359 0.4847286 0
#> [5,] -0.7572746 0.0000000 1.8640966 0
#>
#> Slot "cat":
#> [1] 2 2 2 2 2
#>
#> Slot "poly.mod":
#> An object of class "poly.mod"
#> Slot "model":
#> [1] "drm"
#>
#> Slot "items":
#> $drm
#> [1] 1 2 3 4 5
#>
#>
#>
#> Slot "common":
#> NULL
#>
#> Slot "location":
#> [1] FALSE
#>
#> Slot "groups":
#> [1] 1
#>
#> Slot "dimensions":
#> [1] 2
#>
plot(plinkpars)
plot(mod1, type = 'trace')
cmod <- mirt.model('
F1 = 1,4,5
F2 = 2-4')
model <- mirt(data, cmod)
#>
Iteration: 1, Log-Lik: -2688.760, Max-Change: 0.14173
Iteration: 2, Log-Lik: -2684.008, Max-Change: 0.12454
Iteration: 3, Log-Lik: -2681.310, Max-Change: 0.10331
Iteration: 4, Log-Lik: -2676.246, Max-Change: 0.03288
Iteration: 5, Log-Lik: -2675.922, Max-Change: 0.02425
Iteration: 6, Log-Lik: -2675.741, Max-Change: 0.02320
Iteration: 7, Log-Lik: -2675.347, Max-Change: 0.01978
Iteration: 8, Log-Lik: -2675.302, Max-Change: 0.01833
Iteration: 9, Log-Lik: -2675.265, Max-Change: 0.01694
Iteration: 10, Log-Lik: -2675.110, Max-Change: 0.01201
Iteration: 11, Log-Lik: -2675.093, Max-Change: 0.01168
Iteration: 12, Log-Lik: -2675.078, Max-Change: 0.01132
Iteration: 13, Log-Lik: -2675.008, Max-Change: 0.00942
Iteration: 14, Log-Lik: -2674.999, Max-Change: 0.00918
Iteration: 15, Log-Lik: -2674.992, Max-Change: 0.00860
Iteration: 16, Log-Lik: -2674.955, Max-Change: 0.00705
Iteration: 17, Log-Lik: -2674.950, Max-Change: 0.00691
Iteration: 18, Log-Lik: -2674.946, Max-Change: 0.00589
Iteration: 19, Log-Lik: -2674.932, Max-Change: 0.00597
Iteration: 20, Log-Lik: -2674.928, Max-Change: 0.00513
Iteration: 21, Log-Lik: -2674.926, Max-Change: 0.00512
Iteration: 22, Log-Lik: -2674.911, Max-Change: 0.00501
Iteration: 23, Log-Lik: -2674.909, Max-Change: 0.00425
Iteration: 24, Log-Lik: -2674.908, Max-Change: 0.00425
Iteration: 25, Log-Lik: -2674.898, Max-Change: 0.00385
Iteration: 26, Log-Lik: -2674.897, Max-Change: 0.00347
Iteration: 27, Log-Lik: -2674.896, Max-Change: 0.00353
Iteration: 28, Log-Lik: -2674.890, Max-Change: 0.00313
Iteration: 29, Log-Lik: -2674.889, Max-Change: 0.00304
Iteration: 30, Log-Lik: -2674.888, Max-Change: 0.00304
Iteration: 31, Log-Lik: -2674.884, Max-Change: 0.00240
Iteration: 32, Log-Lik: -2674.883, Max-Change: 0.00240
Iteration: 33, Log-Lik: -2674.883, Max-Change: 0.00209
Iteration: 34, Log-Lik: -2674.881, Max-Change: 0.00254
Iteration: 35, Log-Lik: -2674.881, Max-Change: 0.00204
Iteration: 36, Log-Lik: -2674.880, Max-Change: 0.00200
Iteration: 37, Log-Lik: -2674.879, Max-Change: 0.00227
Iteration: 38, Log-Lik: -2674.879, Max-Change: 0.00192
Iteration: 39, Log-Lik: -2674.879, Max-Change: 0.00186
Iteration: 40, Log-Lik: -2674.878, Max-Change: 0.00214
Iteration: 41, Log-Lik: -2674.877, Max-Change: 0.00178
Iteration: 42, Log-Lik: -2674.877, Max-Change: 0.00174
Iteration: 43, Log-Lik: -2674.877, Max-Change: 0.00184
Iteration: 44, Log-Lik: -2674.876, Max-Change: 0.00175
Iteration: 45, Log-Lik: -2674.876, Max-Change: 0.00168
Iteration: 46, Log-Lik: -2674.875, Max-Change: 0.00177
Iteration: 47, Log-Lik: -2674.875, Max-Change: 0.00293
Iteration: 48, Log-Lik: -2674.875, Max-Change: 0.00304
Iteration: 49, Log-Lik: -2674.873, Max-Change: 0.00214
Iteration: 50, Log-Lik: -2674.873, Max-Change: 0.00151
Iteration: 51, Log-Lik: -2674.873, Max-Change: 0.00132
Iteration: 52, Log-Lik: -2674.872, Max-Change: 0.00124
Iteration: 53, Log-Lik: -2674.872, Max-Change: 0.00124
Iteration: 54, Log-Lik: -2674.872, Max-Change: 0.00124
Iteration: 55, Log-Lik: -2674.871, Max-Change: 0.00119
Iteration: 56, Log-Lik: -2674.871, Max-Change: 0.00114
Iteration: 57, Log-Lik: -2674.871, Max-Change: 0.00112
Iteration: 58, Log-Lik: -2674.870, Max-Change: 0.00047
Iteration: 59, Log-Lik: -2674.870, Max-Change: 0.00039
Iteration: 60, Log-Lik: -2674.870, Max-Change: 0.00036
Iteration: 61, Log-Lik: -2674.870, Max-Change: 0.00037
Iteration: 62, Log-Lik: -2674.870, Max-Change: 0.00174
Iteration: 63, Log-Lik: -2674.870, Max-Change: 0.00217
Iteration: 64, Log-Lik: -2674.870, Max-Change: 0.00075
Iteration: 65, Log-Lik: -2674.870, Max-Change: 0.00045
Iteration: 66, Log-Lik: -2674.870, Max-Change: 0.00034
Iteration: 67, Log-Lik: -2674.870, Max-Change: 0.00032
Iteration: 68, Log-Lik: -2674.870, Max-Change: 0.00160
Iteration: 69, Log-Lik: -2674.870, Max-Change: 0.00167
Iteration: 70, Log-Lik: -2674.869, Max-Change: 0.00146
Iteration: 71, Log-Lik: -2674.869, Max-Change: 0.00038
Iteration: 72, Log-Lik: -2674.869, Max-Change: 0.00027
Iteration: 73, Log-Lik: -2674.869, Max-Change: 0.00029
Iteration: 74, Log-Lik: -2674.869, Max-Change: 0.00138
Iteration: 75, Log-Lik: -2674.869, Max-Change: 0.00030
Iteration: 76, Log-Lik: -2674.869, Max-Change: 0.00144
Iteration: 77, Log-Lik: -2674.869, Max-Change: 0.00071
Iteration: 78, Log-Lik: -2674.869, Max-Change: 0.00041
Iteration: 79, Log-Lik: -2674.869, Max-Change: 0.00028
Iteration: 80, Log-Lik: -2674.869, Max-Change: 0.00131
Iteration: 81, Log-Lik: -2674.869, Max-Change: 0.00040
Iteration: 82, Log-Lik: -2674.869, Max-Change: 0.00028
Iteration: 83, Log-Lik: -2674.869, Max-Change: 0.00125
Iteration: 84, Log-Lik: -2674.869, Max-Change: 0.00077
Iteration: 85, Log-Lik: -2674.869, Max-Change: 0.00028
Iteration: 86, Log-Lik: -2674.869, Max-Change: 0.00024
Iteration: 87, Log-Lik: -2674.869, Max-Change: 0.00131
Iteration: 88, Log-Lik: -2674.869, Max-Change: 0.00069
Iteration: 89, Log-Lik: -2674.869, Max-Change: 0.00042
Iteration: 90, Log-Lik: -2674.869, Max-Change: 0.00026
Iteration: 91, Log-Lik: -2674.869, Max-Change: 0.00121
Iteration: 92, Log-Lik: -2674.868, Max-Change: 0.00025
Iteration: 93, Log-Lik: -2674.868, Max-Change: 0.00116
Iteration: 94, Log-Lik: -2674.868, Max-Change: 0.00058
Iteration: 95, Log-Lik: -2674.868, Max-Change: 0.00036
Iteration: 96, Log-Lik: -2674.868, Max-Change: 0.00025
Iteration: 97, Log-Lik: -2674.868, Max-Change: 0.00115
Iteration: 98, Log-Lik: -2674.868, Max-Change: 0.00112
Iteration: 99, Log-Lik: -2674.868, Max-Change: 0.00023
Iteration: 100, Log-Lik: -2674.868, Max-Change: 0.00113
Iteration: 101, Log-Lik: -2674.868, Max-Change: 0.00039
Iteration: 102, Log-Lik: -2674.868, Max-Change: 0.00024
Iteration: 103, Log-Lik: -2674.868, Max-Change: 0.00108
Iteration: 104, Log-Lik: -2674.868, Max-Change: 0.00019
Iteration: 105, Log-Lik: -2674.868, Max-Change: 0.00022
Iteration: 106, Log-Lik: -2674.868, Max-Change: 0.00100
Iteration: 107, Log-Lik: -2674.868, Max-Change: 0.00049
Iteration: 108, Log-Lik: -2674.868, Max-Change: 0.00031
Iteration: 109, Log-Lik: -2674.868, Max-Change: 0.00098
Iteration: 110, Log-Lik: -2674.868, Max-Change: 0.00044
Iteration: 111, Log-Lik: -2674.868, Max-Change: 0.00028
Iteration: 112, Log-Lik: -2674.868, Max-Change: 0.00096
Iteration: 113, Log-Lik: -2674.868, Max-Change: 0.00035
Iteration: 114, Log-Lik: -2674.868, Max-Change: 0.00022
Iteration: 115, Log-Lik: -2674.868, Max-Change: 0.00095
Iteration: 116, Log-Lik: -2674.868, Max-Change: 0.00020
Iteration: 117, Log-Lik: -2674.868, Max-Change: 0.00089
Iteration: 118, Log-Lik: -2674.868, Max-Change: 0.00068
Iteration: 119, Log-Lik: -2674.868, Max-Change: 0.00042
Iteration: 120, Log-Lik: -2674.868, Max-Change: 0.00026
Iteration: 121, Log-Lik: -2674.868, Max-Change: 0.00093
Iteration: 122, Log-Lik: -2674.868, Max-Change: 0.00035
Iteration: 123, Log-Lik: -2674.868, Max-Change: 0.00021
Iteration: 124, Log-Lik: -2674.868, Max-Change: 0.00089
Iteration: 125, Log-Lik: -2674.868, Max-Change: 0.00021
Iteration: 126, Log-Lik: -2674.868, Max-Change: 0.00018
Iteration: 127, Log-Lik: -2674.868, Max-Change: 0.00085
Iteration: 128, Log-Lik: -2674.868, Max-Change: 0.00018
Iteration: 129, Log-Lik: -2674.868, Max-Change: 0.00016
Iteration: 130, Log-Lik: -2674.868, Max-Change: 0.00085
Iteration: 131, Log-Lik: -2674.868, Max-Change: 0.00029
Iteration: 132, Log-Lik: -2674.868, Max-Change: 0.00018
Iteration: 133, Log-Lik: -2674.868, Max-Change: 0.00082
Iteration: 134, Log-Lik: -2674.868, Max-Change: 0.00015
Iteration: 135, Log-Lik: -2674.868, Max-Change: 0.00017
Iteration: 136, Log-Lik: -2674.868, Max-Change: 0.00076
Iteration: 137, Log-Lik: -2674.868, Max-Change: 0.00037
Iteration: 138, Log-Lik: -2674.868, Max-Change: 0.00023
Iteration: 139, Log-Lik: -2674.868, Max-Change: 0.00075
Iteration: 140, Log-Lik: -2674.868, Max-Change: 0.00034
Iteration: 141, Log-Lik: -2674.868, Max-Change: 0.00021
Iteration: 142, Log-Lik: -2674.868, Max-Change: 0.00073
Iteration: 143, Log-Lik: -2674.868, Max-Change: 0.00028
Iteration: 144, Log-Lik: -2674.868, Max-Change: 0.00018
Iteration: 145, Log-Lik: -2674.868, Max-Change: 0.00072
Iteration: 146, Log-Lik: -2674.868, Max-Change: 0.00017
Iteration: 147, Log-Lik: -2674.868, Max-Change: 0.00013
Iteration: 148, Log-Lik: -2674.868, Max-Change: 0.00072
Iteration: 149, Log-Lik: -2674.868, Max-Change: 0.00013
Iteration: 150, Log-Lik: -2674.868, Max-Change: 0.00015
Iteration: 151, Log-Lik: -2674.868, Max-Change: 0.00067
Iteration: 152, Log-Lik: -2674.867, Max-Change: 0.00025
Iteration: 153, Log-Lik: -2674.867, Max-Change: 0.00016
Iteration: 154, Log-Lik: -2674.867, Max-Change: 0.00067
Iteration: 155, Log-Lik: -2674.867, Max-Change: 0.00015
Iteration: 156, Log-Lik: -2674.867, Max-Change: 0.00012
Iteration: 157, Log-Lik: -2674.867, Max-Change: 0.00067
Iteration: 158, Log-Lik: -2674.867, Max-Change: 0.00017
Iteration: 159, Log-Lik: -2674.867, Max-Change: 0.00014
Iteration: 160, Log-Lik: -2674.867, Max-Change: 0.00064
Iteration: 161, Log-Lik: -2674.867, Max-Change: 0.00013
Iteration: 162, Log-Lik: -2674.867, Max-Change: 0.00060
Iteration: 163, Log-Lik: -2674.867, Max-Change: 0.00037
Iteration: 164, Log-Lik: -2674.867, Max-Change: 0.00023
Iteration: 165, Log-Lik: -2674.867, Max-Change: 0.00014
Iteration: 166, Log-Lik: -2674.867, Max-Change: 0.00061
Iteration: 167, Log-Lik: -2674.867, Max-Change: 0.00012
Iteration: 168, Log-Lik: -2674.867, Max-Change: 0.00013
Iteration: 169, Log-Lik: -2674.867, Max-Change: 0.00058
Iteration: 170, Log-Lik: -2674.867, Max-Change: 0.00016
Iteration: 171, Log-Lik: -2674.867, Max-Change: 0.00011
Iteration: 172, Log-Lik: -2674.867, Max-Change: 0.00058
Iteration: 173, Log-Lik: -2674.867, Max-Change: 0.00055
Iteration: 174, Log-Lik: -2674.867, Max-Change: 0.00019
Iteration: 175, Log-Lik: -2674.867, Max-Change: 0.00012
Iteration: 176, Log-Lik: -2674.867, Max-Change: 0.00053
Iteration: 177, Log-Lik: -2674.867, Max-Change: 0.00033
Iteration: 178, Log-Lik: -2674.867, Max-Change: 0.00012
Iteration: 179, Log-Lik: -2674.867, Max-Change: 0.00010
Iteration: 180, Log-Lik: -2674.867, Max-Change: 0.00055
Iteration: 181, Log-Lik: -2674.867, Max-Change: 0.00028
Iteration: 182, Log-Lik: -2674.867, Max-Change: 0.00017
Iteration: 183, Log-Lik: -2674.867, Max-Change: 0.00011
Iteration: 184, Log-Lik: -2674.867, Max-Change: 0.00052
Iteration: 185, Log-Lik: -2674.867, Max-Change: 0.00011
Iteration: 186, Log-Lik: -2674.867, Max-Change: 0.00049
Iteration: 187, Log-Lik: -2674.867, Max-Change: 0.00024
Iteration: 188, Log-Lik: -2674.867, Max-Change: 0.00015
Iteration: 189, Log-Lik: -2674.867, Max-Change: 0.00011
Iteration: 190, Log-Lik: -2674.867, Max-Change: 0.00049
Iteration: 191, Log-Lik: -2674.867, Max-Change: 0.00048
Iteration: 192, Log-Lik: -2674.867, Max-Change: 0.00048
Iteration: 193, Log-Lik: -2674.867, Max-Change: 0.00078
Iteration: 194, Log-Lik: -2674.867, Max-Change: 0.00048
Iteration: 195, Log-Lik: -2674.867, Max-Change: 0.00030
Iteration: 196, Log-Lik: -2674.867, Max-Change: 0.00010
Iteration: 197, Log-Lik: -2674.867, Max-Change: 0.00009
plot(read.mirt(model))
itemplot(model, 1)
# graded
mod2 <- mirt(Science, 2)
#>
Iteration: 1, Log-Lik: -1634.562, Max-Change: 0.49365
Iteration: 2, Log-Lik: -1608.496, Max-Change: 0.17069
Iteration: 3, Log-Lik: -1604.784, Max-Change: 0.08523
Iteration: 4, Log-Lik: -1604.195, Max-Change: 0.05345
Iteration: 5, Log-Lik: -1603.777, Max-Change: 0.03491
Iteration: 6, Log-Lik: -1603.581, Max-Change: 0.02948
Iteration: 7, Log-Lik: -1603.263, Max-Change: 0.02614
Iteration: 8, Log-Lik: -1603.177, Max-Change: 0.02281
Iteration: 9, Log-Lik: -1603.099, Max-Change: 0.02311
Iteration: 10, Log-Lik: -1602.717, Max-Change: 0.01961
Iteration: 11, Log-Lik: -1602.673, Max-Change: 0.01556
Iteration: 12, Log-Lik: -1602.635, Max-Change: 0.01694
Iteration: 13, Log-Lik: -1602.443, Max-Change: 0.01534
Iteration: 14, Log-Lik: -1602.421, Max-Change: 0.01311
Iteration: 15, Log-Lik: -1602.401, Max-Change: 0.01384
Iteration: 16, Log-Lik: -1602.303, Max-Change: 0.00987
Iteration: 17, Log-Lik: -1602.288, Max-Change: 0.01006
Iteration: 18, Log-Lik: -1602.275, Max-Change: 0.00939
Iteration: 19, Log-Lik: -1602.206, Max-Change: 0.00972
Iteration: 20, Log-Lik: -1602.198, Max-Change: 0.01016
Iteration: 21, Log-Lik: -1602.189, Max-Change: 0.00670
Iteration: 22, Log-Lik: -1602.165, Max-Change: 0.00679
Iteration: 23, Log-Lik: -1602.158, Max-Change: 0.00694
Iteration: 24, Log-Lik: -1602.152, Max-Change: 0.00555
Iteration: 25, Log-Lik: -1602.116, Max-Change: 0.00689
Iteration: 26, Log-Lik: -1602.112, Max-Change: 0.00465
Iteration: 27, Log-Lik: -1602.107, Max-Change: 0.00478
Iteration: 28, Log-Lik: -1602.082, Max-Change: 0.00389
Iteration: 29, Log-Lik: -1602.078, Max-Change: 0.00745
Iteration: 30, Log-Lik: -1602.074, Max-Change: 0.00380
Iteration: 31, Log-Lik: -1602.069, Max-Change: 0.00373
Iteration: 32, Log-Lik: -1602.066, Max-Change: 0.00792
Iteration: 33, Log-Lik: -1602.063, Max-Change: 0.00397
Iteration: 34, Log-Lik: -1602.060, Max-Change: 0.00364
Iteration: 35, Log-Lik: -1602.057, Max-Change: 0.00808
Iteration: 36, Log-Lik: -1602.055, Max-Change: 0.00332
Iteration: 37, Log-Lik: -1602.052, Max-Change: 0.00613
Iteration: 38, Log-Lik: -1602.049, Max-Change: 0.00319
Iteration: 39, Log-Lik: -1602.047, Max-Change: 0.00596
Iteration: 40, Log-Lik: -1602.043, Max-Change: 0.00303
Iteration: 41, Log-Lik: -1602.040, Max-Change: 0.00302
Iteration: 42, Log-Lik: -1602.038, Max-Change: 0.00611
Iteration: 43, Log-Lik: -1602.034, Max-Change: 0.00392
Iteration: 44, Log-Lik: -1602.033, Max-Change: 0.00329
Iteration: 45, Log-Lik: -1602.030, Max-Change: 0.00406
Iteration: 46, Log-Lik: -1602.026, Max-Change: 0.00303
Iteration: 47, Log-Lik: -1602.025, Max-Change: 0.00261
Iteration: 48, Log-Lik: -1602.023, Max-Change: 0.00457
Iteration: 49, Log-Lik: -1602.021, Max-Change: 0.00265
Iteration: 50, Log-Lik: -1602.020, Max-Change: 0.00290
Iteration: 51, Log-Lik: -1602.018, Max-Change: 0.00268
Iteration: 52, Log-Lik: -1602.010, Max-Change: 0.00310
Iteration: 53, Log-Lik: -1602.008, Max-Change: 0.00361
Iteration: 54, Log-Lik: -1602.007, Max-Change: 0.00238
Iteration: 55, Log-Lik: -1602.006, Max-Change: 0.00327
Iteration: 56, Log-Lik: -1602.005, Max-Change: 0.00269
Iteration: 57, Log-Lik: -1602.004, Max-Change: 0.00324
Iteration: 58, Log-Lik: -1602.002, Max-Change: 0.00218
Iteration: 59, Log-Lik: -1602.001, Max-Change: 0.00198
Iteration: 60, Log-Lik: -1602.000, Max-Change: 0.00346
Iteration: 61, Log-Lik: -1601.999, Max-Change: 0.00226
Iteration: 62, Log-Lik: -1601.998, Max-Change: 0.00247
Iteration: 63, Log-Lik: -1601.997, Max-Change: 0.00264
Iteration: 64, Log-Lik: -1601.996, Max-Change: 0.00380
Iteration: 65, Log-Lik: -1601.995, Max-Change: 0.00236
Iteration: 66, Log-Lik: -1601.994, Max-Change: 0.00244
Iteration: 67, Log-Lik: -1601.992, Max-Change: 0.00164
Iteration: 68, Log-Lik: -1601.992, Max-Change: 0.00214
Iteration: 69, Log-Lik: -1601.991, Max-Change: 0.00193
Iteration: 70, Log-Lik: -1601.990, Max-Change: 0.00267
Iteration: 71, Log-Lik: -1601.989, Max-Change: 0.00227
Iteration: 72, Log-Lik: -1601.989, Max-Change: 0.00232
Iteration: 73, Log-Lik: -1601.988, Max-Change: 0.00172
Iteration: 74, Log-Lik: -1601.987, Max-Change: 0.00149
Iteration: 75, Log-Lik: -1601.987, Max-Change: 0.00245
Iteration: 76, Log-Lik: -1601.986, Max-Change: 0.00149
Iteration: 77, Log-Lik: -1601.985, Max-Change: 0.00215
Iteration: 78, Log-Lik: -1601.985, Max-Change: 0.00206
Iteration: 79, Log-Lik: -1601.984, Max-Change: 0.00223
Iteration: 80, Log-Lik: -1601.983, Max-Change: 0.00190
Iteration: 81, Log-Lik: -1601.983, Max-Change: 0.00181
Iteration: 82, Log-Lik: -1601.982, Max-Change: 0.00125
Iteration: 83, Log-Lik: -1601.981, Max-Change: 0.00136
Iteration: 84, Log-Lik: -1601.981, Max-Change: 0.00207
Iteration: 85, Log-Lik: -1601.981, Max-Change: 0.00124
Iteration: 86, Log-Lik: -1601.980, Max-Change: 0.00187
Iteration: 87, Log-Lik: -1601.980, Max-Change: 0.00123
Iteration: 88, Log-Lik: -1601.980, Max-Change: 0.00193
Iteration: 89, Log-Lik: -1601.980, Max-Change: 0.00123
Iteration: 90, Log-Lik: -1601.979, Max-Change: 0.00198
Iteration: 91, Log-Lik: -1601.979, Max-Change: 0.00118
Iteration: 92, Log-Lik: -1601.979, Max-Change: 0.00121
Iteration: 93, Log-Lik: -1601.979, Max-Change: 0.00196
Iteration: 94, Log-Lik: -1601.978, Max-Change: 0.00112
Iteration: 95, Log-Lik: -1601.978, Max-Change: 0.00117
Iteration: 96, Log-Lik: -1601.978, Max-Change: 0.00193
Iteration: 97, Log-Lik: -1601.978, Max-Change: 0.00106
Iteration: 98, Log-Lik: -1601.977, Max-Change: 0.00112
Iteration: 99, Log-Lik: -1601.977, Max-Change: 0.00187
Iteration: 100, Log-Lik: -1601.977, Max-Change: 0.00103
Iteration: 101, Log-Lik: -1601.977, Max-Change: 0.00109
Iteration: 102, Log-Lik: -1601.977, Max-Change: 0.00127
Iteration: 103, Log-Lik: -1601.976, Max-Change: 0.00118
Iteration: 104, Log-Lik: -1601.976, Max-Change: 0.00115
Iteration: 105, Log-Lik: -1601.976, Max-Change: 0.00115
Iteration: 106, Log-Lik: -1601.975, Max-Change: 0.00108
Iteration: 107, Log-Lik: -1601.975, Max-Change: 0.00108
Iteration: 108, Log-Lik: -1601.975, Max-Change: 0.00107
Iteration: 109, Log-Lik: -1601.974, Max-Change: 0.00103
Iteration: 110, Log-Lik: -1601.974, Max-Change: 0.00102
Iteration: 111, Log-Lik: -1601.974, Max-Change: 0.00101
Iteration: 112, Log-Lik: -1601.974, Max-Change: 0.00083
Iteration: 113, Log-Lik: -1601.974, Max-Change: 0.00088
Iteration: 114, Log-Lik: -1601.974, Max-Change: 0.00090
Iteration: 115, Log-Lik: -1601.973, Max-Change: 0.00108
Iteration: 116, Log-Lik: -1601.973, Max-Change: 0.00037
Iteration: 117, Log-Lik: -1601.973, Max-Change: 0.00138
Iteration: 118, Log-Lik: -1601.973, Max-Change: 0.00320
Iteration: 119, Log-Lik: -1601.973, Max-Change: 0.00085
Iteration: 120, Log-Lik: -1601.973, Max-Change: 0.00039
Iteration: 121, Log-Lik: -1601.973, Max-Change: 0.00034
Iteration: 122, Log-Lik: -1601.973, Max-Change: 0.00026
Iteration: 123, Log-Lik: -1601.973, Max-Change: 0.00109
Iteration: 124, Log-Lik: -1601.973, Max-Change: 0.00128
Iteration: 125, Log-Lik: -1601.973, Max-Change: 0.00113
Iteration: 126, Log-Lik: -1601.973, Max-Change: 0.00094
Iteration: 127, Log-Lik: -1601.973, Max-Change: 0.00080
Iteration: 128, Log-Lik: -1601.972, Max-Change: 0.00028
Iteration: 129, Log-Lik: -1601.972, Max-Change: 0.00105
Iteration: 130, Log-Lik: -1601.972, Max-Change: 0.00291
Iteration: 131, Log-Lik: -1601.972, Max-Change: 0.00082
Iteration: 132, Log-Lik: -1601.972, Max-Change: 0.00034
Iteration: 133, Log-Lik: -1601.972, Max-Change: 0.00029
Iteration: 134, Log-Lik: -1601.972, Max-Change: 0.00112
Iteration: 135, Log-Lik: -1601.972, Max-Change: 0.00244
Iteration: 136, Log-Lik: -1601.972, Max-Change: 0.00039
Iteration: 137, Log-Lik: -1601.972, Max-Change: 0.00027
Iteration: 138, Log-Lik: -1601.972, Max-Change: 0.00022
Iteration: 139, Log-Lik: -1601.972, Max-Change: 0.00112
Iteration: 140, Log-Lik: -1601.972, Max-Change: 0.00024
Iteration: 141, Log-Lik: -1601.972, Max-Change: 0.00098
Iteration: 142, Log-Lik: -1601.972, Max-Change: 0.00261
Iteration: 143, Log-Lik: -1601.972, Max-Change: 0.00031
Iteration: 144, Log-Lik: -1601.972, Max-Change: 0.00023
Iteration: 145, Log-Lik: -1601.972, Max-Change: 0.00208
Iteration: 146, Log-Lik: -1601.972, Max-Change: 0.00074
Iteration: 147, Log-Lik: -1601.972, Max-Change: 0.00031
Iteration: 148, Log-Lik: -1601.972, Max-Change: 0.00027
Iteration: 149, Log-Lik: -1601.971, Max-Change: 0.00104
Iteration: 150, Log-Lik: -1601.971, Max-Change: 0.00114
Iteration: 151, Log-Lik: -1601.971, Max-Change: 0.00026
Iteration: 152, Log-Lik: -1601.971, Max-Change: 0.00099
Iteration: 153, Log-Lik: -1601.971, Max-Change: 0.00235
Iteration: 154, Log-Lik: -1601.971, Max-Change: 0.00028
Iteration: 155, Log-Lik: -1601.971, Max-Change: 0.00022
Iteration: 156, Log-Lik: -1601.971, Max-Change: 0.00019
Iteration: 157, Log-Lik: -1601.971, Max-Change: 0.00104
Iteration: 158, Log-Lik: -1601.971, Max-Change: 0.00023
Iteration: 159, Log-Lik: -1601.971, Max-Change: 0.00088
Iteration: 160, Log-Lik: -1601.971, Max-Change: 0.00085
Iteration: 161, Log-Lik: -1601.971, Max-Change: 0.00030
Iteration: 162, Log-Lik: -1601.971, Max-Change: 0.00019
Iteration: 163, Log-Lik: -1601.971, Max-Change: 0.00075
Iteration: 164, Log-Lik: -1601.971, Max-Change: 0.00106
Iteration: 165, Log-Lik: -1601.971, Max-Change: 0.00093
Iteration: 166, Log-Lik: -1601.971, Max-Change: 0.00081
Iteration: 167, Log-Lik: -1601.971, Max-Change: 0.00068
Iteration: 168, Log-Lik: -1601.971, Max-Change: 0.00022
Iteration: 169, Log-Lik: -1601.971, Max-Change: 0.00020
Iteration: 170, Log-Lik: -1601.971, Max-Change: 0.00080
Iteration: 171, Log-Lik: -1601.971, Max-Change: 0.00104
Iteration: 172, Log-Lik: -1601.971, Max-Change: 0.00026
Iteration: 173, Log-Lik: -1601.971, Max-Change: 0.00019
Iteration: 174, Log-Lik: -1601.971, Max-Change: 0.00078
Iteration: 175, Log-Lik: -1601.971, Max-Change: 0.00103
Iteration: 176, Log-Lik: -1601.971, Max-Change: 0.00071
Iteration: 177, Log-Lik: -1601.971, Max-Change: 0.00014
Iteration: 178, Log-Lik: -1601.971, Max-Change: 0.00066
Iteration: 179, Log-Lik: -1601.971, Max-Change: 0.00021
Iteration: 180, Log-Lik: -1601.971, Max-Change: 0.00065
Iteration: 181, Log-Lik: -1601.971, Max-Change: 0.00147
Iteration: 182, Log-Lik: -1601.970, Max-Change: 0.00020
Iteration: 183, Log-Lik: -1601.970, Max-Change: 0.00017
Iteration: 184, Log-Lik: -1601.970, Max-Change: 0.00062
Iteration: 185, Log-Lik: -1601.970, Max-Change: 0.00058
Iteration: 186, Log-Lik: -1601.970, Max-Change: 0.00064
Iteration: 187, Log-Lik: -1601.970, Max-Change: 0.00105
Iteration: 188, Log-Lik: -1601.970, Max-Change: 0.00019
Iteration: 189, Log-Lik: -1601.970, Max-Change: 0.00077
Iteration: 190, Log-Lik: -1601.970, Max-Change: 0.00071
Iteration: 191, Log-Lik: -1601.970, Max-Change: 0.00021
Iteration: 192, Log-Lik: -1601.970, Max-Change: 0.00077
Iteration: 193, Log-Lik: -1601.970, Max-Change: 0.00081
Iteration: 194, Log-Lik: -1601.970, Max-Change: 0.00028
Iteration: 195, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 196, Log-Lik: -1601.970, Max-Change: 0.00069
Iteration: 197, Log-Lik: -1601.970, Max-Change: 0.00063
Iteration: 198, Log-Lik: -1601.970, Max-Change: 0.00020
Iteration: 199, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 200, Log-Lik: -1601.970, Max-Change: 0.00068
Iteration: 201, Log-Lik: -1601.970, Max-Change: 0.00064
Iteration: 202, Log-Lik: -1601.970, Max-Change: 0.00023
Iteration: 203, Log-Lik: -1601.970, Max-Change: 0.00016
Iteration: 204, Log-Lik: -1601.970, Max-Change: 0.00062
Iteration: 205, Log-Lik: -1601.970, Max-Change: 0.00064
Iteration: 206, Log-Lik: -1601.970, Max-Change: 0.00021
Iteration: 207, Log-Lik: -1601.970, Max-Change: 0.00072
Iteration: 208, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 209, Log-Lik: -1601.970, Max-Change: 0.00061
Iteration: 210, Log-Lik: -1601.970, Max-Change: 0.00020
Iteration: 211, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 212, Log-Lik: -1601.970, Max-Change: 0.00067
Iteration: 213, Log-Lik: -1601.970, Max-Change: 0.00065
Iteration: 214, Log-Lik: -1601.970, Max-Change: 0.00021
Iteration: 215, Log-Lik: -1601.970, Max-Change: 0.00014
Iteration: 216, Log-Lik: -1601.970, Max-Change: 0.00058
Iteration: 217, Log-Lik: -1601.970, Max-Change: 0.00060
Iteration: 218, Log-Lik: -1601.970, Max-Change: 0.00020
Iteration: 219, Log-Lik: -1601.970, Max-Change: 0.00068
Iteration: 220, Log-Lik: -1601.970, Max-Change: 0.00019
Iteration: 221, Log-Lik: -1601.970, Max-Change: 0.00059
Iteration: 222, Log-Lik: -1601.970, Max-Change: 0.00019
Iteration: 223, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 224, Log-Lik: -1601.970, Max-Change: 0.00064
Iteration: 225, Log-Lik: -1601.970, Max-Change: 0.00013
Iteration: 226, Log-Lik: -1601.970, Max-Change: 0.00059
Iteration: 227, Log-Lik: -1601.970, Max-Change: 0.00020
Iteration: 228, Log-Lik: -1601.970, Max-Change: 0.00013
Iteration: 229, Log-Lik: -1601.970, Max-Change: 0.00054
Iteration: 230, Log-Lik: -1601.970, Max-Change: 0.00052
Iteration: 231, Log-Lik: -1601.970, Max-Change: 0.00015
Iteration: 232, Log-Lik: -1601.970, Max-Change: 0.00014
Iteration: 233, Log-Lik: -1601.970, Max-Change: 0.00055
Iteration: 234, Log-Lik: -1601.970, Max-Change: 0.00052
Iteration: 235, Log-Lik: -1601.970, Max-Change: 0.00019
Iteration: 236, Log-Lik: -1601.970, Max-Change: 0.00013
Iteration: 237, Log-Lik: -1601.970, Max-Change: 0.00053
Iteration: 238, Log-Lik: -1601.970, Max-Change: 0.00054
Iteration: 239, Log-Lik: -1601.970, Max-Change: 0.00017
Iteration: 240, Log-Lik: -1601.970, Max-Change: 0.00062
Iteration: 241, Log-Lik: -1601.970, Max-Change: 0.00017
Iteration: 242, Log-Lik: -1601.970, Max-Change: 0.00053
Iteration: 243, Log-Lik: -1601.970, Max-Change: 0.00017
Iteration: 244, Log-Lik: -1601.970, Max-Change: 0.00016
Iteration: 245, Log-Lik: -1601.970, Max-Change: 0.00058
Iteration: 246, Log-Lik: -1601.970, Max-Change: 0.00011
Iteration: 247, Log-Lik: -1601.970, Max-Change: 0.00054
Iteration: 248, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 249, Log-Lik: -1601.970, Max-Change: 0.00062
Iteration: 250, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 251, Log-Lik: -1601.970, Max-Change: 0.00054
Iteration: 252, Log-Lik: -1601.970, Max-Change: 0.00017
Iteration: 253, Log-Lik: -1601.970, Max-Change: 0.00016
Iteration: 254, Log-Lik: -1601.970, Max-Change: 0.00058
Iteration: 255, Log-Lik: -1601.970, Max-Change: 0.00012
Iteration: 256, Log-Lik: -1601.970, Max-Change: 0.00055
Iteration: 257, Log-Lik: -1601.970, Max-Change: 0.00018
Iteration: 258, Log-Lik: -1601.970, Max-Change: 0.00012
Iteration: 259, Log-Lik: -1601.970, Max-Change: 0.00050
Iteration: 260, Log-Lik: -1601.969, Max-Change: 0.00048
Iteration: 261, Log-Lik: -1601.969, Max-Change: 0.00014
Iteration: 262, Log-Lik: -1601.969, Max-Change: 0.00013
Iteration: 263, Log-Lik: -1601.969, Max-Change: 0.00051
Iteration: 264, Log-Lik: -1601.969, Max-Change: 0.00049
Iteration: 265, Log-Lik: -1601.969, Max-Change: 0.00017
Iteration: 266, Log-Lik: -1601.969, Max-Change: 0.00059
Iteration: 267, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 268, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 269, Log-Lik: -1601.969, Max-Change: 0.00047
Iteration: 270, Log-Lik: -1601.969, Max-Change: 0.00013
Iteration: 271, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 272, Log-Lik: -1601.969, Max-Change: 0.00049
Iteration: 273, Log-Lik: -1601.969, Max-Change: 0.00047
Iteration: 274, Log-Lik: -1601.969, Max-Change: 0.00016
Iteration: 275, Log-Lik: -1601.969, Max-Change: 0.00058
Iteration: 276, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 277, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 278, Log-Lik: -1601.969, Max-Change: 0.00046
Iteration: 279, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 280, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 281, Log-Lik: -1601.969, Max-Change: 0.00048
Iteration: 282, Log-Lik: -1601.969, Max-Change: 0.00047
Iteration: 283, Log-Lik: -1601.969, Max-Change: 0.00016
Iteration: 284, Log-Lik: -1601.969, Max-Change: 0.00056
Iteration: 285, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 286, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 287, Log-Lik: -1601.969, Max-Change: 0.00045
Iteration: 288, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 289, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 290, Log-Lik: -1601.969, Max-Change: 0.00047
Iteration: 291, Log-Lik: -1601.969, Max-Change: 0.00046
Iteration: 292, Log-Lik: -1601.969, Max-Change: 0.00016
Iteration: 293, Log-Lik: -1601.969, Max-Change: 0.00055
Iteration: 294, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 295, Log-Lik: -1601.969, Max-Change: 0.00010
Iteration: 296, Log-Lik: -1601.969, Max-Change: 0.00044
Iteration: 297, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 298, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 299, Log-Lik: -1601.969, Max-Change: 0.00046
Iteration: 300, Log-Lik: -1601.969, Max-Change: 0.00045
Iteration: 301, Log-Lik: -1601.969, Max-Change: 0.00015
Iteration: 302, Log-Lik: -1601.969, Max-Change: 0.00054
Iteration: 303, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 304, Log-Lik: -1601.969, Max-Change: 0.00010
Iteration: 305, Log-Lik: -1601.969, Max-Change: 0.00043
Iteration: 306, Log-Lik: -1601.969, Max-Change: 0.00012
Iteration: 307, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 308, Log-Lik: -1601.969, Max-Change: 0.00045
Iteration: 309, Log-Lik: -1601.969, Max-Change: 0.00044
Iteration: 310, Log-Lik: -1601.969, Max-Change: 0.00015
Iteration: 311, Log-Lik: -1601.969, Max-Change: 0.00053
Iteration: 312, Log-Lik: -1601.969, Max-Change: 0.00011
Iteration: 313, Log-Lik: -1601.969, Max-Change: 0.00010
plinkpars <- read.mirt(mod2)
plinkpars
#> An object of class "irt.pars"
#> Slot "pars":
#> a1 a2 b1 b2 b3
#> [1,] -1.3350281 0.09676376 5.210669 2.865549 -1.602826
#> [2,] -0.8789508 1.85253465 3.703802 1.153176 -2.904255
#> [3,] -1.4696076 1.16485648 4.663467 1.956626 -1.735796
#> [4,] -1.7220434 0.00000000 3.988789 1.195247 -2.043998
#>
#> Slot "cat":
#> [1] 4 4 4 4
#>
#> Slot "poly.mod":
#> An object of class "poly.mod"
#> Slot "model":
#> [1] "grm"
#>
#> Slot "items":
#> $grm
#> [1] 1 2 3 4
#>
#>
#>
#> Slot "common":
#> NULL
#>
#> Slot "location":
#> [1] FALSE
#>
#> Slot "groups":
#> [1] 1
#>
#> Slot "dimensions":
#> [1] 2
#>
plot(plinkpars)
plot(mod2, type = 'trace')
### multiple group equating example
set.seed(1234)
dat <- expand.table(LSAT7)
group <- sample(c('g1', 'g2'), nrow(dat), TRUE)
dat1 <- dat[group == 'g1', ]
dat2 <- dat[group == 'g2', ]
mod1 <- mirt(dat1, 1)
#>
Iteration: 1, Log-Lik: -1329.923, Max-Change: 0.15128
Iteration: 2, Log-Lik: -1328.209, Max-Change: 0.11669
Iteration: 3, Log-Lik: -1327.382, Max-Change: 0.09194
Iteration: 4, Log-Lik: -1326.251, Max-Change: 0.03589
Iteration: 5, Log-Lik: -1326.130, Max-Change: 0.03273
Iteration: 6, Log-Lik: -1326.043, Max-Change: 0.02967
Iteration: 7, Log-Lik: -1325.810, Max-Change: 0.01307
Iteration: 8, Log-Lik: -1325.804, Max-Change: 0.01068
Iteration: 9, Log-Lik: -1325.800, Max-Change: 0.00885
Iteration: 10, Log-Lik: -1325.793, Max-Change: 0.00306
Iteration: 11, Log-Lik: -1325.792, Max-Change: 0.00301
Iteration: 12, Log-Lik: -1325.792, Max-Change: 0.00277
Iteration: 13, Log-Lik: -1325.791, Max-Change: 0.00148
Iteration: 14, Log-Lik: -1325.791, Max-Change: 0.00138
Iteration: 15, Log-Lik: -1325.790, Max-Change: 0.00028
Iteration: 16, Log-Lik: -1325.790, Max-Change: 0.00027
Iteration: 17, Log-Lik: -1325.790, Max-Change: 0.00112
Iteration: 18, Log-Lik: -1325.790, Max-Change: 0.00088
Iteration: 19, Log-Lik: -1325.790, Max-Change: 0.00025
Iteration: 20, Log-Lik: -1325.790, Max-Change: 0.00089
Iteration: 21, Log-Lik: -1325.790, Max-Change: 0.00080
Iteration: 22, Log-Lik: -1325.790, Max-Change: 0.00034
Iteration: 23, Log-Lik: -1325.790, Max-Change: 0.00015
Iteration: 24, Log-Lik: -1325.790, Max-Change: 0.00063
Iteration: 25, Log-Lik: -1325.790, Max-Change: 0.00058
Iteration: 26, Log-Lik: -1325.790, Max-Change: 0.00015
Iteration: 27, Log-Lik: -1325.790, Max-Change: 0.00011
Iteration: 28, Log-Lik: -1325.790, Max-Change: 0.00047
Iteration: 29, Log-Lik: -1325.790, Max-Change: 0.00039
Iteration: 30, Log-Lik: -1325.790, Max-Change: 0.00011
Iteration: 31, Log-Lik: -1325.790, Max-Change: 0.00010
mod2 <- mirt(dat2, 1)
#>
Iteration: 1, Log-Lik: -1335.081, Max-Change: 0.21322
Iteration: 2, Log-Lik: -1331.357, Max-Change: 0.15695
Iteration: 3, Log-Lik: -1329.775, Max-Change: 0.11402
Iteration: 4, Log-Lik: -1328.640, Max-Change: 0.04541
Iteration: 5, Log-Lik: -1328.561, Max-Change: 0.03568
Iteration: 6, Log-Lik: -1328.511, Max-Change: 0.02841
Iteration: 7, Log-Lik: -1328.423, Max-Change: 0.00706
Iteration: 8, Log-Lik: -1328.420, Max-Change: 0.00634
Iteration: 9, Log-Lik: -1328.418, Max-Change: 0.00534
Iteration: 10, Log-Lik: -1328.415, Max-Change: 0.00141
Iteration: 11, Log-Lik: -1328.415, Max-Change: 0.00120
Iteration: 12, Log-Lik: -1328.415, Max-Change: 0.00027
Iteration: 13, Log-Lik: -1328.415, Max-Change: 0.00025
Iteration: 14, Log-Lik: -1328.415, Max-Change: 0.00104
Iteration: 15, Log-Lik: -1328.415, Max-Change: 0.00084
Iteration: 16, Log-Lik: -1328.415, Max-Change: 0.00021
Iteration: 17, Log-Lik: -1328.415, Max-Change: 0.00015
Iteration: 18, Log-Lik: -1328.415, Max-Change: 0.00063
Iteration: 19, Log-Lik: -1328.415, Max-Change: 0.00056
Iteration: 20, Log-Lik: -1328.415, Max-Change: 0.00014
Iteration: 21, Log-Lik: -1328.415, Max-Change: 0.00010
# convert and combine pars
plinkMG <- read.mirt(list(g1=mod1, g2=mod2))
# equivalently:
# mod <- multipleGroup(dat, 1, group)
# plinkMG <- read.mirt(mod)
combine <- matrix(1:5, 5, 2)
comb <- combine.pars(plinkMG, combine, grp.names=unique(group))
out <- plink(comb, rescale="SL")
equate(out)
#> Maximum iterations reached for true score: 0
#> $tse
#> theta g2 g1
#> 1 -160.732949 0 0.000000
#> 2 -3.233741 1 0.973481
#> 3 -1.910192 2 2.023409
#> 4 -0.999050 3 3.021848
#> 5 0.048359 4 3.973672
#> 6 57.411786 5 5.000000
#>
#> $ose
#> $ose$scores
#> eap.theta.g2 eap.sd.g2 g2 g1
#> 1 -1.896444 0.698843 0 0.0000000
#> 2 -1.474472 0.697610 1 0.9893317
#> 3 -1.002953 0.716007 2 1.9956600
#> 4 -0.455590 0.747730 3 2.9912425
#> 5 0.123384 0.782773 4 3.9901784
#> 6 0.678747 0.812793 5 4.9949444
#>
#> $ose$dist
#> $ose$dist$g2
#> score pop1 pop2 syn
#> [1,] 0 0.03273187 0.03273187 0.03273187
#> [2,] 1 0.17343742 0.17343742 0.17343742
#> [3,] 2 0.48733074 0.48733074 0.48733074
#> [4,] 3 0.99873451 0.99873451 0.99873451
#> [5,] 4 1.61799716 1.61799716 1.61799716
#> [6,] 5 1.56457520 1.56457520 1.56457520
#>
#> $ose$dist$g1
#> score pop1 pop2 syn
#> [1,] 0 0.0338755 0.0338755 0.0338755
#> [2,] 1 0.1748815 0.1748815 0.1748815
#> [3,] 2 0.4863770 0.4863770 0.4863770
#> [4,] 3 1.0132130 1.0132130 1.0132130
#> [5,] 4 1.6175460 1.6175460 1.6175460
#> [6,] 5 1.5489139 1.5489139 1.5489139
#>
#>
#>
equate(out, method = 'OSE')
#> $scores
#> eap.theta.g2 eap.sd.g2 g2 g1
#> 1 -1.896444 0.698843 0 0.0000000
#> 2 -1.474472 0.697610 1 0.9893317
#> 3 -1.002953 0.716007 2 1.9956600
#> 4 -0.455590 0.747730 3 2.9912425
#> 5 0.123384 0.782773 4 3.9901784
#> 6 0.678747 0.812793 5 4.9949444
#>
#> $dist
#> $dist$g2
#> score pop1 pop2 syn
#> [1,] 0 0.03273187 0.03273187 0.03273187
#> [2,] 1 0.17343742 0.17343742 0.17343742
#> [3,] 2 0.48733074 0.48733074 0.48733074
#> [4,] 3 0.99873451 0.99873451 0.99873451
#> [5,] 4 1.61799716 1.61799716 1.61799716
#> [6,] 5 1.56457520 1.56457520 1.56457520
#>
#> $dist$g1
#> score pop1 pop2 syn
#> [1,] 0 0.0338755 0.0338755 0.0338755
#> [2,] 1 0.1748815 0.1748815 0.1748815
#> [3,] 2 0.4863770 0.4863770 0.4863770
#> [4,] 3 1.0132130 1.0132130 1.0132130
#> [5,] 4 1.6175460 1.6175460 1.6175460
#> [6,] 5 1.5489139 1.5489139 1.5489139
#>
#>
# }