Doxygen Source Code Documentation
spm_bf.m
Go to the documentation of this file.00001 function [err, bf] = spm_bf (TR)
00002 %
00003 % [err,] = spm_bf ()
00004 %
00005 %Purpose:
00006 %
00007 %
00008 %
00009 %Input Parameters:
00010 %
00011 %
00012 %
00013 %Output Parameters:
00014 % err : 0 No Problem
00015 % : 1 Problems
00016 %
00017 %
00018 %
00019 %Key Terms:
00020 %
00021 %More Info :
00022 %
00023 %
00024 %
00025 %
00026 % Author : Gang Chen
00027 % Date : Mon Oct 27 17:13:26 EST 2003
00028 % SSCC/NIMH/ National Institutes of Health, Bethesda MD 20892
00029
00030
00031 %Define the function name for easy referencing
00032 FuncName = 'spm_bf.m';
00033
00034 %Debug Flag
00035 DBG = 1;
00036
00037 %initailize return variables
00038 err = 1;
00039
00040
00041 % p(1) - delay of response (relative to onset) 6 -gamd in waver -GAM
00042 % p(2) - delay of undershoot (relative to onset) 16
00043 % p(3) - dispersion of response 1
00044 % p(4) - dispersion of undershoot 1 X
00045 % p(5) - ratio of response to undershoot 6 X
00046 % p(6) - onset (seconds) 0 -gamd
00047 % p(7) - length of kernel (seconds) 32
00048
00049 p = [6 16 1 1 6 0 32]; %default parameters
00050 [bf1 p] = spm_hrf(TR,p); %get gamma hrf from SPM function spm_hrf
00051 dp = 1;
00052 p(6) = p(6) + dp; %for calculating time derivative
00053 bf2 = (bf1 - spm_hrf(TR,p))/dp;
00054 % D = (bf(:,1) - spm_hrf(TR,p))/dp; %time derivative
00055 bf = [bf1 bf2]; %combine the two
00056
00057 err = 0;
00058 return;
00059
