2010/01/14 (Thu) ねこでもできるfMRI解析その３～苦難のろーからいざ編

あー。
画像調整の話追記。

そもそも個人の機能的ローカライザーはしらせてんのにNormaliseかけたら意味ないんじゃないの？とつっこまれた。
…よくよく考えてみたらまったくその通りだー。
うぼあー。

多少ネタバレ（？）すると、今回のわたしの実験は、機能的ローカライザーをはしらせて、その部位がとある課題中にどんなふーに活動するのかなーってかんじのデザインなので、古典的なコントラスト検出用の解析手順したらおかしいわけですわ。
はじめて使うからってうっかりなんでもかんでもふむふむしてる場合とちゃうわー。
いくらねこのひたい程度のあたまとはいえもうちょっと考えて解析せなあきませんね。

…と思ってたんですがSPM5でなぜかNomaliseとばして解析するとエラーはきよる。
なんかorientationがどうとかいうエラー。
Nomalise段階でやってるresliceかなにかが影響してるのか？と思ったけどCoregisterのResliceかましても同じエラーを吐く。
謎。
どうしたらいいんだ。


あとSmoothingてきとーすぎた。
基本はボクセルの2倍だそうです。
なんでもかんでもデフォルトでやるのは手順を覚えるにはいいけど、ちゃんと結果が欲しいときにはかえって遠回りになるなあ。

まあそのへんの反省もこめつつねこ解析（略しすぎ）その３いっくよー。
今日はとうとう、解析っぽいことをしてみるの巻き。
…ただの自分の解析日記じゃないかこれ！
……まあそうやって手順を覚えるのだからいっかー。

今回の目的：個人差込みで（集団解析にしないで）ローカライザーのようすをみてみる。

うらばなし暴露：まあ「はじめてのろーからいざー」の項でほぼばらしてるけど、要するにFFAやらPPAを出したいわけですが。
撮像さしてくれたひとのなかに、みみのうしろのデータがまっくろく抜けてしまうひとがいたのさ…
なんかみみのうしろに空気はいるとこが大きいとそういうことになるらしい。
（なんでそんなとこに空気はいるとこがあるのかは謎。調べたらわかるだろうか？）
で、側頭領域のものは場合によっちゃあ撮れないよ、と。
だからコロナルにきったほうがいいかもねー、とか。
最悪そのひとは分析から除外だねー、とか。
論文ではちらちらそういうケースがあるの知ってたけど実際自分がそう言われるとさびしいね！！！！
まあうまいことくろいあなを避けて活動があればいいんだけど。

というわけで鋭意ろーからいざーの解析だー。

用意するもの。
調整済みデータ。
条件ごとのオンセット時のスキャン枚数なり時間（秒単位）なり。

まずは解析用モデル作り。

１）個人ごとの解析なので、Specify 1st-levelというぼたんをぽちり。
２）Directoryに保存するフォルダを選択
３）Units for designを選択→Scansで数えるかSecondsで数えるか選択（今回はScan)
４）Interscan intervalを選択→TRを秒単位で入力
５）Data & Designを選択→右のNew "Subject/Session"クリック
６）+Subject/Sessionを展開→Scansをクリックして用意したデータ選択
７）Conditionsを選択→New "Conditions"クリック（これが要因になるので、要因数だけ増やす）
８）+Conditionを展開→Nameを選択→右のSpecify Textをクリック→条件名記入
９）Onsetsを選択→右のSpecify Textをクリック→この条件のオンセットをそれぞれ記入（間はスペースまたはタブ）
10)　Durationを選択→右のSpecify Textをクリック→Duration（Scans選んだならスキャン枚数で）記入
11)　でsave→.matを保存でRun！

これでResultボタンがつかえるようになるよ！
結果保存用フォルダにSPM.matというファイルができてるから、これを使っていざコントラスト。

１）Resultをぽっちし
２）ファイル選べっていわれるからさっきできたてほやほやなSPM.matを選択
３）コントラストを作る画面がでてくる。Define New Contrastでコントラストの行列つくろう。
４）さいしょはF-Contrastでまんべんない単位行列つくる。
５）つぎにT-Contrastで見たい効果の行列作る
主効果見たかったら、主効果見たい条件を1、他を0にした行列。
引き算したかったら、引く条件のとこは-1にしたらよい。
ただ合計がプラマイゼロになるような記述。

これがすむと次は仮説検定。
１）mask with other contrastsってのをきかれる（とりあえずNoで）
２）title for comparisonを聞かれる（なまえなんでわかりやすくてきとーに）
３）P value adjustedment for contorolを選択
None：いわゆるUncorrectedってやつ。いちばんゆるい基準
FDR：そこそこの基準。Correctedですね
FWE：結構きびしいめCorrectedな基準
これはお好みで＆データと相談して。
４）有意水準はとりあえずデフォルトで
（Uncorrectedだと勝手にきびしめ水準あたる。すげえ）
５）Extent threshold（ボクセルのかたまりを検出する下限）を聞かれる
まあお好みで。

んでぽちっとなー。
これでなんとか…なればいいけど。

ねむすぎて改訂の改訂がひつようかもしれん。
まあ気づいたことがあったらあとからなおす。



以下愚痴。

某所で「MRIのしくみもろくすっぽわかってないやつがfMRI実験なんかすんな」という文章を見てちょっとへこむ。
…そうですねまったくもってそうですよねー。

とはいえ古くはモル計算からつまづいただめだめの文系っ子に原子核のスピンが励起して、とか言われてもどうしたらいいんだ。
そいやわたしPCがどうして動くのかも理解せずにPC使ってるなあ。

何にもわかってないやつがてきとーにむちゃくちゃなネタフリをしたらぐだぐだになるのはわかってるけど。
あほのこだって実験したいこともあるよー。

じゃあ、あほのこはあほなりに何が理解できて、最低限何を押さえておけば実験ができるのか。
あほのこ目線というのは、頭の良いひとには共有されにくいものだから、わたしはそこんとこをこの恥さらしめもぶろぐでやっとけばいいんじゃないかな、と思った。
あほのこが何もわからないままてきとーなことやってたら怒られてもしかたがない（てゆかそんな実験じゃ論文通らない）けど、あほのこがあほなりに努力することはできるはずだ！
というわけで結局あほのこのまま暴走する路線は変わらず。
だってしょうがないやねー。
あほのこが一朝一夕でかしこくなったりできないものー。
ちょびっとずつ最低ラインを這いずり回るしかないさー。


（追記）
モデル作るときのオンセット入力、しょっぱなは１じゃなくて０だった・・・
考えてみたら、１スキャン後じゃなくて０スキャン目が開始なんだからそーだよなー。
うわあああ。うわああああああああああ。ぜんぶやりなおしかよー。
こういうばかもいるので気をつけましょう。
うっわ。うっわあ・・・
自分のあほさに若干ひくわ。