est_rank {lazy.irt} R Documentation

## Estimation of the Latent Rank of LRT Model

### Description

Estimation of the Latent Rank of LRT Model

### Usage

est_rank(
Uc = NULL,
U = NULL,
V = NULL,
alpha = rep(1, ncol(V)),
vmin = 1e-04,
vmax = 1 - vmin,
rho = NULL,
print = 0,
plot = 0,
title = NULL
)

### Arguments

 Uc Item Response Data in compressed format U Item Response Data V LRT Item Parameter Matrix: nitem x nclass alpha = vector of Dirichle parameters of length nclass or a number A to set alpha=rep(A,nclass) or a negative number -A to set alpha=A*normal pdf vmin minimum value of V vmax maximum value of V rho initial value of probability vector of each latent class. print = 1 to print the result plot = 1 to plot the estimated theta distributions title Title strings

### Details

Note that uLRT returns the V as the nclass x nitem data frame, whereas in this function V is defined as nitem x nclass. (sorry!)

The core part of uLRT, namely, the E-step, is used in this function.

The rank of a person i is defined as the location of the highest posterior probability of H[i,].

### Value

A list of:
rank The estimated rank: nrow(Uc) x 1
H The posterior probability matrix: nrow(Uc) x nclass
rho The prior probability of the class:
alpha The prior of rho:
nclass The number of classes
method EAP

### Examples

#### In the following examples, maxiter is set to 20 which is
#### not large enough to obtain convergence.
####
#
#
#
set.seed(1701)

param=paramB1[c(1:3,7:9,13:15),]
thmin=-2; thmax=2; npoint=5
N=1000

# discrete theta
# theta0=seq(thmin,thmax,length=npoint)
theta0=c(-2, -1, 0, 2, 3)
theta=unlist(lapply( theta0, rep, round(N/npoint) ))
res2 <- gendataIRT( 1, paramB1, theta=theta, compress=1 )
Uc=as.data.frame(res2\$U)
ncat=res2\$ncat
type=res2\$type

nclass=5
res1 <- uLRT( Uc, nclass=nclass, estrho=1, monotone=1, alpha=20
, maxiter=20, plot=1, print=1 )
V1=res1\$V[,seq(2,2*res1\$nitems,2)]
res <- est_rank( Uc=Uc, V=t(V1), print=1, alpha=-4, plot=1, title="Test" )

[Package lazy.irt version 0.1.3 Index]