This function contains some code from a function in psych::psych-package
that's not exported `print.psych.fa.ci`

but useful nonetheless. It
basically takes the outcomes of a factor analysis and extracted the
confidence intervals.

## Arguments

- fa
The object produced by the `psych::fa()`

function from the
psych::psych-package package. It is important that the `n.iter`

argument
of`psych::fa()`

was set to a realistic number, because otherwise, no
confidence intervals will be available.

## Value

A list of dataframes, one for each extracted factor, with in each
dataframe three variables:

- lo
lower bound of the confidence interval

- est
point estimate of the factor loading

- hi
upper bound of
the confidence interval

## Details

THis function extract confidence interval bounds and combines them with
factor loadings using the code from the `print.psych.fa.ci`

in
psych::psych-package.

## Examples

```
if (FALSE) {
### Not run because it takes too long to run to test it,
### and may produce warnings, both because of the bootstrapping
### required to generate the confidence intervals in fa
faConfInt(psych::fa(Thurstone.33, 2, n.iter=100, n.obs=100));
}
```