[Freesurfer] About the pdf_gamma function in fsfast's toolbox

2016-01-08 Thread LiTuX Sol
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


The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.


Re: [Freesurfer] About the pdf_gamma function in fsfast's toolbox

2016-01-11 Thread Douglas N Greve
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

-- 
Douglas N. Greve, Ph.D.
MGH-NMR Center
gr...@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422

Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.