今天推出的是利用Shapley方法，测算DEA中投入产出变量重要性的方法。 Shapley 值法是由 L. S. Shapley 于 1953 年提出的方法。Shapley 值的特点为参与人之间的平等性。可以通过构造联盟、联盟内外平等和联盟中参与人平等性这三个角度体现。基于 Shapley 值进行联盟成员的利益分配体现了各盟员对联盟总目标的贡献程度，避免了分配上的平均主义。通过每个变量的 Shapley 值，可以求出变量在效率中的边际影响，从而合理筛选投入产出指标，构建科学的投入产出指标体系，为后续测算煤炭生态效率奠定基础。 参考文献：《A Shapley value index on the importance of variables in DEA models》Yongjun Li, Liang Liang DEA模型中的效率贡献率(ECR) 1. 引言 Wagner和Shimshak(2007)开发了一种逐步减少变量的程序方法。该方法的核心是计算每个变量对效率评估的影响程度;在他们的论文中,给定变量的影响被定义为包含和不包含该变量的两个DEA模型中所有D ...

今天推出的工具是非平衡面板的malquist指数计算 参考文献：《非平衡面板数据的全要素生产率测算方法》马占新，苏日古嘎 下面是该理论部分的主要内容。 1. 传统 DEA-Malmquist 指数方法及其 RD 分解模型 假设共有 \(n\) 个决策单元，每个决策单元在第 \(t\) 期 \((t=1,2,\cdots,T)\) 均用 \(m\) 种投入获得 \(s\) 种产出。其中，第 \(j\) 个决策单元在第 \(t\) 期的投入产出值为 \((x_j^t, y_j^t) > 0\)。则根据 DEA 方法的基本原理 (Charnes 等, 1978; Banker 等, 1984)，可以给出第 \(t\) 期的生产可能集如下： \[ S_c^t = \{(x^t,y^t)|\sum_{j=1}^n \lambda_jx_j^t \leq x^t, \sum_{j=1}^n \lambda_jy_j^t \geq y^t, \sum_{j=1}^n \lambda_j = \delta, \lambda_j \geq 0, j = 1,2,\cdots, ...

今天推出的是复合系统协调度这个模型的工具。 参考文献:《“干线公路—城市结点”复合系统协调度分析模型》 复合系统整体协调度模型以协同学的序参量原理和役使原理为基础，模型的参变量选择是模型合理性的前提， 重点选择在整个系统发展演变过程中起主导性作用的因素作为模型序参量变量，它决定着系统演化的特征和规律。 模型具体如下： 1、计算子系统有效度 将"干线公路—城市结点"系统表示为 \(S=\{S_1,S_2, S_3,S_4\}\)，其中 \(S_i\) 为第 \(i\) 子系统，\(i=1, 2, 3, 4\)，分别对应干线公路、社会经济、城市交通和生态环境等四个子系统； 设子系统 \(S_i\) 发展过程中的序参量为：\(x_i=(x_{i1},x_{i2},\cdots,x_{in})\)，其中 \(a_{ij} \leq x_{ij} \leq \beta_{ij}\)，\(a_{ij}\)，\(\beta_{ij}\) 为序参量 \(x_{ij}\) 的阈值。 则第 i 个子系统的有序度\(\mu_\mathrm{i}(\mathbf{x}_\mathrm{ ...

今天推出的是存在共享投入的两阶段博弈交叉DEA模型。 网络 DEA 博弈交叉效率模型不仅适用于基本两阶段网络结构，还适用于多种网络结构，本文进一步将模型拓展为存在共享投入的两阶段DEA 博弈交叉效率模型。 存在共享投入的两阶段网络结构中，假设有 n 个 DMU;在第一阶段，\(DMU_j(j=1,2,...,n)\)用 m 种共享比例分别为\(\alpha_{ij}(i=1,2,...,m)\)的外源投入\(x_ij(i=1,2,...,m)\)产生 q 种产出\(z_pj(p=1,2,...,q);\) 在第二阶段，\(DMU_{j}(j=1,2,...,n)\)用 m 种共享比例分别为的外源投入\[(1-\alpha_{i_{i}})\],\(x_{ij}(i=1,2,....,m)\)和第一阶段的 q 种产出\(z_{pj}(p=1,2,...,q)\)产生 S 种产出\(y_{rj}(r=1,2,....,s)\),这 S 种产出离开系统。其中，共享比例\(\alpha_{i_i}(i=1,2,...,m)\)是未知参数，管理者可根据现实情况设置共享比例的下界 ...

今天推出的是EG指数计算工具。 EG 指数法的运用相对成熟，它弥补了早期测度方法的劣势，解决了Gini 系数法失真的原因并结合 了赫芬达尔指数法的使用，同时对能够采集到的数据具有良好的适用性。 其计算公式如下： \[\gamma_{EG}=\frac{G-(1-\sum_{i=1}^{M}x_{i}^{2})H}{(1-\sum_{i=1}^{M}x_{i}^{2})(1-H)}=\frac{\sum_{i=1}^{M}(s_{i}-x_{i})^{2}-(1-\sum_{i=1}^{M}x_{i}^{2})\sum_{j=1}^{N}Z_{j}^{2}}{(1-\sum_{i=1}^{M}x_{i}^{2})(1-\sum_{j=1}^{N}Z_{j}^{2})}\] 其中，N为经济体中某产业企业的个数；M 代表经济体地理区域的划分个数\(;_{\lambda_i}\)表示\(i\)市总就业人数占全国总就业人数的比重；s_i表示\(i\)市旅游产业就业人数占该产业全部就业人数的比重；\(G=\sum_{i=1}^{n}(s_i-x_i)^2\)为空间Gini 系数\( ...

今天推出的是DO指数计算工具。 DO 指数方法由Duranton and Overman（2005）提出，用于对产业空间集聚程度进行度量。 DO 指数的计算分为三步: (1)计算核密度函数。假定行业\(I\)中存在\(n\) 家企业，运用每家企业的经纬度数据，可以计算出\(n( n- 1) /2\)个两两企业对彼此的空间距离，它与实际交通距离存在一定偏差，但对估计结果影响不大，Duranton and Overman(2005)曾详细阐述了这一问题。计算该行业两两企业彼此距离的核密度函数\(\hat{K}(d)\)的公式如下： \[ \hat{K}(d)=\frac{1}{n(n-1)h}\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}f(\frac{d-d_{i,j}}{h}) \] 其中，\(d\) 表示两两企业彼此的空间距离，\(d_{i,j}\)表示企业\(i\) 与企业\(j\) 之间距离值，\(f\)为高斯核密度函数，h 表示窗宽。 因为企业规模的大小会明显地影响企业的分布模式特征，所以本文采用企业就业人数作为规模的权重进入核密度函数 ...

今天推出的是网络SBM模型。 传统的SBM模型处理关于多输入与多输出的dmu相对效率的测量。这些模型的缺点之一是忽略了中间产品或链接活动。在指出将中间产品纳入DEA模型的需求后，作者提出了一个基于松弛的网络DEA模型，称为网络SBM，它可以形式化地处理中间产品。使用该模型，我们可以评估部门效率以及决策单位(dmu)的整体效率。 其结构图如下： image-20240412172944704 变量约束如下： \[ \] \[\begin{aligned} &\mathbf{x}_{o}^{k}=\mathbf{X}^{k}\lambda^{k}+\mathbf{s}^{k-}\quad(k=1,\ldots,K), \\ &\mathbf{y}_{o}^{k}=\mathbf{Y}^{k}\lambda^{k}-\mathbf{s}^{k+}\quad(k=1,\ldots,K), \\ &\mathbf{e}\lambda^{k}=1\quad(k=1,\ldots,K), \\ &\lambda^{k}\geqsl ...

之前做的上市公司团体机构持股比例的代码，重新写了一遍，现在可以直接利用国泰安下载的数据进行运算。 分别执行三个py文件，即可得到团体机构持股比例，最大团团体机构持股比例两个指标，可以直接导入stata进行回归等计算，方便快捷。 需要的同学可以联系微信 canglang12002 往期推文： 空间马尔科夫工具更新 共享投入型关联两阶段生产系统的网络DEA效率测度与分解复现 基尼系数双维度分解计算工具 动态SBM模型（DSBM)复现 共享投入的两阶段DEA模型复现 基于凸分位数回归的影子价格计算工具 Doyle和Green(1994)的对抗性和仁慈性DEA 一种新的两阶段网络DEA模型 适用于面板数据的动态StoNED模型计算工具 零和博弈SBM模型（ZSG-SBM）模型 三阶段动态网络DEA(DNSBM)模型的实现 计算相对资源承载力模型的工具 基于参数化的方向性距离函数（DDF)估算污染物影子价格的工具 基于非期望产出的RAM 碳环境效率模型 使用遗传算法或NSGA2算法解决多式联运问题 莫兰指数计算小工具 Panda_DEA增加至强有 ...

Markov工具已经推出很久了，最近对其做了一些更新，具体如下： 1、可以指定状态数。 之前默认是分4类，可是有些用户有分三类或者五类的需求，现在用户可以自由指定分类状态数。 2、增加了卡方检验 现在很多关于空间markov的文献，增加了卡方检验，这里一并加上。 3、初始分布和稳态分布 结果表中增加了传统markov的初始分布和稳态分布结果，以及空间markov的稳态分布结果。 4、分位数 原先的版本，需要用户手动计算分位数，比较麻烦，本次更新的版本，用户可以指定软件自行计算分位数， 当然，如果有用户需要手动指定分位数，我们也保留了这一功能。 2024年3月14日。 PS：由于某些原因，工具视界官网转到了github，网址为gongju128.github.io 新网站可以评论了哦，欢迎大家留言，提出意见！

今天推出的是一个关于2011年陈凯华和官建成的《共享投入型关联两阶段生产系统的网络DEA效率测度与分解》这篇 文献中建立的DEA模型结果的复现。 该文提出的考虑中间产出的共享投入的两阶段网络DEA模型，在规模报酬不变时，其规划式如下： \[\begin{aligned}&E_k=\max\sum_{p=1}^qW_p^1Z_{pk}+\sum_{r=1}^sU_rY_{rk}\\&\mathrm{s.t.}\begin{cases}\sum_{i=1}^m\pi_i^1X_{ik}+\sum_{i=1}^mV_i^2X_{ik}-\sum_{i=1}^m\pi_i^2X_{ik}+\sum_{p=1}^qW_p^2Z_{pk}=1\\\sum_{i=1}^m\pi_i^1X_{ij}-\sum_{p=1}^qW_p^1Z_{pj}\geq0,j=1,2,\cdots,n\\\sum_{i=1}^mV_i^2X_{ij}-\sum_{i=1}^m\pi_i^2X_{ij}+\sum_{p=1}^qW_p^2Z_{pj}-\sum_{r=1}^sU_rY_{rj}\g ...