Compute the low-dimensional u matrix for simlr

simlr minimizes reconstruction error across related modalities. One crucial component of the reconstruction is the low-dimensional cross-modality basis. This function computes that basis, given a mixing algorithm.

simlrU(
  projections,
  mixingAlgorithm,
  initialW,
  orthogonalize = FALSE,
  connectors = NULL
)

Arguments

projections

A list that contains the low-dimensional projections.
mixingAlgorithm

the elected mixing algorithm. see simlr. can be 'svd', 'ica', 'rrpca-l', 'rrpca-s', 'pca', 'stochastic' or 'avg'.
initialW

initialization matrix size n by k for fastICA.
orthogonalize

boolean
connectors

a list ( length of projections or number of modalities ) that indicates which modalities should be paired with current modality

Value

u matrix for modality i

See also

simlr

Author

BB Avants.

Examples


set.seed(1500)
nsub = 25
npix = c(100,200,133)
nk = 5
outcome = matrix(rnorm( nsub * nk ),ncol=nk)
outcome1 = matrix(rnorm( nsub * nk ),ncol=nk)
outcome2 = matrix(rnorm( nsub * nk ),ncol=nk)
u = simlrU( list( outcome, outcome1, outcome2 ), 2, 'avg' )