Dear Hendrik
As Donald mentions, when you know what you favorit hrf should look like you will need to modify some functions. There is a few ways to do this either by changing the way spm_get_bf calls spm_hrf, or this one which I prefer. I have added a field:
defaults.stats.fmri.hrfpar = [6 16 1 1 6 0 32];
in spm_deaults, although this file should not be edited! :-) But i hope the hack will make in the next release. Next I have changed the way spm_hrf finds its default values, namely using:
p = spm_get_defaults('stats.fmri.hrfpar');
instead of:
p = [6 16 1 1 6 0 32];
You can now change the default hrf to any other mixture of two gamma functions before you specify your fMRI design, e.g. an earlier peak with no poststimulus undershoot:
spm_get_defaults('stats.fmri.hrfpar', [4 16 1 1 inf 0 32])
You can either implement these changes, or replace your files with the ones attached.
Hope this helps
Best
Torben
Torben Ellegaard Lund
Associate Professor, PhD
Center of Functionally Integrative Neuroscience (CFIN)
Aarhus University
Aarhus University Hospital
Building 10G, 5th floor, room 31
Noerrebrogade 44
8000 Aarhus C
Denmark
Phone: +45 7846 4380
Fax: +45 7846 4400
http://www.cfin.au.dk <http://www.cfin.au.dk/>
[log in to unmask] <mailto:[log in to unmask]>
> Den Uge:23 05/06/2015 kl. 16.59 skrev MCLAREN, Donald <[log in to unmask]>:
>
> It's possible. You might need to modify the code. However, you'd need to know what the new basis function is ahead of time. Do you know what the "new" basis function is and what its formula and default parameters are?
>
> Best Regards, Donald McLaren
> =================
> D.G. McLaren, Ph.D.
> Research Fellow, Department of Neurology, Massachusetts General Hospital and
> Harvard Medical School
> Postdoctoral Research Fellow, GRECC, Bedford VA
> Website: http://www.martinos.org/~mclaren <http://www.martinos.org/~mclaren>
> Office: (773) 406-2464
> =====================
> This e-mail contains CONFIDENTIAL INFORMATION which may contain PROTECTED
> HEALTHCARE INFORMATION and may also be LEGALLY PRIVILEGED and which is
> intended only for the use of the individual or entity named above. If the
> reader of the e-mail is not the intended recipient or the employee or agent
> responsible for delivering it to the intended recipient, you are hereby
> notified that you are in possession of confidential and privileged
> information. Any unauthorized use, disclosure, copying or the taking of any
> action in reliance on the contents of this information is strictly
> prohibited and may be unlawful. If you have received this e-mail
> unintentionally, please immediately notify the sender via telephone at (773)
> 406-2464 or email.
>
> On Fri, Jun 5, 2015 at 7:58 AM, Hendrik Backhaus <[log in to unmask] <mailto:[log in to unmask]>> wrote:
> Dear all,
>
> I'm pretty new to SPM (and fMRI in general), so please excuse me for asking questions that may be stupid.
>
> I'm trying to find out how to create my own basis functions. I'm analyzing fMRI data of mice, and therefore I don't want to use the canonical HRF.
> Is it possible to add "new" basis functions to SPM12? If not, is it possible to manipulate one of the existing bf's so it fits my needs?
>
> Thank you!
> Best regards
> Hendrik
>
