谢邀!这个问题不难很基础。我用R语言从下面几个方面来回答这个问题:

一、two step GMM, iterated GMM, continuous updated GMM的具体算法;

二、你给的这个矩条件与LS之间的关系;

三、广义经验似然估计;

四、估计求解;

五、Monte Carlo 模拟研究

一、三种算法

众所周知GMM的矩条件是

,其中

,我们的目的就是估计出

。Hansen(1982)提出一种GMM的two step 迭代估计方法,在R里面gmm包(Pierre Causse,2010)的算法如下:

这就是two step GMM的估计方法,尽管这个估计是一致的且有效,但是也有缺点。在处理内生性问题的时候,如果工具变量选的不好那么这个估计效率很低方差很大;此外,正如很多研究指出的(如Newey&Smith,2004; etc.)这个估计量的小样本性质不太好(有偏)。对此,Hansen et al.(1996)针对这个问题,又提出了改进的两种方法:iterated GMM

以及一种高度非线性的迭代 continuous updated GMM

事实上,这两种方法是等价的,且continuous updated GMM估计量的渐进偏差更小。当然这只是大样本情况下的渐进性质,针对你给的这个问题我觉得小样本情况下需要进行具体分析。

二、与LS 之间的关系

众所周知,GMM是LS的generalization,那么你给的这个GMM下的矩条件与LS之间有什么关系呢?事实上,不妨令矩条件为

,其中

。两边取对数后得到

,注意到当样本容量N趋于无穷的时候,我们可以得到

。将此式子带入上面的等式里,得到

看到这是不是感到很熟悉,没错这不就是一般的线性回归嘛,最左边是Y,系数是

。这样我们就把非线性矩条件转化为了一个线性的方程,我们的目的就是用LS求解出

。注意到这是一个受约束的最小二乘问题:

可以看出我们对参数

施加了一个等式约束来满足你给出的矩条件。利用Langrange乘子法可以求出参数,具体的计算我就不展开了,对乘子

求导后得到

发现,其实这就是矩条件:

!由此可见你给的矩条件下GMM与LS之间的等价关系。

三、广义经验似然估计

GMM固然是个好东西,但是正如第一节所述,小样本下面对弱工具变量等情况下,GMM 估计量往往有偏且无效。对此还有一种广义的经验似然方法(GEL)往往比GMM在统计上更优良,

GEL是经验似然估计方法的推广,有很多优良的性质,只不过没有被更多的计量经济学家或实证研究者关注,具体的就不说了,你可以找OWEN他们的文献来读。后面我将用GMM方法与GEL方法来估计你给的那个参数。

四、估计

你没有r是多少,我就用你给的

数据,且假设

。利用R包gmm就可以直接求解了:

library(gmm)

g2

f

return(f)

}

g2是你给的矩条件,以及

g1

f

return(f)}

g1是LS下的

把数据带入:

x1

x2

x

gmm(g2,x,0,wmatrix = c("optimal"),optfct="nlminb")

gmm(g1,x,0,wmatrix = c("optimal"),optfct="nlminb")

gel(g2,x,0,optfct="nlminb")

gel(g1,x,0,optfct="nlminb")

最后的估计结果是

可以看出四个估计结果都一样:

,GMM与GEL 的计算一样,且你给的矩条件与LS推导下矩条件计算结果一致,证明了上面的推导。但是值得注意的是利用GEL对矩条件g1进行估计得到的估计量方差最小。

五、MC模拟

现在主要针对你给的这个矩条件,研究下GMM与GEL估计的小样本性质。考虑如下DGP:

从均匀分布里面抽出50个随机数,真实的参数值

。采用上述四个计算,得到结果

发现GMM与GEL估计的结果都是0.5,与真实值不差毫厘。说明即使样很小,针对你给的矩条件,GMM与GEL估计的性质都还不错。当然你还可以做一个bootstrap 去看看GMM,GEL估计的抽样分布。哦对了,我用的GMM的计算方法还是two step 的,你可以换做iterated GMM 与CUE试试。

R的gmm包使用指南:CRAN - Package gmm​cran.r-project.orghttps://cran.r-project.org/web/packages/gmm/gmm.pdf​cran.r-project.orghttps://cran.r-project.org/web/packages/gmm/vignettes/gmm_with_R.pdf​cran.r-project.org

当然STATA也可以做,跟R差不多,具体的可以参考https://core.ac.uk/download/pdf/6306768.pdf​core.ac.ukhttps://core.ac.uk/download/pdf/6490419.pdf​core.ac.uk

references

Chaussé P. Computing Generalized Empirical Likelihood and. Generalized Method

of Moments with R[J]. Journal of Statistical Software, 2010,34(11):1-35.

NeweyW K, Smith R J. Higher Order Properties of GMM and Generalized Empirical Likelihood Estimators[J]. Econometrica, 2004, 72(1):219-255.

Anatolyev S. GMM, GEL, Serial Correlation, and Asymptotic Bias[J]. Econometrica, 2005, 73(3):983-1002.

Hansen L P. Large Sample Properties of Generalized Method of Moments Estimators[J]. Econometrica, 1982, 50(4):1029-1054.

Hansen L P, Singleton K J. Generalized Instrumental Variables Estimation of

Nonlinear Rational Expectations Models[J]. Econometrica, 1982,

50(5):1269-1286.

Hansen L P, Heaton J, Yaron A. Finite-Sample Properties of Some Alternative

GMM Estimators[J]. Journal of Business & Economic Statistics, 1996,

14(3):262-280.

python广义矩估计_怎么用软件做广义矩估计GMM的参数估计?相关推荐

  1. 广义典型相关分析_一文教你掌握广义估计方程

    广义估计方程是一种研究纵向数据(比如重复测量数据,面板数据)的方法. 同一测量对象的多次测量数据结果之间很可能有着相关关系,如果不考虑数据之间的相关性会造成信息损失.常见的研究模型(比如线性回归)都要 ...

  2. 广义典型相关分析_数学建模/机器学习:广义加性模型(GAM)及其Python实现

    笔者做过国赛也做过美赛,其中一类典型问题就是分析相关性,从而进行预测或者其他操作.这类问题通常情况下属于比较常规的问题,一般通过matlab或SPSS分析相关性,得到一个较好的数值即可. 然而有的时候 ...

  3. 广义pareto分布_高记录值下广义Pareto分布参数的随机估计及垂直概率密度检验

    The Randomized Estimate and Vertical Density Representation test for the generalized Pareto distribu ...

  4. python拦截广告弹窗_会捆绑软件、弹出广告?Python教你实现纯净中英自动翻译小工具...

    本文原创,未经允许不得转载,本文版权归@itcast @wj14k所有,感谢你的阅读! 翻译工具在安装会捆绑其他软件 翻译工具会跳出你不喜欢的广告 -- 对捆绑.广告说不! Python实现纯净翻译小 ...

  5. python 文件分类清单_时间管理软件的清单如何对文件夹与清单进行分类(如,滴答清单)?...

    滴答清单用了很多年了,先上图. 一一说来 Bank 当然是最重要的,虽然只是起名叫bank,其实是一切和财务相关的提醒,比如房贷.信用卡.理财.股票等等. Health 是与日常健康管理以及习惯有关的 ...

  6. 女生学python工作累吗_数据分析师女生做累么 工资一般多少

    数据分析师女生做累么,工资一般多少,小编整理了相关信息,希望会对大家有所帮助! 数据分析师女生做累不累 女生还是很适合做数据分析的,数据分析师因为敲的代码少,相比起天天敲代码的职业更适合女生一些,没那 ...

  7. python与建筑结构设计_新手进入建筑设计院做结构设计,主要看哪些书籍?

    新手的话,有很多书可以看.我毕业的时候,把跟专业无关的.没什么阅读或收藏价值的书都处理了,剩下的所有的专业书,都带到了设计院.刚参加工作的时候,也一直在买书.后来把这些书从设计院搬回家的时候,真是累煞 ...

  8. cvpr 深度估计_无监督单目视频深度估计中的uncertainty方法(CVPR#x27;20)

    Contribution 对11种使用uncertainty方法的全面评估 深度挖掘uncertainty对depth estimation起到的作用 提出一个新颖的self-teaching方法去m ...

  9. python分词统计词频_基于结巴分词做的全文分词统计词频小脚本

    受朋友之托,写一个小脚本,断断续续做了两天,写一下两天的收获. 起因 有个朋友说专业文档很枯燥难懂,需要一个能把全文的关键词找出来并排序的东西,找不到现成的,问我能不能做一个.我前些天也听车神说有关分 ...

  10. 一阶广义差分模型_向后差分公式法、广义 α 法和龙格-库塔法

    解决方法 瞬态求解器提供三种不同的时间步进方法:隐式向后差分公式.广义 α 方法以及显式龙格-库塔系列方法. "向后差分公式"(BDF) 求解器是一种使用向后差分公式的隐式求解器, ...

最新文章

  1. 从芯片到系统:FPGA加速卡的发展历程与展望
  2. 亚马逊首家“无人超市”系统存在bug?!开业当天,记者中途换装成功骗过摄像头...
  3. Java读取String分行字符串
  4. ICML 2021 | AlphaNet:基于α-散度的超网络训练方法
  5. 死磕 Java 8 的日期处理
  6. Vue菜单栏边框参差不齐,修整完美
  7. apache phoenix的JAVA JDBC操作
  8. kafka权威指南-总结
  9. [Maven实战-许晓斌]-[第二章]-2.5 Eclipse上面安装Maven插件 m2eclipse
  10. 关于chm电子书无法显示网页的解决方
  11. 推荐几个选择时间.日期的免费插件 js
  12. w ndows10图标,win10桌面图标变白怎么解决
  13. 使用XMind编写测试用例
  14. Chrome安装程序遇到错误 0xe0000008解决办法
  15. 数据人故事——【五花肉】7年数据产品的职业成长之路
  16. .net6API使用AutoMapper和DTO
  17. 1D/1D动态规划学习总结
  18. 贝叶斯分类方法-例题
  19. 赔1100万美元!谷歌招聘年龄歧视
  20. Firefox扩展初接触

热门文章

  1. win7关于无线连接的服务器,win7网络连接受限
  2. 华为云 云学院 白嫖
  3. 如何练就DOTA绝世神功宝典(转)
  4. java初学者看什么书比较好
  5. 无产权证的房产能让渡吗?
  6. 串行进位加法器电路图(logisim实现)
  7. Linux Netfilter/NAT的两个典型问题
  8. 概率论 方差公式_概率论学习笔记(6)
  9. Linux修改时间 修改时区 | Linux时间校准
  10. jQuery获取兄弟元素