I want to estimate a model with a fractional outcome variable following Papke & Wooldridge (2002). I run glm dep_var independent_var, family(binomial) link(logit) robust, right?
How can I add a three level fixed effects? My approach until now was to use egen new_var = group(var1 var2 var3)...