This work introduces a novel method to multivariate analysis applied to fused hyperspectral datasets in the field of Cultural Heritage (CH). Hyperspectral Imaging is a well-established approach for the non-invasive examination of artworks, offering insights into their composition and conservation status. In CH field, a combination of hyperspectral techniques is usually employed to reach a comprehensive understanding of the artwork. To deal with hyperspectral data, multivariate statistical methods are essential due to the complexity of the data. The process involves factorizing the data matrix to highlight components and reduce dimensionality, with techniques such as Non-negative Matrix Factorization (NMF) gaining prominence. To maximize the synergies between multimodal datasets, the fusion of hyperspectral datasets can be coupled with multivariate analysis, with potential applications in CH. In this work, I will show examples of this approach with different combinations of datasets, including reflectance and transmittance spectral imaging, Fluorescence Lifetime Imaging and Time-Gated Hyperspectral Imaging, and Raman and fluorescence spectroscopy micro-mapping.