tag:blogger.com,1999:blog-1405968300258104460.post9155244314408951044..comments2024-03-27T01:18:24.341-04:00Comments on Andy's Brain Blog: FSL Tutorials 4-5: Group-Level AnalysisAndrew Jahnhttp://www.blogger.com/profile/16435706598096921650noreply@blogger.comBlogger35125tag:blogger.com,1999:blog-1405968300258104460.post-78605815526132818302020-03-05T12:04:06.388-05:002020-03-05T12:04:06.388-05:00Hi Kirk,
I recently switched over to the CONN too...Hi Kirk,<br /><br />I recently switched over to the CONN toolbox to do these types of analyses; see this writeup for more information: https://andysbrainbook.readthedocs.io/en/latest/FunctionalConnectivity/CONN_Overview.html<br /><br />See, for example, the "Adding Nuisance Covariates" section of the Group-Level Analysis chapter. It's still under construction, but it should help you with what you need. Unfortunately I don't have experience doing that particular analysis in either AFNI or FSL.<br /><br /><br />Best,<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-38606637574628917872020-03-03T03:07:36.788-05:002020-03-03T03:07:36.788-05:00Hi Andy!
Thanks for all your hard work in making ...Hi Andy!<br /><br />Thanks for all your hard work in making all these tutorials. It helps us a tonne and there are almost NO other resources out there which would be able to help us the way you do. <br /><br />Anyway I did have a question I'd be grateful if you could provide advice on. <br /><br />I want to see how resting state connectivity changes with scores on a cognitive test. Is there a way in FSL (or AFNI) to test for this? Obviously one way is to extract the correlations from a specific ROI for each subject and correlate that with their score. But I was wondering if there was a way to include it as a covariate and look at the effect of the covariate on the connectivity (as I've done in AFNI before for univariate analyses). <br /><br />Anyway if you have the time to provide advice on this I'd be super grateful. <br /><br />Hope all is well!<br />KirkAnonymoushttps://www.blogger.com/profile/04786220069195896956noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-25630010645690098322019-02-01T22:27:35.526-05:002019-02-01T22:27:35.526-05:00Hi Ahmet and Andy,
I could find a design matrix d...Hi Ahmet and Andy,<br /><br />I could find a design matrix description for repeated measurement study design in below FSL URL. <br />But, It seems that your guys were following other method..<br />Do you have any reason that?<br /><br />-Su<br /><br />https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/GLM#ANOVA:_2-groups.2C_2-levels_per_subject_.282-way_Mixed_Effect_ANOVA.29Suhttps://www.blogger.com/profile/16735239616316347891noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-70377769424778575782017-12-18T18:56:54.949-05:002017-12-18T18:56:54.949-05:00Hi Andy,
How would I go about implementing third ...Hi Andy,<br /><br />How would I go about implementing third level analysis? I am analysing data from a BART experiment where I am interested in effects related to inflation, last inflation before explosion, cashing out and balloon explosion. I have done first level analysis and then done all in second level analysis (my participant only had one scan each) whilst also regressing out the number of balloon pumps. <br /><br />I now want to do third level analysis and contrast different main effects. Which files should I be inputting in order to do this? I have 4 cope files however fsl tells me these are not registered. I can also input cope.nii.gz files within these files however I feel like this is just re-comparing my 24 participants within the same effect (as I did for level 2 analysis).<br /><br />Could you provide any guidance with my problem or just a run through of third level analysis in general?<br /><br />Many thanks in advance!Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-38798493066372586812017-07-28T11:57:33.928-04:002017-07-28T11:57:33.928-04:00Hi Jim,
I haven't done something like that be...Hi Jim,<br /><br />I haven't done something like that before, but you could try this:<br /><br />1. Open up FEAT and select "Higher-level analysis" from the dropdown menu at the top of the screen;<br /><br />2. In the Data tab, select "Inputs are 3D cope images from FEAT directories"<br /><br />3. Select your individual contrast images from SPM;<br /><br />4. In the Stats tab, selected "Mixed effects: FLAME 1"<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-2985534651721752042017-07-28T11:44:51.281-04:002017-07-28T11:44:51.281-04:00No, I don't. If you're interested in doing...No, I don't. If you're interested in doing FNIRT then the easiest way, as I understand it, is to simply check the "nonlinear" box in the Registration tab during FEAT.<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-3621223697993511362017-07-26T04:36:49.198-04:002017-07-26T04:36:49.198-04:00Hi Andy. Do you have tutorial on how to do FNIRT?Hi Andy. Do you have tutorial on how to do FNIRT?Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-84807570823020246722017-07-10T16:54:55.523-04:002017-07-10T16:54:55.523-04:00Hi Andy,
Thank you for your tutorials !
Because ...Hi Andy,<br /><br />Thank you for your tutorials !<br /><br />Because I am still new to this area, I face a problem about fMRI thresholding that all my friends can't answer. How can I use fsl FLAME1 to do ttest from the results of SPM5, which had already done with preprocessing and 1st level analysis in SPM5? I had looked over all the websites, but still can't find a good method to do so. Thank you if you are willing to answer my question.<br /><br />Thanks,<br />JimAnonymoushttps://www.blogger.com/profile/03929810778169567690noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-67808834790372370772017-04-02T16:12:37.412-04:002017-04-02T16:12:37.412-04:00That made more sense.
Looking at group activation...That made more sense.<br /><br />Looking at group activation results (gfeat) I am unsure which files to be looking at. Do you have any resources on this? Should I be looking at the rendered threth ztats?<br /><br />I have x,y,z coordinates for largest intensity voxel, but I can't seem to use '-' in the cursor tools part of the gui. Any ideas why this is?<br /><br />Cheers,<br />David<br />Davidhttps://www.blogger.com/profile/05783886466213580903noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-50110772376178808802017-03-31T12:13:56.259-04:002017-03-31T12:13:56.259-04:00Hey David,
You're looking at the z-maps, whic...Hey David,<br /><br />You're looking at the z-maps, which can be different from the contrast maps, depending on the variance for that contrast. For example, if there was lower variability in estimating one of the regressors, that regressor will have a higher z-score.<br /><br />Check out the cope images (or pe images, however you generated them for each individual regressor), and see if it holds up.<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-64261565217145887552017-03-31T01:02:18.415-04:002017-03-31T01:02:18.415-04:00Hi Andy,
To check, I chose the strongest activati...Hi Andy,<br /><br />To check, I chose the strongest activating voxel in the subtraction contrast, and the intensity values for A,B,A-B did not add up at all.<br /><br />I overlaid each thresh_zstat image over my example_func from the reg folder and checked voxel (42,49,1).<br /><br />example_func intensity: 692<br />A intensity: 0<br />B intensity: 0<br />A-B intensity: 5.74<br /><br />I have checked my ev's etc, and all seem fine. I can't think why this would occur given the relatively simple analysis. Any ideas?<br /><br />ThanksAnonymousnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-32255825636091915022017-03-24T10:00:49.054-04:002017-03-24T10:00:49.054-04:00Hey there,
It could be that B is negative, which ...Hey there,<br /><br />It could be that B is negative, which is why A-B is larger than A. For example,<br /><br />A = 0.5<br />B = -0.2<br />A - B = (0.5 - (-0.2)) = 0.7<br /><br />Check this out by selecting a voxel, writing down the value for map A and the value for map B, subtract B from A, and then check out the contrast map A-B.<br /><br /><br />Best,<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-44562768964590890042017-03-23T14:43:24.717-04:002017-03-23T14:43:24.717-04:00Having collapsed across runs (12 runs, one subject...Having collapsed across runs (12 runs, one subject), and with a simple block design, I have 4 contrasts (A,B,A-B,B-A). I am seeing a lot more activation for the subtraction contrasts (e.g. A-B) relative to the non-subtractd contrasts (e.g. A). Conteptually I am not sure how this can be. Surely the amount of activation for contrast A alone has to be more than the amount of activation for contrast A minus contrast B? The difference is sizable too. I am getting this result for other subjects too, am I doing something wrong? Or is there a basic explanation here?<br /><br />Thanks for the great resources Andy, very very helpful. I have learned a lot from your website and Youtube channel.Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-75560218702915441722016-08-24T20:46:50.829-04:002016-08-24T20:46:50.829-04:00Hi Andy,
I really can not thank you enough for th...Hi Andy,<br /><br />I really can not thank you enough for those wonderful tutorials on fMRI! Everyone was telling me that it will take ages to figure out fMRI and your videos have really made a difference and got me started! I could not have done it without you!<br /><br />Since you are an expert in the field I just wanted to ask you something... Basically I am using right now the Human Connectome Project task fMRI dataset and I wanted to try the method described here: http://www.nature.com/articles/srep20170#f2 but for now I felt like FSL was easier for me and I have no idea if something like that would be possible to run in FSL. I do not even understand how they ran this in SPM as well so I have nowhere to start to try it out. <br /><br />Ideally I am trying to form variability maps of the hand motor area for those HCP healthy subjects and even though I am not an expert in fMRI I thought this was a robust way to do that. 1) What do you think about that? 2) Do you think this is possible using FSL and could you guide me through if your time allows it? I would be happy to try that in SPM too or anything needed but I know nothing about SPM so far since it seemed more complicated for me to start with that so perhaps it would be harder to communicate the proper way of running the analysis to me. Any guidance would be appreciated. Thanks a lot!<br /><br />In any case thank you so much for all your help and for taking all the time to teach us how to perform neuroimage analysis but especially for making everything look so easy!<br /><br />Thanks!<br /><br />Best regards,<br />Magda Tsintou.Magda Tsintouhttp://www.linkedin.com/in/magda-tsintou-m-d-msc-232b9250noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-15758987458272527322016-08-02T15:37:45.394-04:002016-08-02T15:37:45.394-04:00Hi Andy,
Your tutorials have been of huge help as...Hi Andy,<br /><br />Your tutorials have been of huge help as this area is very new to me! <br /><br />I have a question, I have run dual regression (with unpaired t-test design) between (2) groups differences, and I have used group ICAs derived from the combined group ICAs of my sample (drug exposed, and non-drug exposed). I have got corrected 1-P values larger than 0.95, but when I look at the results in fslview it seems that all of the voxels that are significantly different between groups are not laying within the group networks (for example, if I am overlaying corrected IC 001 over IC-1, the significantly different voxels are not within the "highlighted" group network, but in another region of the brain that does not seem to be part of the functional network). What does this mean? Should I try and use just the non-exposed group ICAs and feed this to the dual regression?<br />Hope I am making sense. <br /><br />Many thanks,<br />NormaNormahttps://www.blogger.com/profile/03869062963098760902noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-56307438637927552692016-06-11T17:26:48.038-04:002016-06-11T17:26:48.038-04:00I forget if I did it at the individual run level (...I forget if I did it at the individual run level (I should check on that), but you should run featquery at the subject level, which is an average over all the runs within that subject. FSL has an odd convention where 1st level refers to individual runs and 2nd level refers to the subject, whereas the other major packages (AFNI, SPM) refer to subject results as 1st level and group results as 2nd level.<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-62374462363642427202016-06-11T00:05:36.437-04:002016-06-11T00:05:36.437-04:00Hi Andy,
Yes, that's what I'm referring t...Hi Andy,<br /><br />Yes, that's what I'm referring to. Thanks so much for your help! I got it to work that way and the report log seems to make sense, so all seems to be well.<br /><br />Speaking of ROIs, I've been reading around about Featquery (your blog, a tutorial by Jeanette Mumford, various forums), and it seems like plugging in first level feat directories into Featquery is a bad bad idea. Is this inherently an fsl issue, or is it generally a bad idea to look at ROI estimates at an individual subject, individual run level?<br /><br />Thanks again!<br />AhmetAhmetnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-57078353978057480012016-06-08T12:35:43.709-04:002016-06-08T12:35:43.709-04:00Hey Ahmet,
Are you referring to the parameter est...Hey Ahmet,<br /><br />Are you referring to the parameter estimates for the early runs and the late runs? In that case, you can do what you proposed, averaging the PE maps for early runs and late runs and then subtracting them from each other (the SPM developers refer to this as a summary statistic approach). If you assume that the variance between the groups is roughly equal, then you could do a contrast between the groups and follow up with ROI stats on those maps.<br /><br /><br />Best,<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-41492712917860556082016-06-06T00:21:29.381-04:002016-06-06T00:21:29.381-04:00Hi Andy,
Thanks for the wonderful tutorials. They...Hi Andy,<br /><br />Thanks for the wonderful tutorials. They've been of great help throughout the years.<br /><br />Being a very novice student of the fMRI school, I have what is probably a very simply group level analysis problem that has still managed to stump me. I have a patient and a control group with 4 subjects each, all undergoing 6 runs. What I'm interested in doing is contrasting late vs. early runs (last 2 runs vs. first 2 runs), and then of course looking at group differences. I understand that this is pretty much a repeated measures ANOVA, but I'm having trouble entering this kind of a model into FEAT. Do I concatenate the early runs and the late runs and then subtract using fslmaths once for the late vs early contrast, and one more time for the patient vs. control contrast? Or is there a way to enter this simple RM ANOVA model into FEAT?<br /><br />Thanks!<br />Ahmet Ahmetnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-21687175451211924072015-12-07T16:36:09.994-05:002015-12-07T16:36:09.994-05:00Hi there,
I don't use FSL that much, but I be...Hi there,<br /><br />I don't use FSL that much, but I believe that the typical FSL pipeline does the same series of steps: Motion correction (via MCFLIRT), followed by using FLIRT to register the functional images to the anatomical, and then warp the anatomical to a standard space and apply those warps to the functional images. FNIRT is the same idea, just carried out slightly differently, allowing for nonlinear warps and transforms.<br /><br />Best,<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-63261559256339405962015-12-06T21:53:52.600-05:002015-12-06T21:53:52.600-05:00Thank you so much for replying. That helps a lot. ...Thank you so much for replying. That helps a lot. I have another question. In SPM you realign, coregister and then normalize. From my limited understanding, in FSL FLIRT is like coregister where we do functional to anatomical registration and FNIRT should be used for aligning functional to standard space (using the anatomical). Is this a correct way for understanding ? Also, does this mean that motion correction should be done separately (or is it done through both of these) ? <br /><br />Thanks,Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-24482806791474325552015-12-06T14:20:44.102-05:002015-12-06T14:20:44.102-05:00Hi there,
Performing a group-level analysis with ...Hi there,<br /><br />Performing a group-level analysis with different numbers of runs is acceptable; you are still getting an average estimate of the parameters for that particular subject, with perhaps less accuracy, but the numbers are still valid.<br /><br /><br />Best,<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-75040782950707337332015-12-06T10:07:15.335-05:002015-12-06T10:07:15.335-05:00Hi Andrew,
I had a question regarding group analy...Hi Andrew,<br /><br />I had a question regarding group analysis. I collected three runs of data for my subjects but on analyzing the data, for one subject one run was incorrectly collected. My question is if I can include this subject during group level analysis where I only take two of his runs, while for others I take three each ? Would that be incorrect ?Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-69850374822685715792015-11-11T15:04:19.304-05:002015-11-11T15:04:19.304-05:00Hey Matthew,
Thanks for the feedback! It is true ...Hey Matthew,<br /><br />Thanks for the feedback! It is true that if you are analyzing one subject, then a fixed effects model is appropriate; you are covering the entire condition space that you are interested in, so there is no need to generalize to other conditions. However, if you recruit additional subjects, then these should be treated as a random factor with a differently modeled error term, since you are attempting to generalize their effects to all subjects.<br /><br />-AndyAndrew Jahnhttps://www.blogger.com/profile/16435706598096921650noreply@blogger.comtag:blogger.com,1999:blog-1405968300258104460.post-63576403132439745012015-11-04T15:26:34.315-05:002015-11-04T15:26:34.315-05:00Thank you for your incredibly useful and straightf...Thank you for your incredibly useful and straightforward tutorials! They're very helpful.<br /><br />Quick comment. This will depend on your purpose, but for repeated scans within the same individual (i.e., multiple runs of a task), it might be worth considering the Fixed Effects model. I don't know whether you get to that later, but it does seem potentially more powerful. Across subjects, of course, you'd want to use mixed effects.Anonymoushttps://www.blogger.com/profile/12970752292344324303noreply@blogger.com