multirec {mand}

Multi components reconstruction

Description

This is a function that returns the weight vector of multiple components obtained by the msma function applied after dimension reduction by the radial basis function to the same dimension as the original image.

Usage

multirec(
  object,
  imagedim,
  B = NULL,
  mask = NULL,
  midx = 1,
  comps = NULL,
  XY = c("X", "Y", "XY")[1],
  signflip = FALSE
)

Arguments

object

an object of class msma. Usually, a result of a call to msma

imagedim

a vector for original dimension.

B

a list or a matrix.

mask

a list or a matrix.

midx

a block number.

comps

a component sequence to be plotted.

XY

a character, indicating "X" or "Y". "XY" for the scatter plots using X and Y scores from msma.

signflip

a logical if the sign in the block is flipped to pose the super as possitive.

Details

multirec requires the output result of msma function.

Examples

data(baseimg)
data(diffimg)
data(mask)
img1 = simbrain(baseimg = baseimg, diffimg = diffimg, mask=mask)
B1 = rbfunc(imagedim=img1$imagedim, seppix=2, hispec=FALSE, mask=img1$brainpos)
SB1 = basisprod(img1$S, B1)
fit111 = msma(SB1, comp=2)
ws = multirec(fit111, imagedim=img1$imagedim, B=B1, mask=img1$brainpos)