R/lattice_pa.R

In cddesja/hemp: Handbook of Educational Measurement and Psychometrics Using R Companion Package

#' @title Parallel analysis plot via lattice
#'
#' @description
#' \code{lattice_pa} returns a parallel analysis based on a reduced correlation matrix using lattice
#'
#' @param data A data.frame with no missing data
#'
#' @import lattice
#' @import parallel
#' @importFrom psych fa
#' @examples
#' library(hemp)
#' cognition <- subset(interest, select = vocab:analyrea)
#' lattice_pa(cognition)
#' @export
#'
#'
lattice_pa <- function(data, ...){
  x <- hemp::fa.parallel.hemp(data)
  dat.tmp <- data.frame(fa.values = x[[1]],
                            fa.sim = x[[5]],
                            nfactors = seq_along(x[[1]]))
  xyplot(fa.values ~ nfactors, data = dat.tmp,
         panel = function(x,y, ...){
           panel.xyplot(x, y, type = c("p", "l"), col = "black")
           panel.lines(x, y = dat.tmp$fa.sim, lty = 3, col = "black")
         },
         xlab = "Number of Factors",
         ylab = "Eigenvalues of reduced correlation matrix",
        )
}