CONNチュートリアル 13: first-level Analysis: ROI-to-ROI & Seed-to-Voxel

(2017.11.26: 追記しました)

Denoisingが終わった後、first-levelの画面に来ます。
first-levelということは、

個人のfunctional connectivityを見る

ということになります。

最初は、下のようなものしか見えません。

ROI-to-ROIとSeed-to-Voxelがひとつ
Voxel-to-VoxelとICA networksが次にきて、
最後にdynamic ICAがきます。

まずは、ROI-to-ROI, Seed-to-Voxelを見ていきましょう。

ROI-to-ROIをクリックします。

ここでは画面の右側でいろいろなことが確認できますので、右側のスクリーンショットで説明します。

上にある “Preview first-level analysis results”
の下に領域を選ぶことができます。

この領域は “Seed” になる領域です。

まず、わかりやすく、デフォルトモードネットワーク(DMN)を確認してみましょう。

atlas.FP r のところをクリックすると領域が多く出てきますので、スクロールして、”networks.DefaultMode.PCC (1,-61,38)”を選択してください。

画像が見えているところの右側にマウスをもっていくとスクロールバーがあらわれますので、適宜スクロールしていただくと、確かにDMNっぽい感じで見えます。

ただ、これは色がつきすぎの感じがしますね。

画面下を見ると threshold 0.25 と見えます。
これは、今の場合Default Mode NetworkのPCCと相関係数0.25以上の関係にある領域を描出しているということになります。

この値を変えてみることで変わります。

そして、一番右側にあるSubjectsを適宜クリックしていくと、各個人のDMNの様子を観察することができます。

続いて、どんな解析をしていくのかを決定します。

画面中央をご覧ください。

ここに

Analysis type
Analysis options

とあり、それぞれ2項目設定する必要があります。

  • Analysis type
  • ここで、何の解析をするかを決定します。デフォルトは、”functional connectivity (Weighted GLM)” となっています。このままで大丈夫です。もし、Psychophysiological interaction (PPI)などをしたい場合は、変更します。

    その下の段のデフォルトは、”ROI-to-ROI analysis only” となっています。ここをクリックすると、”ROI-to-ROI” と “Seed-to-Voxel”を選ぶことができるのがわかります。いろいろな結果を見てみたいので、”ROI-to-ROI and “Seed-to-Voxel”をここでは選んでいます。

  • Analysis options
  • ここでは解析のオプションを指定します。上段の選択肢は4つ、”correlation (bivariate)”, “correlation (semipartial)”, “regression (bivariate)”, “regression (multivariate)” となります。 2変量の相関係数、半偏相関係数、2変量の回帰分析のβ, 多変量回帰分析のβということですね。どんな解析をするかによって異なります。今回は、シンプルにデフォルトの2変量の相関係数とします。

    その下の段のデフォルトは、デフォルトの no weighting で大丈夫です。タスクfMRIの時にここを変える必要が出てくると私は理解しています。

次に、画面左側の Voxel-to-Voxel をクリックします。

そうすると、画面中央が下のような図になるはずです。

  • Analysis type
  • ここでは、主成分分析、独立成分分析、MVPA、ALFFの計算などができます。
    私は、独立成分分析 (ICA) を行ってみたいので、Analysis typeの下をクリックして、”Group ICA” を選択しました。

  • Number of factors
  • これは、独立成分分析の際のFactor数です。デフォルトの20のままにしました。

  • Dimensionality Reduction
  • これは、独立成分分析を行う際に、connectivity matrixの次元を落とすというものです。デフォルトの64にしました。

その下に dyn-ICA circuits (dynamic ICA circuits) がありますが、ここは正直私の手にはおえないので、何も変更しませんでした。

このうえで、左下のDoneをクリックします。

そうすると、first-level Analysisの解析が行われます。

Print Friendly, PDF & Email

36 thoughts on “CONNチュートリアル 13: first-level Analysis: ROI-to-ROI & Seed-to-Voxel

    • いえ、kiyotaka_at_nemotos.netにお願いします。_at_を@に変えてください。

  1. 根本先生
    ご指導ありがとうございます
    1.Between-subject contrast は設定していません。1つのグループでconndition間の違いを検討したいと思いました。ICAは2グループの違いの検討のほうが適していますでしょうか。
    2.ありがとうございます。自分の得たICAの結果と似通っているものはどれになるかというのを知りたいときに使うものだとはわかりませんでした。rois/networks.niiを使って調べてみます。
    お忙しいところ、誠にありがとうございました。今後ともよろしくお願い致します。

    • 失礼しました。1つの群で3つのconditionの違いを確認したいということですね。

      どのように設定したかを示していただけますか?たぶん、コントラストの設定がうまくいっていないと思います。

      • 根本先生
        ご指導ありがとうございます

        メールで資料を添付しましたがうまくいかないようです

        Groupは1つで、conditionが5つです
        Connの最初の設定で5つのconditionを設定しています
        5つのconditionのうちの2つづつを比較しております。
        Close > Open、Open >Restです

        考えてみたら、各conditionのICAが違っていたら、比較できないので、5つのconditionのICAが同じになっているのは当然でしょうか。
        だとすると、必要がないconditionがあるときは、解析の時点で外しておかないとだめでしょうか。

        お忙しいところ誠に申し訳ありませんが、ご指導のほどよろしくお願いいたします

  2. 根本先生、ありがとうございます。大丈夫です。お待ちしてます。お忙しいところ申し訳ありません。どうかよろしくお願いいたします。

  3. いつもお忙しいところ、ご指導誠にありがとうございます。
    日本大学松戸歯学部の石井智浩です

    CONNを用いた独立成分分析についてさらに質問がありメールいたしました。
    基本的なことで申し訳ありませんが、ご指導いただければ幸いです。

    1.1つのグループで、conditionを3つ(A,B,C)にして、独立成分分析を行いました。
    Spatial componentsのなかで、Analysis of group-ICA networks during A、おなじく、Analysis of group-ICA networks during B、Analysis of group-ICA networks during Cとして、また、Condition間の差(A>B) として、独立成分分析を行い。
    SummaryでICA toolsからCompute spatial match to templateを実行すると、Correlation matrixが得られて、それぞれのICAとnetworksやregionとの相関がでるのですが、Condition間、Condition間の差で全く同じCorrelation matrixが得られます。
    これは、解析が間違っているでしょうか。

    2.Compute spatial match to templateを実行する際にSelect template file がでてきますが、ここでは、どのfileを選べばよいでしょうか。

    以上、お忙しいところ誠に申し訳ありませんが、ご指導いただけると幸いです。

    • すみません、今、仕事がものすごく立て込んでいるので、もう数日返信に時間をください。

    • 石井先生

      返信が遅くなりすみません。

      私の環境で確認してみました。

      1. 解析方法について
      Group ICAの画面で、左下の方にある
      Between-subject contrast をどのように設定されましたか?

      そこを設定されていないのではないかなというのがちょっと気になりました。

      2. templateですが、ここは、自分の得たICAの結果と似通っているものはどれになるかというのを知りたいときに使うものです。
       CONNに入っているものであれば、rois/networks.nii を私だったら使います。

      ご参考まで。

  4. 根本先生
    お忙しいところ誠にありがとうございました
    ICAの理解が深まりました
    いつもご丁寧なご指導まことに感謝いたします
    今後ともどうかよろしくお願いいたします

  5. お忙しいところ、ご返事いただきありがとうございます。
    お待ちします。
    誠に申し訳ありませんが何卒よろしくお願いいたします。

  6. いつもご指導誠にありがとうございます。
    日本大学松戸歯学部の石井智浩です

    CONNを用いた独立成分分析について質問がありメールいたしました。
    基本的なことで申し訳ありませんが、ご指導いただければ幸いです。

    独立成分分析を行った際の、spatial componentの結果の部分で、
    1.右側にできる脳の図のカラーバーの数値は何を示していますか?
    2.p-unccrect < 1 デフォルトが1になっています この意味はなんでしょうか?
      デフォルトの1 でよいですか?
    3.条件のコントラストを逆向きにすると、正負が逆転することの意味は?
    4.display resultsにでてくるのはコントラスト間で有意な差がある領域ですか?有意でなければそのspatial componentではコントラスト間で差がないということですか?

    以上、お忙しいところ誠に申し訳ありませんが、ご指導いただけると幸いです

    • 石井先生

      今、諸事情でPCにがっつり触れる環境でないので、お返事に数日ください。よろしくお願いします。

    • 返信が遅くなってすみません。

      独立成分分析を行った際の、spatial componentの結果の部分で、
      1.右側にできる脳の図のカラーバーの数値は何を示していますか?

      →T値を示しています。

      2.p-unccrect < 1 デフォルトが1になっています この意味はなんでしょうか?   デフォルトの1 でよいですか?

      p<1 というのは、すべての結果を表示という意味になります。もし、p<0.05と設定すれば、一般的な統計学的に有意な結果だけ表示するということになります。
      p<1なので、統計学的有意かどうかは関係なく、すべての領域の結果を出しなさいということです。

      3.条件のコントラストを逆向きにすると、正負が逆転することの意味は?

      →帰無仮説と対立仮説の関係です。

      帰無仮説:コネクティビティは0である

      対立仮説:コネクティビティは正である。
      の場合、コントラストは 1 となります。

      コントラストを -1 とすると
      対立仮説:コネクティビティは負である。

      となりますので、もし、そこで正のコネクティビティがある場合、統計値としては負になります。

      4.display resultsにでてくるのはコントラスト間で有意な差がある領域ですか?有意でなければそのspatial componentではコントラスト間で差がないということですか?

      display results の画面の一番上に統計の閾値が記載されています。その統計の閾値を上回るものが表示されています。

      以上、よろしくお願いします。

  7. 脊椎外科 大学院生です
    いつもこちらで勉強させていただいています。
    今マウス痛みモデルのresting state fMRIを麻酔下、覚醒下で行っているのですが、CONNを利用することは可能でしょうか。教えて下さい。よろしくおねがいします。

    • ありむらさん

      CONNを利用すること自体はなんとかできるのですが、

      – 前処理は別に行う必要がある。
      – マウスのROIは別に準備する必要がある。

      など、ハードルはなかなか高いです…。

      まず、SPMでマウスの標準化ができるスキルを身につけるところからでしょうか。
      SPMのマウスのテンプレートは、下記に提供されています。
      https://www.nitrc.org/projects/tpm_mouse

      •  Kiyotaka 様
        ご連絡ありがとうございます。遅くなり申し訳ありません。
        それを聞いて安心いたしました。
        >CONNの前処理は別に行う必要がある
        とのことですが、
        幸い解剖画像を用いて構造画像、機能画像を標準化することまではできました。
        しかし、その次に、Kiyotaka様の作成した手順のどこから進めればいいかわかりません。
        予め前処理しているので、Preproccessingは飛ばして1st level covariatesに進んでいいということでしょうか。それともpreproccessingをするべきなのでしょうか。

        また、マウスのROIも準備できております。ありがとうございます。

        ご返信お待ちしております。
        ありむら

        • すみません、このコメントが、スパムに分類されてしまっていました。

          予め前処理しているので、Preproccessingは飛ばして1st level covariatesに進んでいいということでしょうか。それともpreproccessingをするべきなのでしょうか。

          Preprocessingは飛ばしていいのですが、Settingのところで、GM/WM/CSFを指定したり、ROIを指定したりという作業を丁寧にやっていただければ、その次のプロセスに進めると思います。

          根本清貴

      • 根本先生
        お返事がうまくできておらず、申し訳ありません。再度送らせていただきます。
        幸い、マウスのMEMRIの解析をやっていたので、前処理
        や、ROI、テンプレートの準備はできております。逆にここができれば、あとはチュートリアルを参考にすれば良いでしょうか。
        ご連絡お待ちしております。
        ありむら

        • そうなんですね。
          それならばできると思います。
          細かい課題は出てくるかと思いますが、ほかの研究者の先生とサルについては進められています。

          根本清貴

  8. 根本先生

    ご多忙のところ、失礼いたします。私は昨年度より、MRI研究に初めて従事している者です。先生のVBM書籍ならびにホームページでの解析手法をいつも参考に大変お世話になっております。この度、CONNでの機能結合解析について、ご教授いただきたくご連絡いたしました。

    以下、基本的なご質問となり誠に申し訳ございません。先の方の投稿にもありますように、seed-to-voxelとroi-to-roi解析にて、任意のseed(roi)座標をnii画像ファイルで作成し、そのファイルをインストールされたconnディレクトリの中のroisサブディレクトリに保存・connのSETUPから走らせることで解析で可能であると理解しております。

    大変お恥ずかしい話ではありますが、任意のseed(roi)座標のnii画像ファイルの作成について、どのようなソフトウェア(MarsBar、WFU PickAtlus、SPM Anatomy toolbox)を用い、どのような方法で作成するのでしょうか?

    大変お手数をお掛けいたしますが、お返事をいただけましたら幸いでございます。

    • 山下さん

      ご質問ありがとうございます。

      どのようなROIを作成したいかによるかで使うソフトウェアが変わってきます。
      ROIとしてどのようなものをイメージされていらっしゃいますか?

      多くは以下の2パターンかと思います。

      1. 群間比較で得られた領域
      2. MNI座標で指定して、特定の大きさを持つ球体、もしくは立方体

      根本清貴

      • 根本先生

        ご多忙のところ、ご連絡いただき誠にありがとうございます。

        ROIのイメージとしましては、「VBM解析により、1. 群間比較で得られた領域」となります。

        大変お手数をお掛けいたしますが、ご連絡をいただけましたら幸いでございます。

        山下

      • 根本先生

        ご多忙のところ、ご連絡いただき誠に感謝申し上げます。

        CONNで用いるためのROI作成に関しましては、VBM解析により
        「1. 群間比較で得られた領域」となります。
        もし、同じ内容のメールが複数回に亘り届いておりましたら申し訳ございません。

        大変お手数をお掛けいたしますが、ご教授いただけましたら幸いでございます。

        山下

        • 山下さん

          昨日は連休直前で外来が大変なことになっていたので反応できませんでした。すみません。

          SPMの結果からROIを作成する方法は以下のとおりになります。

          1. SPMで結果を出します。
          2. ROIを作成したい結果にカーソルを移動させます。
          3. 左下のWindowにある “Save…” から、 “Current Cluster” を選択します。
          4. そうするとファイル名を聞かれますので任意の名前をいれます。
          5. そうするとワーキングディレクトリ内にniiファイルが保存されます。

          これでやってもらえたらと思います。

          なお、私のブログはスパム対策で私が承認しないとコメントが表示されないようにしてあります。
          ご了承ください。

          根本清貴

          • 根本先生

            ご多忙のところ、ご教授いただき誠に感謝申し上げます。

            上記の通り、SPM結果上でROI作成を試み、CONNによる一連の解析を
            無事終えることができました。引き続き、本方法による再解析を進めます。

            この度はROI作成からCONNによるFunctional Connectivity解析の手法について
            ご教授いただき誠に感謝申し上げます。

            山下

          • 山下さん

            無事にできたようでよかったです。

            根本清貴

  9. チュートリアルを利用させていただいている学生です。
    ROI-to-ROIとSeed-to-Voxelの解析において任意のSeed(座標?)は指定できるのでしょうか?
    もし方法があるのならばお教えください。よろしくお願いします。

    • はい、可能です。

      まず、Seedを作成する必要があります。

      自分の関心領域だけを1、それ以外は0にしたnifti画像を作成します。
      そして、それを、connフォルダの中のroisファイルに保存します。

      これで、CONNの解析を走らせると、そのseedファイルを解析に使うことができるようになります。

      これがやりたいことでしょうか?

      • ご連絡遅くなりすみません。
        ご指示の通り行ったら無事任意のSeedで解析することができました。
        ありがとうございます。

        • ご丁寧にありがとうございます。
          無事にできたようでよかったです。

  10. 初めて、お便りします。広島共立病院リハビリ科の吉川と申します。非常に基本的な質問で申し訳ありません。14:QA(Quality Assurance)、にやっと到達しましたが、最後のSelect 4 subjects の選択が、subject1,2,16,19 と間が空いているため、マウスで指定できません。クリックを押しながら選択するとなっていますが、具体的に、どうすればよろしいのでしょうか。当方のダウンロードしたconnは17a です。御教示頂けましたら、幸甚と存じます。

    • ご指摘ありがとうございます。これは私の間違いです。
      正しくは、「Ctrlキーを押しながら」です。

      Ctrlキーを押しながら選択することで、とびとびの値を選択できます。
      (ちなみにこれはExcelなどでも使える技ですので、知っていると便利です。)

      ブログの該当部位も修正させていただきました。

      以上、よろしくお願いします。

      • Kiyotaka 様
        早速の御返事、ありがとう御在ます。御指示の通り実施し、上手く行きました。
        感激です。
        今後ともよろしくお願い致します。

        • 無事にできたようでよかったです。続きはそのうち載せますので、よろしくお願いします。

ありむら へ返信するコメントをキャンセル

このサイトはスパムを低減するために Akismet を使っています。コメントデータの処理方法の詳細はこちらをご覧ください