血流動態反応関数(HRF)をPythonでモデリング

fMRI Data Analysis
2018.11.16




血流動態反応関数とは

血流動態反応関数は、HRFと略されます。

編集中

Pythonでモデル化

Pythonで血流動態反応関数をグラフ化・数値化していきます。

PythonのNipyというライブラリを使います。

必要なライブラリ

  1. matplotlib
  2. numpy
  3. nipy

プログラム

01
02
03
04
05
06
07
08
09
10
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
import nipy.modalities.fmri.hrf as hrf
 
t = np.arange(0, 30, 1)
 
ff = hrf.spm_hrf_compat(t)
plt.plot(ff)
plt.savefig('HRF.png')

HRFのグラフ

参考

Course on Practical Neuroimaging in Python — Practical neuroimaging analysis
practical-neuroimaging.github.io
Neuroimaging in Python — NIPY Documentation
nipy.org
fMRI
タイトルとURLをコピーしました