概率论与数理统计9.4 多元线性回归.ppt

三、MATLAB中回归分析的实现 化为多元线性回归求解 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, bint, r, rint, stats]=regress(Y, X, alpha) 　rcoplot(r, rint) 用命令: b=regress(Y, X) 2.求回归系数的点估计和区间估计, 并检验回 归模型, 　用命令: 3.画出残差及其置信区间, 用命令: 符号说明 (1) (2) alpha为显著性水平, (3) bint为回归系数的区间估计; (4) r与rint分别为残差及其置信区间; 　　(5) stats 是用于检验回归模型的统计量, 默认为 0.05; 有三 个数值, FF1-alpha 第一个是相关系数 r2, 其值越接近于 1, 说 明回归方程越显著; 第二个是 F 值, (p,n-p-1)时， 拒绝 H0, F 越大, 说明回归方程越显著; 第三个是与F对应的概率 p, palpha 时拒绝, 模型成立. 回归 102 100 99 98 96 97 98 96 腿长 164 162 160 159 158 157 156 155 身高 95 93 93 92 91 88 85 88 腿长 154 153 150 149 147 146 145 143 身高 例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 残差分析 rcoplot(r,rint) 预测及作图 z=b(1)+b(2)*x plot(x,Y,’k+’,x,z,’r’) 数据比较 残差图形 预测图形 数据比较 帮　　助 程序运行结果 一元多项式回归 　　[p,S]=polyfit(x,y,m) 也可使用命令:polytool(x,y,m) 用命令: 1.确定多项式系数, 结果产生一个交互式的画面, 画面中有拟合曲 2.预测和预测误差估计用命令: 求回归多项式在x处的预测值Y. [Y,DELTA]=polyconf (p, x, S, alpha) 的默认值是 0.05. 一元多项式回归可化为多元线性回归求解. Y=polyval(p, x) 线和 y 的置信区间, 左下方的 Export 可以输出参数. 求回归多项式在 x 处的预测值 Y 以及预测值 的显著性为1-alpha 的置信区间 Y±DELTA, alpha 例2 某件产品每件平均单价Y(元)与批量x(件)之间 的关系的一组数据 x y 30 35 40 50 55 60 65 70 80 90 1.65 1.55 1.48 1.40 1.30 1.26 1.24 1.21 1.20 1.18 25 20 1.81 1.70 解 (XTX)-1 (XTX)-1 得到回归方程 试用一元二次多项式进行回归分析, 输入数据 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’) 回归结果 残差图形 预测图形 帮　　助 程序运行结果 化为多元线性回归 X=[ones(12,1)