Print

---

In spm_get_bf.m, you need to modify this section of code to have a case for
MION (e.g. case {'MION_hrf'}:

% get order and length parameters

%--------------------------------------------------------------------------

switch xBF.name



    case {  'Fourier set','Fourier set (Hanning)',...

        'Gamma functions','Finite Impulse Response'}

    %----------------------------------------------------------------------

    try,    l          = xBF.length;

    catch,  l          = spm_input('window length {secs}',3,'e',32);

        xBF.length = l;

    end

    try,    h          = xBF.order;

    catch,  h          = spm_input('order',4,'e',4);

        xBF.order  = h;

    end

end







% create basis functions

%--------------------------------------------------------------------------

switch xBF.name



    case {'Fourier set','Fourier set (Hanning)'}

    %----------------------------------------------------------------------

    pst   = [0:dt:l]';

    pst   = pst/max(pst);



    % hanning window

    %----------------------------------------------------------------------

    if strcmp(xBF.name,'Fourier set (Hanning)')

        g  = (1 - cos(2*pi*pst))/2;

    else

        g  = ones(size(pst));

    end



    % zeroth and higher Fourier terms

    %----------------------------------------------------------------------

    bf    = g;

    for i = 1:h

        bf = [bf g.*sin(i*2*pi*pst)];

        bf = [bf g.*cos(i*2*pi*pst)];

    end



    case {'Gamma functions'}

    %----------------------------------------------------------------------

    pst   = [0:dt:l]';

    bf    = spm_gamma_bf(pst,h);



    case {'Finite Impulse Response'}

    %----------------------------------------------------------------------

    bin   = l/h;

    bf    = kron(eye(h),ones(round(bin/dt),1));



    case {'NONE'}

    %----------------------------------------------------------------------

    bf = 1;


Then, you need to set xBF.name='MION_hrf' in you code prior to the call to
spm_get_bf.m. Make sure that the bf is a column vector or matrix.
Best Regards, Donald McLaren
=================
D.G. McLaren, Ph.D.
Postdoctoral Research Fellow, GRECC, Bedford VA
Research Fellow, Department of Neurology, Massachusetts General Hospital
and
Harvard Medical School
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 Sun, Dec 4, 2011 at 5:31 PM, Justin Vincent
<[log in to unmask]>wrote:

> How do I use my own hemodynamic impulse response function with SPM8?
>
> I have monkey fMRI data collected with MION contrast agent.  The impulse
> response function for MION is different than BOLD.  I need to use my own
> impulse response function (rather than the canonical HRF) to properly
> model the data.  I didn't see an option for using your own impulse
> response function in the GUI or when I saved the output batch and script
> files.  How do I use my own MION impulse response function to model my
> data? Thanks for any help.
>
> --
> Justin Vincent
>
>
> 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.
>