multinomial choice model with high-dimension fixed effects
I am running a regression with 1 million observation
outcome: occupation choice(5 category)
fixed effects: high school fixed effects(about 1000 groups), college fixed effects(about 800 groups), and neighborhood fixed effects (about 600 groups)