Dear Karl, Thanks very much for the clarification. I'll copy to the list for others following this thread. Best, Andy At 12:02 PM 12/3/2001 +0000, you wrote: >Dear Andy, > >Oh - I see. Your approach seems absolutely fine. You are simply >reporting the parameter estimates about which you made an inference (in >terms of a contrast of these estimates). Just say that they are >relative to the mean over all conditions and report them directly for >each condotion. > >With very best wishes, > >Karl > > >----- Begin Included Message ----- > >Date: Mon, 26 Nov 2001 12:47:50 -0500 >To: [log in to unmask] (Karl Friston) >From: "Andrew J. Saykin" <[log in to unmask]> >Subject: obtaining condition means from beta and contrast image values > >Dear Karl, > I should have provided a bit more background on the rationale for >trying to extract these means. What we are really trying to do is >visualize the regional activation results in hypothesized areas for a 4 >level parametric blocked design task. We have a patient and control group >that we are comparing. I have analyzed the parametric responses as a >contrast ( -6 0 1 0 2 0 3 0 ) for each individual (4 conditions >plus temporal derivatives included in model). I then used a 2nd level >random effects model to compare the above contrast (coding for parametric >increase) between groups. The motivation for my questions about computing >means is that the random effects results don't lend themselves to an >intuitive graph of the parametric activation effects or group >differences. Since we didn't want to enter 4 conditions per subject in a >random effects model because it would violate the 1 scan per subject >restriction, I thought we could get fitted and adjusted means for the 4 >conditions using the betas (1,3,5,7) to examine the parametric >responses. Is there a different way that you would suggest approaching >this analysis and visualization? > >Many thanks, >Andy > > >At 04:36 PM 11/26/2001 +0000, you wrote: > >Dear Andy, > > > > > I know there has been previous correspondence on this topic with > > > regard to calculating percent change and I wanted to check out my > > > understanding. I am interested in obtaining condition means for fmri > > > data on a voxel by voxel basis separately for each of a group of > > > subjects. Each subject was processed separately. We performed a blocked > > > analysis where there were 4 conditions (and also included the 4 temporal > > > derivatives to adjust for the possibility of minor timing > > > errors). Since SPM automatically includes the intercept there are then > > > 9 columns in the design matrix for extracting contrasts > > > (cond1, td1, cond2, td2, cond3, td3, cond4, td4, constant). > > > I think the answers to the following questions are "yes" but would > > > appreciate confirmation of this or guidance if I am missing something. > > > > > >(1) Does the following contrast test for a condition 1 effect versus zero > > >(null hypothesis, no slope)? > > > > > >cond1 td1 cond2 td2 cond3 td3 cond4 td4 constant > > > 1 0 0 0 0 0 0 0 > 0 > > > >Yes it does. > > > > >(2) Do the obtained contrast image values indicate the condition 1 mean, > > >but with the global scaling removed since the intercept was included in > > >the model? > > > >The contrast returns the estimated effects of cond1, discounting any > >effects that can be explained by the other regressors. Because the > >other regressors include the constant, the ensuing contrast is > >'activation' from the mean of this voxel's time-series. Global effects > >are removed before estimation by proportional scaling to 100, in fMRI. > >This means the [adimensional] units are percent of whole brain mean. > > > > >(3) Is the sum of beta1 (ie, cond1) plus beta9 (ie, constant) the actual > > >mean for condition 1? > > > >Not really. The 'mean of condition 1' does not have any meaning. > >Beta9 is the mean over all conditions. Beta1 + Beta9 would be an > >estimate of the fitted response during condition 1 but only because the > >boxcar for cond1 has unit height. The mean over conditions is Beta9. > >Note that if you are comparing Beta1 from different subjects it is not > >necessarily a good idea to divide by Beta9 (implicit in % regional > >activation). This is because you would be confounding any differences > >in cond1 with differences in Beta9 (which could have artefactual or > >technical causes). > > > >I hope this helps - Karl > >----- End Included Message -----