The IRT Procedure

Example 53.3 Multiple-Group Analysis

This example shows how to use the IRT procedure to do multiple-group analysis. The following DATA step creates the data set IrtGroup:

data IrtGroup;
   input item1-item8 GroupVar @@;
   datalines;
1 0 0 0 1 1 2 1 2 1 1 1 1 1 3 3 3 1 0 1 0 0 1 1 1 1 1 1 0 0 1 0 1 2 3
2 0 0 0 0 0 1 1 1 1 1 0 0 1 0 1 3 3 1 0 0 0 0 0 1 1 3 2 0 0 1 0 0 1 2
2 1 0 1 0 0 1 1 1 2 1 0 0 0 0 0 2 2 3 1 0 1 0 1 0 2 3 3 2 0 0 1 0 1 1
2 3 1 1 1 1 1 1 2 2 3 2 0 0 0 0 1 1 2 2 1 1 0 1 1 1 2 3 3 2 0 1 0 0 1
1 2 3 2 1 0 1 1 1 2 3 3 2 0 1 0 1 1 3 2 3 2 1 1 1 0 0 1 3 3 1 1 1 0 0
1 2 3 3 2 0 1 1 1 1 1 2 1 2 1 0 0 1 1 3 1 1 2 1 0 0 0 1 1 1 3 2 0 0 0
0 1 1 3 3 1 0 0 0 0 0 1 1 1 2 1 0 0 0 0 1 3 3 2 1 1 0 1 1 3 1 1 1 1 0
1 1 1 1 3 1 2 1 1 1 1 1 2 3 2 2 0 0 1 0 0 2 2 2 1 0 0 1 0 1 1 2 3 2 1

   ... more lines ...   

1 0 0 2 1 3 2 1 1 1 1 1 1 3 2 1 1 1 1 0 0 3 3 1 2 1 0 0 1 1 3 3 3 2 0
0 1 0 0 1 1 1 2 0 0 0 0 1 3 1 1 2 1 0 0 1 0 1 3 3 2 0 0 1 1 0 2 2 3 2
1 0 0 0 1 1 1 1 1 0 0 0 0 1 1 1 1 2 0 0 1 1 0 1 1 3 2 0 0 0 0 1 1 1 1
2 0 0 1 0 0 1 1 1 1 
; 

To set up a multiple-group IRT model, you need to specify the grouping variable by using the GROUP statement. Very often you also want to specify cross-group equality constraints for different parameters. You can accomplish this by using the EQUALITY statement.

The model that is specified in the following statements is an extension of the model in Example 53.1. The group variable, GroupVar, is specified in the GROUP statement. It has two values, 1 and 2, to indicate group membership. Equality constraints are specified in the EQUALITY statement.

proc irt data=IrtGroup;
   var item1-item8;
   group GroupVar;
   model item1-item4/resfunc=twop,
         item5-item8/resfunc=graded;
   equality item1-item4/parm=[intercept] between_gp=[1 2],
            _allgr_/parm=[slope] within_gp=[1];
run;

Two different sets of equality constraints have been specified. The first entry specifies equality constraints on the intercept parameters for item1 to item4 between group 1 and group 2:

item1-item4/parm=[intercept] between_gp=[1 2]

Notice that 1 and 2 are the actual values for the group variable GroupVar. The second entry specifies equality constraints on the slope parameters for all the graded response items within group 1:

_allgr_/parm=[slope] within_gp=[1]

Output 53.3.1 shows the "Modeling Information" table and the "Group Information" table for this example. For multiple-group analysis, the "Modeling Information" table contains two extra pieces of information: the group variable and the number of groups. The "Group Information" table contains information about the data for each group. There are 272 observations that have been read and used for group 1; this number for group 2 is 328.

Output 53.3.1: Modeling and Group Information

The IRT Procedure

Modeling Information
Data Set WORK.IRTGROUP
Group Variable GroupVar
Link Function Logit
Number of Items 8
Number of Factors 1
Number of Groups 2
Number of Observations Read 600
Number of Observations Used 600
Estimation Method Marginal Maximum Likelihood

Group Information
GroupVar Nobs
Read
Nobs
Used
1 272 272
2 328 328



Because there are two groups in this example, the IRT procedure produces two "Item Information" tables. For this example, these two tables contain the same information. That means that all the items have the same levels for the two groups. It is possible that the same item might have different numbers of levels, or maybe the same number of levels but different values. For example, an item has four levels, from 1 to 4, but one group might observe only levels 1 and 2, and the other group might observe only levels 3 and 4.

Output 53.3.2: Item Information

The IRT Procedure

Item Information
GroupVar = 1
Response
Model
Item Levels Values
TwoP item1 2 0 1
  item2 2 0 1
  item3 2 0 1
  item4 2 0 1
Graded item5 2 0 1
  item6 3 1 2 3
  item7 3 1 2 3
  item8 3 1 2 3

Item Information
GroupVar = 2
Response
Model
Item Levels Values
TwoP item1 2 0 1
  item2 2 0 1
  item3 2 0 1
  item4 2 0 1
Graded item5 2 0 1
  item6 3 1 2 3
  item7 3 1 2 3
  item8 3 1 2 3



Output 53.3.3 includes "Item Parameter Estimates" tables for both groups. You can see that the intercept parameters for item1 are the same for both groups. The same applies to item2 to item4. You can also see that the slope parameters have the same value for item5 to item8 in group 1. These results suggest that equality constraints that are specified in the EQUALITY statement have been fulfilled.

Output 53.3.3: Parameter Estimates

The IRT Procedure

Item Parameter Estimates
GroupVar = 1
Response
Model
Item Parameter Estimate Standard
Error
Pr > |t|
TwoP item1 Difficulty 0.24545 0.07800 0.0008
    Slope 1.89309 0.37484 <.0001
  item2 Difficulty 0.49512 0.10112 <.0001
    Slope 2.00288 0.37954 <.0001
  item3 Difficulty 0.41490 0.10687 <.0001
    Slope 1.40433 0.26977 <.0001
  item4 Difficulty 0.57754 0.16188 0.0002
    Slope 0.90590 0.20213 <.0001
Graded item5 Threshold -0.71312 0.16145 <.0001
    Slope 0.96099 0.10872 <.0001
  item6 Threshold 1 -0.15789 0.14680 0.1411
    Threshold 2 1.42434 0.21483 <.0001
    Slope 0.96099 0.10872 <.0001
  item7 Threshold 1 -0.86494 0.16861 <.0001
    Threshold 2 0.81432 0.17160 <.0001
    Slope 0.96099 0.10872 <.0001
  item8 Threshold 1 -1.24131 0.19158 <.0001
    Threshold 2 0.37597 0.15203 0.0067
    Slope 0.96099 0.10872 <.0001

Item Parameter Estimates
GroupVar = 2
Response
Model
Item Parameter Estimate Standard
Error
Pr > |t|
TwoP item1 Difficulty 0.31862 0.09331 0.0003
    Slope 1.45833 0.28819 <.0001
  item2 Difficulty 0.68599 0.12813 <.0001
    Slope 1.44561 0.28659 <.0001
  item3 Difficulty 0.51327 0.12276 <.0001
    Slope 1.13519 0.22854 <.0001
  item4 Difficulty 0.61629 0.16217 <.0001
    Slope 0.84894 0.18964 <.0001
Graded item5 Threshold -0.90646 0.27324 0.0005
    Slope 0.68749 0.17606 <.0001
  item6 Threshold 1 -0.77880 0.29844 0.0045
    Threshold 2 2.20837 0.57332 <.0001
    Slope 0.54126 0.14327 <.0001
  item7 Threshold 1 -0.84400 0.18217 <.0001
    Threshold 2 0.80547 0.16743 <.0001
    Slope 1.03946 0.18984 <.0001
  item8 Threshold 1 -1.18570 0.20677 <.0001
    Threshold 2 0.19876 0.12012 0.0490
    Slope 1.15609 0.21044 <.0001