[Sebastian Kripfganz]

Arkangel Cordero

There are 3 issues here:
1) When replicating the xtdpdgmm results with ivreg2, you need to specify the untransformed variables.
2) You need to specify all untransformed regressors as endogenous, including time dummies.
3) ivreg2 applies collinearity checks that drop some of the instruments, even though it should not. You need to specify the nocollin option with ivreg2.

webuse abdata
xtdpdgmm L(0/1).n, model(diff) gmm(L.n, l(2 5)) nocons vce(cluster id) teffects nolevel twostep
quietly predict iv*, iv
ivreg2 n (L1.n i.year = iv*), noconstant gmm2s cluster(id) nocollin

[Arkangel Cordero]

Dear Professor @Sebastian Kripfganz,

To provide some more context. I am interested in estimating a "conventional" difference gmm model (as in your 1.d. in your Post # 481) . I have been able to replicate the results of your xtdpdgmm command with with "teffects" and "year-dummies" and xtabond2 with year-dummies as below.

webuse abdata
xtdpdgmm L(0/1).n ,                    model(diff) gmm(L.n , l(2 5))  nocons   vce(cluster id)  teffects            nolevel      twostep overid
estimate store m_xtdpdgmm_teffects
estat overid

xtdpdgmm L(0/1).n  (yr1978 - yr1984),  model(diff) gmm(L.n , l(2 5))  nocons   vce(cluster id)  iv(yr1978 - yr1984) nolevel      twostep overid
estimate store m_xtdpdgmm_timeDum
estat overid

xtabond2 L(0/1).n (yr1978 - yr1984), gmm(L.n , lag(2 5)) iv(yr1978 - yr1984) nocons nolevel cluster(id) twostep ar(6)
estimate store m_xtbound2_timeDum

esttab m_xtdpdgmm_teffects m_xtdpdgmm_timeDum m_xtbound2_timeDum

Results:

------------------------------------------------------------
                      (1)             (2)             (3)  
                        n               n               n  
------------------------------------------------------------
L.n                 0.698**         0.698**         0.698**
                   (2.99)          (2.99)          (2.99)  

1978.year         -0.0103                                  
                  (-0.84)                                  

1979.year         -0.0107                                  
                  (-0.71)                                  

1980.year         -0.0528**                                
                  (-3.19)                                  

1981.year          -0.156***                                
                  (-7.87)                                  

1982.year          -0.171***                                
                  (-4.72)                                  

1983.year          -0.155*                                  
                  (-2.47)                                  

1984.year          -0.153*                                  
                  (-2.21)                                  

yr1978                            -0.0103         -0.0103  
                                  (-0.84)         (-0.84)  

yr1979                            -0.0107         -0.0107  
                                  (-0.71)         (-0.71)  

yr1980                            -0.0528**       -0.0528**
                                  (-3.19)         (-3.19)  

yr1981                             -0.156***       -0.156***
                                  (-7.87)         (-7.87)  

yr1982                             -0.171***       -0.171***
                                  (-4.72)         (-4.72)  

yr1983                             -0.155*         -0.155*  
                                  (-2.47)         (-2.47)  

yr1984                             -0.153*         -0.153*  
                                  (-2.21)         (-2.21)  
------------------------------------------------------------
N                     891             891             751  
------------------------------------------------------------
t statistics in parentheses
* p<0.05, ** p<0.01, *** p<0.001

What I am trying to do is replicate these results with "ivreg2" (as per your 2019 presentation in slides 39 -42) in order to obtain the under-identification and weak-identification tests available in the latter. I have tried the following, without success.

xtdpdgmm L(0/1).n ,   model(diff) gmm(L.n , l(2 5))  nocons   vce(cluster id)  teffects    nolevel   twostep overid
estimate store m_xtdpdgmm_teffects
estat overid
drop iv*
quietly predict iv*, iv

ivreg2 n iv19-iv25 (L1.n   = iv1 - iv18 iv19-iv25),  gmm2s noconstant  first  cluster(id) endog(L1.n)
estimate store m_ivreg2_lev1

ivreg2 n iv19-iv25 (L1.n   = iv1 - iv18), partial(iv19-iv25) gmm2s noconstant  first  cluster(id) endog(L1.n)
estimate store m_ivreg2_lev2

ivreg2 D1.n D1.(iv19-iv25) (D1.L1.n   = D1.(iv1 - iv18 iv19-iv25)),  gmm2s noconstant  first  cluster(id) endog(D1.L1.n)
estimate store m_ivreg2_dif1

ivreg2 D1.n D1.(iv19-iv25) (D1.L1.n   = D1.(iv1 - iv18)), partial(D1.(iv19-iv25)) gmm2s noconstant  first  cluster(id) endog(D1.L1.n)
estimate store m_ivreg2_dif2


esttab  m_ivreg2_lev1  m_ivreg2_dif1  m_ivreg2_lev2  m_ivreg2_dif2

With the following results:

--------------------------------------------------------------------------------------------
                      (1)             (2)             (3)             (4)             (5)  
                        n               n             D.n               n             D.n  
--------------------------------------------------------------------------------------------
L.n                 0.698**         0.953***                        0.953***                
                   (2.99)          (9.75)                          (9.75)                  

LD.n                                                0.172                           0.172  
                                                   (1.88)                          (1.88)  

1978.year         -0.0103                                                                  
                  (-0.84)                                                                  

1979.year         -0.0107                                                                  
                  (-0.71)                                                                  

1980.year         -0.0528**                                                                
                  (-3.19)                                                                  

1981.year          -0.156***                                                                
                  (-7.87)                                                                  

1982.year          -0.171***                                                                
                  (-4.72)                                                                  

1983.year          -0.155*                                                                  
                  (-2.47)                                                                  

1984.year          -0.153*                                                                  
                  (-2.21)                                                                  

iv19                              -0.0488***                                                
                                  (-5.18)                                                  

iv20                              -0.0847***                                                
                                  (-8.13)                                                  

iv21                               -0.137***                                                
                                  (-9.34)                                                  

iv22                               -0.149***                                                
                                  (-7.47)                                                  

iv23                              -0.0684***                                                
                                  (-3.44)                                                  

iv24                             -0.00832                                                  
                                  (-0.46)                                                  

iv25                             0.000355                                                  
                                   (0.03)                                                  

D.iv19                                            -0.0348***                                
                                                  (-6.05)                                  

D.iv20                                            -0.0749***                                
                                                  (-7.22)                                  

D.iv21                                             -0.143***                                
                                                  (-8.58)                                  

D.iv22                                             -0.191***                                
                                                  (-9.15)                                  

D.iv23                                             -0.143***                                
                                                  (-7.71)                                  

D.iv24                                            -0.0736***                                
                                                  (-4.95)                                  

D.iv25                                            -0.0311**                                
                                                  (-3.07)                                  
--------------------------------------------------------------------------------------------
N                     891             891             751             891             751  
--------------------------------------------------------------------------------------------
t statistics in parentheses
* p<0.05, ** p<0.01, *** p<0.001

I have only been able to replicate the results of an "unconventional" difference gmm model ( as in 3. in your Post # 482), i.e. one in which the time dummies are instruments for the untransformed model, as shown below.

xtdpdgmm L(0/1).n ,                    model(diff) gmm(L.n , l(2 5))  nocons   vce(cluster id)  teffects                  twostep overid
estimate store m_xtdpdgmm_teffects
estat overid
drop iv*
quietly predict iv*, iv 

ivreg2 n iv19-iv26 (L1.n   = iv1 - iv18 iv19-iv26),  gmm2s noconstant  first  cluster(id) endog(L1.n) 
estimate store m_ivreg2_lev1

ivreg2 n iv19-iv26 (L1.n   = iv1 - iv18), partial(iv19-iv26) gmm2s noconstant  first  cluster(id) endog(L1.n) 
estimate store m_ivreg2_lev2

esttab m_xtdpdgmm_teffects m_ivreg2_lev1  m_ivreg2_lev2  


------------------------------------------------------------
                      (1)             (2)             (3)  
                        n               n               n  
------------------------------------------------------------
L.n                 0.946***        0.946***        0.946***
                  (20.75)         (37.89)         (37.89)  

1977.year          0.0885                                  
                   (1.63)                                  

1978.year          0.0732                                  
                   (1.45)                                  

1979.year          0.0674                                  
                   (1.35)                                  

1980.year          0.0236                                  
                   (0.47)                                  

1981.year         -0.0704                                  
                  (-1.37)                                  

1982.year         -0.0563                                  
                  (-1.27)                                  

1983.year         -0.0122                                  
                  (-0.33)                                  

1984.year         -0.0172                                  
                  (-0.59)                                  

iv19                               0.0885**                
                                   (2.70)                  

iv20                               0.0732*                  
                                   (2.48)                  

iv21                               0.0674*                  
                                   (2.31)                  

iv22                               0.0236                  
                                   (0.81)                  

iv23                              -0.0704*                  
                                  (-2.33)                  

iv24                              -0.0563*                  
                                  (-2.15)                  

iv25                              -0.0122                  
                                  (-0.57)                  

iv26                              -0.0172                  
                                  (-0.87)                  
------------------------------------------------------------
N                     891             891             891  
------------------------------------------------------------
t statistics in parentheses
* p<0.05, ** p<0.01, *** p<0.001

I wonder if you could please provide some guidance in being able to replicate the results for a "conventional" (as in your 1.d. in your Post # 481, i.e., with the time-dummies as instruments in the transformed--not the level--model) difference gmm model with ivreg2?

Thank you in advance for any help in this regard.

[Arkangel Cordero]

Dear professor Kripfganz,

I am trying to replicate the results of a "conventional" (i.e., one in which the time-dummies are used as instruments in the difference, not the levels, equation) difference gmm model estimated with xtdpdgmm (as per 1.D) in Post # 481) with ivreg2 in order to access the under-identification and weak-instruments tests available through the latter (similar to what you explain in slides 39-42 of your Sept. 2019 presentation). I am able to reproduce the results with an "unconventional" difference gmm model (as 3) in Post # 482). However, I have been unable to do so with the conventional difference gmm. Can you please provide some guidance.

I have tried:

webuse abdata
xtdpdgmm L(0/1).n ,                    model(diff) gmm(L.n , l(2 5))  nocons   vce(cluster id)  teffects                     nolevel      twostep overid
estimate store model_xtdpdgmm_teffects
estat overid
quietly predict iv*, iv 

ivreg2 n iv19-iv25 (L1.n   = iv1 - iv18), partial(iv19-iv25) gmm2s noconstant  first  cluster(id) endog(L1.n) 
estimate store model_ivreg2_lev

ivreg2 D1.n D1.(iv19-iv25) (D1.L1.n   = D1.(iv1 - iv18 iv19-iv25)), gmm2s noconstant  first  cluster(id) endog(D1.L1.n) 
estimate store model_ivreg2_dif


. esttab  model_xtdpdgmm_teffects model_ivreg2_lev model_ivreg2_dif

------------------------------------------------------------
                      (1)             (2)             (3)   
                        n               n             D.n   
------------------------------------------------------------
L.n                 0.698**         0.953***                
                   (2.99)          (9.75)                   

LD.n                                                0.172   
                                                   (1.88)   

1978.year         -0.0103                                   
                  (-0.84)                                   

1979.year         -0.0107                                   
                  (-0.71)                                   

1980.year         -0.0528**                                 
                  (-3.19)                                   

1981.year          -0.156***                                
                  (-7.87)                                   

1982.year          -0.171***                                
                  (-4.72)                                   

1983.year          -0.155*                                  
                  (-2.47)                                   

1984.year          -0.153*                                  
                  (-2.21)                                   

D.iv19                                            -0.0348***
                                                  (-6.05)   

D.iv20                                            -0.0749***
                                                  (-7.22)   

D.iv21                                             -0.143***
                                                  (-8.58)   

D.iv22                                             -0.191***
                                                  (-9.15)   

D.iv23                                             -0.143***
                                                  (-7.71)   

D.iv24                                            -0.0736***
                                                  (-4.95)   

D.iv25                                            -0.0311** 
                                                  (-3.07)   
------------------------------------------------------------
N                     891             891             751   
------------------------------------------------------------
t statistics in parentheses
* p<0.05, ** p<0.01, *** p<0.001

Any guidance is much appreciated.

[Sebastian Kripfganz]

A substantial difference between those two test statistics suggests that the weighting matrix might be imprecisely estimated. A solution might be to use the iterated GMM estimator. But it could also be a symptom of a relatively small sample size or weak instruments.

[Nursena Sagir]

Sebastian Kripfganz Dear Sebastian, I have a question to you about Sargan-Hansen test. How should I interpret diverged 2-step and 3-step weighting matrix results like statistics below?


2-step moment functions, 2-step weighting matrix chi2(2) = 3.1408
Prob > chi2 = 0.2080

2-step moment functions, 3-step weighting matrix chi2(2) = 11.5937
Prob > chi2 = 0.0030

Thanks in advance!

Best regards,
Nursena

[Sebastian Kripfganz]

Also, I struggle to find an initial candidate model that passes the specification tests. I assume this is due to the very large N and that the specification tests can already detect relatively small deviations from the null hypotheses.

This would normally be my answer, yes.

If statistical tests do not seem to be helpful, you might have to resort to economic theory as a guide for model specification / variable classification.

Another possibility might be to consider smaller subsets of the data where the Hansen test has less power. This might sound odd, but if the Hansen test does not detect small deviations from the null hypothesis anymore, it might be helpful to single out models with more severe model specifications.

Forward-orthogonal deviations indeed seem to be appropriate given the unbalanced nature of the panel.

With such a large sample size, the two-step estimator might deliver substantial efficiency gains and is therefore highly recommended.

[Jupp Peters]

Dear Sebastian Kripfganz,

thank you so much for your help. I have one more question that I hope may add to the collection of advice you have already given in this thread.

I struggle to find a model specification. This is a summary of my sample:

• N = ~25.000 households
• T = 12 years
• ~150.000 observations
• Unbalanced panel:
  • Minimum observations per household: 1
  • Maximum observations per household: 12

The goal is to estimate the relationship between consumption of a good (dependent variable) and (1) price of this good, (2) household income, and (3) various household covariates such as household size, employment status of household head, dwelling type, ... + federal state and time dummies. Some of the household covariates are binary or categorical variables. I am also interested in the interaction effect of price and household income.

I tested fixed and random effects models as benchmarks with the following specifications:

L(0/1)y c.x_price##c.x_income X_covariates time_effects

I would like to have the same combinations of regressors in the GMM model to compare the results with the fixed and random effects models.

I started the model selection process. However, it seems impractical in my case as a single estimation run lasts ~20 minutes. Also, I struggle to find an initial candidate model that passes the specification tests. I assume this is due to the very large N and that the specification tests can already detect relatively small deviations from the null hypotheses.

My questions are:

• Do you have any advice on how to configure IV lags and how to specify the variable types (endogenous/predetermined/exogenous) for the different variables?
• Would you agree that I should use forward orthogonal deviations because of the unbalanced panel?
• Would you agree that I should use the two-step estimator because of the large sample size?

I really appreciate any help you can provide.

Jupp

[Sebastian Kripfganz]

Zeenat, Rupali: Apologies; I do not currently have the time to respond to longer queries.

Jupp: If this regressor is not one of your main regressors of interest, you could remove it entirely if all coefficients have high p-values. In principle, you could retain the respective instruments if you believe that they are strong instruments for the other regressors; otherwise, you can remove them as well.

[Jupp Peters]

Dear Sebastian Kripfganz,

I have a question regarding the Sequential model selection process you presented at the 2019 London Stata Conference.

In step 3, you recommend to "remove lags or interaction effects with (very) high p-values in individual or joint significance tests". What if, after removing all lags of one regressor due to their high individual p-values, the remaining contemporaneous effect of that regressor still has a high p-value? Should it then be entirely removed? If yes, what should be done with the respective IVs?

I really appreciate any help you can provide.

[Rupali Vashisht]

Dear Prof. Sebastian Kripfganz,

Thank you very much for your previous response. I am using xtdpdgmm and the gmm technique for the first time and I would really appreciate your help in understanding if I am doing it correctly.

In my model, I consider EPS and ROAA as endogenous. The predetermined variables are - l.ROAA l.EPS TDtoTA SalesGrowth lnTA CFOtoSales CapitalIntensity l.EnvDisclScore. And, the exogenous variables include the time dummies.

I intend to incorporate fixed effects by first differences transformation. And, I wish to take care of the correlation between the error and the lagged EPS by considering their lags as their instruments.

I am using the following code:

xtdpdgmm L(0/1).EPS L(0/1).ROAA TDtoTA SalesGrowth lnTA l.EnvDisclScore i.t, model(diff) collapse gmm(EPS ROAA, lag(2 3)) gmm(EPS ROAA, lag(1 1) diff model(level)) gmm(l.ROAA l.EPS TDtoTA SalesGrowth lnTA CFOtoSales CapitalIntensity l.EnvDisclScore, lag(1 3)) gmm(l.ROAA l.EPS TDtoTA SalesGrowth lnTA CFOtoSales CapitalIntensity l.EnvDisclScore i.t, lag(0 0) diff model(level)) iv(i.t) two vce(cluster id) small overid

And I obtain the following output:

		WC-Robust
EPS	Coefficient	std. err.	t	P>t	[95% conf.	interval]
EPS
L1.	1.009836	0.0749899	13.47	0	0.8620847	1.157588
ROAA	0.2411051	0.1097758	2.2	0.029	0.0248154	0.4573949
L1.ROAA	0.0063018	0.0585874	0.11	0.914	-0.1091322	0.1217357
TDtoTA	0.0151328	0.0461458	0.33	0.743	-0.0757876	0.1060532
SalesGrowth	-0.0158857	0.0163982	-0.97	0.334	-0.0481948	0.0164234
lnTA	0.9796653	1.435696	0.68	0.496	-1.849068	3.808398
EnvDisclScore
L1.	0.0485014	0.0524756	0.92	0.356	-0.0548905	0.1518933
t
2011	0	(empty)
2012	1.378857	1.268998	1.09	0.278	-1.121434	3.879147
2013	-0.3298332	1.159062	-0.28	0.776	-2.613518	1.953851
2014	-1.774592	1.015847	-1.75	0.082	-3.776101	0.2269177
2015	0.5190836	1.136322	0.46	0.648	-1.719796	2.757963
2016	0.6554844	0.9437869	0.69	0.488	-1.204046	2.515015
2017	-0.2927053	0.8384453	-0.35	0.727	-1.944683	1.359272
2018	-2.61427	1.03919	-2.52	0.013	-4.661772	-0.5667682
2019	-0.5594183	0.8990458	-0.62	0.534	-2.330796	1.21196
2020	0	(omitted)
_cons	-11.50282	14.16942	-0.81	0.418	-39.42064	16.415

Instruments corresponding to the linear moment conditions:
1, model(diff):
L2.EPS L3.EPS L2.ROAA L3.ROAA
2, model(level):
L1.D.EPS L1.D.ROAA
3, model(diff):
L3.L.ROAA L3.L.EPS L1.TDtoTA L2.TDtoTA L3.TDtoTA L1.SalesGrowth
L2.SalesGrowth L3.SalesGrowth L1.lnTA L2.lnTA L3.lnTA L1.CFOtoSales
L2.CFOtoSales L3.CFOtoSales L1.CapitalIntensity L2.CapitalIntensity
L3.CapitalIntensity L1.L.EnvDisclScore L2.L.EnvDisclScore
L3.L.EnvDisclScore
4, model(level):
D.TDtoTA D.SalesGrowth D.lnTA D.CFOtoSales D.CapitalIntensity
D.L.EnvDisclScore
5, model(diff):
2013bn.t 2014.t 2015.t 2016.t 2017.t 2018.t 2019.t 2020.t
6, model(level):
_cons

I just want to make sure if the code above does what I want to do. I have the following doubts, I would greatly appreciate your help on these:

1. Does the model above specify the endogenous, predetermined and exogeneous variables correctly?
2. I have introduced ‘CFOtoSales’ and ‘CapitalIntensity’ variables in my set of instruments as external instruments as I wish to instrument for ROAA (which I consider endogenous). Am I correct in doing that?
3. When I use iv(i.t) – am I right in thinking that i.t acts as an instrument for all variables?
4. Is my model first differenced?
5. The text below the table shows that the lagged values of the endogenous variables act as instruments of the differenced model and that the level model has differenced variables as instruments (which I have specified) - is that the correct way to go about it?

Apart from the above questions, Could you please help me understand what is the benefit of specifying instruments for both level and difference models?

Thank you very much in advance!

Kind regards,
Rupali

[Zeenat Murtaza]

Hi Sebastian,

I hope you are doing well. Sebastian, I have few concerns regarding system gmm and marginal plots contrast effects and their interpretation. I used the following model:

D_i;t = β₁D_i;t−1 + β₂D²_i;t−1 + β₃FCR + β₄Crisis + β₅ΔCA_i;t + β₆ΔCA_i;t−1 + β₇CS_i;t + β₈CS_i;t−1 + β₉SI_i;t + β₁₀SI_i;t−1 + β₁₁L_i;t + β₁₂L_i;t−1 + β₁₃B_i;t-1 + β₁₄G_i;t + β₁₅FA + β₁₆EZ + β₁₇(FCR * ΔCA_i;t) + β₁₈(FCR * Crisis * ΔCA_i;t) + β₁₉(FCR * FA* ΔCA_i;t) + β₂₀(FCR * EZ * ΔCA_i;t)+ α_j + d_t +ν_j;t.

Here, FCR, CRISIS, FA & EZ are factor/ binary variables while rest are continuous. Following literature, I treated all financial variables to be potentially endogenous in nature. Furthermore, to resolve endogeneity concerns, I used t-2 to t-5 lagged levels for equation in differences as both the overidentification and serial correlation were satisfied at the 5^th lag. For the equation in levels, I used lagged differences dated t-1 as instruments. Also, both the firm fixed and time effects are included.

System GMM command

xtdpdgmm L(0/1).D D^2 I.FCR I.CRISIS CA CA_L1 CS CS_L1 SI SI_L1 L L_L1 B G FA EZ I.FCR#CA I.FCR#C.CA#I.CRISIS I.FCR#I.FA#C.CA I.FCR#I.EZ#C.CA , model(diff) collapse gmm(D_L1 CA CS CS_L1 SI L B G FCR FA , lag(2 5)) gmm(D^2 EZ I.FCR#CA I.FCR#C.CA#I.CRISIS I.CRISIS I.FCR#I.FA#C.CA I.FCR#I.EZ#C.CA , lag(1 2)) gmm(D_L1 CA CS CS_L1 SI L B G FCR FA , lag(1 1) diff model(level)) gmm(D^2 EZ I.FCR#CA I.FCR#C.CA#I.CRISIS I.CRISIS I.FCR#I.FA#C.CA I.FCR#I.EZ#C.CA, lag(0 1) diff model(level)) teffects two vce(r)

Can you please advise if the system GMM command is in accordance with above description of model.


Second, I have few ambiguities regarding marginal plots and contrast command. First, I ran marginal plots command for an interaction term of two factor (FCR& CRISIS) and 1 continuous variable (CA) which leads to following result:

margins fcr, at (ca = (-2 (0.5) 2) crisis=(0 1))
marginsplot, recast(line) x(ca) by(crisis) yline(0)

1ST CRISIS EFFECTS.gph

Because the confidence intervals of fcr are overlapping, does that mean that there is no significant difference between the two categories of fcr in terms of the effect of ca on D (dependent variable)?? Also, am I right to interpret that for crisis=0, three lines of fcr=1 above the red line are statistically significant and different from zero and the same is true for crisis =1 while for fcr= 1, 6 lines out of 9 show statiscally significant results and same goes for crisis=1. So, there is no effect of crisis?? and there is no statical difference between fcr=0 & 1 as confidence intervals are overlapping?

I also run the contrast command to compare the effects and I took fcr=1 as the reference group: Command is as follows:

margins rb1.fcr, at (ca = (-2 (0.5) 2) crisis=(0 1))
marginsplot, recast(line) x(ca) by(crisis) yline(0)

Results show the following effects:

fcr 1 CRISIS EFFECTS.gph

Further, I also run contrast command with fcr= 0 as the reference category

margins r.fcr, at (ca = (-2 (0.5) 2) crisis=(0 1))
marginsplot, recast(line) x(ca) by(crisis) yline(0)

fcr 0 CRISIS EFFECTS.gph

Can you please advise if I am doing it right and whether the interpretation is correct?

Second, when I ran system GMM, the results show a negative relationship between ca and D but when an interaction term is run between ca#fcr#crisis, the sign is positive. So, given the graphs ( fcr 1 CRISIS EFFECTS.gph fcr 0 CRISIS EFFECTS.gph) how do I interpret these 2 opposite direction FCR graphs?

(fcr & crisis are dummy/factor variables while ca is continuous variables).

Results of GMM are attached for reference: Generalized method of moments estimation.pdf

Please advise on interpretation of results and overlapping CIs .

[Sebastian Kripfganz]

1. It is slightly unusual not to specify the exogenous variables for the first-differenced model as well. I would add:
gmm(TDtoTA SalesGrowth lnTA CFOtoSales CapitalIntensity l.EnvDisclScore, lag(1 3) collapse)

2. The diff option and the model(diff) option do different things. diff first-differences the instruments, while model(diff) specifies that these are instruments for the first-differenced regressors.

3. There is no automatic differencing. You need to specify explicitly what you want the command to do.

[Rupali Vashisht]

Dear Prof. Sebastian Kripfganz ,

I am employing GMM estimation using the xtdpdgmm command. Following is my code:

xtdpdgmm L(0/1).EPS L(0/1).ROAA TDtoTA SalesGrowth lnTA l.EnvDisclScore i.t, model(diff) gmm(L(0/1).EPS L(0/1).ROAA, lag(2 3) collapse) iv(TDtoTA SalesGrowth lnTA CFOtoSales CapitalIntensity l.EnvDisclScore, diff model(level)) iv(i.t) igmm vce(cluster id)

My endogenous variables are L(0/1).EPS and L(0/1).ROAA while the other variables are exogenous. I have the following questions:

1. Do you think my code is correct?
2. I get very different results when I include model(level) suboption in iv() and without it. Could you please help me understand what it means to include 'diff' suboption alongside 'model(level)'suboption? Isn't the first difference transformation automatically done by xtdpdgmm when model (diff) suboption is specified?
3. Finally, are the variables mentioned in gmm() automatically differenced?

Thank you very much for your time.

Kind regards,
Rupali

[Sebastian Kripfganz]

If the dynamic model is the correct model, then any results you obtain for the static model are not reliable because they suffer from an omitted-variable bias. Consequently, the Hausman test for the static model is not reliable. The same is true when some independent variables suffer from endogeneity. You can estimate a dynamic model by GMM in this case, yes, but the Hausman test for the static model has no relevance.

[Sarah Magd]

Thanks Prof. Sebastian Kripfganz ,

Could you please explain your answer more?

I applied the robust standard Hausman test on my static model, and it says that the Random estimator should be applied. However, in the second step, I follow the dynamic specification for my analysis. Also, my independent variables suffer from potential endogeneity due to a reverse causality. Given these two points, can I use the sys-GMM estimator to estimate the dynamic model?

Thanks