R/calc_rel_emp.R

In susanne-frick/MFCblockInfo: compute block information for multidimensional forced-choice questionnaires

#' empirical reliability based on factor scores
#'
#' @param scores matrix of factor scores, rows = persons, columns = traits
#' @param ses matrix of SEs, same as scores
#'
#' @return vector of empirical reliabilities for each trait
#' @export
#'
calc.rel.emp <- function(scores, ses) {
  svar <- apply(scores, 2, var)
  svar/(svar + colMeans(ses^2, na.rm=T))
}