|
此外,我们取ρ=Corr(ui,vi)∈0.5,0.95},和σu=1/P1-ρ(var[uiy2i,zi]=1)。为了表征潜在的仪器弱点,我们通过调整ζ的值来限制内生回归子y2i(y2i)与仪器zito be corr(y2i)之间的相关性,zi)=γ/nλ,在γ=1.5的情况下,我们考虑λ∈0.5的值网格,0.4,0.3,0.2,0.1}。由于DJ检验和标准弱IV检验的性能可能依赖于σz和σv,我们使用以下网格来模拟数据:σz={0.2,0.5,1,5,10}和σv={0.2,0.5,1,5,10}。对于每一个蒙特卡罗试验,我们将样本量取为n=500,5000,10000和考虑n=1000个蒙特卡罗复制之一。在每一个蒙特卡罗设计中,θ=(Eρ,α,β,π,ζ由CUGMM以一个单一的过度程度估计。我们选择仪器函数ai=a(y2i,zi)=(1,y2i,zi,zi,0,0)和bi=b(zi)=(0,0,0,0,1,zi)。DJ检验是按照第3.5节所述的程序进行的。使用5%的信号水平,我们拒绝了弱仪器不符合定理3的无效假设;也就是说,如果jδn>χ0.95(H+1-p),我们拒绝零值,其中在本例中H=6,p=5和χ0.95(H+1-p)=5.99。理论上,DJ检验的假设对应于:λ=0.5和λ<0.5的替代方案。然而,我们注意到,在样本中,很难只有λ决定线索的行为。鉴于此,为了将DJ检验的行为与传统的线性检验进行比较,我们引入了几组标准来评估样本中仪器弱点的潜在影响:线索的行为和相关的Wald统计量的大小失真。具体来说,我们计算偏差,标准偏差(S.D.)和相对均方根误差(rrmse)如下所示(以α为例)来衡量在以下设计下的估计性能:BIAS=bα-α,S.D.=Vuutnnxl=1(bαl-bα),rrmse=Vuutnnxl=1bαl-αα(31)为了计算简单,在蒙特卡罗模拟中,我们采用扰动δn=~ρ/log(log(n)),其中~ρ是每次蒙特卡罗复制中~ρ的CUGMM估计。本程序是3.5节中开发的数据驱动方法的简单版本。我们注意到每个测试的零假设略有不同:dj-h:λ=0.5;SS-Fn<10是一个信息无效假设;SY-用Cragg-Donald统计量,三重{ζ,σv,σz}使得2SLS相对偏差或Wald试验尺寸畸变大于agiven公差;稳健的测试将Nagar偏差超过基准的一小部分视为空。虽然弱仪器的定义对每个测试都是不同的,但它们的零假设在捕捉仪器弱的情况方面是一致的。其中bα=1/npnl=1bαl,bαl代表第L次蒙特卡罗CUGMM估计,α是truevalue。正如3.3和3.4节所证明的,在零的情况下,提示是一致的,而在替代的情况下,估计量将是一致的和渐近正态的,尽管是非标准的。与Stock和Yogo(2005)不同,他们选择2SLS与OLS的相对偏差作为检测弱仪器的一个标准,这里我们考虑偏差,S.D。而rrmse在(31)中被修改,原因如下。对于IV probit模型(30),CUE(和其他通常采用的Dessimation方法)没有封闭的表达。因此,在线性模型中潜在的IV弱下的通常的“biastowards OLS”概念在这种非线性情况下是不成立的,因为IV弱的潜在影响现在被模型的非线性特征复杂化了。在这种情况下,不能保证正负偏差不会相互影响,导致虚假的小总体偏差。
|