Global Linear Model

fMRI=β0+β1BOLD+error
Y=Xβ+ϵ
minϵ 會得到β^=(XtX)1XtY
———————————————————————–
因爲ϵN(0,δ2I)
推出
Y ~N(Xβ,δ2I)
β^ ~N(β,(XtX)1σ2)
βj^ ~N(βj,(XtX)1jjσ2)
βj^βj ~N(0,(XtX)1jjσ2)
所以βjβj^(XTX)1jjσ2 ~N(0,1)
由於σ2 (standard deviation)不能得到,但是可以計算σ^2 (standard error)
所以該問題轉化爲一個服從t分佈的問題:
βjβj^SE ~tnp1 (SE=(XTX)1jjσ^2 )
—————————————————————————————————————-
如果Xi 與Y沒有線性關係,那麼βj 應當爲0,而估計得到的βj^ 只不過是噪聲而已。
——————————————————————————–
計算(βjβj^) 的p-value
P(|(0βj^)|/SE)=P(|βj^/SE|)
如果p-value<0.05,認爲有顯著差異,βj^ 不是一個噪聲,而是反映線性相關性質的參數。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章