统计回归模型举例.pptx
《统计回归模型举例.pptx》由会员分享,可在线阅读,更多相关《统计回归模型举例.pptx(59页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、会计学1统计回归模型举例统计回归模型举例例例1 1 牙膏的销售量牙膏的销售量 问问题题建立牙膏销售量与价格、广告投入之间的模型建立牙膏销售量与价格、广告投入之间的模型 预测在不同价格和广告费用下的牙膏销售量预测在不同价格和广告费用下的牙膏销售量 收集了收集了3030个销售周期本公司牙膏销售量、价格、个销售周期本公司牙膏销售量、价格、广告费用,及同期其它厂家同类牙膏的平均售价广告费用,及同期其它厂家同类牙膏的平均售价 9.260.556.804.253.70307.930.055.803.853.80298.510.256.754.003.7527.38-0.055.503.803.851销售量
2、销售量(百万支百万支)价格差价格差(元)(元)广告费用广告费用(百万元百万元)其它厂家其它厂家价格价格(元元)本公司价本公司价格格(元元)销售销售周期周期第1页/共59页 n n令令令令y y y y表示公司牙膏的销售量,表示公司牙膏的销售量,表示公司牙膏的销售量,表示公司牙膏的销售量,n nx x1 1表示表示表示表示其它厂家与本公司价格差,其它厂家与本公司价格差,其它厂家与本公司价格差,其它厂家与本公司价格差,x x2 2 表示表示表示表示公司广告费用,则数据如下:公司广告费用,则数据如下:公司广告费用,则数据如下:公司广告费用,则数据如下:n n x1=-0.05 0.25 0.6 0
3、0.25 0.2 0.15 x1=-0.05 0.25 0.6 0 0.25 0.2 0.15 x1=-0.05 0.25 0.6 0 0.25 0.2 0.15 x1=-0.05 0.25 0.6 0 0.25 0.2 0.15 0.05-0.15 0.15 0.2 0.1 0.4 0.45 0.35 0.3 0.05-0.15 0.15 0.2 0.1 0.4 0.45 0.35 0.3 0.05-0.15 0.15 0.2 0.1 0.4 0.45 0.35 0.3 0.05-0.15 0.15 0.2 0.1 0.4 0.45 0.35 0.3 0.5 0.5 0.4-0.05-0.0
4、5-0.1 0.2 0.1 0.5 0.5 0.5 0.4-0.05-0.05-0.1 0.2 0.1 0.5 0.5 0.5 0.4-0.05-0.05-0.1 0.2 0.1 0.5 0.5 0.5 0.4-0.05-0.05-0.1 0.2 0.1 0.5 0.6-0.05 0 0.05 0.55;0.6-0.05 0 0.05 0.55;0.6-0.05 0 0.05 0.55;0.6-0.05 0 0.05 0.55;n nx2=5.5 6.75 7.25 5.5 7 6.5 6.75 5.25 x2=5.5 6.75 7.25 5.5 7 6.5 6.75 5.25 x2=5.5
5、6.75 7.25 5.5 7 6.5 6.75 5.25 x2=5.5 6.75 7.25 5.5 7 6.5 6.75 5.25 5.25 6 6.5 6.25 7 6.9 6.8 6.8 7.1 7 6.8 5.25 6 6.5 6.25 7 6.9 6.8 6.8 7.1 7 6.8 5.25 6 6.5 6.25 7 6.9 6.8 6.8 7.1 7 6.8 5.25 6 6.5 6.25 7 6.9 6.8 6.8 7.1 7 6.8 6.5 6.25 6 6.5 7 6.8 6.8 6.5 5.75 5.8 6.5 6.25 6 6.5 7 6.8 6.8 6.5 5.75 5
6、.8 6.5 6.25 6 6.5 7 6.8 6.8 6.5 5.75 5.8 6.5 6.25 6 6.5 7 6.8 6.8 6.5 5.75 5.8 6.8;6.8;6.8;6.8;n n y=7.38 8.51 9.52 7.5 9.33 8.28 8.75 y=7.38 8.51 9.52 7.5 9.33 8.28 8.75 y=7.38 8.51 9.52 7.5 9.33 8.28 8.75 y=7.38 8.51 9.52 7.5 9.33 8.28 8.75 7.87 7.1 8 7.89 8.15 9.1 8.86 8.9 8.87 7.87 7.1 8 7.89 8.
7、15 9.1 8.86 8.9 8.87 7.87 7.1 8 7.89 8.15 9.1 8.86 8.9 8.87 7.87 7.1 8 7.89 8.15 9.1 8.86 8.9 8.87 9.26 9 8.75 7.95 7.65 7.27 8 8.5 8.75 9.26 9 8.75 7.95 7.65 7.27 8 8.5 8.75 9.26 9 8.75 7.95 7.65 7.27 8 8.5 8.75 9.26 9 8.75 7.95 7.65 7.27 8 8.5 8.75 9.21 8.27 7.67 7.93 9.26;9.21 8.27 7.67 7.93 9.26
8、;9.21 8.27 7.67 7.93 9.26;9.21 8.27 7.67 7.93 9.26;第2页/共59页 n n下面探讨下面探讨下面探讨下面探讨y y y y与与与与x x x x1 1 1 1、x x x x2 2 2 2的关的关的关的关系:系:系:系:n n用用用用matlabmatlabmatlabmatlab软件作图:软件作图:软件作图:软件作图:n nplot(xplot(xplot(xplot(x1 1 1 1,y,y,y,y,*););););n nplot(xplot(xplot(xplot(x2 2 2 2,y,y,y,y,*)n n运行得如下图形:运行得如下图
9、形:运行得如下图形:运行得如下图形:x1y从右图看出,从右图看出,y y与与x x1 1成成线性关系,线性关系,y y与与x x2 2成二成二次曲线关系。次曲线关系。x2y第3页/共59页 n n x3=x2.2;x3=x2.2;x3=x2.2;x3=x2.2;n n x=ones(30,1)x1 x2 x3;x=ones(30,1)x1 x2 x3;x=ones(30,1)x1 x2 x3;x=ones(30,1)x1 x2 x3;n n b,bint,r,rint,stats=regress(y,x)b,bint,r,rint,stats=regress(y,x)b,bint,r,rint
10、,stats=regress(y,x)b,bint,r,rint,stats=regress(y,x)n n运行结果:运行结果:运行结果:运行结果:n nb=17.3244b=17.3244b=17.3244b=17.3244,1.30701.30701.30701.3070,-3.6956-3.6956-3.6956-3.6956,0.34860.34860.34860.3486n nbint=5.7282 28.9206bint=5.7282 28.9206bint=5.7282 28.9206bint=5.7282 28.9206n n 0.6829 1.9311 0.6829 1.93
11、11 0.6829 1.9311 0.6829 1.9311 n n -7.4989 0.1077 -7.4989 0.1077 -7.4989 0.1077 -7.4989 0.1077 n n 0.0379 0.6594 0.0379 0.6594 0.0379 0.6594 0.0379 0.6594n nstats=0.9054stats=0.9054stats=0.9054stats=0.9054,82.940982.940982.940982.9409,0.00000.00000.00000.0000,0.04900.04900.04900.0490模型求解模型求解MATLAB 统
12、计工具箱统计工具箱 第4页/共59页结果分析结果分析y的的90.54%可由模型确可由模型确定定 参数参数参数估计值参数估计值置信区间置信区间17.32445.7282 28.92061.30700.6829 1.9311-3.6956-7.4989 0.1077 0.34860.0379 0.6594 R2=0.9054 F=82.9409 p=0.0000 0 1 2 3F远超过远超过F检验的临界检验的临界值值 P x=ones(30,1)x1,x2(x2.2)(x1.*x2);x=ones(30,1)x1,x2(x2.2)(x1.*x2);x=ones(30,1)x1,x2(x2.2)(x
13、1.*x2);x=ones(30,1)x1,x2(x2.2)(x1.*x2);n n b,bint,r,rint,stats=regress(y,x)b,bint,r,rint,stats=regress(y,x)b,bint,r,rint,stats=regress(y,x)b,bint,r,rint,stats=regress(y,x)n nb=29.1133b=29.1133b=29.1133b=29.1133n n 11.1342 11.1342 11.1342 11.1342n n -7.6080 -7.6080 -7.6080 -7.6080n n 0.6712 0.6712 0.
14、6712 0.6712n n -1.4777 -1.4777 -1.4777 -1.4777n nbint=3.7013 44.5252bint=3.7013 44.5252bint=3.7013 44.5252bint=3.7013 44.5252n n 1.9778 20.2906 1.9778 20.2906 1.9778 20.2906 1.9778 20.2906n n -12.6932 -2.5228 -12.6932 -2.5228 -12.6932 -2.5228 -12.6932 -2.5228n n 0.2538 1.0887 0.2538 1.0887 0.2538 1.
15、0887 0.2538 1.0887n n -2.8518 -0.1037 -2.8518 -0.1037 -2.8518 -0.1037 -2.8518 -0.1037n nstats=0.9209stats=0.9209stats=0.9209stats=0.9209,72.777172.777172.777172.7771,0.00000.00000.00000.0000,0.04260.04260.04260.0426第8页/共59页模型比较模型比较x1和和x2对对y的的影响独影响独立立 参数参数参数估计值参数估计值置信区间置信区间17.32445.7282 28.92061.3070
16、0.6829 1.9311-3.6956-7.4989 0.1077 0.34860.0379 0.6594 R2=0.9054 F=82.9409 p=0.0000 0 1 2 3参数参数参数估计值参数估计值置信区间置信区间29.113313.7013 44.525211.13421.9778 20.2906-7.6080-12.6932 -2.5228 0.67120.2538 1.0887-1.4777-2.8518 -0.1037 R2=0.9209 F=72.7771 p=0.0000 3 0 1 2 4x1和和x2对对y的影响的影响有交互作有交互作用用第9页/共59页 n n由于由
17、于R R2 2有所提高,所以模型有所提高,所以模型(*)(*)比模型(比模型(*)有所改进,并)有所改进,并且参数的置信区间不再包含且参数的置信区间不再包含0 0点,点,所以有理由认为模型(所以有理由认为模型(*)比)比模型(模型(*)更符合实际。)更符合实际。n n预测比较:预测比较:n nx1=0.2;x2=6.5;x1=0.2;x2=6.5;n nY=b(1)+b(2)*x1+b(3)*x2+b(4)Y=b(1)+b(2)*x1+b(3)*x2+b(4)*(x2.2)+b(5)*(x1.*x2)*(x2.2)+b(5)*(x1.*x2)n nY=8.3272Y=8.3272第10页/共5
18、9页两模型销售量预测两模型销售量预测比较比较(百万支百万支)区间区间 7.8230,8.7636区间区间 7.8953,8.7592(百万支百万支)控制价格差控制价格差x1=0.2元,投入广告费元,投入广告费x2=6.5百万元百万元预测区间长度更短预测区间长度更短 略有增加略有增加 第11页/共59页完全二次多项式模型完全二次多项式模型 x=x1 x2;x=x1 x2;rstool(x,y,quadratic)rstool(x,y,quadratic)运行结果:运行结果:beta=2.0984beta=2.0984 14.7436 14.7436 -8.6367 -8.6367 -2.1038
19、 -2.1038 1.1074 1.1074 0.7594 0.7594rmse=0.2083rmse=0.2083剩余标准差为剩余标准差为0.2.830.2.83较较小,说明回归模型的显小,说明回归模型的显著性比较好。著性比较好。第12页/共59页n n问题:一家高技术公司人事部门为研究软件开问题:一家高技术公司人事部门为研究软件开问题:一家高技术公司人事部门为研究软件开问题:一家高技术公司人事部门为研究软件开发人员的薪金与他们的资历、管理责任、教育发人员的薪金与他们的资历、管理责任、教育发人员的薪金与他们的资历、管理责任、教育发人员的薪金与他们的资历、管理责任、教育程度等因素之间的关系,要
20、建立一个数学模型,程度等因素之间的关系,要建立一个数学模型,程度等因素之间的关系,要建立一个数学模型,程度等因素之间的关系,要建立一个数学模型,以便分析公司人事策略的合理性,并作为新聘以便分析公司人事策略的合理性,并作为新聘以便分析公司人事策略的合理性,并作为新聘以便分析公司人事策略的合理性,并作为新聘人员的薪金的参考。他们认为目前公司人员的人员的薪金的参考。他们认为目前公司人员的人员的薪金的参考。他们认为目前公司人员的人员的薪金的参考。他们认为目前公司人员的薪金总体上是合理的,可以作为建模的依据。薪金总体上是合理的,可以作为建模的依据。薪金总体上是合理的,可以作为建模的依据。薪金总体上是合理
21、的,可以作为建模的依据。于是调查了于是调查了于是调查了于是调查了46464646名软件开发人员的档案资料,如名软件开发人员的档案资料,如名软件开发人员的档案资料,如名软件开发人员的档案资料,如下表,其中资历一列指从事专业工作的年数,下表,其中资历一列指从事专业工作的年数,下表,其中资历一列指从事专业工作的年数,下表,其中资历一列指从事专业工作的年数,管理一列中:管理一列中:管理一列中:管理一列中:1 1 1 1表示管理人员,表示管理人员,表示管理人员,表示管理人员,0 0 0 0表示非管理人表示非管理人表示非管理人表示非管理人员,教育一列中:员,教育一列中:员,教育一列中:员,教育一列中:1
22、1 1 1表示中学程度,表示中学程度,表示中学程度,表示中学程度,2 2 2 2表示大学表示大学表示大学表示大学程度,程度,程度,程度,3 3 3 3表示更高程度(研究生)。表示更高程度(研究生)。表示更高程度(研究生)。表示更高程度(研究生)。例例2 2 软件开发人员的薪软件开发人员的薪金金第13页/共59页 编号编号薪金薪金资历资历管理管理教育教育编号编号薪金薪金资历资历管理管理教育教育1 113876138761 11 11 1131319800198003 31 13 32 211608116081 10 03 3141411417114174 40 01 13 31870118701
23、1 11 13 3151520263202634 41 13 34 411283112831 10 02 2161613231132314 40 03 35 511767117671 10 03 3171712884128844 40 02 26 620872208722 21 12 2181813245132455 50 02 27 711772117722 20 02 2191913677136775 50 03 38 810535105352 20 01 1202015965159655 51 11 19 912195121952 20 03 3212112366123666 60 01
24、 1101012313123133 30 02 2222221352213526 61 13 3111114975149753 31 11 1232313839138396 60 02 2121221371213713 31 12 2242422884228846 61 12 2第14页/共59页编号编号薪金薪金资历资历管理管理教育教育编号编号薪金薪金资历资历管理管理教育教育252516978169787 71 11 13636168821688212120 02 2262614803148038 80 02 23737241702417012121 13 3272717404174048 8
25、1 11 13838159901599013130 01 1282822184221848 81 13 33939263302633013131 12 2292913548135488 80 01 14040179491794914140 02 23030144671446710100 01 14141256852568515151 13 33131159421594210100 02 24242278372783716161 12 23232231742317410101 13 34343188381883816160 02 23333237802378010101 12 244441748
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 统计 回归 模型 举例
限制150内