The function 'unitaryAlpha()' computes the unitary alpha ((Laurencelle and Cousineau 2023) ). This quantity is a novel way to compute correlation in a matrix where each column is a measure and each line, a subject. This measure is based on Cronbach's alpha (which could be labeled a 'global alpha').

unitaryAlpha( m )

Arguments

m

A data matrix for a group of observations.

Value

A measure of correlation between -1 and +1.

Details

This measure is derived from Cronbach' measure of reliability as shown by Laurencelle and Cousineau (2023) .

References

Laurencelle L, Cousineau D (2023). “Analysis of proportions using arcsine transform with any experimental design.” Frontiers in Psychology, 13, 1045436. doi:10.3389/fpsyg.2022.1045436 .

Examples


# Generate a random matrix (here binary entries)
set.seed(42)
N <- M <- 10
m <- matrix( runif(N*M), N, M)

# compute the unitary alpha from that random matrix
unitaryAlpha(m)
#> [1] 0.03055626