templateICA | R Documentation |
Perform template independent component analysis (ICA) using variational Bayes (VB) or expectation-maximization (EM).
templateICA(
BOLD,
template,
tvar_method = c("non-negative", "unbiased"),
scale = c("template", "global", "local", "none"),
scale_sm_surfL = NULL,
scale_sm_surfR = NULL,
scale_sm_FWHM = "template",
nuisance = NULL,
scrub = NULL,
TR = NULL,
hpf = "template",
GSR = "template",
Q2 = "template",
Q2_max = "template",
brainstructures = "template",
mask = NULL,
time_inds = NULL,
varTol = "template",
spatial_model = NULL,
resamp_res = NULL,
rm_mwall = TRUE,
reduce_dim = FALSE,
method_FC = c("VB1", "VB2", "none"),
maxiter = 100,
miniter = 3,
epsilon = 0.001,
kappa_init = 0.2,
usePar = TRUE,
PW = FALSE,
verbose = TRUE
)
BOLD |
Vector of subject-level fMRI data in one of the following
formats: CIFTI file paths, If multiple BOLD data are provided, they will be independently centered, scaled, detrended (if applicable), and denoised (if applicable). Then they will be concatenated together followed by computing the initial dual regression estimate. |
template |
Template estimates in a format compatible with |
tvar_method |
Which calculation of the template variance to use:
|
scale |
|
scale_sm_surfL , scale_sm_surfR , scale_sm_FWHM |
Only applies if
If If To create a |
nuisance |
(Optional) Signals to regress from the data, given as a
numeric matrix with the same number of rows as there are volumes in the
|
scrub |
(Optional) A numeric vector of integers indicating the indices
of volumes to scrub from the BOLD data. (List the volumes to remove, not the
ones to keep.) If multiple |
TR , hpf |
These arguments control detrending. Note that if multiple |
GSR |
Center BOLD across columns (each image)? This
is equivalent to performing global signal regression. Default:
|
Q2 , Q2_max |
Denoise the BOLD data? Denoising is based on modeling and removing nuisance ICs. It may result in a cleaner estimate for smaller datasets, but it may be unnecessary (and time-consuming) for larger datasets. Set If The defaults for both arguments is |
brainstructures |
Only applies if the entries of |
mask |
Required if and only if the entries of |
time_inds |
Subset of fMRI BOLD volumes to include in analysis, as a
vector of numeric integers. If |
varTol |
Tolerance for variance of each data location. For each scan,
locations which do not meet this threshold are masked out of the analysis.
Default: |
spatial_model |
Should spatial modeling be performed? If If If If |
resamp_res |
Only applies if |
rm_mwall |
Only applies if |
reduce_dim |
Reduce the temporal dimension of the data using PCA?
Default: |
method_FC |
Variational Bayes (VB) method for FC template ICA model:
|
maxiter |
Maximum number of EM or VB iterations. Default: |
miniter |
Minimum number of EM or VB iterations. Default: |
epsilon |
Smallest proportion change between iterations. Default:
|
kappa_init |
Starting value for kappa. Default: |
usePar |
Parallelize the computation? Default: |
PW |
Prewhiten to account for residual autocorrelation? Default: |
verbose |
If |
A (spatial) template ICA object, which is a list containing:
subjICmean
, the V \times L
estimated independent components
S; subjICse
, the standard errors of S; the
mask
of locations without template values due to too many low
variance or missing values; the nuisance
design matrix or matrices if
applicable; and the function params
such as the type of scaling and
detrending performed.
If BOLD
represented CIFTI or NIFTI data, subjICmean
and
subjICse
will be formatted as "xifti"
or "nifti"
objects, respectively.
## Not run:
tm <- estimate_template(cii1_fnames, cii2_fnames, gICA_fname)
templateICA(newcii_fname, tm, spatial_model=TRUE, resamp_res=2000)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.