Hi All,
I used a third party software to obtain t-stats map and now I wan to
calculate the FDR using the following code:
********************************************************
ttologp -logpout logp1 varcope1 cope1 `cat dof`
(this outputs the log(p) image)
fslmaths logp1 -exp p1
(this creates the p-value image)
fdr -i p1 -m ../mask -q 0.05
Probability Threshold is:
0.00339819
(this tells you the FDR threshold - p-values below 0.00339819 are significant)
fslmaths p1 -mul -1 -add 1 -thr .99660181 -mas ../mask thresh_1_minus_p1
(this creates a 1-p image for ease of display, thresholds that at
1-thresh and then remasks)
***************************************************************
However, I don't have the varcope1 and cope1 values since I didn't use
FEAT to generate the t stat maps.
Could you please advise what how can I work it around?
Best,
Leo
|