Hi there!
Hopefully this is an easy question for you. Could you use a func-bucket dataset which contains 64 coefficients (beta coefficients) the same way as a 3d+time dataset in 3dTcorr1D? How do you convert the two where the func-bucket dataset becomes a 3d+time with 64 coefficient timepoints. I want to correlate coefficients instead of timepoints unlike a normal timeseries. I also thought of dumping out the coordinates and 64 coefficients and use Matlab to correlate and use 3dUndump to ‘put’ it back into AFNI. What I would like is to use clusterize program and cluster the correlation coefficient data–would a statistic be necessary for me to do this? Any ideas and comments will help.
Thanks!
-Linda
Hi Linda,
A beta weight dataset should behave the same way as an MRI
time series dataset (though you might have to consider a
demean operation, depending on how the correlations are
computed).
Is there something that you are trying that is not working?
Would you provide some corresponding commands?
- rick
Here is my program where I resample the ROIs into the space of the coefficient file: DeconIMGAMPPI+tlrc. Then I dump out the middle frontal coefficients
into a file called tempmiddle and subsequently transpose the matrix. Finally, I correlate every middle set of coefficients with every set in the Insula giving 1024 Briks but there were I believe 1883 coefficient sets–is 1024 a cutoff? Anyway, it seemed to work giving different correlation values in the Insula ROI which makes me happy. I just need to prove the correlations are correct. I used a pearson correlation. Thanks for your help! -Linda
3dresample -master DeconIMGAMPPI+tlrc -prefix Rmiddlenew -input Rmiddle+tlrc
3dresample -master DeconIMGAMPPI+tlrc -prefix RInsulanew -input RInsula+tlrc
3dmaskdump -mask Rmiddlenew+tlrc -noijk DeconIMGAMPPI+tlrc>tempmiddle.txt
1dtranspose tempmiddle.txt tempmiddle-transpose.txt
3dTcorr1D -pearson -mask RInsulanew+tlrc -prefix middleresults DeconIMGAMPPI+tlrc tempmiddle-transpose.txt
Hi Linda,
To be sure, what is the output of:
1d_tool.py -show_rows_cols -infile tempmiddle.txt
or maybe more directly:
3dROIstats -quiet -nzvoxels -mask Rmiddlenew+tlrc Rmiddlenew+tlrc
However it seems like what you want could more directly
come from 3dAutoTcorrelate. Consider something like:
3dAutoTcorrelate -polort -1 -mask_source RInsulanew+tlrc
-mask Rmiddlenew+tlrc -prefix test.AC DeconIMGAMPPI+tlrc
- rick
Hello Rick,
What I want is correlation between voxels of two different ROIs. What I get out of my analysis is a BRIK for a given voxel in the Rmiddle correlated to each voxel in the Insula. These BRIKs I would like to put each into an FDR giving a list which can be thresholded by q-value leaving the remaining “active” voxels- ones with significant correlations. One can then convert the BRIK values to binary and then one ‘adds’ a given voxel across BRIKs. This results in a ‘count’ for each voxel in the Insula where the larger the count the more significant that voxel is for the correlation of the two brain regions. What is hoped is that it will result in a cluster in the Insula and simiilarily in the Rmiddle frontal. There should be a connection between these regions. Unfortunately, the Rmiddle frontal is huge with 2883 voxels and when I run this there is a limit of 1025 BRIKs to be stored in the bucket. Can this be changed or else I may have to remove parts of the ROI which is the standard one in Talariach space. Although I find a lot of zero voxels. I don’t know if they can be deleted without the other voxels losing there ‘place’ in the list so to speak. Thanks for all your help. -Linda
Hi Linda,
I have not yet found where that 1024 limit is coming from.
But in the meanwhile, it does seem like that 3dAutoTcorrelate
command should do the same thing. Run it and just compare
the first few resulting volumes (hopefully the orders match).
3dAutoTcorrelate -polort -1 -mask_source RInsulanew+tlrc \
-mask Rmiddlenew+tlrc -prefix test.AC DeconIMGAMPPI+tlrc
A related question, why there are so many voxels in each ROI?
What is the resolution of the data?
- rick