Multiple matrices/tensors can be specified and decomposed simultaneously by Probabilistic Latent Tensor Factorisation (PLTF). See the reference section of GitHub README.md <https://github.com/rikenbit/gcTensor>, for details of the method.

Provides complete detailed preprocessing of two-dimensional gas chromatogram (GCxGC) samples. Baseline correction, smoothing, peak detection, and peak alignment. Also provided are some analysis functions, such as finding extracted ion chromatograms, finding mass spectral data, targeted analysis, and nontargeted analysis with either the 'National Institute of Standards and Technology Mass Spectral Library' or with the mass data. There are also several visualization methods provided for each step of the preprocessing and analysis.

Geographical detectors for measuring spatial stratified heterogeneity,
             as described in Jinfeng Wang (2010) <doi:10.1080/13658810802443457> and 
             Jinfeng Wang (2016) <doi:10.1016/j.ecolind.2016.02.052>. Includes the
             optimal discretization of continuous data, four primary functions of 
             geographical detectors, comparison of size effects of spatial unit and 
             the visualizations of results. 
             To use the package and to refer the 
             descriptions of the package, methods and case datasets, please cite 
             Yongze Song (2020) <doi:10.1080/15481603.2020.1760434>.
             The model has been applied in factor exploration of road performance and 
             multi-scale spatial segmentation for network data, as described in 
             Yongze Song (2018) <doi:10.3390/rs10111696> and 
             Yongze Song (2020) <doi:10.1109/TITS.2020.3001193>, respectively.

Datasets used in the book 'Graphical Data Analysis with R' (Antony Unwin, CRC Press 2015).

Processing collections of Earth observation images as on-demand multispectral, multitemporal raster data cubes. Users
    define cubes by spatiotemporal extent, resolution, and spatial reference system and let 'gdalcubes' automatically apply cropping, reprojection, and 
    resampling using the 'Geospatial Data Abstraction Library' ('GDAL'). Implemented functions on data cubes include reduction over space and time, 
    applying arithmetic expressions on pixel band values, moving window aggregates over time, filtering by space, time, bands, and predicates on pixel values, 
    exporting data cubes as 'netCDF' or 'GeoTIFF' files, plotting, and extraction from spatial and or spatiotemporal features.  
    All computational parts are implemented in C++, linking to the 'GDAL', 'netCDF', 'CURL', and 'SQLite' libraries. 
    See Appel and Pebesma (2019) <doi:10.3390/data4030092> for further details.