The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Data » Dataset use in Stata » fixed effects and inflation factors
fixed effects and inflation factors [message #474] Sat, 25 May 2013 14:26 Go to next message
AdrianeW is currently offline  AdrianeW
Messages: 1
Registered: May 2013
Location: California
Member
Hello:
I'm trying to run a fixed-effects logistic regression to look at the association between marital status and hiv testing. I'm looking for advice on what regional variable I can use as a fixed effect? When I use the variable for state (v101), I get the error "Hessian is not negative semidefinite" and when I use the psu (s021), many observations get dropped because the # of women who get tested is limited.

Also, I recoded the marital status variable to be 3 dummies: never marred, currently married, and formerly married. Do I need to use the inflation facto (AWFACTT)? If so, how do I incorporate this?

Thank you for your consideration.
Re: fixed effects and inflation factors [message #477 is a reply to message #474] Sun, 26 May 2013 16:03 Go to previous message
Reduced-For(u)m
Messages: 292
Registered: March 2013
Senior Member


I think it probably depends on how you are trying to interpret your coefficients and what kind of variation you are hoping to use. Disclaimer out of the way, I generally do use V101 (region of residence) as a regional fixed-effect. You say you are losing a lot of observations, but the only way the regional fixed effect is causing that problem is if you have only 1 observation in that region, which seems pretty unlikely. Otherwise, you'd be losing all the observations anyway since the women weren't tested. Am I missing something?

How many women do you have who have been tested. Maybe make a "tested" variable, and "bysort V101: tab tested"... if you have more than 10 or 20 women per region, it's not the regional fixed-effects which are the problem. But I've never worked with the HIV testing data, so maybe there are just so few tested observations that you are not going to be able to run the FE estimator.

Are you doing it by reg Y `X' i.region? or xtreg Y `X', fe (or something like that)? Maybe you've defined your panel/time vars wrong in your xtset somehow? It seems odd to me that V101 wouldn't have enough observations to handle a fixed-effect.



Previous Topic: Panel data
Next Topic: (please) Need help with multi-family´s dwellings
Goto Forum:
  


Current Time: Thu Nov 28 12:31:58 Coordinated Universal Time 2024