模型13：统计模型.pptVIP

数学建模 第十三讲 统计模型 模型一 牙膏的销售量 确定关系： 牙膏销售量——价格、广告投入 内部规律复杂??数据统计分析 常用模型??回归模型?×?数学原理??软件 30个销售周期数据： 销售量、价格、广告费用、同类产品均价 分析 y ~公司牙膏销售量 x1~其它厂家与本公司价格差 x2~公司广告费用 Matlab 统计分析 MATLAB7.0版本 s增加一个统计量: 剩余方差s2 于是 考察 y与x1 y与 x2 有 结果 模型 多元回归模型 结果分析 即： 销售量预测 控制价格差x1=0.2元，投入广告费x2=650万元 得： 模型改进 x1和x2对y的影响有交互作用 比较：销售量预测 控制价格差x1=0.2元，投入广告费x2=6.5百万元 比较：两模型y与x1,x2关系 x2=6.5 讨论：交互作用影响 价格差 x1=0.1 价格差 x1=0.3 完全二次多项式模型 MATLAB 相应面分析 相应面分析 rstool(x(:,2:3),y) all: beta rmse residuals 系数 均方差 残差 评注 回归模型 数据、经验、图形 ??回归变量、函数形式 Matlab求解 统计分析： R2、F、p 回归系数置信区间包含0点??改进 添加二次项、交叉项…… 模型二 软件开发人员的薪金 薪金——资历、岗位、学历 建立模型：分析人事策略的合理性，作为新聘用人员薪金的参考 模型假设 假设： y~ 薪金，x1 ~资历（年） x2 = 1~ 管理人员，0~ 非管理人员 模型求解 x1~资历(年) x2 = 1~ 管理，0~ 非管理 中学：x3=1, x4=0;大学：x3=0, x4=1; 更高：x3=0, x4=0 结果 R2,F, p ? 模型整体上可用 资历增加1年薪金增长546 管理人员薪金多6883 中学程度薪金比更高的少2994 大学程度薪金比更高的多148 a4置信区间包含零点 ?解释不可靠! 结果分析 残差分析法 残差分析 ? 与管理x2—教育x3 、x4的关系 模型改进 增加管理x2与教育x3, x4的交互项 残差分析 消除了不正常现象 异常数据(33号)??去掉 模型改进 去掉异常数据后的结果 残差分析 残差图正常 模型的结果??可以应用 模型应用 制订基础薪金 资历为0 ： x1= 0 管理—教育组合：6种 评注 对定性因素：如管理、教育 可以引入0-1变量处理 0-1变量的个数应比定性因素的水平少1 残差分析：可以发现模型的缺陷 引入交互作用项常常能够改善模型 剔除：异常数据 有助于得到更好的结果 另： 可以直接对6种管理—教育组合引入5个0-1变量 模型三 酶促反应 酶促反应 由酶作为催化剂催化进行的化学反应 生物体内的化学反应绝大多数属于酶促反应 酶促反应中酶作为高效催化剂使得反应以极快的速度（103~1017倍）或在一般情况下无法反应的条件下进行 酶是生物体内进行各种化学反应最重要的因素 酶促反应动力学 问题： 研究酶促反应中——嘌呤霉素——对反应速度与底物（反应物）浓度之间关系的影响 方案 建立数学模型，反映该酶促反应的速度与底物浓度以及经嘌呤霉素处理与否之间的关系 设计了两个实验 酶经过嘌呤霉素处理 酶未经嘌呤霉素处理 实验数据: 分析 酶促反应的基本性质 底物浓度较小时，反应速度大致与浓度成正比； 底物浓度很大、渐进饱和时，反应速度趋于固定值 数据分析 实验数据：散点图 线性化模型 经嘌呤霉素处理后实验数据的估计结果 结果分析 mei4.m Matlab 统计分析 beta的置信区间 Matlab 程序 function y=f1(beta, x) y=beta(1)*x./(beta(2)+x); 结果分析 最终反应速度为 半速度点(达到最终速度一半时的x值 )为 效果 其它输出 给出交互画面 拖动画面的十字线，得y的预测值和预测区间 画面左下方的Export 输出其它统计结果。 剩余标准差s= 10.9337 混合反应模型 在同一模型中考虑嘌呤霉素处理的影响 Matlab 程序 用nlinfit 和 nlintool命令 参数初值：基于对数据的分析 简化混合模型 估计结果和预测 结果分析 一般混合模型与简化混合模型预测比较 评注 酶促反应 END y ~酶促反应速度 x ~底物浓度 x y 0 ?1 经嘌呤霉素处理 x y 未经嘌呤霉素处理 x y 模型 mei.m mei1.m 0.04841 195.8027 估计值 β2=?2 /?1 β1=1/?1 参数 R2=0.8557