 |
 |
Dear all FSL experts I have a statistical question about how to calculate corrected p value (voxel-wised and cluster-wised corrected p value) using FSL command line tool. The following thing is what I want to do and wish some of you could help me to figure out this problem. Rationale : I want to calculate the type I error rate of “one vs many statistical test (one patient vs many healthy controls)” of gray matter volume. I use Z-score to achieve this goal. Step 1: My dataset consist of 50 healthy subjects. I random select 26 subjects from this dataset and perform 100 times resample. For each resample, the one normal control were treated as patient and the other 25 subjects were served as normal controls. I calculate the Z score for every resample using the following formula [(patient GMV - healthy control mean GMV) / standard deviation of healthy controls GMV]. Step 2: After Z-score calculation, I need to set the significant level (p value) of this test to identify which brain area have significant difference between the patient and the normal controls (something like P < 0.05 = Z > 1.73). The area with significant statistical difference was related to the type I error rate. I think I need to consider the multiple comparison correction issue since the origin of this analysis is voxel-wised whole brain analysis , but I don’t know how to calculate the voxel-wised and cluster-wised FWE corrected P value using the FSL command line tool. After quick searching the FSL mail list, I think “fsl_glm” (which calculate the residual images) , “smoothest” (which calculate the data smoothness) may be needed for my demand, is it correct ? If yes, how to use these tools to calculate the voxel / cluster wised FWE corrected P value? Sorry for this long question, I try to describe it more clearly and hope to hear feedbacks from you Thanks Best Paul