数学建模2统计模型.docx

1、数学建模2统计模型数学建模论文题目： 一个医药公司的新药研究部门为了掌握一种新止痛剂的疗效，设计了一个药物试验，给患有同种疾病的病人使用这种新止痛剂的以下4个剂量中的某一个：2 g，5 g，7 g和10 g，并记录每个病人病痛明显减轻的时间（以分钟计）. 为了解新药的疗效与病人性别和血压有什么关系，试验过程中研究人员把病人按性别及血压的低、中、高三档平均分配来进行测试. 通过比较每个病人血压的历史数据，从低到高分成3组，分别记作0.25，0.50和0.75. 实验结束后，公司的记录结果见下表（性别以0表示女，1表示男）.请你为该公司建立一个数学模型，根据病人用药的剂量、性别和血压组别，预测出服

2、药后病痛明显减轻的时间.病人序号病痛减轻时间min用药剂量g性别血压组别135200.25243200.50355200.75447210.25543210.50657210.75726500.25827500.50928500.751029510.251122510.501229510.751319700.251411700.501514700.751623710.251720710.501822710.7519131000.252081000.502131000.7522271010.2523261010.502451010.751、摘要在农某医药公司为了掌握一种新止痛药的疗效，设计了一个

3、药物实验，通过观测病人性别、血压和用药剂量与病痛时间的关系，预测服药后病痛明显减轻的时间。我们运用数学统计工具minitab软件，对用药剂量，性别和血压组别与病痛减轻时间之间的数据进行深层次地处理并加以讨论概率值P（是否0.05）和拟合度R-Sq的值是否更大（越大，说明模型越好）。首先，假设用药剂量、性别和血压组别与病痛减轻时间之间具有线性关系，我们建立了模型。对模型用minitab软件进行回归分析，结果偏差较大，说明不是单纯的线性关系，然后对不同性别分开讨论，增加血压和用药剂量的交叉项，我们在模型的基础上建立了模型，用minitab软件进行回归分析后，用药剂量对病痛减轻时间不显著，于是我们有

4、引进了用药剂量的平方项，改进模型建立了模型，用minitab软件进行回归分析后，结果合理。最终确定了女性病人服药后病痛减轻时间与用药剂量、性别和血压组别的关系模型： Y=31.8-3.49+56.1-9.32+0.26 对模型和模型关于男性病人用minitab软件进行回归分析，结果偏差依然较大，于是改进模型建立了模型，用minitab软件进行回归分析后，结果合理。最终确定了男性病人服药后病痛减轻时间与用药剂量、性别和血压组别的关系模型：Y=32.8-4.02+0.955+0.0.0427 关键词 止痛剂 药剂量 性别 病痛减轻时间2、问题的提出 一个医药公司的新药研究部门为了掌握一种新止痛剂的

5、疗效，设计了一个药物实验，给患有同种病痛的病人使用这种新止痛剂的一下4个剂量中的某一个：2g，5g，7g和10g，并记录每个病人病痛明显减轻的时间（以分钟计）。为了了解新药的疗效与病人性别和血压有什么关系，实验过程中研究人员把病人按性别及血压的低、中、高三档平均分配来进行测试。通过比较给个病人血压的历史数据，从低到高分成三组，分别记作0.25，0.50和0.75.实验结束后，公司的记录结果附录1-1表（性别以0表示，1表示男）。 现在为公司建立一个模型，根据病人用药的剂量、性别和血组别，预测出服药后病痛明显减轻的时间。3、问题的分析 假定每个患该种病的程度相差不大，即病情基本相同，根据现实，用

6、药量与病痛减轻时间会有一定的关系，一般，药用量越高，病痛减轻时间变得越快；而更一般，男性身体素质相对于女性来说比较强壮，病痛减轻的时间也会跟性别有关系，正常而言，身体素质越好，病痛减轻时间越快；另一个，一个人的血压组别的高地也会影响到他的病痛减轻时间的快慢。对1-1表格中的数据进行相关分析如下：相关分析: 用药剂量（g）, 血压组别，知用药剂量（g）和 血压组别 的Pearson相关系数 = 0.000 P 值 = 1.000；由此，可以看出用药剂量与血压组别没有关系，如图1-1所示 1-1图相关分析: 用药剂量（g）, 性别，知用药剂量（g） 和 性别 的Pearson相关系数=0.000

7、P值 = 1.000；由此可以看出用药剂量与性别相互独立。如1-2图所示 1-2图根据所给数据可分别作出病痛减轻时间与用药剂血压组别的散点图量，性别及如下：1-3.1图1-3.2图 1-3.3图4、模型假设与符号假设假设病痛减轻时间只与用药剂量、性别和血压组别有关,不受其他因素的影响，由以上散点图（图1-3.1-图1-3.3）可以作出如下模型假设模型： 符号说明 1、为病痛减轻时间量，单位（min）； 2、表示用药剂量 单位（g）； 3、表示性别 ； 4、 表示血压组别； 5、 S 表示标准差； 6、 R-Sq 表示线性拟合度。5、模型的建立 下面用minitab软件对分别对残差对用药剂量、残

8、差对性别和残差对血压组别进行绘图，到出对应的1-4.1图、1-4.2图和1-4.3图，并对这些图进行分析，分别可以看出残差对用药剂量是正常的、残差对性别是正常的、残差对血压组别正常的。1-4.1图1-4.2图1-4.3图由1-4.11-4.3图分析，可以用药剂量和血压组别的乘积表示对病痛减轻时间的交互式影响，性别对病疼减轻时间有显著影响，因此可以对男性和女性分开讨论，得到如下模型：模型 （1）对女性的进行分析如下：回归分析: 病痛减轻时间（min） 与 用药剂量（g）, 血压组别, 用药剂量及血压组别 回归方程为病痛减轻时间（min） = 23.1 + 0.040 用药剂量（g） + 59.4

9、 血压组 别 - 10.2 用药剂量及血压组别交叉项即 Y=23.1+0.040+59.4-10.2自变量 系数 系数标准误 T P常量 23.096 6.108 3.78 0.005用药剂量（g） 0.0397 0.9767 0.04 0.969血压组别 59.38 11.84 5.02 0.001用药剂量及血压组别 -10.163 2.021 -5.03 0.001S = 3.37051 R-Sq = 96.5% R-Sq（调整） = 95.2%方差分析来源 自由度 SS MS F P回归 3 2486.03 828.68 72.94 0.000残差误差 8 90.88 11.36合计 1

10、1 2576.92来源 自由度 Seq SS用药剂量（g） 1 2184.16血压组别 1 14.52用药剂量及血压组别 1 287.36异常观测值 用药剂 病痛减轻时 拟合值 标准化观测值 量（g） 间（min） 拟合值 标准误 残差 残差 8 7.0 11.000 17.495 1.081 -6.495 -2.03RR 表示此观测值含有大的标准化残差因为用药剂量p值为0.969，所以对病痛减轻时间影响不显著， 不妨引进用药剂量的平方项加以讨论，因此模型进一步改进为:模型 回归分析: 病痛减轻时间（min） 与 用药剂量（g）, 血压组别, 用药剂量及血压组别, 用药剂量的平方 回归方程为：

11、病痛减轻时间（min） = 31.8 - 3.49 用药剂量（g） + 56.1 血压组别 - 9.32 用药剂量及血压组别 + 0.264 用药剂 量的平方即 Y=31.8-3.49+56.1-9.32+0.26自变量 系数 系数标准误 T P常量 31.779 5.755 5.52 0.001用药剂量（g） -3.494 1.558 -2.24 0.060血压组别 56.122 9.141 6.14 0.000用药剂量及血压组别 -9.322 1.579 -5.90 0.001用药剂量的平方 0.2636 0.1020 2.58 0.036S = 2.57789 R-Sq = 98.2%

12、R-Sq（调整） = 97.2%方差分析来源 自由度 SS MS F P回归 4 2530.40 632.60 95.19 0.000残差误差 7 46.52 6.65合计 11 2576.92来源 自由度 Seq SS用药剂量（g） 1 2184.16血压组别 1 14.52用药剂量及血压组别 1 287.36用药剂量的平方 1 44.36由拟合值R-Sq = 98.2%可以确定，该模型比较合理。（2）、对男性用模型进行分析，分析结果如下:回归分析: 病痛减轻时间（min） 与 用药剂量（g）, 血压组别, 用药剂量及血压组别 回归方程为：病痛减轻时间（min） = 31.5 + 0.16

13、用药剂量（g） + 39.0 血压组别 - 7.59 用药剂量及血压组别即 Y=31.5+0.16+39.0-7.59 系数标自变量 系数 准误 T P常量 31.48 13.71 2.30 0.051用药剂量（g） 0.157 2.055 0.08 0.941血压组别 39.03 25.39 1.54 0.163用药剂量及血压组别 -7.588 3.806 -1.99 0.081S = 7.84538 R-Sq = 76.6% R-Sq（调整） = 67.9%方差分析来源 自由度 SS MS F P回归 3 1615.27 538.42 8.75 0.007残差误差 8 492.40 61.

14、55合计 11 2107.67来源 自由度 Seq SS用药剂量（g） 1 1349.42血压组别 1 21.13用药剂量及血压组别 1 244.72因为用药剂量p值为0.941，所以对病痛减轻时间影响不显著， 不妨引进用药剂量的平方项加以讨论，因此可以利用模型进行分析： 回归分析: 病痛减轻时间（min） 与 用药剂量（g）, 血压组别, 用药剂量及血压组别, 用药剂量的平方 回归方程为：病痛减轻时间（min） = 49.8 - 7.84 用药剂量（g） + 39.0 血压组别 - 7.59 用药剂量及血压组别 + 0.667 用药剂量的平方即 Y=49.8-7.84+39.0-7.59+0

15、.667自变量 系数 系数标准误 T P常量 49.81 10.71 4.65 0.002用药剂量（g） -7.843 2.784 -2.82 0.026血压组别 39.03 16.96 2.30 0.055用药剂量及血压组别 -7.588 2.543 -2.98 0.020用药剂量的平方 0.6667 0.2018 3.30 0.013S = 5.24268 R-Sq = 90.9% R-Sq（调整） = 85.7%方差分析来源 自由度 SS MS F P回归 4 1915.27 478.82 17.42 0.001残差误差 7 192.40 27.49合计 11 2107.67来源 自由度

16、 Seq SS用药剂量（g） 1 1349.42血压组别 1 21.13用药剂量及血压组别 1 244.72用药剂量的平方 1 300.00由此，可以看出，在男性方面血压组别的P=0.55，对病痛减轻时间不显著，不妨取消血压组别这个单变量，将模型进一步改进。模型 回归分析: 病痛减轻时间（min） 与 用药剂量（g）, 性别, 用药剂量及血压组别, 用药剂量的平方 * 性别（实质上）是常量 * 性别 已从方程中删除。回归方程为：病痛减轻时间（min） = 32.8 - 4.02 用药剂量（g） + 0.955 用药剂量及血压组别 + 0.00427 用药剂量的平方 Y=32.8-4.02+0.

17、955+0.0.0427自变量 系数 系数标准误 T P常量 32.794 3.437 9.54 0.000用药剂量（g） -4.0229 0.5371 -7.49 0.000用药剂量及血压组别 0.9549 0.8562 1.12 0.297用药剂量的平方 0.0042660 0.0006143 6.94 0.000S = 3.36837 R-Sq = 95.7% R-Sq（调整） = 94.1%方差分析来源 自由度 SS MS F P回归 3 2016.90 672.30 59.25 0.000残差误差 8 90.77 11.35合计 11 2107.67来源 自由度 Seq SS用药剂量

18、（g） 1 1349.42用药剂量及血压组别 1 120.36用药剂量的平方 1 547.12异常观测值 用药剂 病痛减轻时 拟合值 标准化观测值 量（g） 间（min） 拟合值 标准误 残差 残差 12 10.0 5.000 0.793 2.691 4.207 2.08RR 表示此观测值含有大的标准化残差* 注 * 列中的所有值相同。 用药剂量及血压组别的P= 0.297，但是 R-Sq = 95.7% R-Sq（调整） = 94.1%，说明这个模型改进更加合理。6、模型的优缺点与改进方向 通过回归模型的建立及不断改进过程当中，得知该公司的新药的疗效对于男性和女性的作用程度不一样。该模型是针对该公司的新药进行建模，不具有普遍性。7、参考文献1、姜启源，谢金星，叶 俊. 数学模型（第三版）. 高等教育出版社，2003.8（2012重印）2、马 林，何 桢.六西格玛管理 （第二版）.中国人民大学出版社，2007.7（2011.11重印） 3、吴 翊，李永乐，胡庆军.应用数理统计.国防科技大学出版社，1995.8（2010.7重印）8、附录部分