LikertMakeR (version 1.0.0)

makeCorrLoadings: Correlation matrix from item factor loadings

Description

makeCorrLoadings() generates a correlation matrix of inter-item correlations based on item factor loadings as might be seen in Exploratory Factor Analysis (EFA) or a Structural Equation Model (SEM).

Such a correlation matrix can be applied to the makeItems() function to generate synthetic data with those predefined factor structures.

Usage

makeCorrLoadings(
  loadings,
  factorCor = NULL,
  uniquenesses = NULL,
  nearPD = FALSE
)

Value

Correlation matrix of inter-item correlations

Arguments

loadings

(numeric matrix) k (items) by f (factors) matrix of standardised factor loadings. Item names and Factor names are taken from the row_names (items) and the column_names (factors), if present.

@note "Censored" loadings (for example, where loadings less than '0.30' are removed for clarity) tend to severely reduce the accuracy of the makeCorrLoadings() function. For a detailed demonstration, see the file makeCorrLoadings_Validate.pdf in the package website on GitHub.

factorCor

(numeric matrix) f x f factor correlation matrix

uniquenesses

(numeric vector) length k vector of uniquenesses. If NULL, the default, compute from the calculated communalities.

nearPD

(logical) If TRUE, project factorCor and the final correlation matrix onto nearest Positive Definite matrix, if needed.

Examples

# Example loadings
factorLoadings <- matrix(
  c(
    0.05, 0.20, 0.70,
    0.10, 0.05, 0.80,
    0.05, 0.15, 0.85,
    0.20, 0.85, 0.15,
    0.05, 0.85, 0.10,
    0.10, 0.90, 0.05,
    0.90, 0.15, 0.05,
    0.80, 0.10, 0.10
  ),
  nrow = 8, ncol = 3, byrow = TRUE
)

# row and column names
rownames(factorLoadings) <- c("Q1", "Q2", "Q3", "Q4", "Q5", "Q6", "Q7", "Q8")
colnames(factorLoadings) <- c("Factor1", "Factor2", "Factor3")

# Factor correlation matrix
factorCor <- matrix(
  c(
    1.0,  0.7, 0.6,
    0.7,  1.0, 0.4,
    0.6,  0.4, 1.0
  ),
  nrow = 3, byrow = TRUE
)

# Apply the function
itemCorrelations <- makeCorrLoadings(factorLoadings, factorCor)

round(itemCorrelations, 3)