数学建模新止痛剂生效时间预测模型.doc
《数学建模新止痛剂生效时间预测模型.doc》由会员分享,可在线阅读,更多相关《数学建模新止痛剂生效时间预测模型.doc(20页珍藏版)》请在咨信网上搜索。
1、数学建模论文题目:新止痛剂生效时间预测模型 班级:计算机科学与技术1401班 学号:201409824 姓名:马元凯 学号:201409822 姓名:李炳毅 日期:2016年12月15日目录摘要1一、 问题描述1二、 问题分析2三、 模型假设3四、 模型建立3五、 模型求解6六、 模型分析6七、 模型改进9八、 模型评注与推广11九、 参考文献12十、 附录12摘要某医药公司的新药研究部门为了掌握一种新止痛剂的疗效,进行了相应的研究,根据医药公司给出的相关实验数据,通过MATLAB软件进行数据处理工作,根据拟合曲线和散点图建立基本模型(模型(1).用MATLAB软件解出该基本模型,再利用残差图
2、两次剔除数据并且回归,得到最佳模型(模型(2),该模型的拟合度为85.14%.拟合度不是很高,针对不同的性别,引入独立变量的交互作用,对模型进行改进,得到两个拟合度较高的模型(4)和(5),针对男性建立的模型,拟合度为90.87%,针对女性建立的模型,拟合度为97.72%.本文给出的模型较为科学,若实验数据真实可信,则可以进行药物的推广.关键字:MATLAB软件,残差图,拟合度,统计回归模型一、 问题描述 一个医药公司的新药研究部门为了掌握一种新止痛剂的疗效,设计了一个药物试验,给患有同种疾病的病人使用这种新止痛剂的以下4个剂量中的某一个:2 g,5 g,7 g和10 g,并记录每个病人病痛明
3、显减轻的时间(以分钟计). 为了解新药的疗效与病人性别和血压有什么关系,试验过程中研究人员把病人按性别及血压的低、中、高三档平均分配来进行测试. 通过比较每个病人血压的历史数据,从低到高分成3组,分别记作0.25,0.50和0.75. 实验结束后,公司的记录结果见下表1(性别以0表示女,1表示男).请你为该公司建立一个数学模型,根据病人用药的剂量、性别和血压组别,预测出服药后病痛明显减轻的时间.病人序号病痛减轻时间min用药剂量g性别血压组别135200.25243200.50355200.75447210.25543210.50657210.75726500.25827500.5092850
4、0.751029510.251122510.501229510.751319700.251411700.501514700.751623710.251720710.501822710.7519131000.252081000.502131000.7522271010.2523261010.502451010.75表1二、 问题分析 一般来说,药物的疗效可以直观的用服药后病痛明显减轻的时间来衡量.在新药推广中,医药公司的新药研究部门设计了一种药物给患有同种疾病的病人使用后,根据病人的用药剂量、性别和血压组别,预测病痛减轻时间的多少来预测止痛药的疗效,这是一个统计回归问题,针对这个问题,我们需要给
5、出合理的假设,尤其是变量的选取,进而预测出药物的疗效.先将因变量分别与变量进行单独分析,得出两者间的大致函数关系,进一步整合这些关系,得出一个因变量与各个变量间的关系函数模型并进行求解,得出最终结论。 三、 模型假设1、 假设病人只服用了新型止痛药,未服用其它药物.2、 假设题中给出的实验数据真实可信,误差很小.3、 假设24名病人都是在服用新型止痛药的人群中随机选取的.4、 假设病人在实验阶段吃的食物对新型止痛药无影响.5、 假设模型中出现的符号含义如下表2所示.符号含义单位用药剂量性别女-0,男-1血压组别低-0.25,中-0.5,高-0.75病痛减轻时间minp概率值随机误差回归系数置信
6、水平表2四、 模型建立为了大致地分析y与,之间的关系,首先利用表1-1的数据分别作出y对,和的散点图(见图4-1,图4-2和图4-3的圆点).如图1为y对的散点图,图2为y对的散点图,图3为y对的散点图.图1图2图3由上图可知:y对可用二次函数拟合,拟合后如图4所示. 图4根据对散点图图1和图4的分析可得出y对的二次函数模型 ;根据图2可得出y对的线型模型;根据图3可得出y对的线型模型 ;结合上述的3个模型可建立如下(1)的多元线性回归模型: (1)上式右端的,称为回归变量(自变量),给定,时,病痛减轻时间y的平均值为;由表1的数据估计,影响y的其他因素作用都包含在随机误差中,如果模型选的合适
7、,应大致服从均值为0的正态分布.五、 模型求解直接利用MATLAB统计工具箱中的regress求解,进一步求出回归系数估计值及其置信区间(置信水平=0.05)、检验统计量R2,F,p,s2,详细数据见表3所示.参数参数估计值置信区间63.129148.7173,77.5409-10.2706-14.9243,-5.61695.6667-0.0213,11.3546-1.5000-15.4325,12.43250.51110.1319,0.8903R2=0.8275 F=22.7903 p=0.0000 s2=44.3109表3由表3得出的模型为:六、 模型分析表3显示,R2=0.8275 指因
8、变量y(病痛减轻时间)的82.75%可由模型确定,F值远远超过F检验的临界值,p远小于,因而模型(1)从整体来看是可用的.表3的回归系数给出了模型(1)中参数的估计值,的置信区间包含零点,因而对这两个系数的解释是不可靠的,所以需要残差分析.首次回归所得残差图如下图5所示,可以看出,第3个和第24个数据存在异常,剔除,进行第二次回归.图5第二次回归所得残差图如下图6所示,可以看出,第5个数据存在异常,剔除,进行第三次回归.图6第三次回归所得残差图如下图7所示.图7第三次回归的结果如表4所示.参数参数估计值参数置信区间55.812144.3918 67.2324-8.0962-12.0605 -4
9、.13207.13112.5539 11.7083-6.3868-18.5348 5.76120.40250.0869 0.7181R2= 0.8514 F= 22.9153 p= 0.0000 s2= 24.1056表4由图7可知数据没有异常项,因此模型基本可用,此时得出的最佳模型应为: (2)易知因变量y(病痛减轻时间)的85.14%可由模型确定.七、 模型改进从以上的分析可以看出,模型的拟合度最大为85.14%,拟合度不是很高,需要进行模型的改进;因为服药期间要考虑生理反应,而性别对生理反应有直观的影响,故改进的模型需要对男女分开进行讨论;模型(2)中回归变量和对因变量y的影响是相互独立
10、的.根据直觉和经验可以猜想和之间的相互作用会对y有影响,不妨简单地用和的乘积代表它们的交互作用,于是在模型(1)中增加一项,得到 (3)1、 针对男性,利用表1的数据估计模型(3)的系数,利用MATLAB得到如下表5的结果.参数参数估计值参数置信区间49.808824.4805 75.1372-7.8431-14.4259 -1.260439.0294-1.0850 79.1438-7.5882-13.6016 -1.57480.66670.1895 1.1438R2= 0.9087 F= 17.4206 p= 0.0010 s2= 27.4856 表5根据表5,可得出此时的最佳模型为: (4
11、)对表5的数据进行残差分析,得到如下图8的结果.图8从图8可以看出数据无异常项,则由此得出的模型(4)基本可信,病痛减轻时间的90.87%可由模型确定.2、 针对女性,利用表1的数据估计模型(3)的系数,利用MATLAB得到如下表6的结果.参数参数估计值参数置信区间36.939522.9221,50.9570-5.1686-8.8117,-1.525548.323526.1230,70.5240-7.4706-10.7986,-4.14260.35560.0915,0.6196R2=0.9772 F=74.8974 p=0.0000 s2= 8.4184 表6对表6的数据进行残差分析,得到如图
12、9所示的结果,图9从图9可以看出,第8项数据存在异常,但影响不大,所以不用剔除数据,模型基本可信.根据表5得出此时的最佳模型为: (5)病痛减轻时间的97.72%可由模型确定.通过以上分析可以确定,经过对模型的改进,针对不同的性别,得出的两种模型(4)和(5)的拟合程度都达到了百分之九十以上,达到了模型改进的目的.八、 模型评注与推广评注 从这个实例我们可以看出,建立回归模型要根据已知的数据,从常识和经验进行分析,辅以作图,决定最终的函数模型.用MATLAB软件求解后,作统计分析R2的值是对模型的直观评价,决定模型的拟合程度.本次建立的模型,优点在于,由简到繁,先单独考虑变量对结果的影响,再综
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 止痛剂 生效 时间 预测 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【精****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【精****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。