Hello, I'm trying to implement the -logitfe- model shown below. I want to lag my outcome variable Dental_Wave, by 1 time period, but I'm getting a numeric overflow problem. My dataset consists of 3 time periods and Observations: 27,447. I have included my code with the error below, any advice would be helpful.
Code:
sort year
xtset aid year
. logitfe dental_Wave ///
> L1.dental_Wave /// <- 1 to 4 lags of outcome
> c.depressive_symptoms_Wave c.age_Wave i.smoked_marijuana_alt_Wave i.smoked
> _cigarettes_Wave ///
> i.alcohol_Wave i.education_Wave, ///
> analytical ieffects(yes) teffects(yes)
Computing uncorrected fixed effects estimator
note: multiple positive outcomes within groups encountered
note: multiple positive outcomes within time periods encountered
note: 5913 groups (11625 obs) dropped because of all positive or
all zero outcomes
Iteration 0: f(p) = -3957.5092
Iteration 1: f(p) = -2372.2041
Iteration 2: f(p) = -2222.6291
Iteration 3: f(p) = -2196.7665
Iteration 4: f(p) = -2191.5725
Iteration 5: f(p) = -2190.3398
Iteration 6: f(p) = -2189.9916
Iteration 7: f(p) = -2189.9015
Iteration 8: f(p) = -2189.8768
Iteration 9: f(p) = -2189.8702
Iteration 10: f(p) = -2189.8687
Iteration 11: f(p) = -2189.8685
Iteration 12: f(p) = -2189.8685
Computing analytical correction
initial values not feasible
r(1400);
end of do-file
r(1400);
