Dear all
I'm working in Stata 14.2.
I'm doing a two-level logit model with melogit. I'm interested in a cross-level interaction between the variable win (first level) and the variable camp (second level).
I think in this case, I have to take the first level variable in the random part of the model. Is this right?
Because win ist a factor variable i need the R. notation.
In the Stata manuals the .R notation ist normally used with the _all option. But this option takes all cluster together. I don't understand why this should makes sense.
So my question is what ist the difference between the following models? and which one should I take?

Code1
. melogit entscheid i.win##i.wissen i.win##i.prä i.win##i.Wirtschaft age i.educ i.sex i.lang i.win##c.camp || _all: R.win, cov(unstructured) intpoints(30)

Code2 with group identity variable
melogit entscheid i.win##i.wissen i.win##i.prä i.win##i.Wirtschaft i.win##ib2.zukunft age i.educ i.sex i.lang i.win##c.camp || idnr:R.win, covariance(unstructured) intpoints(30)

idnr ist the identifying variable.

And why is there only one var(_cons) in the output? Should not be there a variance for each category of win?

Dear all.

I'm doing my research on Unemployment, and my topic is "The Impact of remittances on unemployment." This is my model:

Unemployment_it = β0 + β1Unemployment_t-1,i + β2(FDI/GDP)_it +β3(GFCF/GDP)_it + β4PopGrowth_it + β5(Remittances/GDP)_it + β6Inflation_ti + β7FemaleLF_ti + β8YouthEmpRatio_ti + β9EnrollTertiary_ti + β10EnrollSecond_ti + β11GovernExp_{ti +} β12Crisis_ti + u<sub>ti


Thirteen countries, 12 variables, and more than 160 observations.

Can you advise me which methodology I should use for analysis? And don't forget please about endogeneity problem and lagged effect of previous years Unemployment to the current year...

Thank you in advance


Best regards,
Farrukh</sub>

Hi,

Im running some regressions to investigate the effects corruption has on enrolment rates. I am using panel data for different countries across different years. One of my control variables is lnGDP, but it was been suggested to me to use dummy variables for High income/Middle income/Low income countries. I have done this and I am getting some strange results. I would have thought that including lnGDP would mean I wouldn't need to include this dummy variable?

Any advice would be appreciated,

Thanks!

Hello,

I am new to Stata and am experiencing problems. There are several threads concerning this problem, for my specific case I did not find a solution.

I am trying to estimate a VAR regression, however Stata keeps returning "no observations r(2000);". In the simplest form, I am just entering "var VARIABLE1 VARIABLE2 VARIABLE3" into the command box.

I checked the data for missing values and string components but there are neither.

With the same variables, I have been able to perform a panel regression, my guess would therefore be that the "no observations" is due to the autoregressive component. Does anyone have an idea what could be the problem?

Thank you in advance
Kai

I'm studying the firm size distribution i need to represent it in a double logarithmic form, my problem is with the inverse cumulative distribution which is in my y axis, in the x axis i have the employment in the logarithmic scale. How can I compute the inverse cumulative distribution(right distribution)?

Reading through the mi append file right now and I fail to understand the difference between the two commands. Could someone please explain the difference to me.

The file says

.

What does this mean? I have 2 waves of data that I combined using append (quite new to stata and this is my first go at trying to do an empirical work) not knowing that mi append exists. In what ways have I messed up the data by using append vs mi append?

Overall I am trying to use the command mi xtset but failing horribly, so I am trying to work my way backwards and find what I have messed up. Thanks for any help!

Hi Members,

This question is in relation to generating variables following Multiple Imputation.

I am using multiple imputation for chained equations (MICE). The data is from a repeat cross-sectional dataset (n=643) with missing data for absolute income (n=345; 53.65%). In my original non-MI model I have used two alternative model specifications, one with log income and the second with the log of equivalised disposable household income.

I am looking to generate the equivalised disposable household variable, following MICE. Log income is entered as an imputed variable and I wanted want to know if it's possible to generate equivalised disposable household income from this variable following MI. Can it be done with mi xeq or mi passive?

For reference, log equivalised disposable household income is generated using the following code:

Many thanks
Becca

Hi statalist,

I am using the nldecompose command to examine intra cohort change and cohort replacement. I haven't been able to find any post estimation commands that would allow me to graph the "Char" (cohort replacement) and "Coef" (intracohort change). I have run the following commands where I get two separate outputs, one where I control for an individual's education and one where I don't:
Can anyone suggest if there is a way to graph these results?

I have tried following the advice in this post: https://www.statalist.org/forums/for...%80%9D-package

However, using the "coefplot" stata command returns the error:
"no coefficients found"

Good afternoon everyone,

I have a replication task, which involves the use of outreg command. In their replication code the authors use outreg in a very unfamiliar way (paper dates back to 2013). Below the code lines, in bold problematic lines :

######

regress voted_10g treat_5-treat_3 if s_r_v == `i', robust
testparm treat_5-treat_7
local secr = `r(p)'
testparm treat_1 treat_2
local sos = `r(p)'
lincom treat_5 + treat_6 - treat_7
test treat_1=treat_2

*Test average of secrecy interventions against SOS civic duty.

lincom (treat_5 + treat_6 + treat_7)/3 - treat_4

*First table does not include the covars.

outreg treat_5-treat_3 using "logs\Table2_OLSTreatmentEffects.out", se ctitle("`collabel'") addstat("F-Test p-value on three Secrecy Treatments",`secr',"F-Test p-value on two SOS Placebos",`sos',"Covariates included?",0) addnote("Full results with covariates in Table A2") bracket adec(6) rdec(6) 3aster `repl'

*Second table includes everything in model.

outreg using "logs\TableA4_OLSTreatmentEffects.out", se ctitle("`collabel'") addstat("F-Test p-value on three Secrecy Treatments",`secr',"F-Test p-value on two SOS Placebos",`sos') bracket adec(6) rdec(6) 3aster `repl'

######

Two errors can be reported :
varlist not allowed
and the fact that addstat is not an option.

The former is because they specify as an argument of outreg not a model name, but variable names. In the current version of outreg, one would rather write keep(treat_5-treat_3) as an option.

A solution anyone? I cannot write again the code (time consuming) and I already tried to use earlier versions of Stata (from 9.0 on). Authors wrote in their code to use version 10.1. But it does not work.

Thank you very much,

Emile T.

Hello dear community!

I am referring to the attached screenshort where you see two columns containing start and end dates with very hetereogeneous formatting. As of now, both columns are stored as strings since "UnKnown" was inserted everytime no date is given. I am now (1) trying to convert these strings to dates and (2) seeking to extract the year from all these dates. Step 2 is fairly straightforward, but I am struggling with the conversion. There are three different formats: e.g. a.) 18.Jul 14; b.) 1995 and c.) Feb 73. Only the last two digits of the years are shown in cases a.) and c.) so that I would need insert "20" whenever the number are [0,19] and a "19" otherwise. I have been experimenting a lot with date() but failed as of now.

Did someone already encounter a similar problem and could provide me with some insight?

Thanks!

Marie

Hello,

I am using Stata version 15.1 and would like to create a file on disk and append the contents of the current dataset in memory to the file on each iteration of a loop:

The scenario is as follows:

The Stata append command appears to be only for appending from disk to the current dataset while I wish to append in the other direction.

If there is no way of doing this I would appreciate any alternative suggestions, such as iterating overall all current rows in the dataset and writing them out one-by-one.

Hi,
I am doing a research paper on panel data taken from BRICS nations as representative of recently industrialized and emerging market developing economies . i am using stata
After declaring my dataset as panel data, I have appointed general pooled ols model to do the regression which gave me R-squared as only 0.7421. intuiting it as not sufficient i went to fixed effect(fe) model, but it further reduced my (overall) R-squared to 0.1251 . then, i did random effect (re) model which gave (overall) R-squared as 0.7421. then, i did hausman test to choose between fe and re, which gave chi-squared as negative to be -205.32 and thus, hausman failed

My questions:
Q.1) how to choose between pooled ols and fixed effect? Please describe command too, if possible.
Q.2) as hausman test failed, how can i choose between fe and re? Please describe command too, if possible.
Q3) does any other model exist that can be used for panel regression other than pooled ols, fe and re in general? does any of them can be used for my study? If yes, how to choose between them?
Q.4) what kind of other test can i employ in my study to check characteristics of my data?
Q5) does things like heteroscedasticity, unit root test, autocorrelation etc creates a problem in panel data? if it does, how can i check them?
Q.6) any other suggestion for my study is helpful and so highly appreciated.
Q7) I am log of variable to signify growth. am i right? if not, when do we take ln of dependent and independent variables?
Q8) does 72 observations sufficient for my study on government spending and unemployment?

thank you for any and all kind of help that canbe provided for the study.

Ps.1) please excuse my poor knowledge about econometrics. please answer as soon as possible as i have a close deadline to submit my work
2) i am extremely sorry if i asked some elementary questions which u may find not worth answering. in such a case, plese tell me the question number and if possible, the URL link where i can read about it.
3) Please excuse length of the post. i know some may find it lengthy but to provide every related information, i choose this structure.

Dear community,

I have a data set enlisting all previous/current occupations of executives (with start and end dates) which further specifies whether this position was at a quoted company or governmental. If sorted by directorid, you have all his/her previous/current positions grouped together. I now want to check whether his governmental career preceeded his company-career. So far I have generated two variables: gen min_gov = min(start) if companytype=="government" and gen min_quoted = min(start) if companytype=="quoted") and would like to assign these value to the WHOLE group and not just to thie specific row. Is there a way to do this?

Thanks in advance!

Marie

Hi,

I would really appreciate any guidance regarding the following:

I am working with HRS data, and have divided my three interest groups by unique healthcare needs and insurance coverage at wave 8. What I'd like to do is track at what point do individuals in each group transition to a different insurance status. For example:

1) Group A during wave 8 are covered by an employer's health insurance, at what point do some or all individuals in Group A transition to Medicare coverage?

I can't seem to figure out what command to run.

Thanks,
Asmaa

Greetings,

I'm running Stata 15.1 on OSX. My ultimate objective is to examine the extent that each US state's political ideology has varied over time. I'm unsure how to go about this though. My current approach--which is computationally demanding, given the size of the dataset (N=392,755)--is to do the following:

And so on and so forth. The obvious issue with this approach (which I'm hoping is surmountable) is that I'm forced to create separate mean-variables for each state. This is both time (and CPU) consuming and also leaves me with dozens of variables. I'm thus wondering if there is a way to store all the annual state means in a single variable. Thanks in advance for your input!

Here is some sample data:

Hi, I want to create a pie chart in the same way like in the link below

https://imaa-institute.org/m-and-a-by-industries/

Kindly help how to do this.

My data is Like this,

Hi, I want to create a pie chart in the same way like in the link below

https://imaa-institute.org/m-and-a-by-industries/

Kindly help how to do this.

My data is Like this,

Hello,

I have some doubts about panel regression in fixed effects:

1) When I run fixed effects should I always test for heteroskedasticity, , autocorrelation and serial correlation independence? Or can I draw my conclusions through the interpretation of fixed effects?

2) In the case the three tests are positive (applied to fixed effects model), how should I solve when N is greater than T and T is greater than N?

Thanks,
Mary

Dear Statalisterss

May you kindly assist me with this question. Through this forum i have learnt how to correct for heteroskedasticity and autocorrelation. The problem is i am getting an error message (Invalid syntax, in some cases it says option vce not allowed) ) and i am failing to fix the error. Before i add the command vce (cluster PanelID) the code runs smoothly without errors, but after adding it gives me an error. I have tried to use vce (Cluster PanelID) or vce(robust) but i get the same error message. I am using Stata/SE 15.0
here is the code i am using (i am pasting the first parts of my code)
Thats the code i am using and below i am pasting the output i am getting

panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe vce cluster (CompCode)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
option vce not allowed
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe cluster (CompCode)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout

. // Section 3
. local s3Varlist Bsize NEDs LnFsize Frsk CAge Lvg CEOwn Gops Block CFOSls InstOwn Indp CoAct ROA K3

. local varlist3 LnNEDelta LnNEDOwn NEDOpsDummy

.
. foreach ret_vari in `varlist3'{
2. qui xtreg res_2a `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe cluster (CompCo
> de)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s3a_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
8. local no=`no'+1
9.
. qui xtreg res_2b `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe cluster (CompCode)
10. local p=Ftail(e(df_m),e(df_r),e(F))
11. qui estat ic
12. matrix S=r(S)
13. local aic=S[1,5]
14. outreg2 using s3b_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
15. local no=`no'+1
16. }
invalid syntax
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=,fe cluster (CompCode)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
invalid syntax
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe vce(robust)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout

. // Section 3
. local s3Varlist Bsize NEDs LnFsize Frsk CAge Lvg CEOwn Gops Block CFOSls InstOwn Indp CoAct ROA K3

. local varlist3 LnNEDelta LnNEDOwn NEDOpsDummy

.
. foreach ret_vari in `varlist3'{
2. qui xtreg res_2a `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(robust)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s3a_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
8. local no=`no'+1
9.
. qui xtreg res_2b `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(robust)
10. local p=Ftail(e(df_m),e(df_r),e(F))
11. qui estat ic
12. matrix S=r(S)
13. local aic=S[1,5]
14. outreg2 using s3b_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
15. local no=`no'+1
16. }
invalid syntax
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe cluster(CompCode)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout

. // Section 3
. local s3Varlist Bsize NEDs LnFsize Frsk CAge Lvg CEOwn Gops Block CFOSls InstOwn Indp CoAct ROA K3

. local varlist3 LnNEDelta LnNEDOwn NEDOpsDummy

.
. foreach ret_vari in `varlist3'{
2. qui xtreg res_2a `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe cluster(CompCode
> )
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s3a_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
8. local no=`no'+1
9.
. qui xtreg res_2b `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe cluster(CompCode)
10. local p=Ftail(e(df_m),e(df_r),e(F))
11. qui estat ic
12. matrix S=r(S)
13. local aic=S[1,5]
14. outreg2 using s3b_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
15. local no=`no'+1
16. }
invalid syntax
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe vce(robust)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout

. // Section 3
. local s3Varlist Bsize NEDs LnFsize Frsk CAge Lvg CEOwn Gops Block CFOSls InstOwn Indp CoAct ROA K3

. local varlist3 LnNEDelta LnNEDOwn NEDOpsDummy

.
. foreach ret_vari in `varlist3'{
2. qui xtreg res_2a `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(robust)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s3a_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
8. local no=`no'+1
9.
. qui xtreg res_2b `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(robust)
10. local p=Ftail(e(df_m),e(df_r),e(F))
11. qui estat ic
12. matrix S=r(S)
13. local aic=S[1,5]
14. outreg2 using s3b_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
15. local no=`no'+1
16. }
invalid syntax
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe vce(cluster CompCode)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout

. // Section 3
. local s3Varlist Bsize NEDs LnFsize Frsk CAge Lvg CEOwn Gops Block CFOSls InstOwn Indp CoAct ROA K3

. local varlist3 LnNEDelta LnNEDOwn NEDOpsDummy

.
. foreach ret_vari in `varlist3'{
2. qui xtreg res_2a `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(cluster Comp
> Code)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s3a_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
8. local no=`no'+1
9.
. qui xtreg res_2b `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(cluster CompCod
> e)
10. local p=Ftail(e(df_m),e(df_r),e(F))
11. qui estat ic
12. matrix S=r(S)
13. local aic=S[1,5]
14. outreg2 using s3b_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
15. local no=`no'+1
16. }
invalid syntax
r(198);

end of do-file

r(198);

. do "C:\Users\01424576\AppData\Local\Temp\STD2158_ 0000 00.tmp"

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Data"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Dat
> a

. use database.dta, clear

. cd "C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Au
> to\Stata output\2. Fixed effects_Robust"
C:\Users\01424576\Desktop\Results\With Hausman Test\Original Model but -InstOwnBlockCont - Hetero+Auto\Sta
> ta output\2. Fixed effects_Robust

. // drop if Year==2010|Year==2011|Year==2012
. tsset CompId Year, yearly
panel variable: CompId (unbalanced)
time variable: Year, 2002 to 2016, but with gaps
delta: 1 year

. set more off

.
. // Section 1
. local no=1

. local varlist_depe OpsTC LnOpsIss LnNEDelta LnNEDOwn LnTC

. local s1_varlist K3 CoAct Time InstOwn CEOwn Gops Bsize LnFsize Lvg NEDs

. foreach vars in `varlist_depe'{
2. qui xtreg `vars' `s1_varlist' Ind1-Ind13 Yr* if s1!=.,fe vce(cluster CompCode)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s1_main.doc, append ctitle("`vars'") addstat( R_sq_within,e(r2), R_sq_between,e(r2_b),
> R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`s1_varlist') label
8. local no=`no'+1
9. }
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout
s1_main.doc
dir : seeout

. // Section 3
. local s3Varlist Bsize NEDs LnFsize Frsk CAge Lvg CEOwn Gops Block CFOSls InstOwn Indp CoAct ROA K3

. local varlist3 LnNEDelta LnNEDOwn NEDOpsDummy

.
. foreach ret_vari in `varlist3'{
2. qui xtreg res_2a `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(cluster Comp
> Code)
3. local p=Ftail(e(df_m),e(df_r),e(F))
4. qui estat ic
5. matrix S=r(S)
6. local aic=S[1,5]
7. outreg2 using s3a_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
8. local no=`no'+1
9.
. qui xtreg res_2b `ret_vari' c.`ret_vari'#c.K3 `s3Varlist' Ind1-Ind13 Yr* if s3!=.,fe vce(cluster CompCod
> e)
10. local p=Ftail(e(df_m),e(df_r),e(F))
11. qui estat ic
12. matrix S=r(S)
13. local aic=S[1,5]
14. outreg2 using s3b_main.doc, append ctitle("`ret_vari'") addstat( R_sq_within,e(r2), R_sq_between,e(
> r2_b), R_sq_overall,e(r2_o),F_test, e(F), Prob>F,`p',AIC,`aic') keep(`ret_vari' c.`ret_vari'#c.K3 `s3Var
> list') label
15. local no=`no'+1
16. }
invalid syntax
r(198);

end of do-file

r(198);

Hi! I am attempting to run a linear regression but keep getting an r(2000) error - no observations. I double checked all variables of interest are floats/int and there are no missing observations (using program/command mdesc). Any thoughts on what else I can do to get the regression to work? Using STATA 15 and there are a total of 133,515 observations.

The regression I am attempting to run is: xtreg std_Fi std_Fj_1 std_SameReg std_C4_3 c.std_Fj_1#c.std_SameReg c.std_Fj_1#c.std_C4_3

Variable Descriptions
std_Fi - standardized number of deals engaged in by Firm i
std_Fj_1 - standardized one year lag of number of deals engaged in by Firm j
std_SameReg -standardized same region for Firm i and Firm j
std_C4_3 - standardized concentration of top 4 firms in industry i