fMRI=β0+β1∗BOLD+error
即
Y=Xβ+ϵ
minϵ
会得到
β^=(XtX)−1XtY
———————————————————————–
因为
ϵ−N(0,δ2I)
推出
Y
~
β^
~
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
~
tn−p−1
(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^
不是一个噪声,而是反映线性相关性质的参数。
Global Linear Model
最新推荐文章于 2022-06-04 23:02:19 发布