Dear all,
I am having problem finding a reliable code for my country fixed effects model.
I have 21 countries in my database and individuals nested within them. I am
running a multilevel analysis first, and then I am using the fixed effects
approach to check the robustness of my findings. At this point, I need to
include in my model some cross-level interaction effects between the country
level factors and individual level variables. However, I am not sure if my r
code is right:
model1<-lm(safety ~ mixed_neigh + ethnic_neigh + age + gndr + eduyrs + domicil
+ partner +
tvpol + income + victim + trust + trustXprison_pop +
trustXforeign_pop + victimXprison_pop +
victimXforeign_pop + mixed_neighXprison_pop +
mixed_neighXforeign_pop + ethnic_neighXprison_pop +
ethnic_neighXforeign_pop + factor(cntry)-1, data=mydata)
As you can see, the interaction are between a individual level variable (such
as trst, victim, mixed_neigh and ethnic_neigh) and two country level variables
(prison population and foreign population). The code runs fine, so I don�t have
warning messages and it is all spelled correctly. I just want to know if you
think that this is the right code to do a fixed effect analysis (only cross
sections and not cross time). Indeed, I I used the lm and not ppm because my
dataset is from only one year. Do you think it is correct?
Thank you very much for your attention,
Cristina
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
