今天推出的是一种可以适用于面板数据的动态StoNED模型的计算工具 静态 StoNED 模型 StoNED（随机非参数数据包络分析法）最早由学者 Kwosmane首先提出，该模型是传统 DEA 模型与 SFA 模型相结合发展而来，在随机非参数数据包络模型中，公式中引入的随机成分的设定方法是运用了 SFA 方法来假设，对于确定部分则继续运用数据包络分析法的非参数处理方式。 StoNED模型最初是对截面数据开展研究，认为生产技术和效率项不随时间变化而变化，属于静态随机非参数数据包络模型。 在静态 StoNED 模型中，x表示生产活动中的投入向量， y表示生产活动的产出向量。y=f(x)表示决策单元的生产技术，f为生产函数。 首先根据 DEA 模型的假设条件，将生产函数看作是没有特定的函数形式的单调递增的凹函数，然后运用 SFA 模型的方法，在模型中引入随机参数,由于存在残差ei，导致观测值yi与真实值f(xi)不相等，残差项ei=vi-ui，其中vi为误差项，ui为非效率项，且ui>0，其具体形式如下： 与SFA模型相比，随机非参数数据包络模型的回归是 ...

今天介绍零和博弈DEA模型中比较常用的一种，基于SBM模型的ZSG-SBM模型。 在产出导向SBM 模型评价的基础上，产出导向ZSG-SBM 模型基于“零和收益”的思想对无效决策单元的非期望产出要素松弛量进行重新分配，以实现所有决策单元到达效率前沿，即实现了系统最优效率条件下对非期望产出的分配，其原理如下图所示： 其规划式如下： 需要的同学，可以联系微信canglang12002 参考文献：《中国“十三五”时期省际碳减排目标的效率分配》郭文，刘小峰，吴孝灵 往期推荐： 三阶段动态网络DEA(DNSBM)模型的实现 计算相对资源承载力模型的工具 基于参数化的方向性距离函数（DDF)估算污染物影子价格的工具 基于非期望产出的RAM 碳环境效率模型 使用遗传算法或NSGA2算法解决多式联运问题 莫兰指数计算小工具 Panda_DEA增加至强有效前沿最近距离-MinDS模型 QLab增加耦合协调度模型 广义SBM模型的matlab代码 QLab1.4正式发布！增加全排列多边形图示指标法 大量数据时的Dagum基尼系数分解工具 Panda-DEA_1 ...

今天要推出的是一个比较新的DEA模型，三阶段动态网络DEA模型。 传统的DEA（数据包络分析）模型处理关于多个输入与多个输出的决策单元（DMU）的相对效率的测量。这些模型的缺点之一是省略了DMU的内部结构。例如，许多公司由几个相互关联的部门组成，这些部门具有特定部门的输入和输出以及与其他部门的链接。为了反映现实世界，后来有学者开发了网络DEA模型，使用链接变量考虑DMU的内部结构。此外，公司的活动通常会持续多个时期。开发动态DEA模型是为了使用结转变量从长期角度评估DMU的性能。 今天提出的三阶段动态DEA(DNSBM)模型，是一种将这两个DEA模型相结合的模型，称之为动态网络DEA。这种组合模型不仅使我们能够获得整个观测周期内DMU的整体效率，而且还可以进行进一步的分析，即观测DMU的周期效率的动态变化。 其结构大致如上图所示，其中Link是不同Division之间的链接变量，Carry-Over是不同时期之间的结转变量。 其约束大致下： 其中链接变量和结转变量分为四种情况，具体在运用该模型的时候，根据实际情况，添加该情况对应的约束即可。 下面是目标函数： ...

今天推出的是计算相对资源承载力模型的工具，文献参考自《相对资源承载力模型的改进及其实证分析》黄常锋，何伦志 1921年人类生态学家Park和Burgess提出承载力（Carrying Capacity）概念，即“某一特定环境条件下（主要指生存空间、营养物质、阳光等生态因子的组合），某种个体存在数量的最高极限” 。该理论最早被引入草原管理应用中，相应提出草地承载力、最大载畜量等相关概念。 2000年国内学者黄宁生等提出了相对资源承载力概念，并以广东省为例分析了该省的相对资源承载力与可持续发展情况。相对资源承载力理论提出后，引起了国内许多学者和专家的研究兴趣并被他们广泛借鉴用来对区域可持续发展情况。 具体如下： AS1 具体的判断标准如下： （1）超载：指实际人口数P大于可承载人口数量Cs，即P - Cs > 0。 （2）富余：指实际人口数P小于可承载人口数量Cs，即P - Cs < 0。 （3）临界：指实际人口数P等于可承载人口数量Cs，即P - Cs = 0。 参考文献对该模型进行改进： 第一：在原模型计算相对土地资源承载力和相对经济资源承 ...

影子价格定价方法是测算污染物边际减排成本的重要方法之一。计算污染物影子价格的思路是将污染物视为生产活动过程中为了生产社会所需要的产品时产生的副产品(undesirable output), 即非合意产出。根据生产过程中合意产出与非合意产出之间的数量替代关系以及合意产出的市场价格, 计算得到污染物的影子价格。 污染物的影子价格定义为减少一个单位非合意产出所需要放弃的合意产出的机会成本。污染物的影子价格是环境政策制定者重点关注的价格信息, 因为它可以为污染排放权交易和环境税提供重要的参考信息。 早期影子价格研究中超越对数形式的谢泼德距离函数占据了很大比例，而方向性距离函数（DDF）则是目前主流的影子价格模型，这是因为DDF函数下期望产出和非期望产出的数量能够向不同方向变化，可同时满足企业期望产出增长和非期望产出减少，更符合碳减排政策和低碳经济的发展要求。根据估计方法的不同DDF可分为参数和非参数两种：参数方法虽不够灵活，但测算结果稳定且易于解释；非参数方法不受距离函数形式限制，因而灵活且便于模型创新。 今天介绍的是使用基于参数法的DDF来计算污染物影子价格的工具。 方向 ...

多式联运在我国是指通过转运衔接两种以上的运输方式，从而将货物从出发点送到目的地的过程。多式联运对多种运输方式的集成，是将相互独立的不同运输方式通过一些科学的算法统筹计划，从而减少整个运输过程的时间与费用，减少中途转运时因为没有计划，产生的等待时间，使用户尽可能得到最好的运输路线，合理的运输成本与运输时间。多式联运同时也是一种以实现货物运输效益达到最优为目标的组织方式。各个分段的分段承运人和作为多式联运运营商的总承运人合作完成货物的运输与中转的过程。它将运输全过程看作一个整体,通过联合各种不同的运输方式,形成连续的、最优的、综合一体化运输体系。 目前主要的多式联运形式有四种，它们分别为 “公铁联运”、“海空联运”、“海陆联运”、“陆桥联运”。 多式联运路径优化问题可以描述为：某次物流作业经过由N 个城市组成的运输网络，需要将一批货品从起始地O 运送至目的地D，在网络中的两个城市之间都有K 种运输的方式可供选择，各种运输方式的运输时间、费用及运输承载量都不相同，每到一个城市都可以选择是否更换运输工具，一种运输方式转运为另一种运输方式时，需要支付该过程所产生的人员的劳务费用和并且 ...

以往的DEA 模型由于角度限制只能将能源作为普通投入要素处理，因此，模型单纯地将能源消耗的增加认为是效率发生了恶化，故而测算出的能源效率忽略了能源替代效应，反映不出不同种类能源内部的替代或能源与劳动力、资源等普通投入要素之间的替代作用。此外，Zhou(2008)也验证了能源“混合效应”的存在性并分析了其作用原理，即当劳动力、资本等普通投入要素、期望产出以及非期望产出均不变时，替代能源作为能源投入的增加往往会使得其他能源的消耗量缩减得更多，因此，某种能源消耗量增加并不一定导致效率恶化，反而会使得效率得到改善。因此，在考察存在碳排放这一非期望产出条件下的碳环境效率时，应该将能源要素单独作为一类投入要素处理，其既可增加又可减少。 因此，定义基于非期望产出的RAM碳环境效率(Carbon Environmental Efficiency，简称为CE)模型如下： 经济效率(PE)核算的假设前提是没有环境管制，此效率核算的是单纯追求经济产出而不顾资源与环境质量的发展模式下的经济效率，而碳环境效率(CE)指凭借能源消耗的持续改进和劳动力、资本等投入的优化配置以及相应的碳减排政策的实施来满 ...

熵权TOPSIS 法是将熵权法的客观赋权与TOPSIS 法的多属性决策排序相结合，在TOPSIS 法构造加权规范化矩阵的过程中，利用熵权法得到的权重加以运算，如此一来，将客观的赋权方法与严谨的多维排序巧妙联系，帮助决策者进行科学评价。 熵权-Topsis的计算步骤主要分为两大步，首先是用熵权法得到指标权重，然后是用Topsis​计算相对贴近度。 熵权法部分可以参考之前的推文，这里​介绍Topsis部分： 1、对原始矩阵进行处理 如果有负向指标，取倒数，做趋同化处理，得到新的决策矩阵，然后使用如下公式得到​规范化决策矩阵。 2、构造加权规范化矩阵，传统的topsis是利用人为给定的权重加以计算，主观因素会影响最后的排序结果，这里采用的是熵权法得到的​权重，避免主观因素。具体操作是将规范化决策矩阵中的每一个向量与该向量对应的指标权重相乘​，从而得到加权规范化矩阵。 3、确定正负理想解 其中，正理想解为正向指标的最大值和负向指标的最小值组成的集合；负理想解为正向指标的最小值和负向指标的最大值组成的集合。 4、计算欧氏距离 分别计算待评价对象到这两个标杆的欧氏​距离， ...

在安装R包的arulesViz包的时候，总是失败，最后的提示就是这个错误。当时没有截图，大家知道我说的是这个就好了。 总结一下我的解决之路： 1、把R和rstudio都卸载，换上最新版本 2、把之前安装下载的R包，找到那个缓存路径，就是C盘中的downloaded_packages，里面内容全部删除。 3、重新安装aruleViz包，中间有个弹窗，选择否！！！ 安装完毕，可以愉快地使用了。

关于全局莫兰指数和局部莫兰指数，网上有很多介绍了，就不再赘述。 目前的教程，基本是基于R，stata，或者ArcGIS，操作不是很友好，大多数需要一个shp文件来，而且有些人又希望自定义空间权重矩阵，就不是很方便。 因此写了这个小工具，主要特点是，依然只需要用户准备好数据，打开小工具，就能得到莫兰指数散点图，以及对应的excel结果文件，包含全局莫兰指数及对应的Z值、P值；局部莫兰指数及对应的Z值、P值；画莫兰指数散点图用的数据，如果用户想进一步美化软件自动生成的散点图，可以根据数据自行画图。 特别地，不管是一年还是多年的数据，该工具都可以一次性计算出莫兰指数，并画出莫兰散点图。 部分结果如下图： 结果和stata的结果比对过，基本一致，可放心使用。 有需要联系微信canglang12002 如果您有DEA模型定制，数值计算，运筹优化等方面的需求，也欢迎联系！ 往期推荐： Panda_DEA增加至强有效前沿最近距离-MinDS模型 QLab增加耦合协调度模型 广义SBM模型的matlab代码 QLab1.4正式发布！增加全排列多边形图示指标法 大量数据时的 ...