在MATLAB中對hrf時間序列進行建模

我試圖對fMRI數據建模，以便檢查實驗設計的功效。我一直在跟隨幾個教程，並有一個問題。我首先需要通過將刺激輸入時間序列與典型血液動力學響應函數（HRF）進行卷積來建模BOLD響應。我查了第一個教程說，一個可以作出HRF是任何幅度，只要在HRF的「形狀」是正確的，使他們創造了以下HRF在MATLAB：在MATLAB中對hrf時間序列進行建模

hrf = [ 0 0 1 5 8 9.2 9 7 4 2 0 -1 -1 -0.8 -0.7 -0.5 -0.3 -0.1 0 ]

然後卷積的HRF與刺激只要使用「CONV」這樣：

hrf_convolved_with_stim_time_series = conv(input,hrf);

這是非常簡單的，但我希望我的模型，最終得到儘可能準確，所以我查了更高級教程，他們做了以下。首先他們創建了20個時間點的矢量，然後使用'gampdf'函數創建HRF。

t = 1:1:20; % MEASUREMENTS 
h = gampdf(t,6) + -.5*gampdf(t,10); % HRF MODEL 
h = h/max(h); % SCALE HRF TO HAVE MAX AMPLITUDE OF 1

這樣做比較簡單嗎？我想我有3個具體問題。

在「gampdf」幫助頁面是超短並且僅說，「6」和「10」中的每個函數調用所代表「A」，這是一個「形狀」參數。這是什麼意思？它沒有提供其他信息。爲什麼第一個電話是6，第二個是10？
這個問題直接關係到上面的問題。該代碼是針對TR = 1且刺激非常短（如1s）的情況編寫的。在我的情況下，我的TR = 2，我的刺激很長（12s）。我試着去適應上面的代碼做，使工作HRF我的情況如下：
```
t = 1:2:40; % 2s timestep with the 40 to try to equate total time to above 
h = gampdf(t,6) + -.5*gampdf(t,10); % HRF MODEL 
h = h/max(h); % SCALE HRF TO HAVE MAX AMPLITUDE OF 1 
```

因爲我不知道什麼「gampdf」參數的含義（或行所做的，在所有的現實）我不知道這給了我什麼，我正在尋找。我基本上得出20個值，其中1-14具有一些數值，但是15-20全是0.我假設在整個12s刺激期（前6個TR值如1-6值）中將會有響應適當的糾正可能是其餘的價值，但我不確定。

最後一個問題。另一個代碼並不會'規模化'HRF的幅度爲1.最終會這樣嗎？

來源

2017-10-18 chainhomelow