DMMD-package: Double Matched Matrix Decomposition

DMMD-packageR Documentation

Double Matched Matrix Decomposition

Description

Extract joint and individual signals from double matched multi-view matrices <doi: 10.1080/10618600.2022.2067860>.

Details

The DESCRIPTION file: This package was not yet installed at build time.
Index: This package was not yet installed at build time.
The goal of DMMD package is to extract joint and individual signals from double-matched multi-view matrices. The main function is DMMD_Fit, which can do double-matched matrix factorization efficiently. DMMD_i is the updated version of DMMD_Fit which also updates the joint structure but is less efficient. The rank selection is based on the method of profile likelihood or edge distribution. There is a function of DoubleDataGen which can generate double-matched matrices that can be used to test functions.

Author(s)

Dongbang Yuan and Irina Gaynanova
Maintainer: Dongbang Yuan <yuandb09@gmail.com>

References

Dongbang Yuan & Irina Gaynanova (2022) Double-Matched Matrix Decomposition for Multi-View Data, Journal of Computational and Graphical Statistics, DOI: 10.1080/10618600.2022.2067860

See Also

Feng, Qing, et al. "Angle-based joint and individual variation explained." Journal of multivariate analysis 166 (2018): 241-265.

Lock, Eric F., et al. "Joint and individual variation explained (JIVE) for integrated analysis of multiple data types." The annals of applied statistics 7.1 (2013): 523.

Zhu, Mu, and Ali Ghodsi. "Automatic dimensionality selection from the scree plot via the use of profile likelihood." Computational Statistics & Data Analysis 51.2 (2006): 918-930.

Examples

data = DoubleDataGen(n = 20, p = 16, rank = c(4, 3), rc = 2, rr = 1, nrep = 1)
X1 = data$X1[[1]]
X2 = data$X2[[1]]

svd_x1 = svd(X1)
svd_x2 = svd(X2)
r1 = ProfileLikCluster(svd_x1$d)$index
r2 = ProfileLikCluster(svd_x2$d)$index

result_DMMD = DMMD_Fit(X1,X2)
result_DMMD$rc
result_DMMDi = DMMD_i(X1,X2)

justicesuker/DMMD_Package documentation built on Aug. 6, 2022, 12:34 p.m.