En este trabajo, presentamos un método para construir una matriz de densidad reducida de primer orden (1-RDM) de una molécula a partir de cálculos de Monte Carlo cuántico variacional (VMC) mediante una función de onda de mapeo correlacionada dada. Dicha función de onda está modelada en una forma de Interacción de Configuración de Espacio Activo Completo de Enlace de Valencia Generalizado y se ajusta mejor a la densidad resultante de la función de onda de Slater-Jastrow de VMC. La precisión del método propuesto ha sido demostrada al comparar la energía cinética resultante con el valor correspondiente de VMC. Este 1-RDM se utiliza para analizar la cantidad de correlación capturada eventualmente en cálculos de Kohn-Sham realizados en un enfoque no restringido (UKS-DFT) y con diferentes funcionales de energía. Realizamos cálculos de prueba en un conjunto seleccionado de moléculas que muestran un carácter multireferencial significativo. En este análisis, comparamos tanto indicadores locales como globales de correlación no dinámica y dinámica. Además, siguiendo la descomposición orbital natural del 1-RDM, también comparamos las temperaturas efectivas de las distribuciones tipo Fermi correspondientes. Aunque hay un acuerdo general entre UKS-DFT y VMC, encontramos la mejor coincidencia con el funcional LC-BLYP.