IRT.RMSD: Root Mean Square Deviation (RMSD) Item Fit Statistic

Description Usage Arguments Details Value References See Also Examples

View source: R/IRT.RMSD.R

Description

Computed the item fit statistics root mean square deviation (RMSD), mean absolute deviation (MAD) and mean deviation (MD). See Oliveri and von Davier (2011) for details.

The RMSD statistics was denoted as the RMSEA statistic in older publications, see itemfit.rmsea.

If multiple groups are defined in the model object, a weighted item fit statistic (WRMSD; Yamamoto, Khorramdel, & von Davier, 2013; von Davier, Weeks, Chen, Allen & van der Velden, 2013) is additionally computed.

Usage

1
2
3
4
5
6
7
IRT.RMSD(object)

## S3 method for class 'IRT.RMSD'
summary(object, file=NULL, digits=3, ...)

## core computation function
IRT_RMSD_calc_rmsd( n.ik, pi.k, probs, eps=1E-30 )

Arguments

object

Object for which the methods IRT.expectedCounts and IRT.irfprob can be applied.

n.ik

Expected counts

pi.k

Probabilities trait distribution

probs

Item response probabilities

eps

Numerical constant avoiding division by zero

digits

Number of digits used for rounding

file

Optional file name for a file in which summary should be sinked.

...

Optional parameters to be passed.

Details

The RMSD and MD statistics are in operational use in PISA studies since PISA 2015. These fit statistics can also be used for investigating uniform and nonuniform differential item functioning.

Value

List with entries

RMSD

Item-wise and group-wise RMSD statistic

RMSD_bc

Item-wise and group-wise RMSD statistic with analytical bias correction

MAD

Item-wise and group-wise MAD statistic

MD

Item-wise and group-wise MD statistic

chisquare_stat

Item-wise and group-wise χ^2 statistic

...

Further values

References

Oliveri, M. E., & von Davier, M. (2011). Investigation of model fit and score scale comparability in international assessments. Psychological Test and Assessment Modeling, 53, 315-333.

von Davier, M., Weeks, J., Chen, H., Allen, J., & van der Velden, R. (2013). Creating simple and complex derived variables and validation of background questionnaire data. In OECD (Eds.). Technical Report of the Survey of Adults Skills (PIAAC) (Ch. 20). Paris: OECD.

Yamamoto, K., Khorramdel, L., & von Davier, M. (2013). Scaling PIAAC cognitive data. In OECD (Eds.). Technical Report of the Survey of Adults Skills (PIAAC) (Ch. 17). Paris: OECD.

See Also

itemfit.rmsea

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
## Not run: 
#############################################################################
# EXAMPLE 1: data.read | 1PL model in TAM
#############################################################################

data(data.read, package="sirt")
dat <- data.read

#*** Model 1: 1PL model
mod1 <- TAM::tam.mml( resp=dat )
summary(mod1)

# item fit statistics
imod1 <- CDM::IRT.RMSD(mod1)
summary(imod1)

#############################################################################
# EXAMPLE 2: data.math| RMSD and MD statistic for assessing DIF
#############################################################################

data(data.math, package="sirt")
dat <- data.math$data
items <- grep("M[A-Z]", colnames(dat), value=TRUE )

#-- fit multiple group Rasch model
mod <- TAM::tam.mml( dat[,items], group=dat$female )
summary(mod)

#-- fit statistics
rmod <- CDM::IRT.RMSD(mod)
summary(rmod)

#############################################################################
# EXAMPLE 3: RMSD statistic DINA model
#############################################################################

data(sim.dina)
data(sim.qmatrix)
dat <- sim.dina
Q <- sim.qmatrix

#-- fit DINA model
mod1 <- CDM::gdina( dat, q.matrix=Q, rule="DINA" )
summary(mod1)

#-- compute RMSD fit statistic
rmod1 <- CDM::IRT.RMSD(mod1)
summary(rmod1)

## End(Not run)

Example output

Loading required package: mvtnorm
**********************************
** CDM 7.0-12 (2018-11-06 16:23:12)      
** Cognitive Diagnostic Models  **
**********************************

....................................................
Processing Data      2019-01-06 08:14:55 
    * Response Data: 328 Persons and  12 Items 
    * Numerical integration with 21 nodes
    * Created Design Matrices   ( 2019-01-06 08:14:55 )
    * Calculated Sufficient Statistics   ( 2019-01-06 08:14:55 )
....................................................
Iteration 1     2019-01-06 08:14:55
E Step
M Step Intercepts   |----
  Deviance = 3955.0021
  Maximum item intercept parameter change: 0.282708
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.001983
....................................................
Iteration 2     2019-01-06 08:14:55
E Step
M Step Intercepts   |---
  Deviance = 3942.7485 | Absolute change: 12.2536 | Relative change: 0.00310788
  Maximum item intercept parameter change: 0.048264
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.038661
....................................................
Iteration 3     2019-01-06 08:14:55
E Step
M Step Intercepts   |---
  Deviance = 3940.8177 | Absolute change: 1.9308 | Relative change: 0.00048995
  Maximum item intercept parameter change: 0.032447
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.033165
....................................................
Iteration 4     2019-01-06 08:14:55
E Step
M Step Intercepts   |---
  Deviance = 3939.8371 | Absolute change: 0.9806 | Relative change: 0.00024889
  Maximum item intercept parameter change: 0.023311
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.027088
....................................................
Iteration 5     2019-01-06 08:14:55
E Step
M Step Intercepts   |---
  Deviance = 3939.3212 | Absolute change: 0.516 | Relative change: 0.00013098
  Maximum item intercept parameter change: 0.016965
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.021496
....................................................
Iteration 6     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3939.0438 | Absolute change: 0.2774 | Relative change: 7.041e-05
  Maximum item intercept parameter change: 0.012453
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.016715
....................................................
Iteration 7     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.8929 | Absolute change: 0.151 | Relative change: 3.832e-05
  Maximum item intercept parameter change: 0.009194
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.012811
....................................................
Iteration 8     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.8101 | Absolute change: 0.0828 | Relative change: 2.101e-05
  Maximum item intercept parameter change: 0.006814
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.009721
....................................................
Iteration 9     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7645 | Absolute change: 0.0456 | Relative change: 1.157e-05
  Maximum item intercept parameter change: 0.005064
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.007324
....................................................
Iteration 10     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7393 | Absolute change: 0.0252 | Relative change: 6.39e-06
  Maximum item intercept parameter change: 0.00377
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.005493
....................................................
Iteration 11     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7254 | Absolute change: 0.0139 | Relative change: 3.54e-06
  Maximum item intercept parameter change: 0.00281
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.004106
....................................................
Iteration 12     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7177 | Absolute change: 0.0077 | Relative change: 1.96e-06
  Maximum item intercept parameter change: 0.002096
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.003063
....................................................
Iteration 13     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7134 | Absolute change: 0.0043 | Relative change: 1.09e-06
  Maximum item intercept parameter change: 0.001564
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.002283
....................................................
Iteration 14     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.711 | Absolute change: 0.0024 | Relative change: 6e-07
  Maximum item intercept parameter change: 0.001168
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.0017
....................................................
Iteration 15     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7097 | Absolute change: 0.0013 | Relative change: 3.4e-07
  Maximum item intercept parameter change: 0.000872
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.001265
....................................................
Iteration 16     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.709 | Absolute change: 7e-04 | Relative change: 1.9e-07
  Maximum item intercept parameter change: 0.000651
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000942
....................................................
Iteration 17     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7085 | Absolute change: 4e-04 | Relative change: 1e-07
  Maximum item intercept parameter change: 0.000486
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000701
....................................................
Iteration 18     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7083 | Absolute change: 2e-04 | Relative change: 6e-08
  Maximum item intercept parameter change: 0.000363
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000522
....................................................
Iteration 19     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7082 | Absolute change: 1e-04 | Relative change: 3e-08
  Maximum item intercept parameter change: 0.000271
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000389
....................................................
Iteration 20     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7081 | Absolute change: 1e-04 | Relative change: 2e-08
  Maximum item intercept parameter change: 0.000202
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000289
....................................................
Iteration 21     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7081 | Absolute change: 0 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000151
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000216
....................................................
Iteration 22     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 3938.7081 | Absolute change: 0 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000113
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.000161
....................................................
Iteration 23     2019-01-06 08:14:55
E Step
M Step Intercepts   |-
  Deviance = 3938.708 | Absolute change: 0 | Relative change: 0
  Maximum item intercept parameter change: 8.4e-05
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 0.00012
....................................................
Iteration 24     2019-01-06 08:14:55
E Step
M Step Intercepts   |-
  Deviance = 3938.708 | Absolute change: 0 | Relative change: 0
  Maximum item intercept parameter change: 6.3e-05
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0
  Maximum variance parameter change: 8.9e-05
....................................................
Item Parameters
   xsi.index xsi.label     est
1          1        A1 -2.1045
2          2        A2 -1.2742
3          3        A3 -0.3379
4          4        A4  0.1953
5          5        B1 -1.1262
6          6        B2 -0.0327
7          7        B3 -2.7346
8          8        B4 -0.9502
9          9        C1 -3.1058
10        10        C2 -1.1262
11        11        C3 -2.3095
12        12        C4 -1.2553
...................................
Regression Coefficients
     [,1]
[1,]    0

Variance:
     [,1]
[1,] 1.19


EAP Reliability:
[1] 0.662

-----------------------------
Start:  2019-01-06 08:14:55
End:  2019-01-06 08:14:55 
Time difference of 0.2555604 secs

------------------------------------------------------------
TAM 2.13-15 (2018-09-30 14:35:48) 
sh: 1: cannot create /dev/null: Permission denied
R version 3.4.4 (2018-03-15) x86_64, linux-gnu | nodename=ip-172-31-1-5 | login=unknown 

Date of Analysis: 2019-01-06 08:14:55 
Time difference of 0.2555604 secs
Computation time: 0.2555604 

Multidimensional Item Response Model in TAM 

IRT Model 1PL

Call:
TAM::tam.mml(resp = dat)

------------------------------------------------------------
Number of iterations= 24 
Numeric integration with 21 integration points

Deviance= 3938.71 
   Log likelihood= -1969.35 
Number of persons= 328 
Number of persons used= 328 
Number of items= 12 
Number of estimated parameters= 13 
    Item threshold parameters= 12 
    Item slope parameters  = 0 
    Regression parameters  = 0 
    (Co)Variance parameters= 1 

AIC = 3965  | penalty= 26    | AIC=-2*LL + 2*p 
AIC3 = 3978  | penalty= 39    | AIC3=-2*LL + 3*p 
BIC = 4014  | penalty= 75.31    | BIC=-2*LL + log(n)*p 
aBIC = 3973  | penalty= 33.91    | aBIC=-2*LL + log((n-2)/24)*p  (adjusted BIC) 
CAIC = 4027  | penalty= 88.31    | CAIC=-2*LL + [log(n)+1]*p  (consistent AIC) 
AICc = 3966  | penalty= 27.16    | AICc=-2*LL + 2*p + 2*p*(p+1)/(n-p-1)  (bias corrected AIC) 

------------------------------------------------------------
EAP Reliability
[1] 0.662
------------------------------------------------------------
Covariances and Variances
     [,1]
[1,] 1.19
------------------------------------------------------------
Correlations and Standard Deviations (in the diagonal)
      [,1]
[1,] 1.091
------------------------------------------------------------
Regression Coefficients
     [,1]
[1,]    0
------------------------------------------------------------
Item Parameters -A*Xsi
   item   N     M xsi.item AXsi_.Cat1 B.Cat1.Dim1
A1   A1 328 0.851   -2.104     -2.104           1
A2   A2 328 0.738   -1.274     -1.274           1
A3   A3 328 0.567   -0.338     -0.338           1
A4   A4 328 0.460    0.195      0.195           1
B1   B1 328 0.713   -1.126     -1.126           1
B2   B2 328 0.506   -0.033     -0.033           1
B3   B3 328 0.909   -2.735     -2.735           1
B4   B4 328 0.683   -0.950     -0.950           1
C1   C1 328 0.933   -3.106     -3.106           1
C2   C2 328 0.713   -1.126     -1.126           1
C3   C3 328 0.872   -2.309     -2.309           1
C4   C4 328 0.735   -1.255     -1.255           1
-----------------------------------------------------------------------------
CDM 7.0-12 (2018-11-06 16:23:12) 

Call:
CDM::IRT.RMSD(object = mod1)

-----------------------------------------------------------------------------
Root Mean Square Deviation (RMSD) 

    Parm     M    SD   Min   Max
1 Group1 0.031 0.012 0.009 0.046

   item Group1
1    A1  0.046
2    A2  0.040
3    A3  0.035
4    A4  0.023
5    B1  0.046
6    B2  0.045
7    B3  0.020
8    B4  0.009
9    C1  0.037
10   C2  0.024
11   C3  0.026
12   C4  0.025
-----------------------------------------------------------------------------
Bias Corrected Root Mean Square Deviation (RMSD) 

    Parm     M   SD    Min   Max
1 Group1 0.019 0.02 -0.022 0.042

   item Group1
1    A1  0.042
2    A2  0.033
3    A3  0.025
4    A4 -0.009
5    B1  0.040
6    B2  0.037
7    B3  0.013
8    B4 -0.022
9    C1  0.035
10   C2  0.007
11   C3  0.020
12   C4  0.011
-----------------------------------------------------------------------------
Mean Absolute Deviation (MAD) 

    Parm     M   SD   Min  Max
1 Group1 0.023 0.01 0.008 0.04

   item Group1
1    A1  0.032
2    A2  0.035
3    A3  0.031
4    A4  0.018
5    B1  0.033
6    B2  0.040
7    B3  0.013
8    B4  0.008
9    C1  0.018
10   C2  0.018
11   C3  0.012
12   C4  0.019
-----------------------------------------------------------------------------
Mean Deviation (MD) 

    Parm M SD Min Max
1 Group1 0  0   0   0

   item Group1
1    A1      0
2    A2      0
3    A3      0
4    A4      0
5    B1      0
6    B2      0
7    B3      0
8    B4      0
9    C1      0
10   C2      0
11   C3      0
12   C4      0
....................................................
Processing Data      2019-01-06 08:14:55 
    * Response Data: 664 Persons and  30 Items 
    * Numerical integration with 21 nodes
    * Created Design Matrices   ( 2019-01-06 08:14:55 )
    * Calculated Sufficient Statistics   ( 2019-01-06 08:14:55 )
....................................................
Iteration 1     2019-01-06 08:14:55
E Step
M Step Intercepts   |----
  Deviance = 22822.8209
  Maximum item intercept parameter change: 0.39373
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.012353
  Maximum variance parameter change: 0.194021
....................................................
Iteration 2     2019-01-06 08:14:55
E Step
M Step Intercepts   |---
  Deviance = 22735.4521 | Absolute change: 87.3688 | Relative change: 0.00384284
  Maximum item intercept parameter change: 0.040335
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.002869
  Maximum variance parameter change: 0.152377
....................................................
Iteration 3     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22729.4836 | Absolute change: 5.9685 | Relative change: 0.00026259
  Maximum item intercept parameter change: 0.016953
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.00124
  Maximum variance parameter change: 0.065705
....................................................
Iteration 4     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.5403 | Absolute change: 0.9434 | Relative change: 4.151e-05
  Maximum item intercept parameter change: 0.007101
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.00083
  Maximum variance parameter change: 0.025931
....................................................
Iteration 5     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3954 | Absolute change: 0.1448 | Relative change: 6.37e-06
  Maximum item intercept parameter change: 0.003092
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000666
  Maximum variance parameter change: 0.00992
....................................................
Iteration 6     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3719 | Absolute change: 0.0235 | Relative change: 1.03e-06
  Maximum item intercept parameter change: 0.001501
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000584
  Maximum variance parameter change: 0.003754
....................................................
Iteration 7     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3668 | Absolute change: 0.0052 | Relative change: 2.3e-07
  Maximum item intercept parameter change: 0.000867
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000535
  Maximum variance parameter change: 0.001414
....................................................
Iteration 8     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3645 | Absolute change: 0.0023 | Relative change: 1e-07
  Maximum item intercept parameter change: 0.000607
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000499
  Maximum variance parameter change: 0.00053
....................................................
Iteration 9     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3628 | Absolute change: 0.0017 | Relative change: 7e-08
  Maximum item intercept parameter change: 0.000492
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.00047
  Maximum variance parameter change: 0.000196
....................................................
Iteration 10     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3613 | Absolute change: 0.0015 | Relative change: 6e-08
  Maximum item intercept parameter change: 0.000434
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000444
  Maximum variance parameter change: 7e-05
....................................................
Iteration 11     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.36 | Absolute change: 0.0013 | Relative change: 6e-08
  Maximum item intercept parameter change: 0.000398
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.00042
  Maximum variance parameter change: 2.6e-05
....................................................
Iteration 12     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3588 | Absolute change: 0.0012 | Relative change: 5e-08
  Maximum item intercept parameter change: 0.000373
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000397
  Maximum variance parameter change: 1e-05
....................................................
Iteration 13     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3578 | Absolute change: 0.001 | Relative change: 5e-08
  Maximum item intercept parameter change: 0.000353
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000375
  Maximum variance parameter change: 4e-06
....................................................
Iteration 14     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3568 | Absolute change: 9e-04 | Relative change: 4e-08
  Maximum item intercept parameter change: 0.000334
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000355
  Maximum variance parameter change: 3e-06
....................................................
Iteration 15     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.356 | Absolute change: 8e-04 | Relative change: 4e-08
  Maximum item intercept parameter change: 0.000316
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000336
  Maximum variance parameter change: 4e-06
....................................................
Iteration 16     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3552 | Absolute change: 8e-04 | Relative change: 3e-08
  Maximum item intercept parameter change: 0.000299
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000318
  Maximum variance parameter change: 4e-06
....................................................
Iteration 17     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3546 | Absolute change: 7e-04 | Relative change: 3e-08
  Maximum item intercept parameter change: 0.000283
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000301
  Maximum variance parameter change: 4e-06
....................................................
Iteration 18     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.354 | Absolute change: 6e-04 | Relative change: 3e-08
  Maximum item intercept parameter change: 0.000268
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000285
  Maximum variance parameter change: 3e-06
....................................................
Iteration 19     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3534 | Absolute change: 5e-04 | Relative change: 2e-08
  Maximum item intercept parameter change: 0.000253
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.00027
  Maximum variance parameter change: 3e-06
....................................................
Iteration 20     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3529 | Absolute change: 5e-04 | Relative change: 2e-08
  Maximum item intercept parameter change: 0.00024
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000255
  Maximum variance parameter change: 3e-06
....................................................
Iteration 21     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3525 | Absolute change: 4e-04 | Relative change: 2e-08
  Maximum item intercept parameter change: 0.000227
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000241
  Maximum variance parameter change: 2e-06
....................................................
Iteration 22     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3521 | Absolute change: 4e-04 | Relative change: 2e-08
  Maximum item intercept parameter change: 0.000215
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000229
  Maximum variance parameter change: 2e-06
....................................................
Iteration 23     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3518 | Absolute change: 3e-04 | Relative change: 2e-08
  Maximum item intercept parameter change: 0.000203
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000216
  Maximum variance parameter change: 2e-06
....................................................
Iteration 24     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3515 | Absolute change: 3e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000192
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000205
  Maximum variance parameter change: 2e-06
....................................................
Iteration 25     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3512 | Absolute change: 3e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000182
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000194
  Maximum variance parameter change: 1e-06
....................................................
Iteration 26     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3509 | Absolute change: 2e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000172
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000183
  Maximum variance parameter change: 1e-06
....................................................
Iteration 27     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3507 | Absolute change: 2e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000163
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000173
  Maximum variance parameter change: 1e-06
....................................................
Iteration 28     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3505 | Absolute change: 2e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000154
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000164
  Maximum variance parameter change: 1e-06
....................................................
Iteration 29     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3503 | Absolute change: 2e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000146
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000155
  Maximum variance parameter change: 1e-06
....................................................
Iteration 30     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3502 | Absolute change: 2e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000138
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000147
  Maximum variance parameter change: 1e-06
....................................................
Iteration 31     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.35 | Absolute change: 1e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000131
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000139
  Maximum variance parameter change: 1e-06
....................................................
Iteration 32     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3499 | Absolute change: 1e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000124
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000132
  Maximum variance parameter change: 1e-06
....................................................
Iteration 33     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3498 | Absolute change: 1e-04 | Relative change: 1e-08
  Maximum item intercept parameter change: 0.000117
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000125
  Maximum variance parameter change: 1e-06
....................................................
Iteration 34     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3497 | Absolute change: 1e-04 | Relative change: 0
  Maximum item intercept parameter change: 0.000111
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000118
  Maximum variance parameter change: 0
....................................................
Iteration 35     2019-01-06 08:14:55
E Step
M Step Intercepts   |--
  Deviance = 22728.3496 | Absolute change: 1e-04 | Relative change: 0
  Maximum item intercept parameter change: 0.000105
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000112
  Maximum variance parameter change: 0
....................................................
Iteration 36     2019-01-06 08:14:56
E Step
M Step Intercepts   |-
  Deviance = 22728.3495 | Absolute change: 1e-04 | Relative change: 0
  Maximum item intercept parameter change: 9.9e-05
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 0.000106
  Maximum variance parameter change: 0
....................................................
Iteration 37     2019-01-06 08:14:56
E Step
M Step Intercepts   |-
  Deviance = 22728.3494 | Absolute change: 1e-04 | Relative change: 0
  Maximum item intercept parameter change: 9.4e-05
  Maximum item slope parameter change: 0
  Maximum regression parameter change: 1e-04
  Maximum variance parameter change: 0
....................................................
Item Parameters
   xsi.index xsi.label     est
1          1       MA1 -2.5225
2          2       MA2 -0.3185
3          3       MA3 -0.1340
4          4       MA4  0.4882
5          5       MB1 -0.6992
6          6       MB2 -0.6019
7          7       MB3  0.7775
8          8       MB4  0.7610
9          9       MC1 -1.4964
10        10       MC2  0.0948
11        11       MC3 -0.8572
12        12       MC4 -0.2799
13        13       MD1 -0.9687
14        14       MD2 -1.3850
15        15       MD3 -0.5380
16        16       MD4 -0.5459
17        17       ME1 -1.0567
18        18       ME2  0.5434
19        19       MF1 -0.3495
20        20       MG1  0.4022
21        21       MG2  0.8189
22        22       MG3  0.3247
23        23       MG4  0.4882
24        24       MH1  0.0110
25        25       MH2  0.6632
26        26       MH3  0.3789
27        27       MH4  0.8692
28        28       MI1  0.8860
29        29       MI2  0.1865
30        30       MI3  1.7326
...................................
Regression Coefficients
         [,1]
[1,]  0.00000
[2,] -0.02654

Variance:
[1] 1.454 1.238


EAP Reliability:
[1] 0.873

-----------------------------
Start:  2019-01-06 08:14:55
End:  2019-01-06 08:14:56 
Time difference of 0.4770472 secs

------------------------------------------------------------
TAM 2.13-15 (2018-09-30 14:35:48) 
sh: 1: cannot create /dev/null: Permission denied
R version 3.4.4 (2018-03-15) x86_64, linux-gnu | nodename=ip-172-31-1-5 | login=unknown 

Date of Analysis: 2019-01-06 08:14:56 
Time difference of 0.4770472 secs
Computation time: 0.4770472 

Multidimensional Item Response Model in TAM 

IRT Model 1PL

Call:
TAM::tam.mml(resp = dat[, items], group = dat$female)

------------------------------------------------------------
Number of iterations= 37 
Numeric integration with 21 integration points

Deviance= 22728.35 
   Log likelihood= -11364.17 
Number of persons= 664 
Number of persons used= 664 
Number of items= 30 
Number of estimated parameters= 33 
    Item threshold parameters= 30 
    Item slope parameters  = 0 
    Regression parameters  = 1 
    (Co)Variance parameters= 2 

AIC = 22794  | penalty= 66    | AIC=-2*LL + 2*p 
AIC3 = 22827  | penalty= 99    | AIC3=-2*LL + 3*p 
BIC = 22943  | penalty= 214.44    | BIC=-2*LL + log(n)*p 
aBIC = 22838  | penalty= 109.47    | aBIC=-2*LL + log((n-2)/24)*p  (adjusted BIC) 
CAIC = 22976  | penalty= 247.44    | CAIC=-2*LL + [log(n)+1]*p  (consistent AIC) 
AICc = 22798  | penalty= 69.56    | AICc=-2*LL + 2*p + 2*p*(p+1)/(n-p-1)  (bias corrected AIC) 

------------------------------------------------------------
EAP Reliability
[1] 0.873
------------------------------------------------------------
Covariances and Variances
Group0 Group1 
 1.454  1.238 
------------------------------------------------------------
Correlations and Standard Deviations (in the diagonal)
Group0 Group1 
 1.206  1.113 
------------------------------------------------------------
Regression Coefficients
         [,1]
[1,]  0.00000
[2,] -0.02654
------------------------------------------------------------
Item Parameters -A*Xsi
    item   N     M xsi.item AXsi_.Cat1 B.Cat1.Dim1
MA1  MA1 664 0.886   -2.522     -2.522           1
MA2  MA2 664 0.562   -0.318     -0.318           1
MA3  MA3 664 0.526   -0.134     -0.134           1
MA4  MA4 664 0.404    0.488      0.488           1
MB1  MB1 664 0.634   -0.699     -0.699           1
MB2  MB2 664 0.616   -0.602     -0.602           1
MB3  MB3 664 0.349    0.778      0.778           1
MB4  MB4 664 0.352    0.761      0.761           1
MC1  MC1 664 0.767   -1.496     -1.496           1
MC2  MC2 664 0.480    0.095      0.095           1
MC3  MC3 664 0.663   -0.857     -0.857           1
MC4  MC4 664 0.554   -0.280     -0.280           1
MD1  MD1 664 0.682   -0.969     -0.969           1
MD2  MD2 664 0.750   -1.385     -1.385           1
MD3  MD3 664 0.604   -0.538     -0.538           1
MD4  MD4 664 0.605   -0.546     -0.546           1
ME1  ME1 664 0.697   -1.057     -1.057           1
ME2  ME2 664 0.393    0.543      0.543           1
MF1  MF1 664 0.568   -0.349     -0.349           1
MG1  MG1 664 0.420    0.402      0.402           1
MG2  MG2 664 0.342    0.819      0.819           1
MG3  MG3 664 0.435    0.325      0.325           1
MG4  MG4 664 0.404    0.488      0.488           1
MH1  MH1 664 0.497    0.011      0.011           1
MH2  MH2 664 0.370    0.663      0.663           1
MH3  MH3 664 0.425    0.379      0.379           1
MH4  MH4 664 0.333    0.869      0.869           1
MI1  MI1 664 0.330    0.886      0.886           1
MI2  MI2 664 0.462    0.187      0.187           1
MI3  MI3 664 0.197    1.733      1.733           1
-----------------------------------------------------------------------------
CDM 7.0-12 (2018-11-06 16:23:12) 

Call:
CDM::IRT.RMSD(object = mod)

-----------------------------------------------------------------------------
Root Mean Square Deviation (RMSD) 

    Parm     M    SD   Min   Max
1 Group1 0.048 0.023 0.021 0.129
2 Group2 0.053 0.015 0.025 0.093
3  WRMSD 0.052 0.015 0.033 0.098

   item Group1 Group2 WRMSD
1   MA1  0.036  0.069 0.055
2   MA2  0.032  0.061 0.048
3   MA3  0.034  0.034 0.034
4   MA4  0.049  0.060 0.055
5   MB1  0.079  0.063 0.072
6   MB2  0.056  0.065 0.060
7   MB3  0.049  0.040 0.045
8   MB4  0.060  0.053 0.056
9   MC1  0.022  0.055 0.041
10  MC2  0.026  0.048 0.038
11  MC3  0.047  0.025 0.038
12  MC4  0.046  0.042 0.044
13  MD1  0.037  0.049 0.043
14  MD2  0.038  0.039 0.039
15  MD3  0.028  0.041 0.035
16  MD4  0.021  0.041 0.033
17  ME1  0.031  0.047 0.040
18  ME2  0.058  0.067 0.063
19  MF1  0.042  0.048 0.045
20  MG1  0.037  0.054 0.046
21  MG2  0.063  0.072 0.067
22  MG3  0.035  0.070 0.055
23  MG4  0.029  0.076 0.057
24  MH1  0.129  0.045 0.098
25  MH2  0.049  0.028 0.040
26  MH3  0.074  0.047 0.062
27  MH4  0.056  0.034 0.047
28  MI1  0.085  0.093 0.089
29  MI2  0.063  0.057 0.060
30  MI3  0.024  0.057 0.043
-----------------------------------------------------------------------------
Bias Corrected Root Mean Square Deviation (RMSD) 

    Parm     M    SD    Min   Max
1 Group1 0.040 0.027 -0.010 0.127
2 Group2 0.046 0.018  0.009 0.090
3  WRMSD 0.046 0.017  0.024 0.095

   item Group1 Group2 WRMSD
1   MA1  0.033  0.067 0.052
2   MA2  0.022  0.056 0.042
3   MA3  0.024  0.024 0.024
4   MA4  0.043  0.055 0.049
5   MB1  0.075  0.058 0.068
6   MB2  0.051  0.060 0.055
7   MB3  0.043  0.032 0.038
8   MB4  0.055  0.047 0.051
9   MC1  0.007  0.051 0.036
10  MC2  0.009  0.041 0.029
11  MC3  0.042  0.009 0.030
12  MC4  0.039  0.034 0.037
13  MD1  0.029  0.043 0.037
14  MD2  0.032  0.033 0.032
15  MD3  0.015  0.033 0.026
16  MD4 -0.010  0.033 0.024
17  ME1  0.022  0.041 0.033
18  ME2  0.053  0.062 0.058
19  MF1  0.034  0.041 0.038
20  MG1  0.028  0.048 0.039
21  MG2  0.058  0.068 0.063
22  MG3  0.026  0.066 0.049
23  MG4  0.018  0.071 0.051
24  MH1  0.127  0.037 0.095
25  MH2  0.044  0.014 0.033
26  MH3  0.070  0.040 0.058
27  MH4  0.051  0.024 0.040
28  MI1  0.081  0.090 0.086
29  MI2  0.058  0.052 0.055
30  MI3  0.013  0.053 0.038
-----------------------------------------------------------------------------
Mean Absolute Deviation (MAD) 

    Parm     M    SD   Min   Max
1 Group1 0.040 0.021 0.016 0.107
2 Group2 0.044 0.014 0.023 0.077

   item Group1 Group2
1   MA1  0.020  0.041
2   MA2  0.027  0.048
3   MA3  0.029  0.024
4   MA4  0.041  0.053
5   MB1  0.073  0.059
6   MB2  0.048  0.050
7   MB3  0.041  0.034
8   MB4  0.050  0.051
9   MC1  0.016  0.044
10  MC2  0.020  0.042
11  MC3  0.037  0.023
12  MC4  0.032  0.036
13  MD1  0.027  0.042
14  MD2  0.034  0.026
15  MD3  0.020  0.032
16  MD4  0.018  0.035
17  ME1  0.027  0.041
18  ME2  0.053  0.053
19  MF1  0.029  0.040
20  MG1  0.029  0.052
21  MG2  0.052  0.059
22  MG3  0.034  0.062
23  MG4  0.022  0.073
24  MH1  0.107  0.035
25  MH2  0.042  0.026
26  MH3  0.066  0.038
27  MH4  0.047  0.029
28  MI1  0.078  0.077
29  MI2  0.052  0.049
30  MI3  0.019  0.047
-----------------------------------------------------------------------------
Mean Deviation (MD) 

    Parm M    SD    Min   Max
1 Group1 0 0.023 -0.056 0.046
2 Group2 0 0.024 -0.049 0.059

   item Group1 Group2
1   MA1 -0.001  0.001
2   MA2  0.001 -0.001
3   MA3  0.016 -0.017
4   MA4 -0.036  0.038
5   MB1 -0.056  0.059
6   MB2 -0.044  0.047
7   MB3  0.012 -0.013
8   MB4  0.015 -0.016
9   MC1  0.006 -0.007
10  MC2 -0.015  0.016
11  MC3  0.007 -0.007
12  MC4 -0.027  0.028
13  MD1 -0.010  0.010
14  MD2  0.017 -0.018
15  MD3 -0.009  0.009
16  MD4  0.008 -0.008
17  ME1  0.010 -0.011
18  ME2  0.036 -0.038
19  MF1 -0.025  0.027
20  MG1  0.021 -0.022
21  MG2  0.046 -0.049
22  MG3 -0.012  0.012
23  MG4  0.008 -0.009
24  MH1  0.021 -0.022
25  MH2  0.003 -0.003
26  MH3 -0.030  0.032
27  MH4  0.014 -0.015
28  MI1  0.014 -0.015
29  MI2  0.006 -0.006
30  MI3  0.006 -0.006
---------------------------------------------------------------------------------
CDM 7.0-12 (2018-11-06 16:23:12) 
GDINA Model 
 Link function: identity 
  ** 2019-01-06 08:14:56 
---------------------------------------------------------------------------------
...........................................................
Iteration 1     2019-01-06 08:14:56 
Deviance = 4390.371
Maximum parameter change: 0.191603 
...........................................................
Iteration 2     2019-01-06 08:14:56 
Deviance = 4120.781 | Deviance change = 269.5906
Maximum parameter change: 0.084701 
...........................................................
Iteration 3     2019-01-06 08:14:56 
Deviance = 4098.962 | Deviance change = 21.81842
Maximum parameter change: 0.088579 
...........................................................
Iteration 4     2019-01-06 08:14:56 
Deviance = 4087.43 | Deviance change = 11.53259
Maximum parameter change: 0.006933 
...........................................................
Iteration 5     2019-01-06 08:14:56 
Deviance = 4086.995 | Deviance change = 0.434914
Maximum parameter change: 0.017085 
...........................................................
Iteration 6     2019-01-06 08:14:56 
Deviance = 4086.637 | Deviance change = 0.3572192
Maximum parameter change: 0.006592 
...........................................................
Iteration 7     2019-01-06 08:14:56 
Deviance = 4086.461 | Deviance change = 0.176217
Maximum parameter change: 0.005834 
...........................................................
Iteration 8     2019-01-06 08:14:56 
Deviance = 4086.377 | Deviance change = 0.0841075
Maximum parameter change: 0.003394 
...........................................................
Iteration 9     2019-01-06 08:14:56 
Deviance = 4086.235 | Deviance change = 0.1420471
Maximum parameter change: 0.014653 
...........................................................
Iteration 10     2019-01-06 08:14:56 
Deviance = 4086.072 | Deviance change = 0.1635339
Maximum parameter change: 0.003625 
...........................................................
Iteration 11     2019-01-06 08:14:56 
Deviance = 4086.953 | Deviance change = -0.8810047
**** Deviances decreases! Check for nonconvergence.   ****
Maximum parameter change: 0.05578 
...........................................................
Iteration 12     2019-01-06 08:14:56 
Deviance = 4085.627 | Deviance change = 1.325659
Maximum parameter change: 0.005208 
...........................................................
Iteration 13     2019-01-06 08:14:56 
Deviance = 4085.623 | Deviance change = 0.003516
Maximum parameter change: 0.015966 
...........................................................
Iteration 14     2019-01-06 08:14:56 
Deviance = 4085.463 | Deviance change = 0.1600612
Maximum parameter change: 0.003566 
...........................................................
Iteration 15     2019-01-06 08:14:56 
Deviance = 4085.42 | Deviance change = 0.0435283
Maximum parameter change: 0.011854 
...........................................................
Iteration 16     2019-01-06 08:14:56 
Deviance = 4085.338 | Deviance change = 0.0820725
Maximum parameter change: 0.002237 
...........................................................
Iteration 17     2019-01-06 08:14:56 
Deviance = 4085.346 | Deviance change = -0.0079826
**** Deviances decreases! Check for nonconvergence.   ****
Maximum parameter change: 0.020078 
...........................................................
Iteration 18     2019-01-06 08:14:56 
Deviance = 4085.206 | Deviance change = 0.1395725
Maximum parameter change: 0.002048 
...........................................................
Iteration 19     2019-01-06 08:14:56 
Deviance = 4085.178 | Deviance change = 0.0281615
Maximum parameter change: 0.002098 
...........................................................
Iteration 20     2019-01-06 08:14:56 
Deviance = 4085.159 | Deviance change = 0.0189789
Maximum parameter change: 0.002161 
...........................................................
Iteration 21     2019-01-06 08:14:56 
Deviance = 4085.144 | Deviance change = 0.0149187
Maximum parameter change: 0.002196 
...........................................................
Iteration 22     2019-01-06 08:14:56 
Deviance = 4085.131 | Deviance change = 0.0129681
Maximum parameter change: 0.00221 
...........................................................
Iteration 23     2019-01-06 08:14:56 
Deviance = 4085.119 | Deviance change = 0.0119553
Maximum parameter change: 0.00221 
...........................................................
Iteration 24     2019-01-06 08:14:56 
Deviance = 4085.108 | Deviance change = 0.0113646
Maximum parameter change: 0.0022 
...........................................................
Iteration 25     2019-01-06 08:14:56 
Deviance = 4085.097 | Deviance change = 0.0109661
Maximum parameter change: 0.002184 
...........................................................
Iteration 26     2019-01-06 08:14:56 
Deviance = 4085.086 | Deviance change = 0.0106571
Maximum parameter change: 0.002163 
...........................................................
Iteration 27     2019-01-06 08:14:56 
Deviance = 4085.076 | Deviance change = 0.0103913
Maximum parameter change: 0.002139 
...........................................................
Iteration 28     2019-01-06 08:14:56 
Deviance = 4085.066 | Deviance change = 0.0101481
Maximum parameter change: 0.002113 
...........................................................
Iteration 29     2019-01-06 08:14:56 
Deviance = 4085.056 | Deviance change = 0.0099183
Maximum parameter change: 0.002087 
...........................................................
Iteration 30     2019-01-06 08:14:56 
Deviance = 4085.046 | Deviance change = 0.0096977
Maximum parameter change: 0.002059 
...........................................................
Iteration 31     2019-01-06 08:14:56 
Deviance = 4085.037 | Deviance change = 0.0094841
Maximum parameter change: 0.002032 
...........................................................
Iteration 32     2019-01-06 08:14:56 
Deviance = 4085.027 | Deviance change = 0.0092768
Maximum parameter change: 0.002004 
...........................................................
Iteration 33     2019-01-06 08:14:56 
Deviance = 4085.018 | Deviance change = 0.0090752
Maximum parameter change: 0.001976 
...........................................................
Iteration 34     2019-01-06 08:14:56 
Deviance = 4085.009 | Deviance change = 0.0088788
Maximum parameter change: 0.001948 
...........................................................
Iteration 35     2019-01-06 08:14:56 
Deviance = 4085.001 | Deviance change = 0.0086876
Maximum parameter change: 0.001921 
...........................................................
Iteration 36     2019-01-06 08:14:56 
Deviance = 4084.992 | Deviance change = 0.0085013
Maximum parameter change: 0.001894 
...........................................................
Iteration 37     2019-01-06 08:14:56 
Deviance = 4084.984 | Deviance change = 0.0083196
Maximum parameter change: 0.001867 
...........................................................
Iteration 38     2019-01-06 08:14:56 
Deviance = 4084.976 | Deviance change = 0.0081424
Maximum parameter change: 0.00184 
...........................................................
Iteration 39     2019-01-06 08:14:56 
Deviance = 4084.968 | Deviance change = 0.0079696
Maximum parameter change: 0.001814 
...........................................................
Iteration 40     2019-01-06 08:14:56 
Deviance = 4084.96 | Deviance change = 0.007801
Maximum parameter change: 0.001787 
...........................................................
Iteration 41     2019-01-06 08:14:56 
Deviance = 4084.952 | Deviance change = 0.0076364
Maximum parameter change: 0.001761 
...........................................................
Iteration 42     2019-01-06 08:14:56 
Deviance = 4084.945 | Deviance change = 0.0074757
Maximum parameter change: 0.001736 
...........................................................
Iteration 43     2019-01-06 08:14:56 
Deviance = 4084.937 | Deviance change = 0.0073188
Maximum parameter change: 0.00171 
...........................................................
Iteration 44     2019-01-06 08:14:56 
Deviance = 4084.93 | Deviance change = 0.0071656
Maximum parameter change: 0.001685 
...........................................................
Iteration 45     2019-01-06 08:14:56 
Deviance = 4084.923 | Deviance change = 0.0070159
Maximum parameter change: 0.001661 
...........................................................
Iteration 46     2019-01-06 08:14:56 
Deviance = 4084.916 | Deviance change = 0.0068697
Maximum parameter change: 0.001636 
...........................................................
Iteration 47     2019-01-06 08:14:56 
Deviance = 4084.91 | Deviance change = 0.0067268
Maximum parameter change: 0.001612 
...........................................................
Iteration 48     2019-01-06 08:14:56 
Deviance = 4084.903 | Deviance change = 0.0065871
Maximum parameter change: 0.001588 
...........................................................
Iteration 49     2019-01-06 08:14:56 
Deviance = 4084.897 | Deviance change = 0.0064506
Maximum parameter change: 0.001564 
...........................................................
Iteration 50     2019-01-06 08:14:56 
Deviance = 4084.89 | Deviance change = 0.0063171
Maximum parameter change: 0.001541 
...........................................................
Iteration 51     2019-01-06 08:14:56 
Deviance = 4084.884 | Deviance change = 0.0061867
Maximum parameter change: 0.001518 
...........................................................
Iteration 52     2019-01-06 08:14:56 
Deviance = 4084.878 | Deviance change = 0.0060591
Maximum parameter change: 0.001495 
...........................................................
Iteration 53     2019-01-06 08:14:56 
Deviance = 4084.872 | Deviance change = 0.0059343
Maximum parameter change: 0.001473 
...........................................................
Iteration 54     2019-01-06 08:14:56 
Deviance = 4084.866 | Deviance change = 0.0058124
Maximum parameter change: 0.00145 
...........................................................
Iteration 55     2019-01-06 08:14:56 
Deviance = 4084.861 | Deviance change = 0.0056931
Maximum parameter change: 0.001428 
...........................................................
Iteration 56     2019-01-06 08:14:56 
Deviance = 4084.855 | Deviance change = 0.0055764
Maximum parameter change: 0.001407 
...........................................................
Iteration 57     2019-01-06 08:14:56 
Deviance = 4084.85 | Deviance change = 0.0054622
Maximum parameter change: 0.001385 
...........................................................
Iteration 58     2019-01-06 08:14:56 
Deviance = 4084.844 | Deviance change = 0.0053506
Maximum parameter change: 0.001364 
...........................................................
Iteration 59     2019-01-06 08:14:56 
Deviance = 4084.839 | Deviance change = 0.0052414
Maximum parameter change: 0.001343 
...........................................................
Iteration 60     2019-01-06 08:14:56 
Deviance = 4084.834 | Deviance change = 0.0051346
Maximum parameter change: 0.001323 
...........................................................
Iteration 61     2019-01-06 08:14:56 
Deviance = 4084.829 | Deviance change = 0.0050301
Maximum parameter change: 0.001302 
...........................................................
Iteration 62     2019-01-06 08:14:56 
Deviance = 4084.824 | Deviance change = 0.0049278
Maximum parameter change: 0.001282 
...........................................................
Iteration 63     2019-01-06 08:14:56 
Deviance = 4084.819 | Deviance change = 0.0048278
Maximum parameter change: 0.001262 
...........................................................
Iteration 64     2019-01-06 08:14:56 
Deviance = 4084.814 | Deviance change = 0.0047299
Maximum parameter change: 0.001243 
...........................................................
Iteration 65     2019-01-06 08:14:56 
Deviance = 4084.81 | Deviance change = 0.0046341
Maximum parameter change: 0.001224 
...........................................................
Iteration 66     2019-01-06 08:14:56 
Deviance = 4084.805 | Deviance change = 0.0045404
Maximum parameter change: 0.001205 
...........................................................
Iteration 67     2019-01-06 08:14:56 
Deviance = 4084.801 | Deviance change = 0.0044488
Maximum parameter change: 0.001186 
...........................................................
Iteration 68     2019-01-06 08:14:56 
Deviance = 4084.796 | Deviance change = 0.004359
Maximum parameter change: 0.001167 
...........................................................
Iteration 69     2019-01-06 08:14:56 
Deviance = 4084.792 | Deviance change = 0.0042713
Maximum parameter change: 0.001149 
...........................................................
Iteration 70     2019-01-06 08:14:56 
Deviance = 4084.788 | Deviance change = 0.0041854
Maximum parameter change: 0.001131 
...........................................................
Iteration 71     2019-01-06 08:14:56 
Deviance = 4084.784 | Deviance change = 0.0041013
Maximum parameter change: 0.001113 
...........................................................
Iteration 72     2019-01-06 08:14:56 
Deviance = 4084.78 | Deviance change = 0.004019
Maximum parameter change: 0.001096 
...........................................................
Iteration 73     2019-01-06 08:14:56 
Deviance = 4084.776 | Deviance change = 0.0039385
Maximum parameter change: 0.001078 
...........................................................
Iteration 74     2019-01-06 08:14:56 
Deviance = 4084.772 | Deviance change = 0.0038597
Maximum parameter change: 0.001061 
...........................................................
Iteration 75     2019-01-06 08:14:56 
Deviance = 4084.768 | Deviance change = 0.0037826
Maximum parameter change: 0.001045 
...........................................................
Iteration 76     2019-01-06 08:14:56 
Deviance = 4084.765 | Deviance change = 0.0037072
Maximum parameter change: 0.001028 
...........................................................
Iteration 77     2019-01-06 08:14:56 
Deviance = 4084.761 | Deviance change = 0.0036333
Maximum parameter change: 0.001012 
...........................................................
Iteration 78     2019-01-06 08:14:56 
Deviance = 4084.757 | Deviance change = 0.003561
Maximum parameter change: 0.000996 
...........................................................
Iteration 79     2019-01-06 08:14:56 
Deviance = 4084.754 | Deviance change = 0.0034903
Maximum parameter change: 0.00098 
...........................................................
Iteration 80     2019-01-06 08:14:56 
Deviance = 4084.75 | Deviance change = 0.0034211
Maximum parameter change: 0.000964 
...........................................................
Iteration 81     2019-01-06 08:14:56 
Deviance = 4084.747 | Deviance change = 0.0033533
Maximum parameter change: 0.000949 
...........................................................
Iteration 82     2019-01-06 08:14:56 
Deviance = 4084.744 | Deviance change = 0.003287
Maximum parameter change: 0.000934 
...........................................................
Iteration 83     2019-01-06 08:14:56 
Deviance = 4084.741 | Deviance change = 0.003222
Maximum parameter change: 0.000919 
...........................................................
Iteration 84     2019-01-06 08:14:56 
Deviance = 4084.737 | Deviance change = 0.0031585
Maximum parameter change: 0.000904 
...........................................................
Iteration 85     2019-01-06 08:14:56 
Deviance = 4084.734 | Deviance change = 0.0030963
Maximum parameter change: 0.00089 
...........................................................
Iteration 86     2019-01-06 08:14:56 
Deviance = 4084.731 | Deviance change = 0.0030354
Maximum parameter change: 0.000875 
...........................................................
Iteration 87     2019-01-06 08:14:56 
Deviance = 4084.728 | Deviance change = 0.0029758
Maximum parameter change: 0.000861 
...........................................................
Iteration 88     2019-01-06 08:14:56 
Deviance = 4084.725 | Deviance change = 0.0029174
Maximum parameter change: 0.000847 
...........................................................
Iteration 89     2019-01-06 08:14:56 
Deviance = 4084.723 | Deviance change = 0.0028603
Maximum parameter change: 0.000834 
...........................................................
Iteration 90     2019-01-06 08:14:56 
Deviance = 4084.72 | Deviance change = 0.0028044
Maximum parameter change: 0.00082 
...........................................................
Iteration 91     2019-01-06 08:14:56 
Deviance = 4084.717 | Deviance change = 0.0027496
Maximum parameter change: 0.000807 
...........................................................
Iteration 92     2019-01-06 08:14:56 
Deviance = 4084.714 | Deviance change = 0.002696
Maximum parameter change: 0.000794 
...........................................................
Iteration 93     2019-01-06 08:14:56 
Deviance = 4084.712 | Deviance change = 0.0026435
Maximum parameter change: 0.000781 
...........................................................
Iteration 94     2019-01-06 08:14:56 
Deviance = 4084.709 | Deviance change = 0.0025921
Maximum parameter change: 0.000768 
...........................................................
Iteration 95     2019-01-06 08:14:56 
Deviance = 4084.706 | Deviance change = 0.0025418
Maximum parameter change: 0.000756 
...........................................................
Iteration 96     2019-01-06 08:14:56 
Deviance = 4084.704 | Deviance change = 0.0024925
Maximum parameter change: 0.000744 
...........................................................
Iteration 97     2019-01-06 08:14:56 
Deviance = 4084.702 | Deviance change = 0.0024443
Maximum parameter change: 0.000732 
...........................................................
Iteration 98     2019-01-06 08:14:56 
Deviance = 4084.699 | Deviance change = 0.0023971
Maximum parameter change: 0.00072 
...........................................................
Iteration 99     2019-01-06 08:14:56 
Deviance = 4084.697 | Deviance change = 0.0023508
Maximum parameter change: 0.000708 
...........................................................
Iteration 100     2019-01-06 08:14:56 
Deviance = 4084.694 | Deviance change = 0.0023055
Maximum parameter change: 0.000696 
...........................................................
Iteration 101     2019-01-06 08:14:56 
Deviance = 4084.692 | Deviance change = 0.0022612
Maximum parameter change: 0.000685 
...........................................................
Iteration 102     2019-01-06 08:14:56 
Deviance = 4084.69 | Deviance change = 0.0022177
Maximum parameter change: 0.000674 
...........................................................
Iteration 103     2019-01-06 08:14:56 
Deviance = 4084.688 | Deviance change = 0.0021752
Maximum parameter change: 0.000663 
...........................................................
Iteration 104     2019-01-06 08:14:56 
Deviance = 4084.686 | Deviance change = 0.0021335
Maximum parameter change: 0.000652 
...........................................................
Iteration 105     2019-01-06 08:14:56 
Deviance = 4084.684 | Deviance change = 0.0020927
Maximum parameter change: 0.000641 
...........................................................
Iteration 106     2019-01-06 08:14:56 
Deviance = 4084.682 | Deviance change = 0.0020527
Maximum parameter change: 0.000631 
...........................................................
Iteration 107     2019-01-06 08:14:56 
Deviance = 4084.68 | Deviance change = 0.0020136
Maximum parameter change: 0.000621 
...........................................................
Iteration 108     2019-01-06 08:14:56 
Deviance = 4084.678 | Deviance change = 0.0019752
Maximum parameter change: 0.00061 
...........................................................
Iteration 109     2019-01-06 08:14:56 
Deviance = 4084.676 | Deviance change = 0.0019377
Maximum parameter change: 6e-04 
...........................................................
Iteration 110     2019-01-06 08:14:56 
Deviance = 4084.674 | Deviance change = 0.0019009
Maximum parameter change: 0.00059 
...........................................................
Iteration 111     2019-01-06 08:14:56 
Deviance = 4084.672 | Deviance change = 0.0018648
Maximum parameter change: 0.000581 
...........................................................
Iteration 112     2019-01-06 08:14:56 
Deviance = 4084.67 | Deviance change = 0.0018295
Maximum parameter change: 0.000571 
...........................................................
Iteration 113     2019-01-06 08:14:56 
Deviance = 4084.668 | Deviance change = 0.0017949
Maximum parameter change: 0.000562 
...........................................................
Iteration 114     2019-01-06 08:14:56 
Deviance = 4084.666 | Deviance change = 0.001761
Maximum parameter change: 0.000553 
...........................................................
Iteration 115     2019-01-06 08:14:56 
Deviance = 4084.665 | Deviance change = 0.0017278
Maximum parameter change: 0.000543 
...........................................................
Iteration 116     2019-01-06 08:14:56 
Deviance = 4084.663 | Deviance change = 0.0016953
Maximum parameter change: 0.000534 
...........................................................
Iteration 117     2019-01-06 08:14:56 
Deviance = 4084.661 | Deviance change = 0.0016634
Maximum parameter change: 0.000526 
...........................................................
Iteration 118     2019-01-06 08:14:56 
Deviance = 4084.66 | Deviance change = 0.0016322
Maximum parameter change: 0.000517 
...........................................................
Iteration 119     2019-01-06 08:14:56 
Deviance = 4084.658 | Deviance change = 0.0016015
Maximum parameter change: 0.000508 
...........................................................
Iteration 120     2019-01-06 08:14:56 
Deviance = 4084.657 | Deviance change = 0.0015715
Maximum parameter change: 5e-04 
...........................................................
Iteration 121     2019-01-06 08:14:56 
Deviance = 4084.655 | Deviance change = 0.0015422
Maximum parameter change: 0.000492 
...........................................................
Iteration 122     2019-01-06 08:14:56 
Deviance = 4084.654 | Deviance change = 0.0015133
Maximum parameter change: 0.000484 
...........................................................
Iteration 123     2019-01-06 08:14:56 
Deviance = 4084.652 | Deviance change = 0.0014851
Maximum parameter change: 0.000476 
...........................................................
Iteration 124     2019-01-06 08:14:56 
Deviance = 4084.651 | Deviance change = 0.0014574
Maximum parameter change: 0.000468 
...........................................................
Iteration 125     2019-01-06 08:14:56 
Deviance = 4084.649 | Deviance change = 0.0014303
Maximum parameter change: 0.00046 
...........................................................
Iteration 126     2019-01-06 08:14:56 
Deviance = 4084.648 | Deviance change = 0.0014038
Maximum parameter change: 0.000452 
...........................................................
Iteration 127     2019-01-06 08:14:56 
Deviance = 4084.646 | Deviance change = 0.0013777
Maximum parameter change: 0.000445 
...........................................................
Iteration 128     2019-01-06 08:14:56 
Deviance = 4084.645 | Deviance change = 0.0013522
Maximum parameter change: 0.000437 
...........................................................
Iteration 129     2019-01-06 08:14:56 
Deviance = 4084.644 | Deviance change = 0.0013271
Maximum parameter change: 0.00043 
...........................................................
Iteration 130     2019-01-06 08:14:56 
Deviance = 4084.642 | Deviance change = 0.0013026
Maximum parameter change: 0.000423 
...........................................................
Iteration 131     2019-01-06 08:14:56 
Deviance = 4084.641 | Deviance change = 0.0012785
Maximum parameter change: 0.000416 
...........................................................
Iteration 132     2019-01-06 08:14:56 
Deviance = 4084.64 | Deviance change = 0.001255
Maximum parameter change: 0.000409 
...........................................................
Iteration 133     2019-01-06 08:14:56 
Deviance = 4084.639 | Deviance change = 0.0012318
Maximum parameter change: 0.000402 
...........................................................
Iteration 134     2019-01-06 08:14:56 
Deviance = 4084.637 | Deviance change = 0.0012092
Maximum parameter change: 0.000396 
...........................................................
Iteration 135     2019-01-06 08:14:56 
Deviance = 4084.636 | Deviance change = 0.001187
Maximum parameter change: 0.000389 
...........................................................
Iteration 136     2019-01-06 08:14:56 
Deviance = 4084.635 | Deviance change = 0.0011652
Maximum parameter change: 0.000383 
...........................................................
Iteration 137     2019-01-06 08:14:56 
Deviance = 4084.634 | Deviance change = 0.0011438
Maximum parameter change: 0.000376 
...........................................................
Iteration 138     2019-01-06 08:14:56 
Deviance = 4084.633 | Deviance change = 0.0011229
Maximum parameter change: 0.00037 
...........................................................
Iteration 139     2019-01-06 08:14:56 
Deviance = 4084.632 | Deviance change = 0.0011023
Maximum parameter change: 0.000364 
...........................................................
Iteration 140     2019-01-06 08:14:56 
Deviance = 4084.631 | Deviance change = 0.0010822
Maximum parameter change: 0.000358 
...........................................................
Iteration 141     2019-01-06 08:14:56 
Deviance = 4084.63 | Deviance change = 0.0010624
Maximum parameter change: 0.000352 
...........................................................
Iteration 142     2019-01-06 08:14:56 
Deviance = 4084.629 | Deviance change = 0.0010431
Maximum parameter change: 0.000346 
...........................................................
Iteration 143     2019-01-06 08:14:56 
Deviance = 4084.627 | Deviance change = 0.0010241
Maximum parameter change: 0.00034 
...........................................................
Iteration 144     2019-01-06 08:14:56 
Deviance = 4084.626 | Deviance change = 0.0010054
Maximum parameter change: 0.000334 
...........................................................
Iteration 145     2019-01-06 08:14:56 
Deviance = 4084.626 | Deviance change = 0.0009872
Maximum parameter change: 0.000329 
...........................................................
Iteration 146     2019-01-06 08:14:56 
Deviance = 4084.625 | Deviance change = 0.0009693
Maximum parameter change: 0.000323 
...........................................................
Iteration 147     2019-01-06 08:14:56 
Deviance = 4084.624 | Deviance change = 0.0009517
Maximum parameter change: 0.000318 
...........................................................
Iteration 148     2019-01-06 08:14:56 
Deviance = 4084.623 | Deviance change = 0.0009344
Maximum parameter change: 0.000313 
...........................................................
Iteration 149     2019-01-06 08:14:56 
Deviance = 4084.622 | Deviance change = 0.0009175
Maximum parameter change: 0.000307 
...........................................................
Iteration 150     2019-01-06 08:14:56 
Deviance = 4084.621 | Deviance change = 0.0009009
Maximum parameter change: 0.000302 
...........................................................
Iteration 151     2019-01-06 08:14:56 
Deviance = 4084.62 | Deviance change = 0.0008847
Maximum parameter change: 0.000297 
...........................................................
Iteration 152     2019-01-06 08:14:56 
Deviance = 4084.619 | Deviance change = 0.0008687
Maximum parameter change: 0.000292 
...........................................................
Iteration 153     2019-01-06 08:14:56 
Deviance = 4084.618 | Deviance change = 0.0008531
Maximum parameter change: 0.000287 
...........................................................
Iteration 154     2019-01-06 08:14:56 
Deviance = 4084.617 | Deviance change = 0.0008377
Maximum parameter change: 0.000283 
...........................................................
Iteration 155     2019-01-06 08:14:56 
Deviance = 4084.617 | Deviance change = 0.0008226
Maximum parameter change: 0.000278 
...........................................................
Iteration 156     2019-01-06 08:14:56 
Deviance = 4084.616 | Deviance change = 0.0008078
Maximum parameter change: 0.000273 
...........................................................
Iteration 157     2019-01-06 08:14:56 
Deviance = 4084.615 | Deviance change = 0.0007933
Maximum parameter change: 0.000269 
...........................................................
Iteration 158     2019-01-06 08:14:56 
Deviance = 4084.614 | Deviance change = 0.0007791
Maximum parameter change: 0.000264 
...........................................................
Iteration 159     2019-01-06 08:14:56 
Deviance = 4084.613 | Deviance change = 0.0007651
Maximum parameter change: 0.00026 
...........................................................
Iteration 160     2019-01-06 08:14:56 
Deviance = 4084.613 | Deviance change = 0.0007514
Maximum parameter change: 0.000255 
...........................................................
Iteration 161     2019-01-06 08:14:56 
Deviance = 4084.612 | Deviance change = 0.000738
Maximum parameter change: 0.000251 
...........................................................
Iteration 162     2019-01-06 08:14:56 
Deviance = 4084.611 | Deviance change = 0.0007248
Maximum parameter change: 0.000247 
...........................................................
Iteration 163     2019-01-06 08:14:56 
Deviance = 4084.61 | Deviance change = 0.0007119
Maximum parameter change: 0.000243 
...........................................................
Iteration 164     2019-01-06 08:14:56 
Deviance = 4084.61 | Deviance change = 0.0006992
Maximum parameter change: 0.000239 
...........................................................
Iteration 165     2019-01-06 08:14:56 
Deviance = 4084.609 | Deviance change = 0.0006867
Maximum parameter change: 0.000235 
...........................................................
Iteration 166     2019-01-06 08:14:56 
Deviance = 4084.608 | Deviance change = 0.0006745
Maximum parameter change: 0.000231 
...........................................................
Iteration 167     2019-01-06 08:14:56 
Deviance = 4084.608 | Deviance change = 0.0006625
Maximum parameter change: 0.000227 
...........................................................
Iteration 168     2019-01-06 08:14:56 
Deviance = 4084.607 | Deviance change = 0.0006507
Maximum parameter change: 0.000223 
...........................................................
Iteration 169     2019-01-06 08:14:56 
Deviance = 4084.606 | Deviance change = 0.0006391
Maximum parameter change: 0.000219 
...........................................................
Iteration 170     2019-01-06 08:14:56 
Deviance = 4084.606 | Deviance change = 0.0006278
Maximum parameter change: 0.000216 
...........................................................
Iteration 171     2019-01-06 08:14:56 
Deviance = 4084.605 | Deviance change = 0.0006166
Maximum parameter change: 0.000212 
...........................................................
Iteration 172     2019-01-06 08:14:56 
Deviance = 4084.605 | Deviance change = 0.0006057
Maximum parameter change: 0.000208 
...........................................................
Iteration 173     2019-01-06 08:14:56 
Deviance = 4084.604 | Deviance change = 0.000595
Maximum parameter change: 0.000205 
...........................................................
Iteration 174     2019-01-06 08:14:56 
Deviance = 4084.603 | Deviance change = 0.0005844
Maximum parameter change: 0.000201 
...........................................................
Iteration 175     2019-01-06 08:14:56 
Deviance = 4084.603 | Deviance change = 0.0005741
Maximum parameter change: 0.000198 
...........................................................
Iteration 176     2019-01-06 08:14:56 
Deviance = 4084.602 | Deviance change = 0.0005639
Maximum parameter change: 0.000195 
...........................................................
Iteration 177     2019-01-06 08:14:56 
Deviance = 4084.602 | Deviance change = 0.000554
Maximum parameter change: 0.000191 
...........................................................
Iteration 178     2019-01-06 08:14:56 
Deviance = 4084.601 | Deviance change = 0.0005442
Maximum parameter change: 0.000188 
...........................................................
Iteration 179     2019-01-06 08:14:56 
Deviance = 4084.601 | Deviance change = 0.0005346
Maximum parameter change: 0.000185 
...........................................................
Iteration 180     2019-01-06 08:14:56 
Deviance = 4084.6 | Deviance change = 0.0005252
Maximum parameter change: 0.000182 
...........................................................
Iteration 181     2019-01-06 08:14:56 
Deviance = 4084.6 | Deviance change = 0.0005159
Maximum parameter change: 0.000179 
...........................................................
Iteration 182     2019-01-06 08:14:56 
Deviance = 4084.599 | Deviance change = 0.0005068
Maximum parameter change: 0.000176 
...........................................................
Iteration 183     2019-01-06 08:14:56 
Deviance = 4084.599 | Deviance change = 0.0004979
Maximum parameter change: 0.000173 
...........................................................
Iteration 184     2019-01-06 08:14:56 
Deviance = 4084.598 | Deviance change = 0.0004892
Maximum parameter change: 0.00017 
...........................................................
Iteration 185     2019-01-06 08:14:56 
Deviance = 4084.598 | Deviance change = 0.0004806
Maximum parameter change: 0.000167 
...........................................................
Iteration 186     2019-01-06 08:14:56 
Deviance = 4084.597 | Deviance change = 0.0004721
Maximum parameter change: 0.000164 
...........................................................
Iteration 187     2019-01-06 08:14:56 
Deviance = 4084.597 | Deviance change = 0.0004638
Maximum parameter change: 0.000162 
...........................................................
Iteration 188     2019-01-06 08:14:56 
Deviance = 4084.596 | Deviance change = 0.0004557
Maximum parameter change: 0.000159 
...........................................................
Iteration 189     2019-01-06 08:14:56 
Deviance = 4084.596 | Deviance change = 0.0004477
Maximum parameter change: 0.000156 
...........................................................
Iteration 190     2019-01-06 08:14:56 
Deviance = 4084.595 | Deviance change = 0.0004399
Maximum parameter change: 0.000154 
...........................................................
Iteration 191     2019-01-06 08:14:56 
Deviance = 4084.595 | Deviance change = 0.0004322
Maximum parameter change: 0.000151 
...........................................................
Iteration 192     2019-01-06 08:14:56 
Deviance = 4084.595 | Deviance change = 0.0004246
Maximum parameter change: 0.000148 
...........................................................
Iteration 193     2019-01-06 08:14:56 
Deviance = 4084.594 | Deviance change = 0.0004172
Maximum parameter change: 0.000146 
...........................................................
Iteration 194     2019-01-06 08:14:56 
Deviance = 4084.594 | Deviance change = 0.0004099
Maximum parameter change: 0.000144 
...........................................................
Iteration 195     2019-01-06 08:14:56 
Deviance = 4084.593 | Deviance change = 0.0004027
Maximum parameter change: 0.000141 
...........................................................
Iteration 196     2019-01-06 08:14:56 
Deviance = 4084.593 | Deviance change = 0.0003957
Maximum parameter change: 0.000139 
...........................................................
Iteration 197     2019-01-06 08:14:56 
Deviance = 4084.593 | Deviance change = 0.0003888
Maximum parameter change: 0.000136 
...........................................................
Iteration 198     2019-01-06 08:14:56 
Deviance = 4084.592 | Deviance change = 0.000382
Maximum parameter change: 0.000134 
...........................................................
Iteration 199     2019-01-06 08:14:56 
Deviance = 4084.592 | Deviance change = 0.0003754
Maximum parameter change: 0.000132 
...........................................................
Iteration 200     2019-01-06 08:14:56 
Deviance = 4084.591 | Deviance change = 0.0003688
Maximum parameter change: 0.00013 
...........................................................
Iteration 201     2019-01-06 08:14:56 
Deviance = 4084.591 | Deviance change = 0.0003624
Maximum parameter change: 0.000127 
...........................................................
Iteration 202     2019-01-06 08:14:56 
Deviance = 4084.591 | Deviance change = 0.0003561
Maximum parameter change: 0.000125 
...........................................................
Iteration 203     2019-01-06 08:14:56 
Deviance = 4084.59 | Deviance change = 0.0003499
Maximum parameter change: 0.000123 
...........................................................
Iteration 204     2019-01-06 08:14:56 
Deviance = 4084.59 | Deviance change = 0.0003438
Maximum parameter change: 0.000121 
...........................................................
Iteration 205     2019-01-06 08:14:56 
Deviance = 4084.59 | Deviance change = 0.0003378
Maximum parameter change: 0.000119 
...........................................................
Iteration 206     2019-01-06 08:14:56 
Deviance = 4084.589 | Deviance change = 0.000332
Maximum parameter change: 0.000117 
...........................................................
Iteration 207     2019-01-06 08:14:56 
Deviance = 4084.589 | Deviance change = 0.0003262
Maximum parameter change: 0.000115 
...........................................................
Iteration 208     2019-01-06 08:14:56 
Deviance = 4084.589 | Deviance change = 0.0003205
Maximum parameter change: 0.000113 
...........................................................
Iteration 209     2019-01-06 08:14:56 
Deviance = 4084.588 | Deviance change = 0.000315
Maximum parameter change: 0.000111 
...........................................................
Iteration 210     2019-01-06 08:14:56 
Deviance = 4084.588 | Deviance change = 0.0003095
Maximum parameter change: 0.000109 
...........................................................
Iteration 211     2019-01-06 08:14:56 
Deviance = 4084.588 | Deviance change = 0.0003042
Maximum parameter change: 0.000108 
...........................................................
Iteration 212     2019-01-06 08:14:56 
Deviance = 4084.587 | Deviance change = 0.0002989
Maximum parameter change: 0.000106 
...........................................................
Iteration 213     2019-01-06 08:14:56 
Deviance = 4084.587 | Deviance change = 0.0002937
Maximum parameter change: 0.000104 
...........................................................
Iteration 214     2019-01-06 08:14:56 
Deviance = 4084.587 | Deviance change = 0.0002886
Maximum parameter change: 0.000102 
...........................................................
Iteration 215     2019-01-06 08:14:56 
Deviance = 4084.587 | Deviance change = 0.0002836
Maximum parameter change: 1e-04 
...........................................................
Iteration 216     2019-01-06 08:14:56 
Deviance = 4084.586 | Deviance change = 0.0002787
Maximum parameter change: 9.9e-05 
---------------------------------------------------------------------------------
Time difference of 0.7033627 secs
---------------------------------------------------------------------------------------------------------- 
CDM 7.0-12 (2018-11-06 16:23:12) 

Call:
CDM::gdina(data = dat, q.matrix = Q, rule = "DINA")

Date of Analysis: 2019-01-06 08:14:57 
Time difference of 0.7033627 secs
Computation Time: 0.7033627 

Generalized DINA Model 

Number of iterations = 216
Iteration with minimal deviance = 216 

Estimation method: WLS
Optimizer: CDM
Monotonicity constraints: FALSE
Number of items at boundary monotonicity constraint: NA

Parameter regularization: FALSE

Deviance = 4084.59  | Log likelihood = -2042.29 

Number of persons = 400 
Number of groups = 1 
Number of items = 9 
Number of estimated parameters = 25 
Number of estimated item parameters = 18 
Number of estimated skill class parameters = 7 ( 8 latent skill classes)

AIC = 4135  | penalty = 50    | AIC = -2*LL + 2*p  
BIC = 4234  | penalty = 149.79    | BIC = -2*LL + log(n)*p   
CAIC = 4259  | penalty = 174.79    | CAIC = -2*LL + [log(n)+1]*p  (consistent AIC)   

----------------------------------------------------------------------------
Used Q-matrix 

      V1 V2 V3
Item1  1  1  0
Item2  0  1  1
Item3  1  0  1
Item4  1  0  0
Item5  0  0  1
Item6  0  1  0
Item7  1  1  1
Item8  0  1  1
Item9  0  1  1

----------------------------------------------------------------------------

Item Parameter Estimates 

       link  item itemno partype rule    est     se partype.attr
1  identity Item1      1       0 DINA 0.0872 0.0166             
2  identity Item1      1     1-2 DINA 0.7035 0.0422        V1-V2
3  identity Item2      2       0 DINA 0.1092 0.0218             
4  identity Item2      2     1-2 DINA 0.6514 0.0375        V2-V3
5  identity Item3      3       0 DINA 0.1294 0.0205             
6  identity Item3      3     1-2 DINA 0.6853 0.0396        V1-V3
7  identity Item4      4       0 DINA 0.2287 0.0263             
8  identity Item4      4       1 DINA 0.5546 0.0431           V1
9  identity Item5      5       0 DINA 0.0058 0.0084             
10 identity Item5      5       1 DINA 0.9942 0.0566           V3
11 identity Item6      6       0 DINA 0.2516 0.0333             
12 identity Item6      6       1 DINA 0.2487 0.0469           V2
13 identity Item7      7       0 DINA 0.2430 0.0248             
14 identity Item7      7   1-2-3 DINA 0.2681 0.0558     V1-V2-V3
15 identity Item8      8       0 DINA 0.2782 0.0313             
16 identity Item8      8     1-2 DINA 0.5971 0.0392        V2-V3
17 identity Item9      9       0 DINA 0.3169 0.0325             
18 identity Item9      9     1-2 DINA 0.6561 0.0345        V2-V3

RMSD (RMSEA) Item Fit Statistics
Item1 Item2 Item3 Item4 Item5 Item6 Item7 Item8 Item9 
0.014 0.034 0.028 0.019 0.000 0.035 0.043 0.113 0.066 

Mean of RMSEA item fit: 0.039 
----------------------------------------------------------------------------
Model Implied Conditional Item Probabilities 

    item rule nessskill itemno skillcomb   prob
1  Item1 DINA     V1-V2      1       A00 0.0872
2  Item1 DINA     V1-V2      1       A10 0.0872
3  Item1 DINA     V1-V2      1       A01 0.0872
4  Item1 DINA     V1-V2      1       A11 0.7907
5  Item2 DINA     V2-V3      2       A00 0.1092
6  Item2 DINA     V2-V3      2       A10 0.1092
7  Item2 DINA     V2-V3      2       A01 0.1092
8  Item2 DINA     V2-V3      2       A11 0.7606
9  Item3 DINA     V1-V3      3       A00 0.1294
10 Item3 DINA     V1-V3      3       A10 0.1294
11 Item3 DINA     V1-V3      3       A01 0.1294
12 Item3 DINA     V1-V3      3       A11 0.8146
13 Item4 DINA        V1      4        A0 0.2287
14 Item4 DINA        V1      4        A1 0.7833
15 Item5 DINA        V3      5        A0 0.0058
16 Item5 DINA        V3      5        A1 1.0000
17 Item6 DINA        V2      6        A0 0.2516
18 Item6 DINA        V2      6        A1 0.5003
19 Item7 DINA  V1-V2-V3      7      A000 0.2430
20 Item7 DINA  V1-V2-V3      7      A100 0.2430
21 Item7 DINA  V1-V2-V3      7      A010 0.2430
22 Item7 DINA  V1-V2-V3      7      A001 0.2430
23 Item7 DINA  V1-V2-V3      7      A110 0.2430
24 Item7 DINA  V1-V2-V3      7      A101 0.2430
25 Item7 DINA  V1-V2-V3      7      A011 0.2430
26 Item7 DINA  V1-V2-V3      7      A111 0.5111
27 Item8 DINA     V2-V3      8       A00 0.2782
28 Item8 DINA     V2-V3      8       A10 0.2782
29 Item8 DINA     V2-V3      8       A01 0.2782
30 Item8 DINA     V2-V3      8       A11 0.8753
31 Item9 DINA     V2-V3      9       A00 0.3169
32 Item9 DINA     V2-V3      9       A10 0.3169
33 Item9 DINA     V2-V3      9       A01 0.3169
34 Item9 DINA     V2-V3      9       A11 0.9730
----------------------------------------------------------------------------

Skill Probabilities 

   skill.prob0 skill.prob1
V1      0.6370      0.3630
V2      0.4234      0.5766
V3      0.2012      0.7988
----------------------------------------------------------------------------

Polychoric Correlations 

Group 1
      V1    V2    V3
V1 1.000 0.436 0.395
V2 0.436 1.000 0.249
V3 0.395 0.249 1.000

----------------------------------------------------------------------------

Skill Pattern Probabilities 

   000    100    010    110    001    101    011    111 
0.1035 0.0097 0.0639 0.0240 0.2309 0.0792 0.2386 0.2501 
-----------------------------------------------------------------------------
CDM 7.0-12 (2018-11-06 16:23:12) 

Call:
CDM::IRT.RMSD(object = mod1)

-----------------------------------------------------------------------------
Root Mean Square Deviation (RMSD) 

    Parm     M    SD Min  Max
1 Group1 0.028 0.024   0 0.08

   item Group1
1 Item1  0.010
2 Item2  0.024
3 Item3  0.020
4 Item4  0.014
5 Item5  0.000
6 Item6  0.025
7 Item7  0.031
8 Item8  0.080
9 Item9  0.047
-----------------------------------------------------------------------------
Bias Corrected Root Mean Square Deviation (RMSD) 

    Parm     M    SD    Min   Max
1 Group1 0.016 0.029 -0.016 0.077

   item Group1
1 Item1 -0.013
2 Item2  0.016
3 Item3  0.009
4 Item4 -0.016
5 Item5 -0.002
6 Item6  0.008
7 Item7  0.021
8 Item8  0.077
9 Item9  0.043
-----------------------------------------------------------------------------
Mean Absolute Deviation (MAD) 

    Parm     M   SD Min   Max
1 Group1 0.021 0.02   0 0.068

   item Group1
1 Item1  0.007
2 Item2  0.022
3 Item3  0.016
4 Item4  0.006
5 Item5  0.000
6 Item6  0.019
7 Item7  0.018
8 Item8  0.068
9 Item9  0.030
-----------------------------------------------------------------------------
Mean Deviation (MD) 

    Parm M SD Min Max
1 Group1 0  0   0   0

   item Group1
1 Item1      0
2 Item2      0
3 Item3      0
4 Item4      0
5 Item5      0
6 Item6      0
7 Item7      0
8 Item8      0
9 Item9      0

CDM documentation built on March 19, 2019, 1:05 a.m.