数学建模医院投入产出

陕西医院投入产出效率评估摘要随着人民群众健康保健需求的日益增长，为了实现效益的最大化，一些医院投入大量资金，增加病床、扩建医院。但这是否能够提高医院的整体效率仍然是一个问题。因此需要建立数学模型客观评价、分析医院的效率和运行特征，给出不同规模下医院的投入产出预测模型和最优投入产出分配方案，为政府相关部门制定医疗资源分配政策提供参考。对于问题一，首先运用聚类分析法借助SPSS软件根据相关系数这一指标对具有相同量纲的投入值进行分类，剔除错误数据。随后建立DEA的CCR效率评价模型，利用MAXDEA软件求解出每个医院每年对应的效率值，最后将四年效率值进行平均求值，得出每个医院总的效率值。对于问题二，在问题一的基础上，借助RTS系数得出在CCR模型中所有最优解的情况进而判断DMU所处的状态，即,说明该DMU处于IRS状态；，说明该DMU处于DRS状态；，说明该DMU处于CRS状态。同时利用投导向入包络模型直接求解出上规模弹性和下规模弹性。根据的散点图分析得出医院的运行特征。对于问题三，在问题一和问题二的基础上，以规模效益将医院分为收益递增、收益不变、收益递减三类。利用近4年各项投入的平均增长率预测下一年的各项投入，同时利用SPSS软件对产出建立逐步回归模型。在各类医院中选取具有代表性的医院进行投入产出预测。最后，评价模型的优缺点，并对模型作出进一步改进与推广，指出DEA模型中的不足之处，以及DEA模型在其他领域的应用。关键字:DEASPSS软件逐步回归RTS系数CCR效率评价模型1一、问题重述随着人民群众健康保健需求的日益增长，医院作为从事医疗服务的独立经济实体，在满足社会对医疗保健事业的需求中获得经营效益的同时，也要促进自身的发展。因此医院必须按照社会对医疗保健事业的需求变化，有计划地安排和投入自身的人力资源和物质资源。但是盲目的扩容增建不一定会提高医院的整体效率，有时甚至会起到相反的作用，增大政府对公立医院总体投入负担、浪费宝贵的社会资源。故在某种程度上能客观评价、分析医院的现状，不仅能为医院的管理层制定相关经营策略提供循证依据，也能够为政府相关部门制定医疗资源分配政策提供参考。问题一，结合附件给定的700余家较大医院1997-2000年的投入产出数据，通过对相关投入产出指标进行分析，建立医院的效率评价模型，并对附件中所涉医院的工作效率进行评估。根据附件中数据及问题一中给出的效率评价模型，对陕西省700余家医院进行规模分类，并在其基础上分析不同规模下医院的运行特征。在问题一和问题二的基础上，建立不同规模下医院的投入产出预测模型，并给出医院发展的最优投入产出方案，辅助行政决策部门进行财政调配。同时，试从不同规模医院中各选取一个，给出具体的投入产出方案。二、模型假设1题中所给数据为可信数据2其他医疗人员指除医师，护士外的其他工作人员3假设医院可以接收所有患者三、符号约定符号代表含义DMU决策单元第j个DMU对第i种输入的投入量，第j个DMU对第r种输出的投入量，第i种输入的一种度量第r种输出的一种度量产出投入比线性组合系数2四、模型的分析、建立与求解4.1问题一4.1.1问题一分析在1997~2000年700余家医院的投入产出数据中，剔除错误数据，根据各个医院近几年发展的情况补充不完整的数据。汇总每年各个医院的投入产出数据，鉴于数据信息量的庞大，运用SPSS将投入量进行聚类，使具有相同量纲的数据可以叠加，减少数据类别，运用CCR模型，根据效率=分别求解出每一年对应的效率值。将每个医院四年的所有效率值加起来平均即可得出各个医院的效率值。最后进行大小比较对每个医院的效率值进行评估。4.1.2问题一模型建立：聚类分析：聚类分析即就是根据事物的相似程度进行分类。首先定义一个度量事物相似程度的指标，这个指标称为聚类统计量，然后按照相似程度的大小，把对象逐一归类，关系密切的聚集到一个小的分类单位，关系疏远的聚合到一个大的分类单位里，直到所有的对象都归类完毕，把不同类型都一一划分出来，形成一个亲戚关系谱系图。谱系图直观地显示分类对象的差异和联系。假设有n样品，每个样品观测了p个指标，用表示第i个样品的第j个指标，则有表1所列数据阵。表1样本数据库样品指标均值对变量进行分类的方法称为R型聚类法，所有的聚类统计量称为相似系数。在此问题中，采用相关系数这一相似系数，如下：记变量：，则可以用两变量的相关系数作为它们的相似系数。3利用系统聚类法由粗到细的指出多种分类情况，最后的系统聚类结果由一个谱系图展示出来。具体的基本步骤如下：设样品集为S，则（1）计算n个样品两两之间的距离,记为矩阵D；（2）首先构造n个类，每一个类中只包含一个样品，每一个类的平台高度均为零；（3）合并距离最近的两类为新类，并且以这两类间的聚合指数作为谱系图中的平台高度；（4）计算新类与当前类的距离，若类的个数已经等于1，转入下一个步骤；否则，回到步骤3；（5）画谱系图；（6）决定类的个数和类。多指标评价的DEA模型CCR设有n个决策单元，每个决策单元都有m种“输入”（表示该决策单元对“资源”的耗费），以及s种“输出”（表示该决策单元消耗了“资源”之后，表明“成效”的数量），由图1给出。其中12jn——第j个DMU对第i种输入的投入量，——第j个DMU对第r种输出的投入量，——第i种输入的一种度量——第r种输出的一种度量i1,2，，m；j1,2，，n;r1,2，，s,为已知数据，可以由题目附件得到；v=和u=为变量。对于第j个决策单元有相应的效率评价指数：，j1，2，我们总可以适当地选取权系数v和u，使得满足：，2，现在，对第个决策单元进行相对效率评价：以权系数v和u为变量；以第个决策单元的效率指数为目标；以所有的决策单元的效率为约束，即，2，于是构成了评价的最优化模型：1m4其中，v,有表示对于r=1，2,，s，0,并且存在有。根据公式计算产出投入比：（4—1—1）4.1.3问题一的模型求解：运用SPSS软件根据相关系数这一指标对投入值进行分类，得到如下所示的树状图：图1平均联接的树状图根据上述树状图所得的数据分类，对相关系数接近的数据进行加和，利用上述公式一计算出各个公司每年的效率值，最后将每个公司四年的效率取平均值，得到公司整个的效率值。（见附录表1）4.2问题二4.2.1问题二的分析：根据生产技术的收益要先后经历规模收益递增、规模收益不变和规模收益递减三个阶段，我们通过RTS的判断方法得出在CCR模型中所有最优解的情况，进而判断DMU所处的状态。DMU规模收益状态的判断方法：（1）在CCR模型中的所有最优解中，,则说明该DMU处于IRS状态。（2）如果在CCR模型的所有最优解中，，则说明该DMU处于DRS状态。（3）在CCR模型的所有最优解中，只要其中一个解，则说明该DMU处于CRS状态。4.2.2模型建立与求解在实际求解时，求解DEA模型的最优解并非易事，确定以上条件并不需要求解所有最优解。在实际计算中，只需要在求解CCR模型之后，再额外求解两个模型，分别5获得的最大值和最小值即可。以投入导向的模型为例，求解最大值（最小值）的线性规划为：max(min)s.t.，，λi=1，2，；j（4—2—1）在上述规划式中，如果max,则说明在所有最优解中均小于1；如果min均大于1。规模弹性：规模弹性就是对规模收益状态的量化，各项投入指标的等比例变动会引起各项产出指标的等比例变动，规模弹性是指产出指标的等比例变动的比例与投入指标等比例变动的比例之比值。用公式表示为：E通过投导向入包络模型直接求解规模弹性：通过投入导向模型计算规模弹性的原理与利用产出导向模型求解得原理相似，区别在于，通过投入导向模型，是先等比例微量增加（或减少）产出，然后求解投入增加（或减少）的比例。计算步骤简述如下：1）求解模型（2），计算所有DMU目标值。mins.t.=i=1，2，（4—2—2）投影点的投入和产出向量分别记为：然后，在DEA模型的数据中，用目标值替换原始数据值。2）计算上方的规模弹性。求解模型（4—2—3），获得在产出微量等比例（）增加后，投入增加的比。mins.t.；j=1，2，，n上方的规模弹性为=（4—2—3）3）求解模型（4—2—4），获得在产出微量等比例减少之后，投入减少的比例。maxs.t.6；j=1，2，，n下方的规模弹性为：=（4—2—4）根据CCR模型运用公式1，解出每一年各个公司的值，根据之间的关系为划分指标，即=1为同一规模的公司，1为同一规模的公司，1为同一规模的公司。得到其散点图如下图所示：图2散点图同样。由公式（2）、（3）、（4）可以求出其上方规模弹性和下方规模弹性，由上述所求结果可以得出：对700余家医院将其规模可以分为三类，即收益递增（IRS）、规模收益不变（CRS）和规模收益递减（DRS）三种规模。通过图表分析，处于规模收益递减和规模收益不变的医院占据大多数，仅有少数医院处于规模收益递增状态。4.3问题三4.3.1问题分析根据问题一和问题二建立效益规模分类，利用年平均增长率，逐步回归分析法确立回归方程，同时给出下一年投入产出预测方案。4.3.2问题三模型建立采用多元线性回归分析，对收益递增型、收益递减型、收益不变型这三种规模的住院人数、门诊病房、出院病人数进行回归分析，得到线性回归方程，便于投入产出的预测。多元线性回归：多个自变量和一个因变量的关系多元线性回归一般模型：Y（4—3—1）Y：因变量，，，：自变量7，，，：是未知参数，称为回归系数：随机因素的影响，通常假定服从期望值为零方差为的正态分布。将n组观测数据代入（4—3—1）式，可得到样本回归模型：矩阵形式：Y=X其中：Y=，X=，，（1）解释变量是确定性变量，不是随机变量，且要求矩阵X中列不相关,样本容量个数应大于解释变量的个数。（2）随机误差项具有零均值和同方差，即Ei1,2，，ncov（，）=i,j=1，2，（3）正态分布的假设条件：~N（0，），i=1，2，参数估计：使得估计值与观测值y之间的残差在所有样本上达到最小。求，，，,使得==min由多元函数求极值点的方法可求得回归系数的最小二乘估计值为：未知参数的一个无偏估计为：==显著性t检验：1）提出假设：：=0；：0（i=1，2，，p）2）t检验的计算公式3）给定显著性水平，确定临界值(np)4）若(np)，则拒绝，接受备择假设，即总体回归系数0。下图是多元线性回归变量的定义表：表2多元线性回归变量定义表变量类别变量名称定义变量名称定义变量名称定义因变量住院病人数出院病人数门诊治疗数Y3自变量可变投1—其他可变投1—其他可变投1-其他可变投2—保养维护可变投2—保养维护可变投2-保养维护可变投3—基础维护可变投3—基础维护可变投3-基础维护可变投4—设备周边补给可变投4—设备周边补给可变投4-设备周边补给可变投5—设备可变投5—设备可变投5-设备可变投6—药品可变投6—药品可变投6-药品准固投1—其他医疗人员数准固投1—其他医疗人员数准固投1-其他医疗人员数准固投2—动能消耗准固投2—动能消耗准固投2-动能消耗8准固投3—学术培养准固投3—学术培养准固投3-学术培养准固投4—病床数准固投4—病床数准固投4-病床数准固投5—护士数准固投5—护士数准固投5-护士数准固投6—医师数准固投6—医师数准固投6—医师数准固投7—办公消耗准固投7—办公消耗准固投7-办公消耗4.3.3模型的求解以下是对收益递增型的住院病人数，出院病人数，门诊治疗，给的回归分析：住院病人数：通过逐步回归法对变量进行筛选，对选入的变量进行F检验和t检验进行剔除，计算回归预测方程的回归系数R方，R方为0.597，说明拟合效果较好。具体计算在SPSS统计软件上进行，具体结果见表由系数写回归方程：=24088.293+724.5220.593+12345.322+1093.612+491.179（4—3—2）‘表三各类系数相关性住院病人数住院病人数与其他相关性Person1.000可变投1-其他—0.031可变投2-保养维护0.705可变投3-基础维护0.478可变投4-设备