Rdimtools (version 1.0.0)

do.fa: Exploratory Factor Analysis

Description

do.fa is an optimization-based implementation of a popular technique for Exploratory Data Analysis. It is closely related to principal component analysis.

Usage

do.fa(
  X,
  ndim = 2,
  preprocess = c("center", "scale", "cscale", "decorrelate", "whiten"),
  maxiter = 1000,
  tolerance = 1e-06
)

Arguments

X

an \((n\times p)\) matrix or data frame whose rows are observations and columns represent independent variables.

ndim

an integer-valued number of loading variables, or target dimension.

preprocess

an additional option for preprocessing the data. Default is "center". See also aux.preprocess for more details.

maxiter

maximum number of iterations for updating.

tolerance

stopping criterion in a Frobenius norm.

Value

a named list containing

Y

an \((n\times ndim)\) matrix whose rows are embedded observations.

trfinfo

a list containing information for out-of-sample prediction.

projection

a \((p\times ndim)\) whose columns are basis for projection.

loadings

a \((p\times ndim)\) matrix whose rows are extracted loading factors.

noise

a length-\(p\) vector of estimated noise.

References

spearman_general_1904Rdimtools

Examples

# NOT RUN {
## use iris data
data(iris)
X = as.matrix(iris[,1:4])
lab = as.factor(iris[,5])

## 1. use centered data
output1 <- do.fa(X,ndim=2)

## 2. use decorrelated data
output2 <- do.fa(X,ndim=2,preprocess="decorrelate")

## 3. use whitened data
output3 <- do.fa(X,ndim=2,preprocess="whiten")

## Visualize three different projections
opar <- par(no.readonly=TRUE)
par(mfrow=c(3,1))
plot(output1$Y, col=lab, main="FA::centered")
plot(output2$Y, col=lab, main="FA::decorrelated")
plot(output3$Y, col=lab, main="FA::whitened")
par(opar)
# }
# NOT RUN {
# }