第二节 回归分析
2.1 regress 命令
b = regress(y,X)
[b,bint] = regress(y,X)
[b,bint,r] = regress(y,X)
[b,bint,r,rint,stats] = regress(y,X)
[b,bint,r,rint,stats] = regress(y,X,alpha)
=
=
pnn
p
p
n
xx
xx
xx
X
y
y
y
y
L
MMMM
L
L
M
2
222
121
2
1
1
1
1
=
p
b
b
b
b
M
2
1
bint是置信度为 100(1- alpha)%的 b的区间估计。它是一个 p×2的矩阵。
缺失情况下,alpha 为 0.05。r 是残差,即实际值与估计值之差,它一个 n×1 的矩阵。rint
置信度为 100(1- alpha)%的 r 的区间估计。如果 rint 的第 i 个区间不包括 0,那么(x
i
,y
i
)
为野值。
YXXXb ')'(
1−
=
bint 的算法:第 i 行的区间为:
个对角元的第( iXXSEpntb
12
2
)'()
−
ו−±
α
r 的计算:
Xbyyye −=−=
∧
rint 的计算,参考 2.2 节
i
i
i
i
i
i
i
i
studres
r
leverage
ileveragepntr
=−
−−±
∧
∧∧
1
s2_i2.2(1)(
2
σ
σσ
α
)
节的个数据点误差的估计。为除掉第,
stats 包括四项
1、 可决系数