数学建模用回归分析探究止痛剂疗效模型

用回归分析探究止痛剂疗效模型摘要：本文通过结合实际情况，对病人疼痛减轻时间和影响它的因素：用药量，性别，血压作了较为合理的分析。首先，我们将用药量分为四组，求其每组用药量的平均值，利用Minitab工具，画出散点图。我们通过散点图分析，得到了病痛减轻时间和平均用药量的数学模型：二次模型。然后，我们假设病人疼痛减轻时间分别与性别，血压呈线性关系。综合上述，我们就第一次得到了较为粗略的数学模型。结合数据，通过MATLAB回归分析检验，我们得到了可以接受的数学模型。同时，我们也发现其存在问题，需要优化，所以，我们提出对模型的改进：同样，将用药量分为四组，将血压分为三组，分别求其平均值。考虑不同性别的病痛患者疼痛减轻时间与平均用药量的关系，不同性别的病痛患者疼痛减轻时间与平均用药量和平均血压共同作用的关系；不同性别的病痛患者病痛减轻时间与平均用药量和不同性别的共同作用的关系。同理，我们利用Minitab工具，画出不同性别患者病痛减轻时间和分好组的平均用药量的散点图；不同性别病痛患者病痛减轻时间和分好组的平均用药量和平均血压共同作用的散点图；不同性别病痛患者病痛减轻时间和分好组的平均用药量和性别共同作用的散点图。这样我们就得到了个改进后的二次模型。我们通过MATLAB工具，经过回归分析检验我们发现改进后的数学模型较以前的模型更加优化，更加接近实际。因此，我们接受改进后的模型。这样我们就得到了病人疼痛减轻时间与用药量和性别以及血压间的关系，即最终数学模型。此数学模型可以作为我们根据患者的用药量，性别，血压来预测其病痛减轻时间。关键字：回归方程矩阵置信区间估计值问题重述一个医药公司的新药研究部门为了掌握一种新止痛药的疗效。设计了一种药物试验，给患者有同种病痛的病人使用这种新止痛剂的一下4个剂量中的某一个：2g，5g，7g和10g，并记录每个病人病痛明显减轻的时间（以分钟计）。为了解新药的疗效与病人性别和血压又什么关系，实验过程中研究人员把病人按性别及血压的低、中、高三档平均分配来进行测试。通过比较每个病人血压的历史数据，从低到高分成三组，分别记作0.25,0.50,和0.75.实验结束后，公司的记录结果见下表（性别以0表示女，1表示男）。请你为公司建立一个模型，根据病人用药的剂量、性别和血压组别，预测出服药后病痛明显减轻的时间。问题分析及建立模型根据常识，我们知道病人的病痛减轻时间与用药量有直接的关系，除此之外，一般的话，病人病痛减轻时间与性别，血压也有明显的关系。我们假设病痛减轻时间为y,用药剂量为1x,性别为2x，血压组别为3x。分析不同用药剂量下的病痛减轻时间，根据所给数据用minitab作散点图：(分别作出线性、二次、三次关系拟合)1109876543216050403020100用药剂量（g）病痛减轻时间（min）S7.94724R-Sq71.5%R-Sq（调整）70.2%拟合线图病痛减轻时间（min）=51.16-4.137用药剂量（g）109876543216050403020100用药剂量（g）病痛减轻时间（min）S7.19838R-Sq78.8%R-Sq（调整）75.6%拟合线图病痛减轻时间（min）=63.47-8.98用药剂量（g）+0.261用药剂量（g）**2+0.0139用药剂量（g）**32109876543216050403020100用药剂量（g）病痛减轻时间（min）S7.19838R-Sq78.8%R-Sq（调整）75.6%拟合线图病痛减轻时间（min）=63.47-8.98用药剂量（g）+0.261用药剂量（g）**2+0.0139用药剂量（g）**3以上三幅图的比较中，我们发现，随着用药剂量的增加，y有向上弯曲增加的趋势，图中的曲线是用二次函数模型可知：病痛减轻时间与用药剂量成二次函数关系，则可建立模型：y=0+11x+212x+分别作病痛减轻时间与性别和血压组别的散点图：31.00.80.60.40.20.06050403020100性别病痛减轻时间（min）病痛减轻时间（min）与性别的散点图0.80.70.60.50.40.30.26050403020100血压组别病痛减轻时间（min）病痛减轻时间（min）与血压组别的散点图在病痛减轻时间与性别的散点图中，性别的不同并未导致病痛减轻时间很大的变化，而在病痛减轻时间与血压组别的散点图中，血压分为三个组别，不易看出两者的关系。于是我们假设病人疼痛减轻时间分别与性别、血压呈线性关系，可建立模型：210xy及310xy4由以上分析可建立综合模型：2143322110xxxxy（1）模型求解直接利用MATLAB统计工具箱中的命令regress求解：x=1.00002.000000.25004.00001.00002.000000.50004.00001.00002.000000.75004.00001.00002.00001.00000.25004.00001.00002.00001.00000.50004.00001.00002.00001.00000.75004.00001.00005.000000.250025.00001.00005.000000.500025.00001.00005.000000.750025.00001.00005.00001.00000.250025.00001.00005.00001.00000.500025.00001.00005.00001.00000.750025.00001.00007.000000.250049.00001.00007.000000.500049.00001.00007.000000.750049.00001.00007.00001.00000.250049.00001.00007.00001.00000.500049.00001.00007.00001.00000.750049.00001.000010.000000.2500100.00001.000010.000000.5000100.00001.000010.000000.7500100.00001.000010.00001.00000.2500100.00001.000010.00001.00000.5000100.00001.000010.00001.00000.7500100.0000[b,bint,r,rint,stats]=regress(y,x)b=63.1291-10.27065.6667-1.50000.51115bint=48.717377.5409-14.9243-5.6169-0.021311.3546-15.432512.43250.13190.8903stats=0.827522.79030.000044.310b表示的各个参数估计值，bint表示其各个参数估计值的置信区间。Stats下的各个数据依次表示回归方程的决定系数，即相关系数，F统计量值，F统计量对应的概率值p。结果分析：从上表我们可以看出，2R=0.8275指因变量y的82.75﹪可由模型确定。F=22.7903远远超过F检验的临界值。P=0.0000远远小于a(置信水平)，所以，该模型大体上还是可以应用的。上表的回归系数给出了的估计值，即0=63.1291，1=-10.2706，2=5.6667，3=-1.5000，4=0.5111.根据其各自的置信区间，我们可以检查发现：2、3的置信区间包含了零点。它表明了回归变量2x，3x还不是太显著。病痛减轻时间预测：将回归系数值代入简化后的模型。即可预测病痛减轻时间y，预测值记为y，得到预测方程：y63.1291-10.27061x+5.66672x-1.50003x+0.511121x我们只需知道用药量1x，性别2x，血压3x，我们就能计算预测值y。如对男性，血压组别为0.5，用药剂量为5时的疼痛时间进行预测y63.1291-10.2706*5+5.6667*1-1.5000*0.5+0.5111*25=29.47（四舍五入）如对女性，血压组别为0.5，用药剂量为5时的疼痛时间进行预测y63.1291-10.2706*5+5.6667*0-1.5000*0.5+0.5111*25=23.80（四舍五入）模型改进由于的置信区间包含了零点，这就表明了此模型还有待改进。模型（1）中6分析回归变量1x、2x、3x对y的影响都是相互独立的。下面讨论自变量之间的相互作用对y的影响：为探究用药剂量与性别的共同作用对疼痛减轻时间的影响，作不同性别时，用药剂量与疼痛减轻时间的关系图：109876543215040302010C1Y数据C2C3变量C2,C3与C1的散点图其中，虚线表示男性用药剂量与疼痛减轻时间分析图，实线表示女性用药剂量与疼痛减轻时间分析图。由图可以明显看出来性别与用药剂量的相互作用对疼痛减轻时间有影响，可用1x与2x的乘积代表他们的相互作用。为探究用药剂量与血压组别的共同作用对疼痛减轻时间的影响，作不同血压组别时，用药剂量与疼痛减轻时间的关系图：7109876543216050403020100C1Y数据C3C6C9变量20.021.027.541.017.015.524.543.04.018.028.556.0C3,C6,C9与C1的散点图其中，C3、C6、C9分别是血压组别为0.25、0.50、0.75时用药剂量与疼痛减轻时间的关系。由图可以看出来用药剂量与血压组别的相互作用对疼痛减轻时间有影响，可用1x与3x的乘积代表他们的相互作用。由表中所给数据可看出性别与血压组别的相互作用对疼痛减轻时间的影响不大。于是将模型（1）中增加两项1x3x和1x2x，得到模型：3162152143322110xxxxxxxxy（2）利用MATLAB的统计工具得到如下结果：其中矩阵x的第六列表示1x2x，第七列表示1x3x。x=1.00002.000000.25004.000000.50001.00002.000000.50004.000001.00001.00002.000000.75004.000001.50001.00002.00001.00000.25004.00002.00000.50001.00002.00001.00000.50004.00002.00001.00001.00002.00001.00000.75004.00002.00001.50001.00005.000000.250025.000001.25001.00005.000000.500025.000002.50001.00005.000000.750025.000003.75001.00005.00001.00000.250025.00005.00001.250081.00005.00001.00000.500025.00005.00002.50001.00005.00001.00000.750025.00005.00003.75001.00007.000000.250049.000001.75001.00007.000000.500049.000003.50001.00007.000000.750049.000005.25001.00007.00001.00000.250049.00007.00001.75001.00007.00001.00000.500049.00007.00003.50001.00007.00001.00000.750049.00007.00005.25001.000010.000000.2500100.000002.50001.000010.000000.5000100.000005.00001.000010.000000.7500100.000007.50001.000010.00001.00000.2500100.000010.00002.50001.000010.00001