非线性回归
非线性回归可用命令nlinfit,nlintool,nlparci,nlpredci来实现。
进行非线性回归时可使用nlinfit指令,其语法如下:
beta = nlinfit(X,y,fun,beta0)
[beta,r,J] = nlinfit(X,y,fun,beta0)
[...] = nlinfit(X, y, fun, beta0, options)
回归:
回归可用命令
[beta,r,J]=nlinfit(x,y,model,beta0) 或者
nlintool(x,y,model,beta0,alpha)
来实现。
其中:
命令 [beta,r,J]=nlinfit(x,y,model,beta0)的作用为确定回归系数;
而命令 nlintool(x,y,model,beta0,alpha)
产生一个交互式的画面,画面中有拟合曲线和y的置信区间。
通过左下方的Export下拉式菜单,可以输出回归系数等。
这里的输入数据:
x、y分别为 m×n 矩阵和n维列向量,对一元非线性回归,x为n维列向量;
model是事先用m-文件定义的非线性函数;
beta0是回归系数的初值。
Beta是估计出的回归系数,
r(残差)、 J(Jacobian矩阵)是估计预测误差需要的数据。
alpha为显著性水平,缺省时为0.05。
预测和预测误差估计:
预测和预测误差估计的命令格式为
betaci=nlparci(beta,r,J)
其用途为记算回归系数的置信区间;
[Y,DELTA]=nlpredci(model,x,beta,r,J)
其用途为求nlinfit或nlintool所得的回归函数在x处的预测值Y及预测值的显著性为