Back to PLS Help

seed PLS with fMRI resting
mbahri
Posted on 02/14/13 12:14:54
Number of posts: 26
mbahri posts:

Dear PLS Experts,

I would like to investigate the brain functional connectivity (FC) of the default network (DN) by performing a See-PLS analysis on two sets of resting state fMRI data (two runs: Rest1 and Rest2). the two runs are separated by a task. I would like to assess the spatial and temporal aspects of the FC. So, as I understood from the PLS User's Guide: my data should be divided into blocks of equal volume number and used in the "regular Behav PLS" PLS option with the seed data as behavioral input.

The extraction of the seed data need to be done with the "Multiple Voxel Extraction" option on the PLS results. So, my question is which option should be used for this step?

How I can assess the temporal aspect of the FC?

Any suggestion or/and technical detail will be of a reat help,

Many thanks in advance,

Mohamed

Replies:

Untitled Post

I'm Online
jshen
Posted on 02/14/13 12:30:30
Number of posts: 291
jshen replies:

PLS User's Guide did not say that your data need to be divided into blocks of equal volume number. In order to run Seed PLS, extraction of the seed data need to be done with the "Multiple Voxel Extraction" option on the PLS results. Then, using the extracted seed data in Regular Behavior PLS.

So the step is:

  1. Running Mean-Centering PLS, and find seed locations from activation map (Brain LV).
  2. Using above seed locations to extract voxel values with the "Multiple Voxel Extraction" option on the PLS results.
  3. Running Regular Behavior PLS. The behavior data will be the voxel values that you extracted above.

I don't know how to assess the temporal aspect of the FC.




Untitled Post
rmcintosh
Posted on 02/14/13 14:06:44
Number of posts: 394
rmcintosh replies:

Hi Mohamed, Can you provide a bit more detail on exactly what you want to do? There are many ways to estimate FC in the literature. For example, you can take a voxel time series and correlate that time series with every other voxel in the brain. That gives you a subject-specific correlation map that you can then further analyse. Another option, which is what is done in PLS, is to extract a voxel for each subject and correlate the individual differences in voxel amplitude with voxel amplitudes across the rest of the brain (and in fMRI across multiple time points/TRs). The first method is not easily done in PLS, though some users have done it on the command-line. We are planning to incorporate this in an upcoming release. So, please clarify what you want to analyze and I can give you more assistance. cheers Randy


Untitled Post
mbahri
Posted on 02/15/13 11:06:49
Number of posts: 26
mbahri replies:

Hi Randy,

Many thanks for your kind answer.
Here is a small description of our data and what we want to do. So, we have two sets of data including two sessions (day1 and day2). Each session includes two functional runs (rest1 & rest2) separated by a memory task. We would like to explore the effect of the memory task on the functional connectivity (FC) of the default network (DN). An example of contrast would be "day1_rest2 - day_rest1". So, we oriented our choice to the seed-voxel analysis based on its possibility to select only one network by using a seed belonging to the desired network (eg. PCC for the DN). As I understood, Seed-PLS could provide a statistical map depicting the difference between correlation maps (correlation maps: seed with rest1 & seed with rest2).
The idea of dividing the time series into blocks mentioned in my previous message is borrowed from the paper of  "Grigg O. and Grady Cheryl L., PloS one, Vol 5, issue 10, 2010". The time series resampling was done in order to asses the temporal aspect of FC of DN as each block represents a time point on the time series data. The used method was the Seed-PLS but I am not sure that they process all blocks in one step.

Hope that this could help us to select the suitable method.

Regarding the two options (methods) that you proposed in your answer, I am not sure to be able to say which one I should use as I am beginner in this field.

Many Thanks,

Mohamed


Untitled Post
mbahri
Posted on 02/21/13 05:49:20
Number of posts: 26
mbahri replies:

Dear PLS experts,

Can anybody have any suggestion? (see mail above).

Many thanks in advance..

Best regards,

Mohamed


Untitled Post

I'm Online
nlobaugh
Posted on 02/21/13 08:17:41
Number of posts: 229
nlobaugh replies:

Hi Mohamed..

in the Grigg paper, the full resting state run was divided into 10-sec blocks (5TRs).. All blocks were entered into the analysis.. Look more carefully at Fig 1. The brain pattern identifies the regions that show mostly positive correlations with the seed across those 10-s blocks for both RS runs. Fig2 shows the second pattern, where the correlations with the seed differed between RS1 and RS2.

So, yes, you can do this kind of analysis using PLS. You would run a Block fMRI Seed PLS, specifying block-lengths smaller than the length of the RS run. If you only want to look at run-to-run differences, you would could derive a metric outside of PLS, as Randy mentioned, or use "properly" preprocessed data and enter those into PLS, using a block length equal to the number of TRs in your study - it's hard to be more helpful than this without really knowing what the question is (Run1-Run2 is not specific enough).

Note: PLS, depending on the structure of the dataset entered, is unlikely to give you "one network" - a region can and most likely will, correlate with activity in multiple sets of regions. - see Fig4 in Grigg- only a subset of regions in LV1 and LV2 showed overlap - i.e., more than one functional network was found. cheers, nancy


Untitled Post
mbahri
Posted on 02/22/13 06:40:07
Number of posts: 26
mbahri replies:

Hi Nancy,

Thank you very much for your answer. Just one more question, how can I get the temporal structure (Figure 1B in this paper). Are these results saved somewhere in the .mat results file or should I run an individual seed-PLS for each block?

Best regards,

Mohamed


Untitled Post

I'm Online
nlobaugh
Posted on 02/22/13 07:02:08
Number of posts: 229
nlobaugh replies:

Hi Mohamed

Those are part of the results.. the "singular image" is derived from the "brainlv" (u) and the correlations for each block are derived from the "designlv" (v).

[u s v]=svd(data');

This is explained more fully in the 2004 review paper - you may want to re-read that.
n


Functional connectivity comparaison
mbahri
Posted on 11/12/13 09:37:21
Number of posts: 26
mbahri replies:

Dear PLS experts,

I have a set of data including the resting state data (fMRI) of two groups (healthy aged volunteers (CTRL) and Alzheimer patients (AD)). I have also the memory score (The memory awareness rating scale “MARS”) for all participants. I would like to answer the following questions: 1) do the Alzheimer patients differ from the healthy controls in terms of functional connectivity of the default mode network (DMN)? 2) Is there a significant correlation between the score MARS and the functional connectivity of the DMN for the Alzheimer patients?

So, I am wondering if PLS could be used to answer to these questions? If yes, it would be very kind if you could suggest a paper talking about the same topic (with PLS)?

Many thanks,

Mohamed



Untitled Post
rmcintosh
Posted on 11/12/13 10:03:51
Number of posts: 394
rmcintosh replies:

quote:

Dear PLS experts,

I have a set of data including the resting state data (fMRI) of two groups (healthy aged volunteers (CTRL) and Alzheimer patients (AD)). I have also the memory score (The memory awareness rating scale “MARS”) for all participants. I would like to answer the following questions: 1) do the Alzheimer patients differ from the healthy controls in terms of functional connectivity of the default mode network (DMN)? 2) Is there a significant correlation between the score MARS and the functional connectivity of the DMN for the Alzheimer patients?

So, I am wondering if PLS could be used to answer to these questions? If yes, it would be very kind if you could suggest a paper talking about the same topic (with PLS)?

Many thanks,

Mohamed

This is definitely something can be done with PLS.  Cheryl Grady has answered similar questions to the ones you have:

 

http://www.ncbi.nlm.nih.gov/pubmed/19789183

http://www.ncbi.nlm.nih.gov/pubmed/15989937

http://www.ncbi.nlm.nih.gov/pubmed/12921348



Untitled Post
mbahri
Posted on 11/17/13 06:53:14
Number of posts: 26
mbahri replies:

Many thanks for your answer!



fmri preprocessing
mbahri
Posted on 11/17/13 07:23:58
Number of posts: 26
mbahri replies:

Dear PLS experts,

When using univariate method for connectivity analysis (SPM), white matter and CSF signals are usually regressed out from fMRI data prior analysis. fMRI data also are bandpass filtered to keep only the neuronal component. I am wondering if these preprocessing steps also are needed when using PLS?

Any suggestion/explanation will be very welcome,

Many thanks



Untitled Post

I'm Online
nlobaugh
Posted on 11/17/13 15:50:17
Number of posts: 229
nlobaugh replies:

See new thread.



PLS for non-specialists?
mbahri
Posted on 11/20/13 05:07:16
Number of posts: 26
mbahri replies:

Dear PLS experts,

I am giving a talk in my lab in two days about the use of PLS as a tool for neuroimaging data analysis and I am wondering if anyone have a simple way (schems, example, or ...) to explain the SVD and the resulted U and V matrix and their relation to the data, to a non-specialist audiance?

Many thanks in advance,

mohamed



Untitled Post

I'm Online
jshen
Posted on 11/20/13 09:17:08
Number of posts: 291
jshen replies:

There are too many papers, articles, and books covering that topic. Here is one of them:

http://www.rotman-baycrest.on.ca/index.php?section=101

 



Y and X preprocessing steps
mbahri
Posted on 11/22/13 10:51:08
Number of posts: 26
mbahri replies:

Dear PLS experts,

Are there any documentation describing the preprocessing steps (and the methodological reasons) done on the Y (behavioral/design) and X (brain activity) matrices for each PLS variations (Task PLS, behavioral PLS, Mean-Centering PLS, Multi-bloc PLS)?

Many thanks,

Mohamed



Untitled Post

I'm Online
nlobaugh
Posted on 11/22/13 11:05:18
Number of posts: 229
nlobaugh replies:

PLS is "blind" to any choices made by the investigator with respect to preprocessing.  What to do is  your decision.

nancy



Untitled Post
mbahri
Posted on 11/22/13 11:33:07
Number of posts: 26
mbahri replies:

I want to speak about the PLS transformations (mean centering, ....) done on the data matrices (Y and X) beforme combining them on a one "correlation" matrix and not the spatial preprocessing of the data.




Login to reply to this topic.

  • Keep in touch

Enter your email above to receive electronic messages from Baycrest, including invitations to programs and events, newsletters, updates and other communications.
You can unsubscribe at any time.
Please refer to our Privacy Policy or contact us for more details.

  • Follow us on social
  • Facebook
  • Instagram
  • Linkedin
  • Pinterest
  • Twitter
  • YouTube

Contact Us:

3560 Bathurst Street
Toronto, Ontario
Canada M6A 2E1
Phone: (416) 785-2500

Baycrest is an academic health sciences centre fully affiliated with the University of Toronto