{smcl} {* *! version 1.0.0 07sep2013}{...} {hline} help for {hi: weakivtest} {hline} {title: Weak instrument test for TSLS and LIML} {title:Syntax} {p 8 14 2} {cmdab:weakivtest} [{cmd:,} level({it:#}) eps({it:#}) n2({it:#})] {title:Description} {p 4 4 2}{cmdab:weakivtest} implements the weak instrument test of Montiel Olea and Pflueger (2013). It is a postestimation command for {cmd:ivreg2} and {cmd:ivregress}. {p 4 4 2}{cmd:weakivtest} tests the null hypothesis of weak instruments for both Two-Stage Least Squares (TSLS) and Limited Information Maximum Likelihood (LIML) with one single endogenous regressor. The test rejects the null hypothesis when the effective F statistic exceeds a critical value, which depends on the estimator (TSLS or LIML), the significance level, and the desired weak instrument threshold tau.{cmd: weakivtest} extends the Stock and Yogo (2005) weak instrument tests available in {cmd:ivreg2} and in the {cmd: ivregress} postestimation command {cmd:estat firststage}. {p 4 4 2}Note: You must install {cmd:avar} by typing "ssc install avar" before running {cmd:weakivtest}. {title:Options} {p 4 4 2}{cmd:level(}{it:#}{cmd:)} specifies the confidence level. The default is {cmd:level(0.05)}. {p 4 4 2}{cmd:eps(}{it:#}{cmd:)} specifies the input parameter for the Nelder-Mead optimization technique; default is {cmd:eps(1e-3)}. {p 4 4 2}{cmd:n2(}{it:#}{cmd:)} specifies the denominator degrees of freedom for the inverse non-central F distribution. Set n2(#) to a large positive number to approximate an inverse non-central chi-squared distribution; default is {cmd:n2(1e7)}. {p 4 4 2}{cmd:weakivtest} estimates the variance-covariance matrix of errors as specified in the preceding {cmd:ivreg2} or {cmd:ivregress} command. The following options are supported: {p 8 4 2}{cmdab: robust} estimates an Eicker-Huber-White heteroskedasticity robust variance-covariance matrix. {p 8 4 2}{cmdab:cl:uster}{cmd:(}{it:varname}{cmd:)} estimates a variance-covariance matrix clustered by the specified variable. {p 8 4 2}{cmd:robust bw(}{it:#}{cmd:)} (for {cmd: ivreg2}) estimates a heteroskedasticity and autocorrelation-consistent variance-covariance matrix computed with a Bartlett (Newey-West) kernel with {it:#} bandwidth. {p 8 4 2}{cmd: bw(}{it:#}{cmd:)} without the {cmd:robust} option (for {cmd: ivreg2}) requests estimates that are autocorrelation-consistent but not heteroskedasticity-consistent. {p 8 4 2}{cmd:vce (hac nw} {it:#}{cmd:)} (for {cmd: ivregress}) estimates a heteroskedasticity and autocorrelation-consistent variance-covariance matrix computed with a Newey-West kernel with number of lags {it:#}. {marker s_examples}{title:Examples} {pstd} Change webuse URL {p_end} {phang2}{stata "webuse set http://fmwww.bc.edu/repec/bocode/d" : . webuse set http://fmwww.bc.edu/repec/bocode/d}{p_end} {pstd}Load Yogo (2004) data{p_end} {phang2}{stata "webuse Data_USAQ.dta, clear" : . webuse Data_USAQ.dta, clear}{p_end} {pstd} Reset URL to the default {p_end} {phang2}{stata "webuse set" : . webuse set, clear}{p_end} {phang2}{stata "gsort date" : . gsort date}{p_end} {phang2}{stata "tsset quarter" : . tsset quarter}{p_end} {pstd}Baseline example - robust Bartlett (Newey-West) kernel with bandwidth 7 .{p_end} {phang2}{stata "ivreg2 dc (rrf=z1 z2 z3 z4), robust bw(7) " : . ivreg2 dc (rrf=z1 z2 z3 z4), robust bw(7)}{p_end} {phang2}{stata "weakivtest" : . weakivtest }{p_end} {pstd}Implement {cmd:weakivtest} as a postestimation command for {cmd:ivregress}.{p_end} {phang2}{stata "ivregress 2sls dc (rrf=z1 z2 z3 z4), vce(hac nw 6)" : . ivregress 2sls dc (rrf=z1 z2 z3 z4), vce(hac nw 6)}{p_end} {phang2}{stata "weakivtest" : . weakivtest }{p_end} {pstd}Load National Longitudinal Survey data{p_end} {phang2}{stata "webuse nlswork, clear" : . webuse nlswork, clear}{p_end} {pstd}Implement {cmd:weakivtest} as a postestimation command for {cmd:ivreg2} in IV regressions with fixed effects.{p_end} {phang2}{stata "ivreg2 ln_w (ttl=age) grade collgrad i.year , ffirst cluster(idcode) " : . ivreg2 ln_w (ttl=age) grade collgrad i.year , ffirst cluster(idcode)}{p_end} {phang2}{stata "weakivtest" : . weakivtest }{p_end} {pstd}Implement {cmd:weakivtest} as a postestimation command for {cmd:reghdfe} in IV regressions with fixed effects.{p_end} {phang2}{stata "reghdfe ln_w grade collgrad (ttl=age) , ffirst absorb(year) cluster(idcode) old stage(first) " : . reghdfe ln_w grade collgrad (ttl=age) , ffirst absorb(year) cluster(idcode) old stage(first)}{p_end} {phang2}{stata "weakivtest" : . weakivtest }{p_end} {pstd}Implement {cmd:weakivtest} as a postestimation command for {cmd:ivreghdfe} in IV regressions with fixed effects.{p_end} {phang2}{stata "ivreghdfe ln_w (ttl=age) grade collgrad, first absorb(year) cluster(idcode) " : . ivreghdfe ln_w (ttl=age) grade collgrad, first absorb(year) cluster(idcode)}{p_end} {phang2}{stata "weakivtest" : . weakivtest }{p_end} {title:Saved results} {p 4 4 2}{cmd:weakivtest} saves the following results in {cmd:r()}: {p 4 4 2}Scalars {p_end} {col 4} {cmd:r(N)}{col 25} Number of Observations {col 4} {cmd:r(K)} {col 25} Number of Instruments {col 4} {cmd:r(n2)} {col 25} Denominator degrees of freedom non-central F {col 4} {cmd:r(level)} {col 25} Test Significance Level {col 4} {cmd:r(eps)} {col 25} Optimization Parameter {col 4} {cmd:r(F_eff)} {col 25} Effective F Statistic {col 4} {cmd:r(c_TSLS_5)} {col 25} TSLS Critical Value for tau{it:=}5{it:%} {col 4} {cmd:r(c_TSLS_10)} {col 25} TSLS Critical Value for tau{it:=}10{it:%} {col 4} {cmd:r(c_TSLS_20)} {col 25} TSLS Critical Value for tau{it:=}20{it:%} {col 4} {cmd:r(c_TSLS_30)} {col 25} TSLS Critical Value for tau{it:=}30{it:%} {col 4} {cmd:r(c_LIML_5)} {col 25} LIML Critical Value for tau{it:=}5{it:%} {col 4} {cmd:r(c_LIML_10)} {col 25} LIML Critical Value for tau{it:=}10{it:%} {col 4} {cmd:r(c_LIML_20)} {col 25} LIML Critical Value for tau{it:=}20{it:%} {col 4} {cmd:r(c_LIML_30)} {col 25} LIML Critical Value for tau{it:=}30{it:%} {col 4} {cmd:r(c_simp_5)} {col 25} TSLS Simplified Conservative Critical Value for tau{it:=}5{it:%} {col 4} {cmd:r(c_simp_10)} {col 25} TSLS Simplified Conservative Critical Value for tau{it:=}10{it:%} {col 4} {cmd:r(c_simp_20)} {col 25} TSLS Simplified Conservative Critical Value for tau{it:=}20{it:%} {col 4} {cmd:r(c_simp_30)} {col 25} TSLS Simplified Conservative Critical Value for tau{it:=}30{it:%} {col 4} {cmd:r(x_TSLS_5)} {col 25} TSLS Non-Centrality Parameter for tau{it:=}5{it:%} {col 4} {cmd:r(x_TSLS_10)} {col 25} TSLS Non-Centrality Parameter for tau{it:=}10{it:%} {col 4} {cmd:r(x_TSLS_20)} {col 25} TSLS Non-Centrality Parameter for tau{it:=}20{it:%} {col 4} {cmd:r(x_TSLS_30)} {col 25} TSLS Non-Centrality Parameter for tau{it:=}30{it:%} {col 4} {cmd:r(K_eff_TSLS_5)} {col 25} TSLS Effective Degrees of Freedom for tau{it:=}5{it:%} {col 4} {cmd:r(K_eff_TSLS_10)} {col 25} TSLS Effective Degrees of Freedom for tau{it:=}10{it:%} {col 4} {cmd:r(K_eff_TSLS_20)} {col 25} TSLS Effective Degrees of Freedom for tau{it:=}20{it:%} {col 4} {cmd:r(K_eff_TSLS_30)} {col 25} TSLS Effective Degrees of Freedom for tau{it:=}30{it:%} {col 4} {cmd:r(x_LIML_5)} {col 25} LIML Non-Centrality Parameter for tau{it:=}5{it:%} {col 4} {cmd:r(x_LIML_10)} {col 25} LIML Non-Centrality Parameter for tau{it:=}10{it:%} {col 4} {cmd:r(x_LIML_20)} {col 25} LIML Non-Centrality Parameter for tau{it:=}20{it:%} {col 4} {cmd:r(x_LIML_30)} {col 25} LIML Non-Centrality Parameter for tau{it:=}30{it:%} {col 4} {cmd:r(K_eff_LIML_5)} {col 25} LIML Effective Degrees of Freedom for tau{it:=}5{it:%} {col 4} {cmd:r(K_eff_LIML_10)} {col 25} LIML Effective Degrees of Freedom for tau{it:=}10{it:%} {col 4} {cmd:r(K_eff_LIML_20)} {col 25} LIML Effective Degrees of Freedom for tau{it:=}20{it:%} {col 4} {cmd:r(K_eff_LIML_30)} {col 25} LIML Effective Degrees of Freedom for tau{it:=}30{it:%} {col 4} {cmd:r(x_simp_5)} {col 25} TSLS Simplified Non-Centrality Parameter for tau{it:=}5{it:%} {col 4} {cmd:r(x_simp_10)} {col 25} TSLS Simplified Non-Centrality Parameter for tau{it:=}10{it:%} {col 4} {cmd:r(x_simp_20)} {col 25} TSLS Simplified Non-Centrality Parameter for tau{it:=}20{it:%} {col 4} {cmd:r(x_simp_30)} {col 25} TSLS Simplified Non-Centrality Parameter for tau{it:=}30{it:%} {col 4} {cmd:r(K_eff_simp_5)} {col 25} TSLS Simplified Effective Degrees of Freedom for tau{it:=}5{it:%} {col 4} {cmd:r(K_eff_simp_10)} {col 25} TSLS Simplified Effective Degrees of Freedom for tau{it:=}10{it:%} {col 4} {cmd:r(K_eff_simp_20)} {col 25} TSLS Simplified Effective Degrees of Freedom for tau{it:=}20{it:%} {col 4} {cmd:r(K_eff_simp_30)} {col 25} TSLS Simplified Effective Degrees of Freedom for tau{it:=}30{it:%} {marker references}{...} {title:References} {marker Baum2007}{...} {phang} Baum, C. F., M. E. Schaffer, and S. Stillman. 2007. Enhanced routines for instrumental variables/generalized method of moments and testing. {it:Stata Journal} 7:465-506. {marker Baum2010}{...} {phang} Baum, C. F., M. E. Schaffer, and S. Stillman. 2010. IVREG2: Stata module for extended instrumental variables/2SLS and GMM estimation. {browse "http://ideas.repec.org/c/boc/bocode/s425401.html"}. {marker C2003}{...} {phang} Campbell, J. Y. 2003. Consumption-based asset pricing. {it:Handbook of the Economics of Finance, Vol. 1} 803-887. {marker MOP2013}{...} {phang} Montiel Olea, J. L. and C. E. Pflueger. 2013. A robust test for weak instruments. {it:Journal of Business and Economic Statistics} 31:358-369. {marker MOPW2013}{...} {phang} Pflueger, C. E. and Su Wang. 2013. A robust test for weak instruments in Stata. {browse "http://papers.ssrn.com/sol3/papers.cfm?abstract_id=2323012"}. {marker NW1987}{...} {phang} Newey, W. and K. D. West. 1987. A simple, positive semi-definite, heteroskedasticity and autocorrelation consistent covariance matrix. {it:Econometrica} 55:703-708. {marker StockYogo2005}{...} {phang} Stock, J. and M. Yogo. 2005. Testing for weak instruments in linear IV regression. {it: In Identification and Inference for Econometric Models: Essays in Honor of Thomas Rothenberg} Chapter 5 80-108. {marker Yogo2004}{...} {phang} Yogo, M. 2004. Estimating the elasticity of intertemporal substitution when instruments are weak. {it:Review of Economics and Statistics} 86:797-810. {p_end} {marker authors}{...} {title:Authors} Carolin E. Pflueger, University of Chicago, Chicago IL, 60637, cpflueger@uchicago.edu Su Wang, ShanghaiTech University, Shanghai, 201210, China, wangsu@shanghaitech.edu.cn The package is currently maintained by Carolin Pflueger, Su Wang and Luis Yepez (luisyepezsa@uchicago.edu). We thank Peter Hull for pointing out and fixing an error with the weights. We thank Johannes Stroebel and Drew Johnston for help with reghdfe. All remaining errors are our own. {marker also}{...} {title:Also see} {p 4 4 2}{help ivregress}, {help ivreg}, {help ivregress_postestimation}, {help ivreg2}, {help reghdfe}, {help ivreghdfe} (if installed), {help avar} (if installed) {p_end}