R/thres5.R

In InDisc: Obtaining and Estimating Unidimensional and Multidimensional IRT Dual Models

thres5 <- function(T){

  n <- size(T)[2]

  b <- matrix(0,4,1)

  term1 <- 0
  term2 <- 0
  term3 <- 0
  term4 <- 0

  for (i in 1:n){
    if (T[i] == 1){
      term1 <- term1 + 1
    }
    if (T[i] == 2){
      term2 <- term2 + 1
    }
    if (T[i] == 3){
      term3 <- term3 + 1
    }
    if (T[i] == 4){
      term4 <- term4 + 1
    }
  }

  if (term1 == 0){ term1 <- .5; n <- n + 1}
  if (term2 == 0){ term2 <- .5; n <- n + 1}
  if (term3 == 0){ term3 <- .5; n <- n + 1}
  if (term4 == 0){ term4 <- .5; n <- n + 1}
  if (term1+term2+term3+term4 == n){ term4 = term4 - 0.5; n <- n + 1}

  term1 <- term1 / n
  term2 <- term2 / n
  term3 <- term3 / n
  term4 <- term4 / n

  b[1,1] <- qnorm(term1)
  b[2,1] <- qnorm(term1+term2)
  b[3,1] <- qnorm(term1+term2+term3)
  b[4,1] <- qnorm(term1+term2+term3+term4)

  return(b)

}