x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]'
X=[ones(16,1) x]
Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]'
[b,bint,r,rint,stats]=regress(Y,X)
rcoplot(r,rint)
用polyfit(X,Y,1)得到的拟合函数只能得到a,b,但不能得到线性相关系数R^2。如想要得到其线性相关系数,可以用regress(y,X),其使用格式[b,bint,r,rint,stats]
=
regress(y,X)
b——拟合系数
bint——b的置信区间
r——残差值
rint——r的置信区间
stats——检验统计量,第一个就是相关系数
例如:
x=[。。。]y=[。。。]
X=[x
ones(n,1)]
%x的行数(列数)
[b,bint,r,rint,stats]
=
regress(y,X)