Comments (11)
Is this going to be supported in the near future? Is there an alternate way to mimic this with current BROCCOLI version?
from broccoli.
Yes another user is also requesting it.
2015-08-17 15:09 GMT+02:00 ImmanuelSamuel [email protected]:
Is this going to be supported in the near future?
—
Reply to this email directly or view it on GitHub
#13 (comment).
Anders Eklund, PhD
from broccoli.
Any idea on when you will be able to release it? I need to meet my deadline in a month and so i might start using FSL if this is going to take a long time.
from broccoli.
I will try to include it this week.
2015-08-17 15:43 GMT+02:00 ImmanuelSamuel [email protected]:
Any idea on when you will be able to release it? I need to meet my
deadline in a month and so i might start using FSL if this is going to take
a long time.—
Reply to this email directly or view it on GitHub
#13 (comment).
Anders Eklund, PhD
from broccoli.
I have now added support for several runs, and hope that I have found all
the new bugs.
To get the update, do a "git pull" in a linux terminal, while being in a
BROCCOLI directory. If you can't use the precompiled functions, go to
BROCCOLI/code/Bash_Wrapper and run recompile.sh
For each run, 4 detrending regressors are added that are only "ON" for the
duration of that run. Each run can have a separate number of volumes. An
arbitrary number of runs is supported. The general syntax is
FirstLevelAnalysis -runs 3 run1.nii run2.nii run3.nii T1.nii MNI.nii
regressors1.txt regressors2.txt regressors3.txt contrasts.txt
(if you only have 1 run you use the old syntax FirstLevelAnalysis run1.nii
T1.nii MNI.nii regressors.txt contrasts.txt)
If there is one or several runs for which the subject did not do a certain
task, you can simply add a dummy regressor text file, e.g. if you have 3
runs and the subject performed 2 tasks, but only performed one task for the
last run, the regressor files can look like
regressors1.txt:
NumRegressors 2
taskA_run1.txt
taskB_run1.txt
regressors2.txt:
NumRegressors 2
taskA_run2.txt
taskB_run2.txt
regressors3.txt:
NumRegressors 2
dummy.txt
taskB_run3.txt
where dummy.txt only contains "NumEvents 0"
Make use of -savedesignmatrix and -saveoriginaldesignmatrix to check if the
design matrix looks as expected (e.g. 1dplot total_designmatrix.txt, where
1dplot is an AFNI function)
2015-08-17 16:01 GMT+02:00 Anders Eklund [email protected]:
I will try to include it this week.
2015-08-17 15:43 GMT+02:00 ImmanuelSamuel [email protected]:
Any idea on when you will be able to release it? I need to meet my
deadline in a month and so i might start using FSL if this is going to take
a long time.—
Reply to this email directly or view it on GitHub
#13 (comment).Anders Eklund, PhD
Anders Eklund, PhD
from broccoli.
Would it be possible to add an ouput directory argument where the results can get saved?
from broccoli.
That option exists for other BROCCOLI function, but for some reason not yet
for FirstLevelAnalysis. I will add it.
2015-08-20 0:21 GMT+02:00 ImmanuelSamuel [email protected]:
Would it be possible to add an ouput directory argument where the results
can get saved?—
Reply to this email directly or view it on GitHub
#13 (comment).
Anders Eklund, PhD
from broccoli.
I have now added -output option for FirstLevelAnalysis, if you for example run
FirstLevelAnalysis -runs 3 bold1.nii.gz bold2.nii.gz bold3.nii.gz brain.nii.gz MNI152_T1_2mm_brain.nii.gz regressors1.txt regressors2.txt regressors3.txt contrasts.txt -output /home/andek/Data/subject1
all the results will be saved in the directory /home/andek/Data , and all the files will start with subject1
from broccoli.
That works well. How do you do higher level analysis after this. FSL asks for registration before running higher level. Since the directory structure is different it complains.
from broccoli.
I don't know if you can run the group analysis through the FSL GUI, the
registration has already been done in the FirstLevelAnalysis (all results
are in MNI space).
You can use randomise in FSL to run a permutation test. First merge the
beta/cope image for each subject into a single nifti file
fslmerge -t allcopes.nii.gz subject*_cope1.nii.gz
then run randomise, e.g. for simple group mean
randomise -i allcopes.nii.gz -o permtest -1 -c 2.3 -m
/usr/local/fsl/data/standard/MNI152_T1_2mm_brain_mask.nii.gz
In BROCCOLI you can use the function RandomiseGroupLevel, which also
performs a permutation test (like randomise in FSL).
RandomiseGroupLevel allcopes.nii.gz -groupmean -mask
/usr/local/fsl/data/standard/MNI152_T1_2mm_brain_mask.nii.gz -cdt 2.5
Both randomise and RandomiseGroupLevel will generate new nifti files which
contain p-values for each voxel/cluster. Note that these p-values are
actually 1-p, so a value of 0.95 corresponds to p = 0.05 (corrected for
multiple comparisons).
2015-08-20 17:12 GMT+02:00 ImmanuelSamuel [email protected]:
That works well. How do you do higher level analysis after this. FSL asks
for registration before running higher level. Since the directory structure
is different it complains.—
Reply to this email directly or view it on GitHub
#13 (comment).
Anders Eklund, PhD
from broccoli.
See the following paper on permutation tests in neuroimaging, I prefer
non-parametric statistics as a permutation test as it is based on a lower
number of assumptions compared to parametric statistical methods
http://www.fil.ion.ucl.ac.uk/spm/doc/papers/NicholsHolmes.pdf
2015-08-20 17:30 GMT+02:00 Anders Eklund [email protected]:
I don't know if you can run the group analysis through the FSL GUI, the
registration has already been done in the FirstLevelAnalysis (all results
are in MNI space).You can use randomise in FSL to run a permutation test. First merge the
beta/cope image for each subject into a single nifti filefslmerge -t allcopes.nii.gz subject*_cope1.nii.gz
then run randomise, e.g. for simple group mean
randomise -i allcopes.nii.gz -o permtest -1 -c 2.3 -m
/usr/local/fsl/data/standard/MNI152_T1_2mm_brain_mask.nii.gz
In BROCCOLI you can use the function RandomiseGroupLevel, which also
performs a permutation test (like randomise in FSL).RandomiseGroupLevel allcopes.nii.gz -groupmean -mask
/usr/local/fsl/data/standard/MNI152_T1_2mm_brain_mask.nii.gz -cdt 2.5
Both randomise and RandomiseGroupLevel will generate new nifti files which
contain p-values for each voxel/cluster. Note that these p-values are
actually 1-p, so a value of 0.95 corresponds to p = 0.05 (corrected for
multiple comparisons).2015-08-20 17:12 GMT+02:00 ImmanuelSamuel [email protected]:
That works well. How do you do higher level analysis after this. FSL asks
for registration before running higher level. Since the directory structure
is different it complains.—
Reply to this email directly or view it on GitHub
#13 (comment).Anders Eklund, PhD
Anders Eklund, PhD
from broccoli.
Related Issues (20)
- BROCOLLI crashes without error
- TFCE is currently turned off error in RandomiseGroupLevel HOT 4
- permutations tmaps different from fsl randomize output by a factor of ~470 HOT 4
- GLM: segmentation fault HOT 2
- permutations turned off HOT 1
- Memory issue HOT 6
- Registration Failure HOT 1
- output uncorrected results for RandomiseGroupLevel HOT 2
- TFCE support? HOT 2
- AMD gpu and Docker failed HOT 4
- Roadmap? HOT 9
- Error building kernelBayesian.cpp for GPU on macOS HOT 2
- Elastic net regularization for GLM? HOT 2
- RandomiseGroupLevel for repeated measures HOT 2
- RandomiseGroupLevel crashing HOT 5
- python wrapper compilation HOT 1
- Segmentation fault - Multiple runs per participant still not fixed
- Compiling the matlab wrapper on Mac HOT 1
- nipype: BROCCOLI nodes do not produce any output within workflows HOT 1
- Broccoli in WSL2 - Ubuntu 22.04 - HOT 14
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from broccoli.