# NOT RUN {
varfac <- factor.model.stat(retmat, zero=TRUE, output="factor")
# perhaps modify loadings and uniquenesses for missing values
varian2 <- fitted(varfac) # get variance matrix from factor model
specif <- fitted(varfac, output="specific") # diagonal matrix
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab