# 1 Functionalities

The equateMultiple package computes:

• Equating coefficients between multiple forms.
• Synthetic item parameters (sort of mean of the item parameter estimates from different forms).
• Standard errors of the equating coefficients and the synthetic item parameters.

# 2 Data preparation

Data preparation follows the same steps of the equateIRT package.

Load the package equateMultiple and the data

library("equateMultiple")
## Loading required package: equateIRT
data("data2pl", package = "equateIRT")

Estimate a two parameter logistic model for 5 data sets with the R package ltm

library("ltm")
m1 <- ltm(data2pl[[1]] ~ z1)
m2 <- ltm(data2pl[[2]] ~ z1)
m3 <- ltm(data2pl[[3]] ~ z1)
m4 <- ltm(data2pl[[4]] ~ z1)
m5 <- ltm(data2pl[[5]] ~ z1)

Extract the item parameter estimates and the covariance matrices

estm1 <- import.ltm(m1, display = FALSE)
estm2 <- import.ltm(m2, display = FALSE)
estm3 <- import.ltm(m3, display = FALSE)
estm4 <- import.ltm(m4, display = FALSE)
estm5 <- import.ltm(m5, display = FALSE)
estm1$coef[1:3, ] ## (Intercept) z1 ## I1 -0.06213808 1.076155 ## I2 -0.03090993 1.122379 ## I3 -0.07939847 1.091369 estm1$var[1:3, 1:3]
##              [,1]         [,2]         [,3]
## [1,] 0.0012285184 0.0002460322 0.0002391000
## [2,] 0.0002460322 0.0012628923 0.0002495126
## [3,] 0.0002391000 0.0002495126 0.0012407430

Create a list of coefficients and covariance matrices

estc <- list(estm1$coef, estm2$coef, estm3$coef, estm4$coef, estm5$coef) estv <- list(estm1$var, estm2$var, estm3$var, estm4$var, estm5$var)
test <- paste("test", 1:5, sep = "")

Create an object of class modIRT

mods <- modIRT(coef = estc, var = estv, names = test, display = FALSE)
coef(mods$test1)[1:5] ## Dffclt.I1 Dffclt.I2 Dffclt.I3 Dffclt.I4 Dffclt.I5 ## 0.05774085 0.02753964 0.07275128 0.41568210 -0.00716265 The linkage plan lplan<-linkp(coef = estc) lplan ## [,1] [,2] [,3] [,4] [,5] ## [1,] 20 10 0 0 10 ## [2,] 10 20 10 0 0 ## [3,] 0 10 20 10 0 ## [4,] 0 0 10 20 10 ## [5,] 10 0 0 10 20 # 3 Multiple equating coefficients Estimation of the equating coefficients using the multiple mean-mean method. Form 1 is the base form. eqMM <- multiec(mods = mods, base = 1, method = "mean-mean") ## Computation of equating coefficients . . . . ## Computation of standard errors . . . . summary(eqMM) ## Equating coefficients: ## EQ Form Estimate StdErr ## A test1 1.00000 0.000000 ## A test2 0.84051 0.018648 ## A test3 0.84347 0.021334 ## A test4 0.83937 0.020694 ## A test5 1.02343 0.021524 ## B test1 0.00000 0.000000 ## B test2 0.10692 0.022427 ## B test3 0.20236 0.024039 ## B test4 0.36774 0.024111 ## B test5 0.50207 0.024026 Estimation of the equating coefficients using the multiple mean-geometric mean method. eqMGM <- multiec(mods = mods, base = 1, method = "mean-gmean") ## Computation of equating coefficients . . . . ## Computation of standard errors . . . . summary(eqMGM) ## Equating coefficients: ## EQ Form Estimate StdErr ## A test1 1.00000 0.000000 ## A test2 0.83860 0.018695 ## A test3 0.84045 0.021383 ## A test4 0.83633 0.020746 ## A test5 1.02135 0.021591 ## B test1 0.00000 0.000000 ## B test2 0.10695 0.022410 ## B test3 0.20277 0.023938 ## B test4 0.36764 0.024044 ## B test5 0.50188 0.024001 Estimation of the equating coefficients using the multiple item response function method. eqIRF<-multiec(mods = mods, base = 1, method = "irf") ## Computation of equating coefficients . . . . ## Computation of standard errors . . . . summary(eqIRF) ## Equating coefficients: ## EQ Form Estimate StdErr ## A test1 1.00000 0.000000 ## A test2 0.83635 0.018353 ## A test3 0.83610 0.020919 ## A test4 0.82919 0.020173 ## A test5 1.01249 0.021184 ## B test1 0.00000 0.000000 ## B test2 0.10802 0.021770 ## B test3 0.20935 0.023031 ## B test4 0.37199 0.023090 ## B test5 0.49709 0.023555 Estimation of the equating coefficients using the multiple item response function method. The initial values are the estimates obtained with the multiple mean-geometric mean method. eqMGM <- multiec(mods = mods, base = 1, method = "mean-gmean", se = FALSE) ## Computation of equating coefficients . . . . eqIRF<-multiec(mods = mods, base = 1, method = "irf", start = eqMGM) ## Computation of equating coefficients . . . . ## Computation of standard errors . . . . summary(eqIRF) ## Equating coefficients: ## EQ Form Estimate StdErr ## A test1 1.00000 0.000000 ## A test2 0.83635 0.018353 ## A test3 0.83610 0.020919 ## A test4 0.82919 0.020173 ## A test5 1.01249 0.021184 ## B test1 0.00000 0.000000 ## B test2 0.10802 0.021770 ## B test3 0.20935 0.023031 ## B test4 0.37199 0.023090 ## B test5 0.49709 0.023555 Estimation of the equating coefficients using the multiple test response function method. eqTRF<-multiec(mods = mods, base = 1, method = "trf") ## Computation of equating coefficients . . . . ## Computation of standard errors . . . . summary(eqTRF) ## Equating coefficients: ## EQ Form Estimate StdErr ## A test1 1.00000 0.000000 ## A test2 0.83681 0.018421 ## A test3 0.83744 0.021048 ## A test4 0.83154 0.020298 ## A test5 1.01646 0.021209 ## B test1 0.00000 0.000000 ## B test2 0.10639 0.021817 ## B test3 0.20580 0.023118 ## B test4 0.36873 0.023154 ## B test5 0.49499 0.023600 # 4 Synthetic item parameters The synthetic item parameters and their standard errors can be extracted as follows (using the multiple item response function method). eqIRF$as
##  Dscrmn.I1 Dscrmn.I10 Dscrmn.I11 Dscrmn.I12 Dscrmn.I13 Dscrmn.I14 Dscrmn.I15
##  1.0263300  1.3186514  1.0628256  1.1077121  1.3780255  1.1849060  1.0797961
## Dscrmn.I16 Dscrmn.I17 Dscrmn.I18 Dscrmn.I19  Dscrmn.I2 Dscrmn.I20 Dscrmn.I21
##  1.3717438  1.2509322  1.1473431  1.3149445  1.1374079  1.0070348  1.1943277
## Dscrmn.I22 Dscrmn.I23 Dscrmn.I24 Dscrmn.I25 Dscrmn.I26 Dscrmn.I27 Dscrmn.I28
##  1.2399568  1.2972537  1.0413250  1.2803009  1.3294265  1.3129234  1.1536342
## Dscrmn.I29  Dscrmn.I3 Dscrmn.I30 Dscrmn.I31 Dscrmn.I32 Dscrmn.I33 Dscrmn.I34
##  1.1056767  1.0782685  1.3051024  1.4679605  1.3688577  1.3333632  1.4005527
## Dscrmn.I35 Dscrmn.I36 Dscrmn.I37 Dscrmn.I38 Dscrmn.I39  Dscrmn.I4 Dscrmn.I40
##  1.2974183  1.0511445  1.2689048  1.4780647  1.3871910  1.3145141  1.4580333
## Dscrmn.I41 Dscrmn.I42 Dscrmn.I43 Dscrmn.I44 Dscrmn.I45 Dscrmn.I46 Dscrmn.I47
##  1.1895968  1.1991245  1.1820406  1.1438774  1.3954929  1.3527279  1.2618251
## Dscrmn.I48 Dscrmn.I49  Dscrmn.I5 Dscrmn.I50  Dscrmn.I6  Dscrmn.I7  Dscrmn.I8
##  1.1551069  0.9800754  1.0270970  1.2328251  0.9410529  1.0052578  1.1833536
##  Dscrmn.I9
##  1.0040658
eqIRF$bs ## Dffclt.I1 Dffclt.I10 Dffclt.I11 Dffclt.I12 Dffclt.I13 Dffclt.I14 ## 0.04680159 0.67020856 0.93109634 0.79043223 0.43683497 0.75463195 ## Dffclt.I15 Dffclt.I16 Dffclt.I17 Dffclt.I18 Dffclt.I19 Dffclt.I2 ## 0.82951627 0.11667937 0.59893566 0.57922077 0.89350937 0.01160208 ## Dffclt.I20 Dffclt.I21 Dffclt.I22 Dffclt.I23 Dffclt.I24 Dffclt.I25 ## 0.07364461 0.23644361 -0.09992649 -0.41581481 0.62576552 0.39477998 ## Dffclt.I26 Dffclt.I27 Dffclt.I28 Dffclt.I29 Dffclt.I3 Dffclt.I30 ## 0.28054816 0.74507324 0.38431305 -0.41217210 0.04639957 0.81854107 ## Dffclt.I31 Dffclt.I32 Dffclt.I33 Dffclt.I34 Dffclt.I35 Dffclt.I36 ## 0.52952744 0.80600281 0.55941089 -0.36603676 0.45840659 0.85483709 ## Dffclt.I37 Dffclt.I38 Dffclt.I39 Dffclt.I4 Dffclt.I40 Dffclt.I41 ## 0.48752328 -0.12828412 0.51856681 0.38850704 0.85860534 -0.58246660 ## Dffclt.I42 Dffclt.I43 Dffclt.I44 Dffclt.I45 Dffclt.I46 Dffclt.I47 ## -0.28920713 1.26110283 1.12164448 0.56075671 -0.31084350 0.47861241 ## Dffclt.I48 Dffclt.I49 Dffclt.I5 Dffclt.I50 Dffclt.I6 Dffclt.I7 ## 0.20196743 0.27988322 0.01748403 0.34034716 -0.77150306 0.14004614 ## Dffclt.I8 Dffclt.I9 ## 0.33614495 -0.17864329 eqIRF$se.as
##  Dscrmn.I1 Dscrmn.I10 Dscrmn.I11 Dscrmn.I12 Dscrmn.I13 Dscrmn.I14 Dscrmn.I15
## 0.03495588 0.04187786 0.04221454 0.04287726 0.04899610 0.04464139 0.04231166
## Dscrmn.I16 Dscrmn.I17 Dscrmn.I18 Dscrmn.I19  Dscrmn.I2 Dscrmn.I20 Dscrmn.I21
## 0.04855524 0.04598656 0.04329678 0.04869024 0.03717496 0.03979567 0.04483082
## Dscrmn.I22 Dscrmn.I23 Dscrmn.I24 Dscrmn.I25 Dscrmn.I26 Dscrmn.I27 Dscrmn.I28
## 0.04645440 0.04913040 0.04132664 0.04703964 0.04825891 0.04844480 0.04380390
## Dscrmn.I29  Dscrmn.I3 Dscrmn.I30 Dscrmn.I31 Dscrmn.I32 Dscrmn.I33 Dscrmn.I34
## 0.04393780 0.03598386 0.04846504 0.04237596 0.04101176 0.03931653 0.04122265
## Dscrmn.I35 Dscrmn.I36 Dscrmn.I37 Dscrmn.I38 Dscrmn.I39  Dscrmn.I4 Dscrmn.I40
## 0.03831698 0.03406463 0.03771167 0.04231999 0.04048852 0.04112994 0.04347242
## Dscrmn.I41 Dscrmn.I42 Dscrmn.I43 Dscrmn.I44 Dscrmn.I45 Dscrmn.I46 Dscrmn.I47
## 0.04359194 0.04247471 0.04266961 0.04113018 0.04616402 0.04670080 0.04278476
## Dscrmn.I48 Dscrmn.I49  Dscrmn.I5 Dscrmn.I50  Dscrmn.I6  Dscrmn.I7  Dscrmn.I8
## 0.04032199 0.03641452 0.03503836 0.04209305 0.03478125 0.03458554 0.03819397
##  Dscrmn.I9
## 0.03477342
eqIRF\$se.bs
##  Dffclt.I1 Dffclt.I10 Dffclt.I11 Dffclt.I12 Dffclt.I13 Dffclt.I14 Dffclt.I15
## 0.02628883 0.02738288 0.03792320 0.03400229 0.02585817 0.03222779 0.03526458
## Dffclt.I16 Dffclt.I17 Dffclt.I18 Dffclt.I19  Dffclt.I2 Dffclt.I20 Dffclt.I21
## 0.02458570 0.02885806 0.02965635 0.03344091 0.02495266 0.02826440 0.02658532
## Dffclt.I22 Dffclt.I23 Dffclt.I24 Dffclt.I25 Dffclt.I26 Dffclt.I27 Dffclt.I28
## 0.02838793 0.03314244 0.03086282 0.02615282 0.02546005 0.02955263 0.02725851
## Dffclt.I29  Dffclt.I3 Dffclt.I30 Dffclt.I31 Dffclt.I32 Dffclt.I33 Dffclt.I34
## 0.03552543 0.02561236 0.03089467 0.02378732 0.02750672 0.02504117 0.02565134
## Dffclt.I35 Dffclt.I36 Dffclt.I37 Dffclt.I38 Dffclt.I39  Dffclt.I4 Dffclt.I40
## 0.02466081 0.03195958 0.02499000 0.02321337 0.02443617 0.02454001 0.02702879
## Dffclt.I41 Dffclt.I42 Dffclt.I43 Dffclt.I44 Dffclt.I45 Dffclt.I46 Dffclt.I47
## 0.03788098 0.03163075 0.03764630 0.03511477 0.02537526 0.03036798 0.02589473
## Dffclt.I48 Dffclt.I49  Dffclt.I5 Dffclt.I50  Dffclt.I6  Dffclt.I7  Dffclt.I8
## 0.02680342 0.02885080 0.02631625 0.02583940 0.03742903 0.02669245 0.02535449
##  Dffclt.I9
## 0.02745309

# 5 Equated scores

Equated scores with the true score equating method

scTSE<-score(eqIRF)
## The following scores are not attainable: 0
## The following scores are not attainable: 0
## The following scores are not attainable: 0
## The following scores are not attainable: 0
round(scTSE,3)
##     theta test1 test2.as.test1 StdErr_test2.as.test1 test3.as.test1
## 1  -2.345     1          1.073                 0.027          0.784
## 2  -1.662     2          2.072                 0.034          1.652
## 3  -1.243     3          3.041                 0.039          2.551
## 4  -0.930     4          3.992                 0.041          3.470
## 5  -0.673     5          4.933                 0.042          4.403
## 6  -0.450     6          5.870                 0.042          5.347
## 7  -0.248     7          6.806                 0.041          6.301
## 8  -0.061     8          7.742                 0.040          7.264
## 9   0.118     9          8.682                 0.040          8.236
## 10  0.293    10          9.626                 0.041          9.217
## 11  0.466    11         10.576                 0.043         10.208
## 12  0.642    12         11.534                 0.047         11.208
## 13  0.824    13         12.501                 0.052         12.221
## 14  1.018    14         13.481                 0.058         13.248
## 15  1.230    15         14.476                 0.062         14.293
## 16  1.472    16         15.491                 0.065         15.358
## 17  1.764    17         16.533                 0.065         16.450
## 18  2.152    18         17.611                 0.060         17.577
## 19  2.782    19         18.746                 0.044         18.752
## 20 35.219    20         20.000                 0.000         20.000
##    StdErr_test3.as.test1 test4.as.test1 StdErr_test4.as.test1 test5.as.test1
## 1                  0.041          0.933                 0.050          0.751
## 2                  0.060          1.968                 0.071          1.634
## 3                  0.070          3.017                 0.080          2.563
## 4                  0.074          4.062                 0.082          3.517
## 5                  0.074          5.096                 0.080          4.487
## 6                  0.071          6.117                 0.076          5.467
## 7                  0.067          7.124                 0.070          6.456
## 8                  0.063          8.119                 0.065          7.453
## 9                  0.060          9.102                 0.061          8.456
## 10                 0.060         10.075                 0.059          9.466
## 11                 0.063         11.041                 0.059         10.483
## 12                 0.068         12.002                 0.061         11.508
## 13                 0.074         12.960                 0.065         12.540
## 14                 0.081         13.920                 0.069         13.580
## 15                 0.087         14.884                 0.073         14.629
## 16                 0.091         15.858                 0.074         15.688
## 17                 0.089         16.847                 0.073         16.756
## 18                 0.081         17.861                 0.066         17.836
## 19                 0.059         18.909                 0.049         18.926
## 20                 0.000         20.000                 0.000         20.000
##    StdErr_test5.as.test1
## 1                  0.034
## 2                  0.047
## 3                  0.052
## 4                  0.052
## 5                  0.050
## 6                  0.047
## 7                  0.045
## 8                  0.043
## 9                  0.041
## 10                 0.040
## 11                 0.041
## 12                 0.042
## 13                 0.044
## 14                 0.047
## 15                 0.050
## 16                 0.051
## 17                 0.051
## 18                 0.047
## 19                 0.036
## 20                 0.000

Equated scores with the observed score equating method, avoiding computation of standard errors

scOSE<-score(eqIRF, method = "OSE", se = FALSE)
round(scOSE,3)
##    test1 test2.as.test1 test3.as.test1 test4.as.test1 test5.as.test1
## 1      0          0.031         -0.164         -0.024         -0.158
## 2      1          1.032          0.714          0.984          0.725
## 3      2          2.014          1.612          2.014          1.639
## 4      3          2.983          2.525          3.049          2.577
## 5      4          3.941          3.445          4.079          3.532
## 6      5          4.893          4.380          5.102          4.499
## 7      6          5.841          5.327          6.115          5.475
## 8      7          6.788          6.285          7.118          6.460
## 9      8          7.735          7.252          8.110          7.453
## 10     9          8.683          8.227          9.094          8.453
## 11    10          9.635          9.211         10.070          9.458
## 12    11         10.590         10.202         11.041         10.470
## 13    12         11.551         11.203         12.007         11.489
## 14    13         12.519         12.214         12.970         12.515
## 15    14         13.496         13.235         13.935         13.547
## 16    15         14.487         14.269         14.902         14.588
## 17    16         15.489         15.317         15.875         15.636
## 18    17         16.507         16.378         16.857         16.692
## 19    18         17.551         17.451         17.853         17.756
## 20    19         18.626         18.547         18.868         18.828
## 21    20         19.742         19.703         19.906         19.906