@线性拟合函数:regress()**
调用格式:b=regress(y,X)
[b,bint,r,rint,stats]= regress(y,X)
[b,bint,r,rint,stats]= regress(y,X,alpha)
该函数求解线性模型:y=Xβ+ε
β是p1的参数向量;ε是服从标准正态分布的随机干扰的n1的向量;y为n1的因变量向量;X为np自变量矩阵。
bint返回β的95%的置信区间。r中为形状残差,rint中返回每一个残差的95%置信区间。Stats向量包含R2统计量、回归的F值和p值。
例 设y的值为给定的x的线性函数加服从标准正态分布的随机干扰值得到。即y=10+x+ε ;求线性拟合方程系数。
x=[ones(10,1) (1:10)’]
y=x*[10;1]+normrnd(0,0.1,10,1)
[b,bint, r,rint,stats]=regress(y,x,0.05)
rcoplot(r,rint)