mfbox_hrf.m

toolbox for spm 5 for data, model free analysis

function hrf=mfbox_hrf(t)
%
% Copyright by Peter Gruber and Fabian J. Theis
% Signal Processing & Information Theory group
% Institute of Biophysics, University of Regensburg, Germany
% Homepage: http://research.fabian.theis.name
%           http://www-aglang.uni-regensburg.de
%
% This file is free software, subject to the 
% GNU GENERAL PUBLIC LICENSE, see gpl.txt

if (nargin<1), t = 1; end

hrf = [0.0000 0.0001 0.0007 0.0027 0.0076 0.0171 0.0334 0.0589 0.0960 0.1471 ...
0.2144 0.2997 0.4047 0.5304 0.6777 0.8467 1.0374 1.2492 1.4812 1.7321 2.0002 ...
2.2839 2.5809 2.8890 3.2060 3.5294 3.8568 4.1857 4.5138 4.8387 5.1582 5.4703 ...
5.7730 6.0645 6.3433 6.6078 6.8568 7.0892 7.3042 7.5010 7.6790 7.8380 7.9775 ...
8.0976 8.1982 8.2796 8.3421 8.3860 8.4118 8.4201 8.4115 8.3867 8.3465 8.2917 ...
8.2230 8.1413 8.0475 7.9424 7.8269 7.7018 7.5680 7.4264 7.2777 7.1228 6.9624 ...
6.7973 6.6282 6.4556 6.2804 6.1031 5.9243 5.7445 5.5642 5.3839 5.2041 5.0251 ...
4.8473 4.6710 4.4965 4.3242 4.1542 3.9869 3.8223 3.6606 3.5020 3.3466 3.1946 ...
3.0459 2.9008 2.7591 2.6210 2.4864 2.3555 2.2281 2.1043 1.9841 1.8674 1.7542 ...
1.6444 1.5381 1.4350 1.3353 1.2389 1.1456 1.0554 0.9683 0.8842 0.8030 0.7246 ...
0.6490 0.5762 0.5060 0.4384 0.3733 0.3107 0.2505 0.1926 0.1370 0.0836 0.0324 ...
-0.0167 -0.0637 -0.1088 -0.1519 -0.1930 -0.2324 -0.2699 -0.3056 -0.3397 -0.3721 ...
-0.4028 -0.4319 -0.4595 -0.4856 -0.5102 -0.5334 -0.5552 -0.5756 -0.5946 -0.6124 ...
-0.6289 -0.6442 -0.6584 -0.6713 -0.6831 -0.6939 -0.7035 -0.7122 -0.7198 -0.7265 ...
-0.7322 -0.7370 -0.7410 -0.7441 -0.7463 -0.7478 -0.7485 -0.7485 -0.7478 -0.7464 ...
-0.7444 -0.7418 -0.7385 -0.7347 -0.7303 -0.7255 -0.7201 -0.7143 -0.7080 -0.7014 ...
-0.6943 -0.6869 -0.6791 -0.6710 -0.6626 -0.6540 -0.6451 -0.6359 -0.6266 -0.6170 ...
-0.6073 -0.5974 -0.5874 -0.5772 -0.5670 -0.5567 -0.5462 -0.5358 -0.5253 -0.5147 ...
-0.5042 -0.4936 -0.4831 -0.4726 -0.4621 -0.4516 -0.4413 -0.4309 -0.4207 -0.4105 ...
-0.4004 -0.3904 -0.3805 -0.3708 -0.3611 -0.3516 -0.3422 -0.3329 -0.3237 -0.3147 ...
-0.3059 -0.2972 -0.2886 -0.2802 -0.2719 -0.2638 -0.2559 -0.2481 -0.2405 -0.2330 ...
-0.2257 -0.2185 -0.2115 -0.2047 -0.1981 -0.1916 -0.1852 -0.1790 -0.1730 -0.1671 ...
-0.1614 -0.1558 -0.1504 -0.1451 -0.1400 -0.1350 -0.1302 -0.1255 -0.1209 -0.1165 ...
-0.1122 -0.1080 -0.1040 -0.1000 -0.0963 -0.0926 -0.0890 -0.0856 -0.0823 -0.0791 ...
-0.0760 -0.0730 -0.0701 -0.0673 -0.0645 -0.0619 -0.0594 -0.0570 -0.0547 -0.0524 ...
-0.0502 -0.0481 -0.0461 -0.0442 -0.0423 -0.0405 -0.0388 -0.0371 -0.0355 -0.0340 ...
-0.0325 -0.0311 -0.0297 -0.0284 -0.0271 -0.0259 -0.0248 -0.0236 -0.0226 -0.0216 ...
-0.0206 -0.0196 -0.0187 -0.0179 -0.0171 -0.0163 -0.0155 -0.0148 -0.0141 -0.0134 ...
-0.0128 -0.0122 -0.0116 -0.0111 -0.0105 -0.0100 -0.0095 -0.0091 -0.0086 -0.0082 ...
-0.0078 -0.0074 -0.0071 -0.0067 -0.0064 -0.0061 -0.0058 -0.0055 -0.0052 -0.0049 ...
-0.0047 -0.0045 -0.0042 -0.0040 -0.0038 -0.0036 -0.0034 -0.0033 -0.0031 ...
-0.0029]; ...

hrf = interp1([0:319]/8,hrf,[0:t:29]);

return