多元回归分析-2.ppt

化为多元线性回归求解 x1=[1000,600,1200,500,300,400,1300,1100,1300,300]; x2=[5,7,6,6,8,7,5,4,3,9]; y=[100,75,80,70,50,65,90,100,110,60]; X=[ones(10,1) x1 x2 (x1.^2) (x2.^2)]; [b,bint,r,rint,stats]=regress(y,X) 回归系数的点估计以及区间估计 残差及其置信区间 检验回归模型的统计量 一、多元线性回归的数学模型 二、数学模型的分析与求解 三、MATLAB中回归分析的实现 四、小结 　多元线性回归 一、多元线性回归的数学模型 多元线性回归模型 用最大似然估计法估计参数. 达到最小. 二、数学模型的分析与求解 化简可得 正规方程组 引入矩阵 正规方程组的矩阵形式 最大似然估计值 P元经验线性回归方程 多元线性回归 　　1.确定回归系数的点估计值,用命令: 　　　b=regress(Y,X) 　　2.求回归系数的点估计和区间估计,并检验回 归模型,用命令: 　　　[b,bint,r,rint,stats]=regress(Y,X,alpha) 　　3.画出残差及其置信区间,用命令: 　　　rcoplot(r,rint) 三、MATLAB中回归分析的实现 符号说明 (1) (2) alpha为显著性水平, 默认为 0.05; (3) bint为回归系数的区间估计; (4) r与rint分别为残差及其置信区间; 　　(5) stats 是用于检验回归模型的统计量, 有三个 数值, 第一个是相关系数 r2, 其值越接近于 1, 说明回 归方程越显著; 第二个是 F 值, FF1-alpha(p,n-p-1) 时 拒绝 H0, F 越大, 说明回归方程越显著; 第三个是与 F对应的概率 p, palpha 时拒绝, 回归模型成立. 身高 143 145 146 147 149 150 153 154 腿长 88 85 88 91 92 93 93 95 身高 155 156 157 158 159 160 162 164 腿长 96 98 97 96 98 99 100 102 例1　测得16组的身高和腿长如下(单位:cm): 试研究这些数据之间的关系. 输入数据 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); b,bint,stats 一元多项式回归 　　1.确定多项式系数,用命令: 　　　[p,S]=polyfit(x,y,m) 也可使用命令:polytool(x,y,m) 　　结果产生一个交互式的画面, 画面中有拟合曲 线和 y 的置信区间, 左下方的 Export 可以输出参数. 2.预测和预测误差估计用命令: 求回归多项式在x处的预测值Y. [Y,DELTA]=polyconf(p,x,S,alpha) 　　求回归多项式在 x 处的预测值 Y 以及预测值的 显著性为1-alpha 的置信区间 Y±DELTA,alpha 的默 认值是 0.05. 一元多项式回归可化为多元线性回归求解. Y=polyval(p,x) 例2　下面给出了某种产品每件平均单价 Y(元) 与 批量 x (件) 之间的关系的一组数据 . x 20 25 30 35 40 50 y 1.81 1.70 1.65 1.55 1.48 1.40 x 60 65 70 75 80 90 y 1.30 1.26 1.24 1.21 1.20 1.18 试用一元二次多项式进行回归分析. 输入数据 x=[20,25,30,35,40,50,60,65,70,75,80,90]; y=[1.81,1.70,1.65,1.55,1.48,1.40,1.30,1.26,1.24,1.21, 1.20,1.18]; 作二次多项式回归 [p,S]=polyfit(x,y,2) 预测及作图 Y=polyconf(p,x,y) plot(x,y,’b+’,x,Y,’r’) 作二次多项式回归 [p,S]=polyfit(x,y,2) p=polyfit(x,y,2) p = 0.0001 -0.0225 2.1983 Y=polyval(p,x) Y = 1.7978 1.7134 1.6352 1