Factor density estimation using bootstrapping
res = fac_density_boot1(Z,meth,n_samples)
* Z = a (nobs x k) matrix, collecting the variables
* meth = 'stan' if the user wants to perform the analysis on standardized variables
* n_samples = an integer, the number of resamples
* res = a results tlist, with:
- res('meth') = 'factor density via bootstrapping'
- res('y') = the (nobs x k) matrix collecting the variables
- res('standardization') = a boolean, indicating whether variables have been standardized
- res('# of draws') = an integrer, the # of draws
- res('factors') = a (nobs x nb_draws) matrix, the factor density
- res('factors percentiles') = a (nobs x 5) matrix, the percentiles (2.5, 25, 50, 75 and 97.5)