arf3DS4-package: Activated Region Fitting fMRI data analysis
Description
Activated Region Fitting (ARF) is an analysis method for fMRI data. The ARF method uses Gaussian shaped functions to model active brain regions.
An entire fMRI volume can be described by multiple parameterized Gaussian regions of activation. This parameterization allows for hypotheses on the
location of an active region, the spatial extent of an active region, and the amplitude of an active region. ARF can be invoked after standard GLM analysis,
and uses the beta-values or t-values from this analysis. In addition ARF can be used to estimate functional connectivity, for this the raw time-series
(on which the GLM was performed) are also needed.
The arf3DS4 package can read and write fMRI data in Nifti format (fmri.data
). It is build to be compatible with fMRI analysis packages (e.g. FSL or SPM). The package
uses a predefined directory- and file-structure to store its data. All objects (S4-class) and fMRI data files are stored in this directory- and file-structure.Quick start
The main directory- and file-structure (termed an experiment) is handled by an object of class experiment
. To make an empty experiment structure
call makeExpDirs
. This will create a directory- and file-structure to hold the fMRI data (which will have to be copied manually to the appropriate directories, see data
for an overview of the data-directories ).
Once this is done an experiment is loaded by calling loadExp
.
To fit the ARF models (class model
) to the data, first create a model using newModel
, adjust the options of the model via the options
-object using loadOptions
and saveOptions
,
and finally fit the model using fitModel
. To perform hypothesis tests call varcov
and wald
to calculate the parameter (co)variance matrix and Wald statistics respectively.
In addition to the hypothesis test connectivity between brain regions can be estimated (call fitConnectivity
to perform this analysis). This function requires that the
raw time-series of the fMRI data are available and that single-trial estimates (call makeSingleTrialEvents
) are calculated.Example data
The package comes with an example dataset (load it using data("arf-example-data")
). By then calling makeExample
, a directory-structure containing
the example data is then created and loaded. Call getExp
to view the data within the experiment.Details
ll{
Package: arf3DS4
Type: Package
Version: 2.5-4
Date: 2011-10-28
Depends: R (>= 2.12.0), methods, graphics, tcltk
Repository/R-Forge/Revision: 203
Publication: 2011-10-28 12:38:50
Packaged: 2011-10-28 12:38:50; rforge
Built: R 2.12.0; universal-apple-darwin9.8.0; 2011-10-28 12:38:50; unix
LazyLoad: yes
}References
l{
Wouter D. Weeda, Lourens J. Waldorp, Ingrid Christoffels and Hilde M. Huizenga (2009). Activated Region Fitting: A Robust High-Power Method for fMRI Analysis Using Parameterized Regions of Activation. Human Brain Mapping, 30(8), 2595-2605.
Wouter D. Weeda, Lourens J. Waldorp, Raoul P.P.P. Grasman, Simon van Gaal and Hilde M. Huizenga (2011). Functional Connectivity Analysis of fMRI Data Using Parameterized Regions of Interest. NeuroImage, 54(1), 410-416.
Wouter D. Weeda, Frank de Vos, Lourens J. Waldorp, Raoul P.P.P. Grasman, Hilde M. Huizenga (2011). arf3DS4: An Integrated Framework for Localization and Connectivity Analysis of fMRI Data. Journal of Statistical Software, 44(14), 1-33. URL http://www.jstatsoft.org/v44/i14/
}