Hi, I’m working on some fMRI data processing these days. In order to construct a HRF, I find this Matlab function in fsfast. But after some careful comparison, it seems that the gamma pdf in this function is not consistent with FSL’s version, and differs from the standard gamma pdf. So perhaps it is a typo in your code? In this function, the gamma pdf is written as:
pdfx = (b.^2) .* (…);
But comparing to FSL’s feat_model.cc, I think this term should be
pdfx = (b.^a) .*(…);
Which will be consistent with the original gamma pdf from Wikipedia. Thanks for your nice software and thanks for your attention. 😊 Su, ---
From Beijing Normal University.
No, it is not. Have you tried using -fslhrf option in mkanalysis-sess?
On 01/09/2016 12:13 AM, LiTuX Sol wrote:
Hi,
I’m working on some fMRI data processing these days. In order to construct a HRF, I find this Matlab function in fsfast. But after some careful comparison, it seems that the gamma pdf in this function is not consistent with FSL’s version, and differs from the standard gamma pdf. So perhaps it is a typo in your code?
In this function, the gamma pdf is written as:
pdfx = (b.^2) .* (…);
But comparing to FSL’s feat_model.cc, I think this term should be
pdfx = (b.^a) .*(…);
Which will be consistent with the original gamma pdf from Wikipedia.
Thanks for your nice software and thanks for your attention.😊
Su,
From Beijing Normal University.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu
-
Douglas N Greve -
LiTuX Sol