hi appu....
i agree with christian...just an add on...we've also tried the
approach similar to what christian is suggesting (running separate ica
estimations)... except not in fsl...see: calhoun et al, neuroimage 14(5)
2001 pp 1080-1088. hope this helps.
cheers,
vince
Christian Beckmann wrote:
> Date: Fri, 20 Dec 2002 17:28:40 +0000
> From: Christian Beckmann <[log in to unmask]>
> Subject: Re: [FSL] ICA
> To: [log in to unmask]
>
> Hi Appu,
>
> Appu Mohanty wrote:
>> I came accross two ways of doing ICA with fMRI data. In one method, ICA is
>> performed directly on the fMRI time series data resulting in components.
>> The ICA derived components are then correlated with the design matrix to
>> identify the task-related components. In the other method the fMRI times
>> series data is cross correlated with the design matrix and the ICA is
>> performed on the result of the cross correlation to yield eigenimages. I
>> think both methods could be implemented in FEAT. Are there any suggestions
>> about which might be a better method? Is there any other method that is
>> commonly used?
>>
>
> ICA on fmri data is normally either done in the temporal or the spatial
> domain, the later being what is more common and what is offered with
> melodic. Calculating correlation scores first and performing ICA on the
> output strikes me as not very sensible: when you correlate against a
> pre-defined time-course I'd expect you to have an associated 'meaning'
> for each regressor. Those regressor time-courses jointly form a basis
> set that encodes your believe in the temporal characteristics of signal
> in the data and it is not at all obvious why if you're willing to
> strictly constrain your signal representation in the first place (by
> pre-specification of regressors) and then mess everything up again by
> looking for a different basis representation that now asks for spatial
> independence. It seems you loose exactly what can make ICA a powerful
> technique (the exploratory waythat a basis representation is generated)
> by forcing the result to be within the span of your original regressors.
>
>
>> My other question is about how to do an ICA on a group of subjects.
>> The following is a well established (though not very common) procedure
>> in the event-related brain potential literature. Has a similar procedure
>> been used in the fMRI literature? Is it (or: how much of it is) implemented
>> in FSL? If this procedure is not used then what is usually recommended for
>> doing ICA on a group of subjects?
>>
>> Assume that the data set for the entire experiment consists of
>> multiple subjects, multiple conditions per subject, multiple electrodes,
>> and multiple time-point observations for each electrode. (The electrode
>> can be seen as corresponding to MR voxels, and the time points to the MR
>> time series across TRs.)
>> Put all of the data into a single matrix, with columns = time
>> points and rows = everything else, not distinguishing subjects,
>> conditions, or electrodes. Conduct a PCA or ICA, producing factors and
>> factor loadings as a function of time.
>> For the time-point vector in each subject x condition x electrode
>> cell, apply the factor loading to the vector, producing a factor
>> score. The factor scores are suitable dependent variables in an ANOVA
>> (or whatever). The ANOVA could have condition and electrode as factors,
>> with subjects as the within-cell replications to estimate
>> error. Rejecting the null hypothesis means that the factor scores vary
>> systematically as a function of conditions and/or electrodes.
>>
>> All of the above can also be done by exchanging electrodes and time
>> points, so that with columns = electrodes and rows = everything
>> else. The result is that the factor loadings are a function of electrode
>> location rather than time.
>>
>>
>
>
> Group analysis is a difficult thing; yes, you can always concatenate in
> space or time or whatever to form a new single data set - you still need
> to think hard about what domain you then enforce independence and if
> that really makes sense. If you concatenate subjects/runs in time you'll
> end up with an analysis that asks for a single spatial map but will
> allow for different temporal characteristics for each subject/run and
> vice versa if you concatenate in space. You can do both using avwmerge
> on the data but either way there are many problems with these
> approaches: If you concatenate in space you'll enforce a single
> time-course that is the same for each subject. Even in the case of very
> simple experiments the temporal characteristics of the hdrf can vary
> between subjects though. Conversely, when you concatenate in time you'll
> enforce a spatial maps for any one of the subjects. This can potentially
> be very restrictive wrt processes that do occur in some of the subjects
> but not others. An alternative has been suggested in
>
> Calhoun V, Adali T, Pearlson G, Pekar J. A Method for Making Group
> Inferences from Functional MRI Data Using Independent Component
> Analysis. Hum Brain Map 2001; 14: 140-151.
>
> This, however, has the added complication of generating an analysis
> where at the first PCA stage you allow for a completely different
> separation into what will be the noise and what will be the signal+
> noise space retained for further analysis.
>
> With melodic 2 I'd personally do the following
> - run each subject/ run separately
> - find the relevant spatial maps (if there's more than one you're in
> trouble, a problem which is hidden using any of the approaches above at
> different expenses)
> - coregister all low-res brains to the same reference image and apply
> the transformation to each of the relevant probability map (these are
> saved using the --Ostats option at the command line or the 'Output full
> stats folder' using the GUI)
> - at this point you have a set of probability maps of that you can
> average and overlay ontop of the reference image
>
> hope this all makes sense
> happy christmas everyone
> Christian
>
>
> --
> Christian F. Beckmann
> Address: Oxford University Centre for Functional
> Magnetic Resonance Imaging of the Brain,
> John Radcliffe Hospital, Headington, Oxford OX3 9DU, UK
> Email: [log in to unmask] http://www.fmrib.ox.ac.uk/~beckmann/
> Phone: +44(0)1865 222782 Fax: +44(0)1865 222717 Mob: +44(0)7811 189123
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
Vince D. Calhoun, Ph.D.
Director, Medical Image Analysis Laboratory
Olin Neuropsychiatry Research Center, Institute of Living
Assistant Clinical Professor
Department of Psychiatry
Yale University
Assistant Professor (part time)
Department of Psychiatry
Johns Hopkins University
MAILING ADDRESS:
Institute of Living
200 Retreat Avenue
Center Building, 2nd Floor
Hartford, CT 06106
Tel: (860) 545-7768, Fax: (860) 545-7510
Email: [log in to unmask]
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
"An author [or researcher] should never conceive himself as bringing into
existence beauty or wisdom which did not exist before, but simply and
solely as trying to embody in terms of his own art some reflection of
eternal Beauty and Wisdom...And always, of every idea and of every
method the Christian will ask not, 'Is it mine?' but 'Is it good?'"
----CS Lewis----
|