多元线性回归预测法在服装制造中的应用

  • 格式:docx
  • 大小:358.36 KB
  • 文档页数:4

多元线性回归预测法在服装制造中的应用乔亮亮(上海海事大学物流工程学院, 上海200135)[ 摘要] 预测是指对复杂变化的事物进行大量调查研究, 应用系统分析的方法找出事物发生变化的固有规律, 揭示事物未来的状况和面貌。

影响服装制造厂物流需求的因素很多, 可以用多元线性回归法分析其主要因素。

用多元线性回归法进行预测, 求出其数学模型, 并对其进行了误差和精度分析。

[ 关键词] 多元线性回归; 物流需求; 预测; 模型; 回归系数; SPSS 软件[ 中图分类号] O211 [ 文献标识码] A [ 文章编号] 167321409 (2010) 012N145203根据中国服装协会对部分服装产业集群抽样调查显示, 每年受到人力成本、原料成本、能源成本、政策成本上升的影响, 中小企业面临生存威胁, 特别是人数在100 人以内的小企业关、停现象普遍。

在这样的条件下, 如何提高需求预测的精度是企业所面临的问题之一。

在物流需求预测中, 物流需求的多少受到多种因素的影响, 可以通过在各相关影响因素间建立回归预测模型来实现对物流量的预测。

回归就是研究自变量与因变量之间的关系的分析方法。

已知某青年女装1993~2005 年13a 中每个季度的销售额及其相关影响因素的数据, 见文献[ 1 ] 。

下面, 笔者用多元线性回归法进行预测, 求出其数学模型, 并进行误差及精度检验。

1 估计参数如果不能确定哪些自变量应包括在变量内, 可以利用所考虑的所有变量建立一个相关矩阵[ 2 ] , 保留因变量与自变量高度相关的因素, 而把能引起多重共线性的自变量删去或替换。

考虑年份、男女比例、购买力、人均购买数、文化程度、生产状况、人口数, 销售额等因素, 由文献[ 1 ] 中数据计算其相关矩阵如表1 所示。

表1 相关矩阵年份男女比例购买力人均购买数文化程度生产状况人口数销售额1 01 22 01 042 01 44 01 17 01 37 01 36 01 7791 01 11 01 38 01 359 01 28 01 31 01 321 01 14 01 168 01 36 01 13 01 5311 01 28 01 47 01 149 01 241 01 356 01 27 01 0421 01 08 01 0091 01 131由表1 数据可知, 年份和购买力与销售额的相关系数较大, 其他因素不存在多重共线问题, 故选用年份和购买力为数学模型的参数, 分别记做X1 , X2 。

2 模型建立分别以年份和购买力为X 、Y 轴, 销售额为Z 轴绘制散点图, 如图1 。

图1 相关因素的散点图[ 收稿日期] 2009212225[ 作者简介] 乔亮亮( 19842) , 女, 2008 年大学毕业, 硕士生, 现主要从事物流系统规划方面的教学与研究工作。

1 22R 2 = 1 - 14146398 31 S 由图 1 可以看出 , 这些散点大致呈直线型 , 与 前 面假 设相 符 , 所 以 可将 该模 型设 为 二元 线性 回 归[2 ] , 销售额的预测值 :^Y = β0 +β1 X 1 +β2 X 23 系数求解回归系数的确定用自小二乘法求得 。

通过最小二乘法[ 3 ] 可得如下方程组 :ΣY = N β0 +β1ΣX 1 +β2ΣX 2 ΣX 1 Y = β0ΣX 1 +β1ΣX 2+β2ΣX 1 X 2 ΣX 2 Y = β0ΣX 2 +β1ΣX 1 X 2 +β2ΣX 2将表 1 中的数据带入上式中 , 得到 : β0 = 23321 94405 β1 = 1611 348691 β2 = - 01 0531439 所以其数学模型为 :Y = 23321 94405 + 1611 348691 X 1 - 01 0531439 X 2 依据以上计算 , 可以得销售额的预测值 , 如表 2 所示 。

4 误差及精度分析41 1 拟合程度评价因变量 y 的各个观察值点聚集在回归直线周围的紧密程度 , 称为回归直线对样本数据点的拟合程 度[ 2 ] 。

通常用可决系数 R 2 来衡量 。

它取值于 0 和 1 之间 , 并取决于回归模型所解释的 y 方差的百分比 。

可决系数 R 2 的公式为 :R 2= 1 -Σ ( Y - ^Y ) Σ ( Y - Y ) 2其中 , Σ ( y - ^y ) 2 称为残差平方和 ( 剩余平方和) ; Σ ( y - y ) 2称为离差平方和 。

显然残差平方和占离差平方和的比重越小 , 可决系数 R 2 越大 , 回归直线的拟合程度越强 。

可决系数 R 2 的取值区间为 [ 0 , 1 ] , 实际上 , 可决系数 R 2 是线性相关关系 r 的平方 , | R | 越接近于 1 , 则因 变量与自变量的线性相关关系越密切 , 回归直线拟合程度越高[ 3 ] 。

带入数据得 :1 29191984= 01 5154 R = 01 7179由此可以看出此回归模型解释了服装销售量变差的 511 54 % 。

41 2 标准误差标准误差又称剩余标准差 , 是评价回归直线代表性大小或实际值与估计值的标准误差大小的综合指 标 , 也是计算置信区间估计值和其他拟合优度的基础指标 。

计算公式如下 :S E =将数据带入得 :S E =1 682141 3 回归系数的显著性检验回归系数的显著性检验是用 t 参数检验的 。

t 服从自由度为 n - 3 的 t 分布 , 取显著性水平α= 01 05 , 查表得 t α = 21 021 , 若 | t β | > t α 则说明回归系数显著 。

在该模型中 :t β = βi S βi = 1 , 2 i βi i带入数据得 :1 1 = 61 88Sβ= 231 44 tβ 5601 6821=1611 348691231 44- 01 0531439Sβ2 =224271 284= 01 025 tβ2=01 025= - 21 31表2 预测值及其它参数年份销售额Y 预测值^Y Y - Y Y - ^Y^Y- Y X 1- X 1X 2- X 2 19931 30171 6 27981 39 - 6721 88 2191 21 - 8921 1 - 51 5 77471 72 30431 54 28001 675 - 6461 94 2421 865 - 8891 8 - 51 5 20261 93 20941 35 28001 09 - 15961 1 - 7051 74 - 8901 4 - 515 31381 34 28091 84 28001 834 - 8801 64 91 00565 - 8891 7 - 51 5 17651 7 19941 32741 8 29621 289 - 4151 68 3121 511 - 72812 - 41 5 16011 72 31631 28 29551 54 - 5271 2 2071 74 - 7341 9 - 41 5 278963 21141 31 29541 211 - 15761 2 - 8391 9 - 7361 3 - 41 5 368724 30241 57 29541 955 - 6651 91 691 6145 - 7351 5 - 41 5 31691 19951 33271 48 31151 507 - 363 2111 973 - 575 - 31 5 372572 34931 48 31181 961 - 197 3741 519 - 5711 5 - 31 5 163893 24391 93 31321 619 - 12501 6 - 6921 69 - 5571 9 - 31 5 166364 34901 79 31211 087 - 1991 69 3691 703 - 5691 4 - 315 77471 7 19961 36851 08 32821 383 - 51 4048 4021 697 - 4081 1 - 21 5 79241 72 36611 23 32801 682 - 291 255 3801 548 - 4091 8 - 21 5 146463 23781 43 32811 904 - 13121 1 - 9031 47 - 4081 6 - 215 96081 14 34591 55 32801 682 - 2301 93 1781 868 - 4091 8 - 21 5 14646 19971 38491 63 34441 316 1591 145 4051 314 - 24612 - 11 5 60871 22 37011 18 34461 335 101 6952 2541 845 - 2441 1 - 11 5 16011 73 26421 38 34411 127 - 10481 1 - 7981 75 - 24914 - 115 190504 35851 52 34381 204 - 1041 96 1471 316 - 2521 3 - 115 37257 19981 40781 66 36171 728 3881 175 4601 932 - 721 76 - 01 5 221952 39071 06 36131 211 2161 575 2931 849 - 771 27 - 01 5 40931 33 28281 46 36131 158 - 8621 02 - 7841 7 - 771 33 - 01 5 39661 44 408915 36131 105 3991 015 4761 395 - 771 38 - 01 5 38411 4 19991 43391 61 37791 927 6491 125 5591 683 891 443 01 5 272182 41481 6 37791 927 4581 115 3681 673 891 443 01 5 272183 29161 45 37791 396 - 7741 03 - 8621 79 881 752 01 5 230984 40841 64 37791 396 3941 155 3051 244 881 911 01 5 24019 20001 42421 42 39401 904 5511 935 3011 516 2501 42 11 5 249572 39971 58 39401 001 3071 095 571 5795 2491 52 11 5 198753 28811 01 39351 855 - 8091 47 - 10541 8 2451 37 11 5 39661 44 40361 23 39271 299 3451 745 1081 931 2361 81 11 5 96081 1 20011 43601 33 40911 464 6691 845 2681 866 4001 98 21 5 20261 92 43601 53 40901 88 6701 045 2691 65 40014 215 31381 33 31721 18 41071 567 - 5281 3 - 9351 39 4171 08 21 5 665534 42231 76 41021 625 5331 275 1211 135 4121 14 21 5 27218 20021 46901 48 42611 688 9991 995 4281 792 5711 2 31 5 148792 46941 48 42631 123 1004 4311 357 5721 64 315 221953 33421 35 42631 92 - 3481 13 - 9211 57 5731 44 315 268894 45771 63 42631 92 8871 145 3131 657 5731 49 31 5 27218 20031 49651 46 44181 148 12741 98 5471 312 7271 66 41 5 8981 752 50261 05 44181 52 13351 57 6071 53 7281 03 41 5 13671 53 34701 24 44171 191 - 2201 24 - 9461 95 7261 71 41 5 1431 54 45251 94 44171 191 8351 455 1081 749 7261 71 41 5 14316 20041 52581 71 45791 496 15681 23 6791 214 8891 01 51 5 8981 752 51891 58 45731 225 14991 1 6161 355 8821 74 51 5 77471 73 35961 76 45701 568 - 931 725 - 9731 81 8801 08 51 5 190504 38811 6 45731 172 1911 115 - 6911 57 8821 69 51 5 79241 7(下转第160 页)2A 样品所得的回归方程为 : l g C = 201 001 t + 11 9779 , 求得 K 1 = 21 303 ×10 - 3。