LMDI(对数平均迪氏指数法）分解法的stata实现案例+code

zjying2000

29808

收藏 2021-10-12

LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应的资料、步骤还是比较少。本人在写论文的时候，通过搜集资料、研究，找到了stata的实现方法，来自于
Kerry Du老师写的一个模块ldmi，网址为：
LMDI: Stata module to compute Logarithmic Mean Divisia Index (LMDI) Decomposition (repec.org)
大家可以自己去研究一下。对应的参考文献就是Ang, B.W., 2005. The LMDI approach todecomposition analysis: a practical guide. Energy Policy 33, 867–871.

目前LMDI分解基本都以这篇文献为参考，论坛里也有，大家可以下载去深入研究一下，总体并不难。

但是对于LMDI方法的解读以及stata实现的步骤、过程，目前基本没有具体的中文资料。

本文在写论文的过程中，进行了总结，现上传上来供大家参考，里面包括stata的程序文件（ado、pkg、sthlp）、案例数据、中文的使用文档：

基本上看完这些资料，半小时内大家就能掌握LMDI的原理，以及stata操作实现的过程，若有问题可以联系本人。

资料如下：

LMDI composition.rar
大小:(742.11 KB)

只需: RMB 15元马上下载

本附件包括：

LMDI方法与案例解析与stata操作过程.pdf
lmdi.ado
lmdi.pkg
lmdi.sthlp
lmdi1.xlsx
lmdi2.xlsx

有些同学运行的过程中出现了：

ERROR：The specifiedvariables can not form an identity的问题，我查了一下资料，进行了一些改进，具体见下面这个附件，还有问题的话大家多沟通！！

关于ERROR The specified variables can not form an identity的解决.docx
大小:(57.67 KB)

马上下载

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

fukongwei4

2021-10-22 10:18:26

zjying2000 发表于 2021-10-12 14:09
LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应 ...

你好，请问我用stata做LMDI模型时，出现ERROR:The specified variables can not form an identity，是怎么回事呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2021-10-22 14:04:36

根据提示：The specified variables can not form an identity
大概意思就是说指定的变量不能形成标示，我感觉你的变量没有定义：区域、年份、部门等等
因为LMDI要从部门去分解的，你可以多看看我发的那个文件里数据的样式

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fukongwei4

2021-10-24 15:08:52

zjying2000 发表于 2021-10-22 14:04
根据提示：The specified variables can not form an identity
大概意思就是说指定的变量不能形成标示，我 ...

好的，万分感谢！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

希达wonderful

2021-11-13 10:51:40

看看

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

希达wonderful

2021-11-14 18:04:27

fukongwei4 发表于 2021-10-22 10:18
你好，请问我用stata做LMDI模型时，出现ERROR:The specified variables can not form an identity，是怎么 ...

啊，我也出现这个问题，请问您解决了吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

fredomom

2021-11-18 14:52:52

我也出现了ERROR:The specified variables can not form an identity的问题，而且我拿文件中的原始数据试了一下，只有分为八个部门时才能正常运行，分为其他数量的部门就不能运行，请问这是怎么回事啊

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fredomom

2021-11-18 15:00:34

你好我也出现了ERROR:The specified variables can not form an identity，而且我用文件中的原始数据测了一下，发现只有分为八个部门才能运行出结果，其他数量的部分不行，请问这是为什么呀

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2021-11-18 19:00:16

fredomom 发表于 2021-11-18 15:00
你好我也出现了ERROR:The specified variables can not form an identity，而且我用文件中的原始数据测了 ...

多谢提醒，我也试了试，也发现了这个问题，那说明原作者就是只设定了8部门才可以的吧？那这个就是原作者的一个问题了。
后来通过help ldmi，发现了一个默认的参数是用来控制这个问题的：
tolerance(real) specifies the tolerance; default is 0.01. The value is first used to check whether the specified variables can form an identity. Then it is used to check the difference between the real change and the total of the decomposed effects.
通过看了ldmi的源代码文件以及自己试了试，默认是tol（0.01），但是把它改大一些就可以，如下
lmdi E = I S Y, t(year) over(sector) tol(1)
具体可以看我又上传的一个文件。
我感觉是原作者貌似把部门标准控制的有些严格了。。。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fredomom

2021-11-22 15:49:13

zjying2000 发表于 2021-11-18 19:00
多谢提醒，我也试了试，也发现了这个问题，那说明原作者就是只设定了8部门才可以的吧？那这个就是原作者的 ...

谢谢解答！这个方法确实有用。但是还有一个问题就是，影响因素只能分为I S Y，可以减少但是不能增加，请问有解决的办法吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2021-11-24 22:14:04

fredomom 发表于 2021-11-22 15:49
谢谢解答！这个方法确实有用。但是还有一个问题就是，影响因素只能分为I S Y，可以减少但是不能增加，请问 ...

他的help说了，只能分解除出最多5个：
LMDI attributes the change of E into 5 components: effect of energy composition (change in F_{ijk}), intensity effect (change in I_{ij}), structural effect at the level 2 (change in S2_{ij}), structural effect at the level 1 (change in S1_{i}), and the activity effect (change in Y).
也就是 lmdi E= F I S2 S1 Y, t(t) over(sector subsector type)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

嘻嘻嘻茜茜

2021-12-25 21:18:58

我改了tol(1)但是还是会出现the specified variables can not form an identity ,是不是因为部门分的太多了呀我分了30多个，还有就是会不会是有零值的缘故？求大佬们解答

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

嘻嘻嘻茜茜

2021-12-25 21:22:54

我改了tol(1)也还是会出现the specified variables can not form an identity 怎么办是不是因为我部门分的太多了我分了三十多个，还有就是会不会是有零值的缘故？或者分解成五个因素了的缘故？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2021-12-25 22:01:15

嘻嘻嘻茜茜发表于 2021-12-25 21:22
我改了tol(1)也还是会出现the specified variables can not form an identity 怎么办是不是因为我部门分的 ...

你可以再改大一些试试？比如tol（10），这里tol（）括号里是个阈值，来控制数据计算算过程的一些条件的。。。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

嘻嘻嘻茜茜

2021-12-26 10:29:33

zjying2000 发表于 2021-12-25 22:01
你可以再改大一些试试？比如tol（10），这里tol（）括号里是个阈值，来控制数据计算算过程的一些条件的。 ...

唉我改了可大还是不行

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-22 15:13:20

zjying2000 发表于 2021-10-12 14:09
LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应 ...

不会算公式，这个文件可以解决吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-22 15:15:14

zjying2000 发表于 2021-10-12 14:09
LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应 ...

我不会算LMDI的公式，这个文件有具体步骤吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2022-1-22 19:04:34

不是山谷gzh 发表于 2022-1-22 15:15
我不会算LMDI的公式，这个文件有具体步骤吗

有的，有步骤，有代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-22 19:26:48

zjying2000 发表于 2022-1-22 19:04
有的，有步骤，有代码

谢谢谢谢！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-23 21:16:14

zjying2000 发表于 2021-10-12 14:09
LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应 ...

不是啊，我只是想简单的算个碳排放的影响因素，只需要算个公式就可以了，可是我数学不好不会算，但是用不到stata什么的，这个文件我不看不懂啊

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-23 21:32:34

zjying2000 发表于 2021-10-12 14:09
LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应 ...

能不能把文件的钱退给我呀，求求啦，这个什么state代码什么窝看不懂

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2022-1-24 13:00:37

不是山谷gzh 发表于 2022-1-23 21:32
能不能把文件的钱退给我呀，求求啦，这个什么state代码什么窝看不懂

哪有你这样的人？第一：我的帖子里有详细介绍、截图，相关的讨论也涉及了比较复杂的内容，这些你都看不懂吗，就想着下载？
第二：你已经下载了，也看过了，不说stata代码，就说LMDI的理论、步骤你也看了呀，看完说退了，再说也不贵，万一你以后要用或者给别人呢，又不可能完全从你的记忆里删掉

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-24 17:45:01

zjying2000 发表于 2022-1-24 13:00
哪有你这样的人？第一：我的帖子里有详细介绍、截图，相关的讨论也涉及了比较复杂的内容，这些你都看不懂 ...

LMDI我本来就会，我都说了我不会公式计算，我还问你了，你说有步骤，有的是那什么state步骤，

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-24 17:48:31

zjying2000 发表于 2022-1-24 13:00
哪有你这样的人？第一：我的帖子里有详细介绍、截图，相关的讨论也涉及了比较复杂的内容，这些你都看不懂 ...

我买东西不满意还可以退款呢，你那截图有不完整，回答我的也不完整，我只是不会计算公式，你文件里一堆什么state，根本就没有LMDI公式的计算步骤。你觉的十五块钱便宜，我可不觉的，都两顿饭钱了，你不给我退我就举报

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-24 17:57:27

zjying2000 发表于 2021-10-12 14:09
LMDI(对数平均迪氏指数法）目前在能源消耗、碳排放等领域应用很多，总体来说并不是一个很难的方法，但是相应 ...

你觉的不贵，哪你为啥不肯把钱退给我，我觉得贵啊，都两顿饭钱了，里面只有state的计算步骤，根本就没有LMDI分解的具体计算步骤。而且你文件里关于LMDI网上随处可见。你也别担心我以后会用到咋了的，我学地理的，根本就用不到什么stata软件

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-24 18:12:41

zjying2000 发表于 2022-1-24 13:00
哪有你这样的人？第一：我的帖子里有详细介绍、截图，相关的讨论也涉及了比较复杂的内容，这些你都看不懂 ...

我明确问你的说我不会算LMDI，而不是我不会LMDI，如果是在我没问你的情况下下载了文件，最后发现文件对我没用，那我认。现在是我明确问你了，你是有具体步骤，可是文件里只有stata步骤，并没有LMDI自己的计算步骤，那我为什么不能要求退钱

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zjying2000

2022-1-24 22:04:40

不是山谷gzh 发表于 2022-1-24 18:12
我明确问你的说我不会算LMDI，而不是我不会LMDI，如果是在我没问你的情况下下载了文件，最后发现文件对我 ...

那也不能说没有计算呀，明明有（你看下边是什么），难道这不是计算过程吗？而且LMDi也不难，你想想要什么计算呀，只是你自己没说清楚，也不能造谣呀，说里边只有stata代码，那下边这个是什么，为什么这样呀那你就举报吧，大不了我就不发了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

不是山谷gzh

2022-1-25 18:00:47

zjying2000 发表于 2022-1-24 22:04
那也不能说没有计算呀，明明有（你看下边是什么），难道这不是计算过程吗？而且LMDi也不难，你想想要什 ...

你这计算过程网上随便一搜就有了，算什么具体步骤。而且已经有好心的小姐姐给了我计算的步骤，你那钱我不要也罢，挣你的黑心钱去吧。还有需要LMDI计算具体步骤的可以联系我:2433619981，我免费分享！！！大家别被楼主骗了！！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zzz150594

2022-3-5 23:53:50

嘻嘻嘻茜茜发表于 2021-12-26 10:29
唉我改了可大还是不行

请问你的问题解决了吗，我也出现了这个问题

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

李佳佳2333

2022-3-12 11:13:40

楼上精彩的对话把我看呆了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群