berkeley-stat159 / project-eta

BSD 3-Clause "New" or "Revised" License
3 stars 4 forks source link

New brain images! #85

Open jonmigueljara opened 8 years ago

jonmigueljara commented 8 years ago

Hey Matt. We updated our design matrix so that we are convolving the gain and loss regressors as well like you suggested.

screen shot 2015-12-12 at 3 47 04 am

screen shot 2015-12-12 at 3 45 04 am

Can you help us interpeting these images? We notice that there are some small clusters of activiation, but what else can we do to analyze and interpet these images?

matthew-brett commented 8 years ago

Are these images from analysis of the smoothed data? Preprocessed data? What else is in the design matrix? Could you also post an image of the design matrix?

jonmigueljara commented 8 years ago

yeah this is the data we smoothed but its not the preprocessed data that we got towards the end of the class. The design matrix has the convolved gain,loss, and onset values.

These are the values. design_plot_gain_loss

design_plot_task_onset

matthew-brett commented 8 years ago

What other regressors did you include?

Please post an 'imshow' of your design matrix.

What smoothing did you use?

jonmigueljara commented 8 years ago

the 3 regressors are. the general task onset times, gain, and loss. All 3 convolved.

we smoothed using a Gaussian Filter by 2 SD's in all 3 spatial dimensions.

oh right, here it is:

screen shot 2015-12-12 at 1 40 42 pm

matthew-brett commented 8 years ago

Is the first column a column of 1?

Did you try including the linear drift term? The squared linear drift term? The PCA components that looked artefactual? The reason I ask is because your image looks noisy and seems to be picking up the edge of the brain, so might be picking up some movement effects.

I also recommend you subtract the mean from your convolved regressors. This makes the betas easier to interpret, because you can interpret the 'gain' and 'loss' betas as the increase / decrease in activity having allowed for activity due to the task event (your 'onset' regressor).

jonmigueljara commented 8 years ago

yes the first column is 1's for the intercept. We were not including the drift terms so I attempted to use them but I found that it still looks "edgey". Is this still possibly still due to noise? (I still need to subtract the mean from the regressors.)

image

matthew-brett commented 8 years ago

I agree - it does look very edgey. I would try subtracting the mean from the regressors, and I would try also looking at the PCA of the data to see if you can find edgey components to remove by regression.