Print

Print


Hi Pablo,

Although it is certainly possible to use subtractions to later study group and timepoint differences, as well as interactions, it becomes easily complicated for 3 timepoints (and 3 groups). I think for these cases, it is more straightforward to use the full design matrix, and deal with the repeated measurements using exchangeability blocks.

Please, see at the link below a mock-up for designs and contrasts (note that there are two sets) to test both within- and between-subject effects.
https://dl.dropboxusercontent.com/u/2785709/outbox/mailinglist/design_pablo.ods

For design1, the hypotheses are for within-subject effects of time (C1 and C2), and interactions group vs. time (C3-C6). To run these in randomise, specify a file with the exchangeability blocks (EBs, with the option "-e"), so that permutations will happen within subject. This will assume compound symmetry, that is, the correlations between t1 and t2, t2 and t3, and t1 and t3 are all the same.

For design2, the hypotheses are for between-subject effects, i.e., group differences. Use the same EB file, but in randomise, include the option "--permuteBlocks", so that permutations will happen of these blocks as a whole.

The above is using randomise. If you'd like to correct for all these contrasts, it's possible to use PALM: you'd enter both designs (with "-d"), both contrast files ("-t"), both EB files ("-eb"), the options "-within" and "-whole", so that permutations will happen within-block and also whole-block, and "-corrcon", to correct across all contrasts (in these cases, not all permutations serve to test both designs simultaneously, but this is fine).

Hope this helps.

All the best,

Anderson






--
Anderson M. Winkler
FMRIB / Analysis Group
Blog | Twitter ]


On 3 September 2015 at 16:29, Pablo Ripolles <[log in to unmask]> wrote:
Hi everyone,
We have a study with 3 patient groups (who underwent different clinical interventions) and 3 time-points and we would like to test several effects using TBSS. I have been reading the mailing list and the GLM guide, but I still have some doubts about how to create the design and contrasts with Glm in order to use them later with randomise.

First, after TBSS processing, I have created a set of FA images subtracting one time point minus each other for each subject: T3-T2, T3-T1, T2-T1. Then to calculate the interaction between any pair of times and Group, I have created the following design and contrasts (assuming I only have 2 subjects per group and that we will later use the T3-T1 FA images):

Design
Group A   Group B    Group C
1                 0             0
1                 0             0
0                 1             0
0                 1             0
0                 0             1
0                 0             1

Contrast
1    -1   0 --->This would calculate: Group A (T3-T1) > Group B (T3-T1)
1     0  -1 --->This would calculate: Group A (T3-T1) > Group C (T3-T1)
0    1   -1 --->This would calculate: Group B (T3-T1) > Group B (T3-T1)

Is this correct so far?

Reading the GLM guide (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/GLM ; in 1 level 4 factors) I think that if I select contrasts number 2 and 3 to do an F-test, I will get a contrast that will test for any significant effect in any of my interactions, i.e., an overall F-test. Is this right?

Then if I want to check for the effect of time in each group, i would specify these contrasts:

1      0       0----> Significant areas in which FA has increased in T3 from T1 in Group A
0      1       0----> Significant areas in which FA has increased in T3 from T1 in Group B
0      0       1----> Significant areas in which FA has increased in T3 from T1 in Group C

And if I want to see the overall effect of time (for T3-T1) in all groups then I would use these three contrasts in an F-Test. Can I do this?

Does this make any sense?
Many thanks in advance.