log_reg: A logistic regression DIF test for MSTs

In RDebelak/mstDIF: A Collection of DIF Tests for Multistage Tests

A logistic regression DIF test for MSTs

Description

This function allows the detection of itemwise DIF for Multistage Tests. It is based on the comparison of three logistic regression models for each item. The first logistic regression model (Model 1) predicts the positiveness of each response solely on the estimated ability parameters. The second logistic regression model (Model 2) predicts the positiveness based on the ability parameters and the membership to the focal and reference group as additive predictor variables. The third model (Model 3) uses the same predictors as Model 2 to predict the positiveness of the responses, but also includes an interaction effect. Three model comparisons are carried out (Models 1/2, Models 1/3, Models 2/3) based on two criteria: The comparison of the Nagelkerke R squared values, and the p-values of a likelihood ratio test.

Usage

log_reg(resp, DIF_covariate, theta = NULL)

Arguments

resp	A data frame containing the response matrix. Rows correspond to respondents, columns to items.
DIF_covariate	A factor indicating the membership to the reference and focal groups.
theta	A vector of ability estimates for each respondent.

Details

Author: Sebastian Appelbaum, with minor changes by Rudolf Debelak and Dries Debeer

Value

A list with four elements. The first element is the response matrix, the second element is the name of the DIF covariate, and the third element is the name of the test. The fourth element is a data frame where each row corresponds to an item. The columns of this data frame correspond to the following entries:

N

The number of responses observed for this item.

overall_chi_sq

The chi squared statistic of the likelihood ratio test comparing Model 1 and Model 3.

overall_p_value

The p-values of the likelihood ratio test comparing Model 1 and Model 3 as an indicator for the overall DIF effect.

Delta_NagelkerkeR2

The difference of the Nagelkerke R squared values for Model 1 and Model 3.

UDIF_chi_sq

The chi squared statistic of the likelihood ratio test comparing Model 1 and Model 2.

UDIF_p_value

The p-values of the likelihood ratio test comparing Model 1 and Model 2.

UDIF_Delta_NagelkerkeR2

The difference of the Nagelkerke R squared values for Model 1 and Model 2.

CDIF_chi_sq

The chi squared statistic of the likelihood ratio test comparing Model 2 and Model 3.

CDIF_p_value

The p-values of the likelihood ratio test comparing Model 2 and Model 3.

CDIF_Delta_NagelkerkeR2

The difference of the Nagelkerke R squared values for Model 2 and Model 3.

Examples

data("toydata")
resp <- toydata$resp
group_categ <- toydata$group_categ
theta_est <- toydata$theta_est
log_reg(resp, DIF_covariate = factor(group_categ), theta = theta_est)

RDebelak/mstDIF documentation built on Dec. 8, 2022, 7:33 a.m.