- 博客(45)
- 收藏
- 关注
RSS订阅原创 三阶段DEA之SFA
R中的frotier包，这个包是对frotier4.1的改进，修复了刚才提到的非平衡面板问题，在这个基础上，我们做了一个界面化的工具，用户只需要输入数据，就可以得到结果，最棒的是，可以直接获得调整后的投入，无需各种繁琐的替换等操作，简直不要太方便！frontier由于年代久远，操作并不方便，很多人按照教程操作，却出不来结果，除此之外，这个软件在非平衡面板数据的时候，效率值会存在一些问题，excel模板对大多数人来说也不够方便。第三个阶段，利用调整后的投入变量，重做一次DEA模型。大家如果需要，请联系微信。
2024-06-28 10:28:47
852
原创 存在共享投入的两阶段博弈交叉DEA模型
网络 DEA 博弈交叉效率模型不仅适用于基本两阶段网络结构，还适用于多种网络结构，本文进一步将模型拓展为存在共享投入的两阶段DEA 博弈交叉效率模型。关于该模型，使用matlab和julia分别完成了模型建模，均能得到和论文相同的结果（指总效率，分阶段效率有差异，但是排名与论文中的结果近似。两阶段博弈DEA模型，由于需要大量计算，如有N个DMU，那么需要进行G*N^2次的线性规划，为了加快运行速度，在计算式采用文献中提供的新的算法设计。今天推出的是存在共享投入的两阶段博弈交叉DEA模型。
2024-05-10 18:38:26
890
原创 EG指数计算工具
官方网址：https://gongju128.github.io/2024/04/22/%E8%AE%A1%E9%87%8F%E7%BB%8F%E6%B5%8E/EG%E6%8C%87%E6%95%B0%E6%8E%A8%E6%96%87/EG 指数法的运用相对成熟，它弥补了早期测度方法的劣势，解决了Gini 系数法失真的原因并结合。Elilsion & Glaeser指出集中度指标分为三个区间：若。了赫芬达尔指数法的使用，同时对能够采集到的数据具有良好的适用性。,表明该区域内的产业分布比较均匀；
2024-04-22 22:30:38
857
原创 DO指数计算工具
的空间分布模式，在计算该行业两两企业彼此距离的核密度函数后，再将其与无约束条件下企业随机分布假定的核密度函数进行比较。然而因为自然地理环境以及土地利用政策的局限，很多地区无法建造工业企业，因此，要利用目前约束条件下企业随机分布的核密度函数。个两两企业对彼此的空间距离，它与实际交通距离存在一定偏差，但对估计结果影响不大，Duranton and Overman(2005)曾详细阐述了这一问题。的 1000 次反事实的核密度按升序进行排序，并选择 5%与 95%作为置信区间的下限。
2024-04-22 22:26:20
362
原创 DO指数计算工具
根据固定距离,本文将获得的行业的 1000 次反事实的核密度按升序进行排序，并选择 5%与 95%作为置信区间的下限与上限,伴随着距离的增加可得到置信带，该置信带一般被称为局部置信区间。假定行业中存在 家企业，运用每家企业的经纬度数据，可以计算出个两两企业对彼此的空间距离，它与实际交通距离存在一定偏差，但对估计结果影响不大，Duranton and Overman(2005)曾详细阐述了这一问题。其中，和分别表示企业与企业的就业人数。今天推出的是DO指数计算工具。
2024-04-22 22:25:45
832
原创 网络SBM模型（NSBM)复现
传统的SBM模型处理关于多输入与多输出的dmu相对效率的测量。这些模型的缺点之一是忽略了中间产品或链接活动。在指出将中间产品纳入DEA模型的需求后，作者提出了一个基于松弛的网络DEA模型，称为网络SBM，它可以形式化地处理中间产品。使用该模型，我们可以评估部门效率以及决策单位(dmu)的整体效率。“固定”链接值的情况。可以看到，结果与论文中的完全一致。如果需要做该模型，请联系微信。今天推出的是网络SBM模型。
2024-04-12 18:10:11
341
原创 空间马尔科夫工具更新
PS：由于某些原因，工具视界官网转到了github，网址为gongju128.github.io。原先的版本，需要用户手动计算分位数，比较麻烦，本次更新的版本，用户可以指定软件自行计算分位数，结果表中增加了传统markov的初始分布和稳态分布结果，以及空间markov的稳态分布结果。之前默认是分4类，可是有些用户有分三类或者五类的需求，现在用户可以自由指定分类状态数。现在很多关于空间markov的文献，增加了卡方检验，这里一并加上。当然，如果有用户需要手动指定分位数，我们也保留了这一功能。
2024-03-14 22:30:17
230
原创 共享投入型关联两阶段生产系统的网络DEA效率测度与分解复现
值的一提的是，作者给出的CRS结果中，DMU10的结果和DMU9的结果相同，这里应该是处理结果时出了问题，根据原始数据，和DMU10数据相近的是DMU11，在VRS的结果中也正是如此，DMU10和DMU11的结果非常接近。今天推出的是一个关于2011年陈凯华和官建成的《共享投入型关联两阶段生产系统的网络DEA效率测度与分解》这篇。的最优解组合后，可进一步通过下述公式计算内部第一个子过程和第二个生产过程的技术效率值.的借助规划 获得最优解组合后，可通过下述公式 计算被测评的生产系统 DMU。
2024-01-26 00:56:53
889
原创 基尼系数双维度分解计算工具
今天推出的是一个基尼系数双维度分解工具，之前的基尼系数分解，要么是基不同的子群，或者基于不同的分项，但是有时候我们需要同时基于分组和分项来研究，如GDP，既需要根据区域分组研究区域差异，也需要揭示总体差异的产业来源，这时候常用的泰尔指数、Dagum基尼系数分解就不再适用了。上式表明基尼系数被分解成 q 个部分，每个部分就是该指标各分项指标对于总体差异的贡献。我们这个工具特点之一，就是可以在结果中指明对应分组的名称，而非用数字代替，简洁明了。基尼系数可以用下面的公式来计算，其中，个子群，每个子群可以表示为。
2024-01-19 12:16:58
488
原创 动态SBM模型（DSBM)复现
今天推出的是动态SBM模型，参考文献是：《Dynamic DEA: A slacks-based measure approach》，文献doi是：10.1016/j.omega.2009.07.003。（1） 理想的（好的）链接，这表示需要结转。在我们的模型中，期望的链接被视为输出，并且链接值被限制为不小于观察到的值。在我们的模型中，不希望的链接被视为输入，其值被限制为不大于观察到的值。与自由链接类似，固定链接通过两个术语之间的连续性条件间接影响效率得分。它的值可以从观察到的值增加或减少。
2024-01-02 13:13:29
1000
原创 共享投入的两阶段DEA模型复现
今天推出的是含有共享投入的两阶段DEA模型,参考文献是：Yao Chen的《DEA model with shared resources and efficiency decomposition》，doi是：10.1016/j.ejor.2010.03.031。可以看到，总效率值，和最大化第一阶段和第二阶段的效率值是基本相同的，但是由于不同的优化求解器，得到的w1和w2不同，就导致最大化其中一个阶段效率后，求得的另外一个效率值，和论文结果是有差异的，这点注意即可。如果需要该模型，欢迎联系微信。
2023-12-29 10:50:06
446
原创 基于凸分位数回归的影子价格计算工具
另外一个就是模型求解后不同分位数的取舍判断，以及计算MRT和MP的部分。3、这里tau有10个不同的取值，求出使残差最大的非负的tau的取值，记为tau*，将MRT和MP计算为，tau*和tau*+0.1的情况下，以残差绝对值进行加权的加权平均值。今天介绍给大家一种新的影子价格计算工具，基于凸分位数回归来计算污染物的影子价格，弥补了以DEA模型计算影子价格所不具备的一些缺陷，是一种不错的新工具。有鉴于此，我们已经做出了相应的工具，可以轻松完成所有流程，直接得到影子价格的估计结果，需要的同学可以联系微信。
2023-09-26 11:50:10
134
2
转载 Doyle和Green(1994)的对抗性和仁慈性DEA
DEA 方法基于 自评思想 ，表现在权重由决策单元自身决定 ，最大化自身效率而尽量最小化他人效率．因此 ，此权重往往对其他决策单元的分配极为悬殊 ，对自己有利的指标赋权很大，对自己不利的指标赋权很小 ，甚至赋予零权重，这种权重下计算出的效率值并不能完全反映出 DMU的优劣。在保证当前DMU得到最大效率值的前提下，使其他DMU的效率值之和尽可能小，这是对抗性DEA，相反，使其他DMU的效率值之和尽可能大，这是仁慈性DEA。第一步，求解CCR模型，计算出各个决策单元的自评效率值。下面是该工具得到的两种效率值。
2023-09-15 17:07:59
185
转载 一种新的两阶段网络DEA模型
目前，主流的两阶段DEA方法，包括四种，第一种是独立的两阶段DEA，它将标准的DEA方法，独立应用于两个不同的阶段；第二种是连通的二阶段方法，这种方法在计算整体效率的时候，考虑了相互作用；第三种是关系型两阶段DEA方法，它构建了整个效率和每个子过程效率之间的关系；越来越多的研究人员通过考虑共享输入来改进传统的两阶段模型。然而，很少有研究同时考虑共享输入、额外中间输入和自由中间输出，这篇文献即提出了一个同时考虑共享输入、额外中间输入和自由中间输出的网络结构进行建模，以进一步发展DEA方法。
2023-08-18 11:17:24
321
转载 适用于面板数据的动态StoNED模型计算工具
与SFA模型相比，随机非参数数据包络模型的回归是运用了凹面非参数最小二乘法（CNLS），而SFA模型是运用的改进最小二乘法(MOLS),其具体步骤为：1、对残差ei进行最小二乘估计，2、基于残差分布的二阶中心矩和三阶中心距，对方差进行估计，3、基于条件分布，对非效率项进行估计。在静态 StoNED 模型中，x表示生产活动中的投入向量， y表示生产活动的产出向量。StoNED模型最初是对截面数据开展研究，认为生产技术和效率项不随时间变化而变化，属于静态随机非参数数据包络模型。
2023-08-13 14:45:48
334
转载 零和博弈SBM模型（ZSG-SBM）模型
今天介绍零和博弈DEA模型中比较常用的一种，基于SBM模型的ZSG-SBM模型。需要的同学，可以联系微信。
2023-07-30 18:06:06
183
原创 三阶段动态网络DEA(DNSBM)模型的实现
传统的DEA（数据包络分析）模型处理关于多个输入与多个输出的决策单元（DMU）的相对效率的测量。今天提出的三阶段动态DEA(DNSBM)模型，是一种将这两个DEA模型相结合的模型，称之为动态网络DEA。这种组合模型不仅使我们能够获得整个观测周期内DMU的整体效率，而且还可以进行进一步的分析，即观测DMU的周期效率的动态变化。其中链接变量和结转变量分为四种情况，具体在运用该模型的时候，根据实际情况，添加该情况对应的约束即可。今天要推出的是一个比较新的DEA模型，三阶段动态网络DEA模型。
2023-07-22 21:45:47
676
转载 计算相对资源承载力模型的工具
1921年人类生态学家Park和Burgess提出承载力（Carrying Capacity）概念，即“某一特定环境条件下（主要指生存空间、营养物质、阳光等生态因子的组合），某种个体存在数量的最高极限”。2000年国内学者黄宁生等提出了相对资源承载力概念，并以广东省为例分析了该省的相对资源承载力与可持续发展情况。改进4：在上文所提出的承载状态标准基础上，进一步提出各状态下承载状态度的计算标准，这一步并不是我们关注的重点，这里略过，大家可以直接看文献。今天推出的是计算相对资源承载力模型的工具，文献参考自《
2023-07-12 15:26:09
68
原创 基于参数化的方向性距离函数（DDF)估算污染物影子价格的工具
早期影子价格研究中超越对数形式的谢泼德距离函数占据了很大比例，而方向性距离函数（DDF）则是目前主流的影子价格模型，这是因为DDF函数下期望产出和非期望产出的数量能够向不同方向变化，可同时满足企业期望产出增长和非期望产出减少，更符合碳减排政策和低碳经济的发展要求。污染物的影子价格定义为减少一个单位非合意产出所需要放弃的合意产出的机会成本。为了求解碳排放的影子价格，在产出距离函数与收益函数之间建立联系，根据方向性距离函数的平移性质，引入方向向量g ＝（gy，–gb），则收益函数可以转化为。
2023-05-10 22:15:42
648
转载 使用遗传算法或NSGA2算法解决多式联运问题
多式联运路径优化问题可以描述为：某次物流作业经过由N 个城市组成的运输网络，需要将一批货品从起始地O 运送至目的地D，在网络中的两个城市之间都有K 种运输的方式可供选择，各种运输方式的运输时间、费用及运输承载量都不相同，每到一个城市都可以选择是否更换运输工具，一种运输方式转运为另一种运输方式时，需要支付该过程所产生的人员的劳务费用和并且占用一定的时间，因此问题就是怎么选择良好的运输路线与相应的运输方式，使得总物流费用最低。一次运输作业，总会有它的限定时间，能否在规定的时间内将货物送达，非常关键。
2023-04-25 17:06:22
418
转载 QLab增加熵权-topsis模块
2、构造加权规范化矩阵，传统的topsis是利用人为给定的权重加以计算，主观因素会影响最后的排序结果，这里采用的是熵权法得到的权重，避免主观因素。熵权TOPSIS 法是将熵权法的客观赋权与TOPSIS 法的多属性决策排序相结合，在TOPSIS 法构造加权规范化矩阵的过程中，利用熵权法得到的权重加以运算，如此一来，将客观的赋权方法与严谨的多维排序巧妙联系，帮助决策者进行科学评价。如果有负向指标，取倒数，做趋同化处理，得到新的决策矩阵，然后使用如下公式得到规范化决策矩阵。5、最后计算相对贴近度。
2023-03-16 18:28:13
130
转载 Panda_DEA增加至强有效前沿最近距离-MinDS模型
其中,Q（即软件中的M值）是一个比较大的正数，对于不同的求解器，M的值可能有所不同，过大或者过小可能会影响到求解的结果，软件中默认该值=1000000，用户也可以根据自己的实际情况来更改。MinDS 模型在SBM 模型的基础上增加了一组混合整数线性约束，使得被评价DMU 的参考标杆位于同一个超平面内，进而通过求解混合整数线性规划来获得MinDS 模型的效率值。3、增加了至强有效前沿最近距离-MinDS模型，并为此支持外接Gurobi和GLPK两种求解器，显著提高求解速度。1、增加了两阶段网络DEA模型；
2023-02-24 21:28:11
343
转载 VSCode中使用armadillo+openblas的详细步骤
百度了下，基本上都是在VS2019中使用armadillo和openblas的教程，我大多数时候写的都是一些小项目，用VS未免有些杀鸡用牛刀的感觉，还是VSCode，好看又好用。armadillo是C++中非常好用的一个线性代数运算库，有着和matlab非常相似的语法，可以非常方便的将matlab代码转移到C++中。第二步：下载armadillo，只需要其中的include文件夹即可，然后放在项目文件夹下，用vscode打开该项目文件夹。，解压到任意盘，最好不要有中文路径，避免出错。
2023-02-13 21:21:44
946
转载 QLab增加耦合协调度模型
可以计算二元及以上子系统的协调耦合度，在计算协调发展度D的时候，提供平均权重和自定义权重两种方式供用户选择，操作还是一如既往的简单，只需准备好数据即可。耦合模型是从物理学中演变而来，它所体现的是两个系统之间通过彼此的相互作用互相影响，耦合度是对不同系统耦合程度的定量描述。QLab新版本增加了耦合协调度模型，同时对界面做了一定程度的更新，模型被分为经济模型和评价模型两个类别。
2023-02-11 12:39:08
640
转载 三行R代码搞定山脊图绘图
这里的limits是对刻度的范围进行限制，breaks这个数组里，可以写入自定义的数组。非常简单，一列是数据，一列是年份。有时候我们对比几组数据的分布情形，需要把几组数据的分布叠加在一起，这就用到山脊图了。今天带来一个最简单的山脊图绘图教程。
2023-02-05 20:56:10
183
转载 Panda_DEA增加RAM-DEA和BAM-DEA模型
由于RAM模型中的参数是由投入和产出的极差构成，会出现无效率值过小，因此导致不同的决策单元效率极为接近，难以有效区分的情形，因此Cooper et al(2011)在RAM模型的基础上提出了一种同时具有非径向性、非角度性、可加性以及更高辨识度等优点的BAM模型。RAM-DEA是1998年由Aida和Cooper等提出的一种非径向DEA模型，由于该模型构建对投入产出指标的设计考虑了指标之间的相互影响效应，并对这些效应进行了定性分析，因此得到了广泛的研究与应用。
2023-01-04 12:44:45
425
转载 QLab1.4正式发布！增加全排列多边形图示指标法
全排列多边形图示指标法定义为：设共有n 个指标(标准化后的值),以这些指标的上限值为半径构成一个中心n 边形,各指标值的连线构成一个不规则中心n 边形,这个不规则中心n 边形的顶点是n 个指标的首尾相接的全排列,n 个指标总共可以构成(n-1)!用户可以自行选择数据文件，以及对应的保存目录，数据文件不需要使用特定的文件名，运行结果的文件名自动加上时间戳，多次运行，也不用关闭之前的结果文件。其中U为指标x的上限，取最大值，L为指标x的下限，取最小值，T为指标x的临界值，取平均值。
2022-12-16 21:55:27
382
转载 大量数据时的Dagum基尼系数分解工具
去年写过的Dagum基尼系数分解工具，本来已经足够使用，但是前天有一位用户，数据量有近15w行，我从未设想过会出现如此大的数据，程序足足运行了将近5个小时才跑出结果，然而用户还有很多类似的数据，不能总是如此。于是在原先的工具基础上，做了优化，优化后的运行速度得到极大提升，同样的数据，在优化后，只用了不到30分钟就运行出结果，时间缩短到未优化前的1/10。
2022-10-21 16:45:10
579
转载 按收入来源分解基尼系数的工具
Wf为第f种收入在所有收入额中的占比，gf为第f种来源的人均值，g为全部收入的人均值，Cf表示f收入来源的集中率。Qfk是第f种收入从第1个县城累积到k个县城的合计数在全部收入中的占比，Wfk是第k个县的第f种收入在全部县f来源总和中的份额，mfk是第k个县f来源的人均值。其中，n表示全省一共n个县城，Pk表示第k个县城的人口在全省人口总数中的份额，Wk为第k个县城收入在全省县级收入总额中的比重，Qk是第一个县累积到第k个县的收入比重之和。支持分解任意多种的收入来源，一年还是多年数据都可以一键生成结果。
2022-10-01 10:52:15
184
转载 Panda-DEA_1.0正式发布！
经过一段时间的努力工作，Panda-DEA_1.0版本终于正式发布了！软件的模型选项界面如下：距离函数导向规模报酬模型面板模型此外，一些比较前沿的DEA模型论文的复现和常见的基础模型，Panda-DEA提供了快捷选择模型模块，只需准备好相应的数据，可以一键复现论文的结果，参考理论，复现结果，最大减轻大家写论文的苦恼。
2022-09-29 21:56:43
358
1
转载 含有非期望产出的ZSG-DEA模型
在ZSG模型中，一次调整一般并不能达到零和DEA有效，这里作者参考林坦和宁俊飞提出的迭代法进行迭代求解，直到各个地区的ZSG效率值均达到1为止。其中yn、en、un和Xk,n分别表示第n个地区的期望产出(GDP），能源投入、非期望产出（CO2)，和非能源投入。根据零和博弈DEA的思想，当一个地区的用能权和碳排放权增加时，其他地区必须减少相同数量的用能权和碳排放权，以保持总量不变。分别为能源投入和CO2排放效率值的权重，这里作者设置的权重各为0.5，也可以根据研究主题自行设置。即是零和DEA效率。
2022-09-08 17:31:49
707
原创 两阶段网络DEA及其计算
在两阶段网络DEA模型中，第一阶段的产出，作为第二阶段的投入。不同于之前的研究将整体的生产过程和两个子过程视作独立的，Kao这篇文献在测量效率时，考虑了两个子过程之间的串联关系，并且整体效率是两个子过程效率的乘积，因此这种关系两阶段DEA方法计算得到的效率相对于从独立两阶段DEA方法计算的效率更有意义。但是在乘数模型中，获得的最优解可能并不唯一，这使得不同DMU之间的效率比较可能会出现问题。然后根据Ek=Ek1*Ek2这一关系，计算Ek2的值。今天推出的是两阶段网络DEA模型。可以看到，结果大体保持一致。.
2022-08-25 19:00:24
2487
1
原创 Panda-DEA，一款新的DEA模型软件
经过数十年的发展，DEA模型及理论的发展十分迅速，国内外期刊，不断有人提出崭新的DEA模型，如零和收益（ZSG）DEA、网络DEA模型、含有非期望产出的SBM模型等等。大家在运用DEA模型的时候，碰到的一大难题就是如何求出DEA模型的结果，那一行行的约束公式，往往令人望而生畏。虽然市场上已经出现了为数众多的DEA模型软件，但是大多只能计算一些常见的DEA模型，对于一些前沿模型，则无能为力。2、结果以论文为导向，如零和收益ZSG-DEA模型，会直接进行多次迭代，并将迭代结果写入文件，用户可以直接放进论文；.
2022-07-26 22:02:00
2586
2
原创 DKM_1.3收敛交叉映射
绿色曲线表示，用时间序列A的影子流形预测时间序列B，随着库长的增加，相关系数仍然很低，这说明A对B的驱动比B对A的驱动更强，也就是说，A是B的原因。如，我们有一个模型来预测时间序列Y，如果在模型中移除时间序列X，则模型对Y的预测能力大大降低，此时可以说X是Y的一个原因。此前嵌套泰尔指数分解是单独的脚本，这次一并合并到DKM_1.3中。计算估计值和实际值的相关系数，即CCM相关系数，在不同的时间序列长度中重复计算CCM相关系数，绘制CCM相关系数随长度L的变化情形，判断相关系数的收敛性。分别估计Y和X的值。.
2022-07-16 14:59:43
649
转载 用R语言做分组核密度图
多密度图是表示多个组的密度图。它允许比较它们的分布。这种图表的问题在于它很容易变得混乱：组相互重叠并且图形变得难以阅读。一个简单的解决方法是使用透明度。但是，它不能完全解决问题，通常最好考虑本文档中进一步建议的示例。不透明的分组核密度图：详细代码：用R语言做分组核密度图 - 工具视界...
2022-07-14 10:30:14
396
转载 最基本的桑基图
桑基图（Sankey diagram），即桑基能量分流图，也叫桑基能量平衡图。它是一种特定类型的流程图，右图中延伸的分支的宽度对应数据流量的大小，通常应用于能源、材料成分、金融等数据的可视化分析。因1898年Matthew Henry Phineas Riall Sankey绘制的“蒸汽机的能源效率图”而闻名，此后便以其名字命名为“桑基图”。代码参见工具视界网站：最基本的桑基图 - 工具视界...
2022-07-13 14:06:05
190
转载 使用ggplot2添加文本标签
这里演示如何使用将文本添加为标记。它的工作原理与geom_point()几乎相同，但添加的是文本而不是圆圈。代码如下：详细代码参见：使用ggplot2添加文本标签 - 工具视界
2022-07-12 16:07:36
1381
原创 dagum基尼系数分解工具
dagum基尼系数分解工具相比于传统的基尼系数而言，Dagum 基尼系数能够将其分解为地区内差距、地区间差距以及超变密度。Dagum基尼系数的相关计算公式如下：1、总体基尼系数：2.子群内部基尼系数3.子群之间基尼系数 4.子群内差异对总体基尼系数贡献5.子群间差异对总体基尼系数的贡献 6.超变密度的贡献 网上目前计算dagum基尼系数分解的方法，一般是matlab代码，该方法要求首先要把数据按子群划分为不同的列，同时要把列按照平均值从大到小的顺序依次排列，且每次只能计算一年的数据，有些则对分组有限制，
2022-07-11 15:03:36
3281
转载 利用R语言绘制小提琴图加箱线图
小提琴图用于显示数据的分布状态和概率密度，它同时具有箱线图和密度图的特征，用于显示数据的分布形状。 代码：利用R语言绘制小提琴图加箱线图 - 工具视界
2022-07-11 12:17:39
578
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人