dual_reg {templateICAr} | R Documentation |
Dual Regression
Description
Dual Regression
Usage
dual_reg(
BOLD,
GICA,
scale = c("global", "local", "none"),
scale_sm_xifti = NULL,
scale_sm_FWHM = 2,
detrend_DCT = 0,
center_Bcols = FALSE,
normA = FALSE
)
Arguments
BOLD |
Subject-level fMRI data matrix ( |
GICA |
Group-level independent components ( |
scale |
|
scale_sm_xifti , scale_sm_FWHM |
Only applies if |
detrend_DCT |
Detrend the data? This is an integer number of DCT bases
to use for detrending. If |
center_Bcols |
Center BOLD across columns (each image)? This
is equivalent to performing global signal regression. Default:
|
normA |
Scale each IC timeseries (column of |
Value
A list containing
the subject-level independent components S (V \times Q
),
and subject-level mixing matrix A (TxQ
).
Examples
nT <- 30
nV <- 400
nQ <- 7
mU <- matrix(rnorm(nV*nQ), nrow=nV)
mS <- mU %*% diag(seq(nQ, 1)) %*% matrix(rnorm(nQ*nT), nrow=nQ)
BOLD <- mS + rnorm(nV*nT, sd=.05)
GICA <- mU
dual_reg(BOLD=BOLD, GICA=mU, scale="local")