Sunday, 21 December 2014

PPI and gPPI

PPI and gPPI

A PPI analysis starts with an ROI and a design matrix. It's a way of searching among all other voxels in the brain (outside the seed ROI) for regions that are highly connected to that seed. One of the most straightforward ways of doing connectivity analyses would be to start with one ROI and simply measure the correlation of all other voxels in the brain to that voxel's timeseries, looking for high correlation values. As Friston and other pointed out a while ago, though, it's not quite as interesting if the correlation between two regions is totally static across the experiment - or if it's driven by the fact that they're both totally non-active during rest conditions, say. What might be more interesting is if the connection strength between a voxel and your seed ROI varied with the experiment - i.e., there was a much tighter connection during condition A between these regions than there was during condition B. That may tell you something about how connectivity influences your actual task (and vice versa).
PPIs are relatively simple to perform; you extract the timeseries from a seed voxel or ROI/VOI and convolve it with a vector representing a contrast in your design matrix (say, A vs. B). You then put this new PPI regressor into a general linear model analysis, along with the timeseries itself and the vector representing your contrast; you'll use those to soak up the variance from the main effects, which you'll ignore in favor of the PPI interaction term. When you estimate the parameters of this new GLM, the voxels where the PPI regressor has a very high parameter are those who showed a signficant change in connectivity with your experimental manipulation.
PPIs are good to do if you have one ROI of interest and want to see what's connected with it. They're tricky to interpret, and they can take a really long time to re-estimate if you have several ROIs to explore and many subjects.


Jung D, Sul S and Kim H (2013) Dissociable neural processes underlying risky decisions for self versus other. Front. Neurosci. 7:15. doi: 10.3389/fnins.2013.00015

http://mindhive.mit.edu/book/export/html/58
 

Friday, 31 October 2014

BROCCOLI

BROCCOLI is an open-source parallel processing software package for fMRI analysis. The goal is to highlight the advantages and limitations of parallel CPU and GPU-based processing. BROCCOLI is written in OpenCL (Open Computing Language).

---------------------------

Motion correction


---------------------------

Normalization, normalize


Normalization, segment


---------------------------

First level analysis, through a batch script (modified from http://www.fil.ion.ucl.ac.uk/spm/data/face_rep/face_rep_spm5_batch.m)


---------------------------

For further details, see our github repository. Also please feel free to try out our software (see the beta testing folder) and comment on that as well.

I am testing it now on resting state data and the author of BROCCOLI

Anders Eklund is very helpfull!!

Friday, 14 March 2014

AlphaSim

AlphaSim parameters

AlphaSim is a Monte Carlo simulation. 

Schematic representation

1. FWHM: 4
2. rmm: 5;
rmm is used to define connection cretiria. For example, when voxel size = 3*3*3, rmm = 4: face connection; rmm = 5: edge connection (SPM use this); rmm = 6 corner connection (FSL use this).
3. threshold: 0.01
4. iterations: 1000
5. mask:image

Ref: http://afni.nimh.nih.gov/pub/dist/doc/manual/AlphaSim.pdf

Wednesday, 5 March 2014

MR imaging: basics and physics

Video to prepare children for MRI scanning:
http://www.youtube.com/watch?v=LaAjrPbahBA

MR imaging
http://www.youtube.com/watch?v=jWRIKNeCXjI
http://www.youtube.com/watch?v=MiL0wCZr0Mw
http://www.youtube.com/watch?v=Ok9ILIYzmaY


fMRI course by Prof. Geoffrey Aguirre
1. http://www.youtube.com/watch?v=vGLd-bUwVXg
2. http://www.youtube.com/watch?v=J1XYcIj86TI
3. http://www.youtube.com/watch?v=_Qo44isGcxw
4. http://www.youtube.com/watch?v=vR8jwPMTick










Events 2014

2014

Brain Connectivity Workshop 2014
http://sfb936.net/index.php/events/brain-connectivity-workshop-2014

OHBM 2014, Hamburg

Resting state/Brain connectivity conference 2014

CAOS Rovereto 2014

LNA conference 2014

2015

OHBM 2015








fMRI preprocessing and functional connectivity


I. fMRI preprocessing

a) Acquisition and Quality control
b) What do we measure?
c) General preprocessing (RS-fMRI, taks fMRI)
d) RS-fMRI additional steps
e) GLM for task-fMRI


II. Functional connectivity

a) seed-based connectivity
b) ROI-based connectivity
c) ICA

************

I. fMRI preprocessing

Toolboxes for preprocessing: SPM, FSL, AFNI

Toolboxes for functional connecitivity: REST, GIFT, CONN.


SPM is widely used, a lot of documentation, wikibook, tutorials, datasets etc.




SPM is based on Matlab scripts and compatible with Windows, Mac and Linux.
User inteface: GUI or command line (also Batch).


1. How to start?
Install Matlab and SPM8 toolbox.


Reference:


2. Imaging data provided in DICOM format

● SPM uses the NIFTI (.nii) format

● conversion to NIFTI from DICOM


Matlab has a function:
SPM -> Import dicom
MRIcro http://www.mccauslandcenter.sc.edu/mricro (dcm2niigui)

3. Read header information


repetition time (TR)


number of slice in a volume
Order of the scanning sequence (ascending, descending, interleaved)

4. fMRI analysis is performed in 3 steps:
  • Preprocessing
  • 1st level analysis
  • 2nd level analysis
Preprocessing:

slice timing 
realign

coregister
segment
normalize
smooth

Slice timing - 

Realigment = motion correction


Coregister two modalities: T1 and T2

Segmentation (AC/PC), GM, WM


Normalization to MNI space

Smoothing -

5. Resting state fMRI 
  • detrending
  • filtering
Detrending - scanner drift removal.

Filtering [0.01-0.1]

6. First level analysis (GLM)

Block design (3 blocks)


Even-related design (3 events)




Tuesday, 18 February 2014

Quality check

BOLD effects of interest are small, so temporal stability during functional acquisition is important. In order to accurately measure such small signal changes, an MR system must have intrinsic image time series fluctuation levels much lower than these expected signal changes. Quality checks allow you to assess if your data are worth being analyzed.

Acquisition (artifacts)

  • SNR - signal to noise - single image (contrast to noise, dropout/susceptibility; gray/white matter contrast in structurals)
  • SNR - temporal (e.g., maps of mean / std. deviation across time)
  • ghosting - image wrap-around artifact*

  • image intensity - high enough to avoid information loss? (e.g., >1000)
  • distortions - spatial distortions in acquisition

    Motion artifacts mainly propagate in the phase-encode direction. This is due to movement of the spins between 2 excitations or between phase-encoding and signal reading: in the first case, the spins will not be recorded at the same position between excitations, in the second case, their phase-encoding will not be correct. As a result, the phase-encoding of these voxels is corrupted and this will be responsible for artifacts in the phase-encode direction.
    On the other hand, signal sampling and spatial-encoding in the frequency-encode direction are done so fast that physiological motion will only produce a small amount of spatial blurring in that direction.

  • Spikes - gradient artifacts and bad images *

  • drift - large signal drift over time (look at FFT or timeseries)
  • periodic noise - low-frequency noise artifacts

    When the movements are periodical (cardiac beats, arterial or CSF pulsations, respiration), they can produce ghost images, propagated in the phase-encode direction, even outside the anatomic limits. The intensity of these ghost images becomes more extreme with the intensity of the moving structure and with the amplitude of movement. These ghost images can show up as an increase or decrease of the true image signal.
    The spacing between ghost images varies with the direction of the movement, its amplitude and its periodicity relative to the phase-sampling interval (TR).

  • streaks/striping in images (RF room leaks)

Processing

  • orientation (reconstruction/header problems)
  • skull stripping or segmentation failures (missing brain)
  • coregistration (anatomical to functional overlap)
  • normalization (warping was ok? bad alignment/distortion?)
  • movement estimates (reasonable? too much movement?)

Modeling/statistics

  • colinearity and predictor variance in model
  • mask (which voxels are analyzed?)
  • image registration (consistent across all images in analysis?)
  • contrast scaling (consistent across all images?)
  • indicators of task-correlated artifacts: global shift, large std across contrast
  • outliers (very unusual images)

    Reference: http://wagerlab.colorado.edu/wiki/
    http://www.imaios.com/en/e-Courses/e-MRI/Image-quality-and-artifacts/image-quality