【初心者向け】基礎&実践プログラミング

初心者がつまづきやすいところ、最短で実力が身につく方法をお伝えします。

【Python】カッパ(k)係数の計算

目的 コーエンの重み付きカッパ(k)係数(Cohen's kappa)をPythonで計算 scikit-learnのインストール カッパ係数の計算には、scikit-learnを使うと簡単にできます。scikit-learnをインストールしていない場合は、以下のコマンドを実行してインストールして…

【FSL】ROI(VOI)を用いた画像計測

目的 関心領域ROI(VOI)内の画像の定量値を計測 計測コマンド FSLのfslstatsを使う。 各オプションは、以下の通り。 -K: 関心領域ROI -M: ROI内の平均値の算出 fslstats -K <ROI> <input> -M 複数人被験者がいる場合 複数人の被験者がいる場合、以下のようなフォルダ構</roi>…

【FSL】Affine変換(平行移動+線形変換)

目的 異なる画像(NIFITI形式)の位置合わせをAffine変換(平行移動+線形変換)を用いて実行する。 FSLを用いて実装します。 必要なファイル T1WI: t1.nii.gz QSMの元画像: qsm_raw.nii.gz QSMの定量マップ: qsm_map.nii.gz 大脳基底核のROI: basalganglia_…

【Linux】ファイルやフォルダを小分けにする方法

目的 並列でプログラムを実行したい場合、データセットを小分けにしたい場合がある。小分けにする数が多いほど手間がかかる。 そこで、以下を目的とする。 シェルスクリプトを使ってデータセット(ファイルやフォルダ)を小分けにする フォルダ構造 フォルダ…

【Linux】DICOMのソート処理

目的 様々なSeriesDescriptionで撮像されたDICOMデータから特定のSeriesDescriptionのDICOMデータのみを取り出す。 シェルスクリプトで実行 dcmdump dcmdumpを使うことでDICOMのヘッダー情報を取得することができる。 dcmdumpのインストール $ apt-get insta…

【FSL】大脳基底核(Basal Ganglia)のセグメントと体積(Volume)計測

目的 FSLを用いて大脳基底核(Basal Ganglia)のセグメントと体積(Volume)計測 フォルダ構造 dataフォルダを作成し、その中に各被験者のフォルダ用意。大脳基底核のセグメントに必要な3D-T1WIは各被験者フォルダに<フォルダ名>_vol.nii(.gzも可)というファイル…

【Python】画像の類似度の計算(Dice係数)

目的 NIFTI形式の画像からPythonを使ってDice係数を計算する。 準備 必要となるパッケージは以下の通りです。必要に応じてインストールしてください。 pip3 install nibabel pip3 install numpy pip3 install glob pip3 install pandas フォルダ構造 フォル…

【SPM】画像の解像度の変更(Coregistration)

目的 画像の位置は同じで解像度が異なる2枚の画像の解像度を合わせる SPMのCoregisterを用いる。 方法 1. MATLABからSPMを開く。 spm 2. 上段(spatial pre-processing)からconergister(EST/RES)を選択。 3. Reference Imageに最終的に合わせたい解像度…

線形近似を用いた交絡因子の排除

はじめに データを扱う上で、年齢や測定器、環境などの交絡因子が結果にどのような影響を及ぼすのか。 それらの交絡因子を取り除くにはどうすべきかを、線形近似を用いて解決する。 間違えていたら、ご指摘いただければ幸いです。 従来のz-score 本来のz-sco…

【MATLAB】ComBatを用いた交絡因子の排除 (実践編)

問題 グループ間の差を解析する際、計測した値は、年齢、性別、計測器等の交絡因子の影響を含む。 目的 年齢、性別、計測器等の交絡因子の影響を排除するように計測値を補正する。 詳細は、こちらのManualをみること。 準備 ComBatのDownload データの補正に…

【Python】相関行列 (Correlation Matrix)とヒートマップ (Heat Map)の作り方

目的 Pythonで相関行列 (Correlation Matrix)を作成。 ヒートマップ (Heat Map)で図示。 Better Heatmaps and Correlation Matrix Plots in Pythonを参考にしました。 シンプルなヒートマップ まずは、基本的な相関行列のヒートマップを作成します。 ライブ…

【Python】縦断データの可視化(プロット・箱ひげ図・バイオリン図・信頼区間・ヒストグラム)

目的 Pythonを用いて、縦断データ(beforeとafter)の可視化 主にやることは、 プロット (dot plot) 箱ひげ図 (box plot) バイオリン図 (violin plot) 信頼区間 (CI bar) ヒストグラム (histogram) 最終的には、このような図を作成する。 準備 open-visualizat…

【PyTorch】サンプル⑨ 〜 動的グラフ 〜

目的 PyTorchの特徴の一つである動的グラフに挑戦する。 前準備 PyTorchのインストールはこちらから。 初めて、Google Colaboratoryを使いたい方は、こちらをご覧ください。 コマンドラインの「>>>」の行がPythonで実行するコマンドです。 それ以外の行は、…

Pythonを使った株価の自動収集

はじめに 最近、株価を予想するAIを作って欲しいということで、全くいい精度がでる期待もなかったが作ることにした。 過去にkerasやtensorflowを使ってAIを作ったことがあったが、株価を集めることはしたことがなかったのでやり方をまとめておく。 ビックデ…

【PyTorch】サンプル⑧ 〜 複雑なモデルの構築方法 〜

目的 このチュートリアルに至るまでは、ニューラルネットワークモデルの定義を積み木を積み重ねるように単純なシーケンスtorch.nn.Sequentialで構築していました。 このtorch.nn.Sequentialを用いた方法は、モデルの定義が簡単である反面、ネットワーク構造…