Hi Anderson,A few other questions also came to my mind.1. Orthogonality means that there is no interaction between group and NP?
2. In the setup that you suggested, we are looking at positive interaction [0 0 0 1 0 ...] and negative interaction [0 0 0 -1 0 ...]. Is it possible that I get significance for both, in different areas though?
3. If I do find significance for NP interaction, then what's the next step?
4. For: "If not significant, you can investigate the effect of group and NP with additional contrasts then.". Then I treat NP like any other co-variate (of no interest) like age and sex. I guess I will just need to look at the contrasts:[0 1 0 0 ...] and [0 -1 0 0 ...] for group ( cont > pat and pat > cont respectively- considering EV2 as the group). )[0 0 1 0 ...] and [ 0 0 -1 0 ...] for NP (Considering EV3 as NP).
5. How do I handle time since injury? It is something which is 0 for all controls. I am assuming that I am demeaning the age, gender and NP test performance scores. Can I demean this time variable too?
Also what contrast do I construct for this?
Can I do an interaction term for this as well?
6. For age and sex, do I simply quote that these were included as covariates while performing t-test or do I still need to perform a significance test for these and state that there were no significant difference for these?
Will the contrast [0 0 0 0 1 ...] indicate that controls mean age was significantly higher than patients'? (Considering EV5 as the age).
This is my first time performing DTI analyses and I greatly appreciate your help.Eyesha
Date: Fri, 29 Apr 2016 10:35:21 +0100
From: [log in to unmask]
Subject: Re: [FSL] Statistical analyses- which mean values to use
To: [log in to unmask]Hi Eyesha,Please, see below:On 29 April 2016 at 03:23, Eyesha <[log in to unmask]> wrote:Dear experts,
I have a question about statistical analyses of DTI data.
I am comparing DTI parameters between two groups (controls and patients) and eventually I also want to look at the correlation between performance on neuropsychological tests and DTI parameters. This is my plan of action (after performing the basic DTI processing) after much thought:
A- Run the randomise tool to obtain t-statistics and p-values for voxel-wise comparison between the two groups, corrected for multiple comparisons,
B- Threshold the output at the desired p-value. I am choosing 0.05 and further dividing it by 4 to account for multiple analyses (4 DTI parameters: FA, MD, RD and L1).
C- Use the thresholded image to find tracks with significantly different values.
D- Get the mean value in each of the significant track (here I use only that part of the tract where significant differences were found, i.e. only the part of the tract that was included in the thresholded image and not the entire track).
E- Use the mean values found in D to investigate the correlation with NP test performance in the patient group.
F- Use regression to account for variables like age, sex and disease stage.
If the whole sample (patients and controls) were used for this follow up test, this would be a circular analysis, unless NP were perfectly orthogonal to group. Using only one of the groups, this isn't really an issue.However, it's still problematic in other grounds: suppose there is a correlation in that particular area found in A above in the patient group. How can one know if that is a noteworthy result if one hasn't looked into the control group? What if there is similar correlation between the DTI measures and NP in the control group in that particular region? Then there is nothing interesting to be reported as far as NP goes, and although the correlation may be significant in patients, it's still a false positive. Conversely, if there is a different correlation between NP and DTI measures according to group, then the difference between groups in terms of DTI varies according to NP, and the initial analysis of group differences is likewise misleading.To deal with this, perhaps be best thing is to drop these multiple models/stages and run a single one, that has:EV1: interceptEV2: group (+1/-1)EV3: NPEV4: NP*group (interaction)EV5, EV6, etc: age, sex, disease stage, etc.Then test for the interaction with contrasts [0 0 0 1 0 ...] and [0 0 0 -1 0 ...]. If significant, then it's something interesting. If not significant, you can investigate the effect of group and NP with additional contrasts then.The point where I got stuck was if I should use only the tracts with significant difference for further analyses or all the tracts. And if I should just use the mean value in the entire tract from the skeletonized parametric image, or un-skeletonized mean value in the entire tract or the skeletonized value from the significantly different region of the tract? I see people doing all of these things but I am not convinced for all.Here I would use the whole brain (voxelwise).Hope this helps.All the best,Anderson
Thank you for your expert opinion and help.
Eyesha