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]