投入與產(chǎn)出效率評價(jià)

效率=）產(chǎn)出：outputOutput=f（ip,hp,op）ip=inpatient住院病人數(shù)hp=healthypatient出院病人數(shù)op=outpatient門診治療數(shù)采用包絡(luò)DEAHP模型對醫(yī)院的投入和產(chǎn)出指標(biāo)進(jìn)行效率分析DEA,即數(shù)據(jù)包絡(luò)分析，是一種常用的數(shù)據(jù)分析方法，主要用來評價(jià)決策單元的投入產(chǎn)出效率。醫(yī)院一共有3個(gè)產(chǎn)出變量，6個(gè)可變投入變量和7個(gè)準(zhǔn)固定投入變量。此時(shí)醫(yī)院為決策單元（DMU）常用DEA模型：假定存在n個(gè)同類型的DMU（決策單元）且每個(gè)DMU含有m個(gè)投入指標(biāo)和s個(gè)產(chǎn)出指標(biāo)。應(yīng)用傳統(tǒng)CCR來計(jì)算它們的DEA效率則決策單元的投入指標(biāo)數(shù)據(jù)可用矩陣X=x（（i=1，2,...，m；j=1，2,...，口）表示xxj表示第j個(gè)DMU的第i個(gè)投入指標(biāo)；同理，決策單元的產(chǎn)出指標(biāo)數(shù)據(jù)可用矩陣Y=yoj(o=1,2,...,s；j=1,2,...，口)表示yoj表示第j個(gè)DMU的第o個(gè)投入指標(biāo)；設(shè)vi=(v1，v2，...，vm)和uo=(u1，u2，...，us)分別為輸入和輸出指標(biāo)的權(quán)向量，規(guī)定用產(chǎn)出指標(biāo)的線性加權(quán)組合與投入指標(biāo)的線性加權(quán)組合比來表示DMq的效率EjOEC^E.=(E)/()n=size(X',1);m=size(X,1);s=size(Y,1);epsilon=10A-10;%定義非阿基米德無窮小f=[zeros(1,n)-epsilon*(1,m+s)1];A=zeros(1,n+m+s+1);b=0;LB=zeros(n+m+s+1,1);UB=[];LB(n+m+s+1)=-Inf;fori=1:n;Aeq=[Xeye(m)zeros(m,s)-X(:,i);Yzeros(s,m)-eye(s)zeros(s,1)];beq=[zeros(m,1);Y(:,i)];w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);%解線性規(guī)劃，得DMU：的最佳權(quán)向量w;endw%輸出最佳權(quán)向量lambda=w(1:n,:)%輸出最優(yōu)值s_minus=w(n+1:n+m,:)%輸出s*-s_plus=w(n+m+1:n+m+s,:)%輸出s*=theta=w(n+m+s+1,:)%輸出.*Aeq=[Xeye(m)zeros(m,s)-X(:,i) Yzeros(s,m)-eye(s)zeros(s,1)];clc,clear;x=xlsread('111.xls','C13:F32');%多目標(biāo)輸入y=xlsread('111.xls','G13:G32');%多目標(biāo)產(chǎn)出X=x';Y=y';n=size(X',1);m=size(X,1);s=size(Y,1);A=[-X'Y'];b=zeros(n,1);LB=zeros(m+s,1);UB=[];fori=1:n;f=[zeros(1,m)-Y(:,i)'];Aeq=[X(:,i)'zeros(1,s)];beq=1;w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);%解線性規(guī)劃，得到DMUi的最佳權(quán)向量wiE(i,i)=Y(:,i)'*w(m+1:m+s,i);%求出DMUi的相對效率值Eiiendfori=1:20xiaolv(i)=E(i,i);endw;%輸出最佳權(quán)向量E;%輸出相對效率值xiaolvomega=w(1:m,: