sirt (version 0.36-30)

plausible.value.imputation.raschtype: Plausible Value Imputation in Generalized Logistic Item Response Model

Description

This function performs unidimensional plausible value imputation (Adams & Wu, 2007; Mislevy, 1991).

Usage

plausible.value.imputation.raschtype(data=NULL, f.yi.qk=NULL, X, 
   Z=NULL, beta0=rep(0, ncol(X)), sig0=1, b=rep(1, ncol(X)), 
   a=rep(1, length(b)), c=rep(0, length(b)), d=1+0*b, 
   alpha1=0, alpha2=0, theta.list=seq(-5, 5, len=50), 
   cluster=NULL, iter, burnin, nplausible=1, printprogress=TRUE)

Arguments

data

An $N \times I$ data frame of dichotomous responses

f.yi.qk

An optional matrix which contains the individual likelihood. This matrix is produced by rasch.mml2 or rasch.copula2. The use of this argument

X

A matrix of individual covariates for the latent regression of $\theta$ on $X$

Z

A matrix of individual covariates for the regression of individual residual variances on $Z$

beta0

Initial vector of regression coefficients

sig0

Initial vector of coefficients for the variance heterogeneity model

b

Vector of item difficulties. It must not be provided if the individual likelihood f.yi.qk is specified.

a

Optional vector of item slopes

c

Optional vector of lower item asymptotes

d

Optional vector of upper item asymptotes

alpha1

Parameter $\alpha_1$ in generalized item response model

alpha2

Parameter $\alpha_2$ in generalized item response model

theta.list

Vector of theta values at which the ability distribution should be evaluated

cluster

Cluster identifier (e.g. schools or classes) for including theta means in the plausible imputation.

iter

Number of iterations

burnin

Number of burn-in iterations for plausible value imputation

nplausible

Number of plausible values

printprogress

Print progress?

Value

• A list with following entries:
• coefs.X
  Sampled regression coefficients for covariates $X$
• coefs.Z
  Sampled coefficients for modeling variance heterogeneity for covariates $Z$
• pvdraws
  Matrix with drawn plausible values
• posterior
  Posterior distribution from last iteration
• EAP
  Individual EAP estimate
• SE.EAP
  Standard error of EAP estimate
• pv.indexes
  Index of iterations when plausible values were drawn

Details

Plausible values are drawn from the latent regression model with heterogeneous variances: $$\theta_p = X_p \beta + \epsilon_p \quad , \quad \epsilon_p \sim N( 0 , \sigma_p^2 ) \quad , \quad \log( \sigma_p ) = Z_p \gamma + \nu_p$$

References

Adams, R., & Wu. M. (2007). The mixed-coefficients multinomial logit model: A generalized form of the Rasch model. In M. von Davier & C. H. Carstensen: Multivariate and Mixture Distribution Rasch Models: Extensions and Applications (pp. 57-76). New York: Springer. Mislevy, R. J. (1991). Randomization-based inference about latent variables from complex samples. Psychometrika, 56, 177-196.

See Also

For estimating the latent regression model see latent.regression.em.raschtype.

Examples

############################################################
### SIMULATED EXAMPLE 1
set.seed(899)
############################################################

I <- 21     # number of items
b <- seq(-2,2, len=I)   # item difficulties
n <- 2000       # number of students

# simulate theta and covariates
theta <- rnorm( n )
x <- .7 * theta + rnorm( n , .5 )
y <- .2 * x+ .3*theta + rnorm( n , .4 )
dfr <- data.frame( theta , 1 , x , y )

# simulate Rasch model
dat1 <- sim.raschtype( theta = theta , b = b )

# Plausible value draws
pv1 <- plausible.value.imputation.raschtype(data=dat1 , X=dfr[,-1] , b = b ,
    nplausible=3 , iter=10 , burnin=5)
# estimate linear regression based on first plausible value
mod1 <- lm( pv1$pvdraws[,1] ~ x+y )
summary(mod1)
##               Estimate Std. Error t value Pr(>|t|)    
##   (Intercept) -0.27755    0.02121  -13.09   <2e-16 ***
##   x            0.40483    0.01640   24.69   <2e-16 ***
##   y            0.20307    0.01822   11.15   <2e-16 ***