Quantcast
Channel: NITRC CONN : functional connectivity toolbox Forum: help
Viewing all 6869 articles
Browse latest View live

ICA analysis

$
0
0
Hi!

I am trying to analyze resting state fMRI data with conn.
I had run ICA, and got the results.
To find detailed description of ICA in conn (ie informax, GICA etc.), I tried to find study estimating ICA with conn.
But I couldn't.
I need it because I have to write detailed description in methods section in my paper.
Do you offer a detailed description of ICA implemented in conn?

Thank you,
Bryan

MVPA explained variance

$
0
0
Dear Alfonso,

There is some option to obtain the percentage of explained variance for each principal component in MVPA analysis?

Thanks for your help!

Roger

Questions about using the Conn software

$
0
0
Hello,

My colleague and I are learning how to use the Conn software to analyze task-fMRI data. We have used the fsl to analyze the task-fMRI data before, and we are wondering if we could import the data that have been defined in fsl (e.g., no. of condition and onset/length of the conditions of interest etc.) for the Conn analysis? Also we are interested in examining differences between task runs (i.e., five task fMRI runs). Is there a way to statistically test the differences?

Any advice would be appreciated. Thank you in advance.

Uploading FreeSurfer-processed volumes

$
0
0
Hello,

First of all, I wanted to thank you for creating such an efficient and elegant software!

I have been trying to upload to the latest CONN distribution (v.16.a) a few subjects that I have processed with FreeSurfer to analyze their resting-state fMRIs. However, I have been running into an issue that I cannot seem to solve:
After I upload the "brainmask.mgz" volume under the "Session-specific structural data," I get a message that tells me that the Freesurfer aseg.mgz segmentation files were found and whether I would like to import them. After I click yes and "OK" on the "1 files assigned to 1 subjects/sessions" message, an error message appears saying:

ERROR DESCRIPTION:

Undefined function 'isnan' for input arguments of type 'struct'.
Error in conn_menu (line 527)
title(isnan(title))=0;
Error in conn (line 1525)
conn_menu('updateimage',CONN_h.menus.m_setup_00{5},vol);
Error in conn_filesearchtool (line 150)
if length(h.callback)>1, feval(h.callback{1},h.callback{2:end},names); else, feval(h.callback{1},names); end
CONN v.16.a
SPM12 + DEM FieldMap MEEGtools
Matlab v.2014b

I was wondering whether you had any feedback on how to solve this. I apologize if it is a trivial issue and/or has been addressed before!

Thank you,
Panos

RE: Adding covariates using does not work

$
0
0
It works like a dream, thank you so much for your help (and for fixing my silly typos)!!

Temporal filtering question

$
0
0
Hi Alfonso,

I have a question about the temporal filtering processs in Conn. I read on the forum that Conn uses Fast Fourier Transformation (FFT) to bandpass filter the data. Is this correct? Also, is it the same as the original method from Biswal 1995 and Uddin 2009? I believe they used this, "temporal bandpass filtering (highpass temporal filtering: Gaussian-weighted least-squares straight line fitting, with sigma=100.0s; Gaussian lowpass temporal filtering HWHM 2.8s)"

Thanks for your help.

Best wishes,
Mike

RE: ROI problem (csf)

$
0
0
I would like to verify the WM and CSF ROIs created when using the batch.Setup.cwthreshold.  When the erosion flags are set an "e" is appended to the masks. This allows me to easily check the masks created after erosion.   However, I am unable to find a similar file when just the threshold is set. 

How can I check the masks created if the erosion step is skipped?

RE: Setup

$
0
0
Once again I am confused.  You are talking about smoothed and smoothed data.  If I wanted to do all the preprocess outside of Conn (normalization, slice timing correction, distortion correction, etc) what would I need to enter into CONN for the Setup.  My current understanding is that I would need to enter in the normalized fMRI data, the structural image, and the tissue masks (gm, wm, and csf). I would also need to include the motion parameters and the scrubbing matrix.  

Is there anything I am missing? 

Bob

2nd level analysis - ROIs lack symmetry

$
0
0
Dear Alfonso,

I was hoping you might be able to help me with a query regarding second-level analysis in CONN.

I have conducted a gPPI analysis (bivariate regression). I have a patient group, and a control group, and a block design task. I am aiming to look at differences in connectivity between the patient and control group during the task.

My question is this:

When I select ROI-A i obtain a beta and a p value for the connectivity with ROI-B. However, when I select ROI-B the values are different. I have checked the forum and one explanation I read explains why FDR corrected p values might be different. However I am seeing differences in the uncorrected values.

Any explanation would be most helpful!

Many Thanks,

Rob

RE: Errror in dynamic FC - again

$
0
0
[color=#000000]Dear Alfonso,[/color]
[color=#000000]I ran into the same issue as mentioned by Eugenio.[/color]
[color=#000000]Has there been an update on the patch?[/color]
I am using CONN 16.1 in Mac OS 10.11.4 and MATLAB R2016a.

[color=#000000]error:[/color]
[color=#000000]"Error using eig
Input to EIG must not contain NaN or Inf.[/color]

[color=#000000]thanks[/color]

[color=#000000]j.[/color]

[i]Originally posted by Alfonso Nieto-Castanon:[/i][quote][color=#000000]Hi Eugenio,[/color]

[color=#000000]Sorry about that, this error arises when using only a portion of your entire timeseries data for dynamic connectivity analyses. I am working on a patch for this and will send you that asap (hopefully later tonight)[/color]

[color=#000000]Best[/color]
[color=#000000]Alfonso[/color]
[color=#000000] 
[/color][i]Originally posted by Eugenio Abela:[/i][quote]Dear all,
 
I'm still running into following error when using the dynamic FC option - any ideas?
 
Thanks a lot,
 
Eugenio
 
ERROR DESCRIPTION:
Error using eig
Input to EIG must not contain NaN or Inf.
Error in conn_ica (line 47)
[Q,D]=eig(X*X'/(Ns-1));
Error in conn_invPPI (line 164)
[S,W]=conn_ica(reshape(B1,[Nr*Nr Nk1])',[],'dodisp',DOPLOT,'rndseed',RNDSEED);
Error in conn_process (line 2871)
[H,B,H0,B0]=conn_invPPI(Y,Ncomponents,X,Xfilter,1);
Error in conn_process (line 40)
case 'analyses_gui_dyn',disp(['CONN: RUNNING DYNAMIC CONNECTIVITY STEP']); conn_process([13.5],varargin{:});
Error in conn (line 4469)
conn_process('analyses_gui_dyn');
Error in conn_menumanager (line 119)
feval(CONN_MM.MENU{n0}.callback{n1}{1},CONN_MM.MENU{n0}.callback{n1}{2:end});
CONN v.15.h
SPM8 + Anatomy Beamforming DEM FieldMap LI MEEGtools Masking TFCE marsbar psc rwls vbm8 xjview
Matlab v.2014a
storage: 599.4Gb available  [/quote][/quote]

RE: Voxel degree centrality maps

$
0
0
Hello Arkan !

I am also trying to plot voxel degree centrality graphs using CONN. I am wondering whether CONN really has some feature which enables us to plot a degree centrality graph directly ?

I have spent a lot of time researching on this but still have not stumbled upon any reliable source to know CONN specs and the fact that whether it could support this or not.

It would be great if you could share your experience and knowledge, like how you overcame these difficulties ?

Regards

Vinay

Custom mask.surface.brainmask.nii

$
0
0
Hello CONN forum,

I have been trying to modify the default surface files (/conn/utils/surf/) to map ROI-to-ROI results onto a cat brain.  We have a surface model made in FreeSurfer and I have been simply replacing the appropriate default *.surf files with the new ones.  Results explorer correctly shows the cat brains around the Connectome Ring, however the 3D Display does not work.  After some debugging, I think the issue is that I need to also change the mask.surface.brainmask.nii file; however, I have no idea what this file actually represents.  I replaced the mask.volume.brainmask.nii file with the appropriate cat brain mask, but the mask.surface.brainmask.nii file is hard to decipher.

Anyone have experience with using custom surface files?  Anyone have any idea what the mask.surface.brainmask.nii file contains?

While solving this part will get me a bit closer to a 3D display of the cat brain, does anyone have any further advice on some additional files I would need to replace?  

Thank you in advance,
Dan

RE: MVPA explained variance

$
0
0
[color=#000000]Hi Roger,[/color]

[color=#000000]That makes sense. I would suggest, instead of analyzing each component separately (which would require a somewhat conservative multiple comparison correction to account for all those individual analyses), to simply test the first-few components (e.g. the first three or four components; just testing the first two components would be a bit biased by the results that you have already observed). [/color]

[color=#000000]To clarify, typically in MVPA analyses you want to test the first-N components (and use an eye(N) between-measures contrast) to evaluate potential between-condition and/or between-subject differences in connectivity. Yet, the choice of exactly how many components N to best include in your MVPA second-level analyses depends on the actual data. If the differences in connectivity between your groups are somewhat strong/salient (compared to the between-subjects residual differences in connectivity) they will typically be captured by the first few MVPA components, while if they are relatively weak they will typically be captured by higher-order MVPA components. So, ideally, you would want to include as many components as possible to make sure that those between-group differences in connectivity are well captured/represented by your group-agnostic MVPA components. Yet, practically, the power of your second-level analyses to detect those between-group differences in connectivity drops as you include more and more MVPA components in those analyses, mainly from the associated drop in degrees of freedom when simultaneously testing more measures. The number of components extracted by the first-level MVPA analyses (in your case 10 components, this is defined by default by attempting to keep approximately a 1:5 ratio between the number of components extracted and the number of subjects in your study) represents roughly a higher-bound on the number of components you might possibly want to test before the power drops dramatically, but it does not necessarily represent the "optimal" number of components that you would want to test (particularly if the number of subjects in your study is not too low). Typically you need to define a priori how many MVPA components you would like to include in your second-level analyses. To play it conservatively (if you do not have pilot or previous experimental data to base this descision on), it is not a bad idea to plan a couple of N values and associated second-level analyses (e.g. one using only a few components, e.g. 3 or 4, and another using more components, e.g. close to the 1:5 ratio above). In this way you would only need to apply a multiple comparison correction across these two analyses (compared for example to testing each component separately which would require a multiple comparison correction across all those individual components). In my experience 3/4 components is becoming my default choice for many MVPA analyses unless I have considerably large datasets and I know a priori that the effects that I am looking for a relatively subtle/weak (and of course unless I have very few subjects and I cannot even reasonably afford 3/4 components, in which case I use the default 1:5 ratio value), but your experience might vary. Last, again just to clarify since I have not explicitly mentioned this, because MVPA components are obtained using PCA, when one extracts M components in the first-level MVPA analysis step, and then only includes N (N<=M) components in the second-level analyses, that is exactly equivalent to having just extracted N components in the original first-level analyses and then including all of them in the second-level analyses, so the choice of N (how many components to actually include in your analyses) can be delayed to the second-level analysis step as long as you have extracted [i]at least[/i] that many components in your first-level MVPA analyses (which is why the default in CONN is set to the 1:5 ratio rule which roughly represents this higher-bound on the number of components one would typically want to include in their second-level analysis step)[/color]

[color=#000000]Hope this helps[/color]
[color=#000000]Alfonso [/color]
[i]Originally posted by Roger Mateu:[/i][quote]Thanks Alfonso! It helps so much.

I have another important question about MVPA

I'm using Conn MVPA method to extract 10 between-groups ([1 -1]) components (I have 52 subjects), which clusters must be used as a seed in a second-level analysis. I have some questions about it:

As I read in this forum, the proper way to do it is selecting all components (eye(10)) on voxel-to-voxel between-measures second-level analysis step (and [1 -1] between-group selection) , and use each resulting significant cluster as a seed. Is this correct?. This second level Seed-to-voxel analysis showed few low-signification results.

In contrast, using each significant cluster of each component independently, offered high-significant and meaningful results, specially clusters from the two first components. As you explained in this forum, individual MVPA components do not have simple interpretations, but our results suggest that first and second component clusters contain areas with divergent iFC between patients and controls (and high correlation with structural analysis differences). How can I explain this situation in order to statistically validate our results with this last method?

Thanks again

Roger[/quote]

RE: conn gui error in 2nd level results display

$
0
0
[color=#000000]Dear Yong,[/color]
[color=#000000]The latest CONN release (16b) includes a few fixes to try to better recover from this sort of scenario. Please give that a try and let me know if you still run into any issues.[/color]
[color=#000000]Best[/color]
[color=#000000]Alfonso[/color]
[i]Originally posted by Yong Li:[/i][quote]Dear Alfonso,

I encountered an error by accident. I rerun a finished voxel-to-voxel analysis by optimizing some parameters but accidentally moved the /results/secondlevel/VoxeltoVoxel/ folder to another place. Then I got:

1. the warning window-'Data not prepared for analyses: Condition (rest ) have not been processed yet...';

2. error message:
'ERROR DESCRIPTION:
Attempted to access CONN_h.menus.m_results.shownconditions(0); index must be a positive integer or logical.
Error in conn (line 5288)
nconditions=CONN_h.menus.m_results.shownconditions(nconditions);
Error in conn_menumanager (line 119)
feval(CONN_MM.MENU{n0}.callback{n1}{1},CONN_MM.MENU{n0}.callback{n1}{2:end});

CONN v.16.a
SPM12b + DEM FieldMap MEEGtools suit vbm8
Matlab v.2013a
storage: 150.5Gb available'

3. in the conn gui second level results, only the 'rest' condition of the 'Conditions' column was gone and all contents in the 'Subject effects' and Voxel-to-Voxel Measures are listed correctly.

This persists after moving back the misplaced folders and even I rerun the first level analysis and overwrite all subjects' data.

My guess is that there might exist a bug in the conn gui menu during showing second level results? Any fix? Or if it is not a bug, any idea why and how to avoid it? Thank you in advance!

Best regards,

Yong[/quote]

best way to close a conn project?

$
0
0
Dear Alfonso

Sorry, I have a somewhat silly question... When I work on a conn project using the GUI, and want to upload or start a different project, what's the right way to do it so that my project will not be run over or loose data? Just click the "load"  or "new" buttons?  There's no "exit" or "close" button that I could find in order to first close a file before uploading another one. 

Thanks!!

L.

RE: Uploading FreeSurfer-processed volumes

$
0
0
[color=#000000]Hi Panos,[/color]

[color=#000000]This issue has been fixed in the latest release (currently 16b). Please give that a try and let me know if you still run into any issues[/color]
[color=#000000]Best[/color]
[color=#000000]Alfonso[/color]
[i]Originally posted by Panagiotis Fotiadis:[/i][quote]Hello,

First of all, I wanted to thank you for creating such an efficient and elegant software!

I have been trying to upload to the latest CONN distribution (v.16.a) a few subjects that I have processed with FreeSurfer to analyze their resting-state fMRIs. However, I have been running into an issue that I cannot seem to solve:
After I upload the "brainmask.mgz" volume under the "Session-specific structural data," I get a message that tells me that the Freesurfer aseg.mgz segmentation files were found and whether I would like to import them. After I click yes and "OK" on the "1 files assigned to 1 subjects/sessions" message, an error message appears saying:

ERROR DESCRIPTION:

Undefined function 'isnan' for input arguments of type 'struct'.
Error in conn_menu (line 527)
title(isnan(title))=0;
Error in conn (line 1525)
conn_menu('updateimage',CONN_h.menus.m_setup_00{5},vol);
Error in conn_filesearchtool (line 150)
if length(h.callback)>1, feval(h.callback{1},h.callback{2:end},names); else, feval(h.callback{1},names); end
CONN v.16.a
SPM12 + DEM FieldMap MEEGtools
Matlab v.2014b

I was wondering whether you had any feedback on how to solve this. I apologize if it is a trivial issue and/or has been addressed before!

Thank you,
Panos[/quote]

RE: Denoising Residual Degrees of Freedom

$
0
0
Hi Alfonso,

Thank you for all the helpful information you have already provided on this subject, however I am very inexperienced with MRI preprocessing so still have a couple of queries. I have just run ART based scrubbing (at default conservative thresholds) on my 30 resting state control subjects that we scanned for 5 minutes and have a highly variable number of outliers (most commonly between 0 - 7  outliers per subject but at it's most extreme I have 60 outliers and 30 outliers per subjects (I figure I should probably exclude these)), however in my denoising step looking at the histograms some subjects, even those that seem to have very few outliers will have low or even 0 residual degrees of freedom. Does this mean my thresholds were too harsh for the sample and at what point should I exlude subjects? I have attached my histogram and .mat file.

Thank you,
Isabella

Conn Conjunction analysis

$
0
0
This is a simple question, but I am looking for a way to find statistically significant connections that occur in both of my two groups.  I don't want any statistically significant connection that is in one group, but not the other.  

Any help would be much appreciated!  

Thanks,
Jane

Error loading functional data

$
0
0
Hi,


I preprocessed my data in SPM and tried to select them in Conn.
However, every time I select the 197 preprocessed files for one subject I receive following error:
I already tried selecting different preprocessed images (preprocessed with the use of fsl), but receive the same error.
Does anybody know the solution to this?

(Note: I don't get this error when I only select a single file)


ERROR DESCRIPTION:

Error using flipud (line 19)
X must be a 2-D matrix.
Error in conn_menu (line 673)
title=fliplr(flipud(data.x1(:,:,data.n,:)));
Error in conn (line 1288)
conn_menu('updateimage',CONN_h.menus.m_setup_00{5},CONN_x.Setup.functional{nsub}{nses}{3});
Error in conn_filesearchtool (line 150)
if length(h.callback)>1, feval(h.callback{1},h.callback{2:end},names); else, feval(h.callback{1},names); end
CONN v.16.b
SPM8 + Beamforming DEM FieldMap MEEGtools
Matlab v.2012b

RE: Second level random effects analysis

$
0
0
Hi Alfonso Nieto-Castanon,

I am not able to find the "launch REX gui". I have updated my CONN toolbox to the most recent version I believe (16b) and would like to extract the connectivity values (transformed correlation) for a particular region, per subject and subject effect.

Please let me know how to proceed. Thanks again, love this awesome toolbox!

Kevin
Viewing all 6869 articles
Browse latest View live