Hi statalist community

For cross-section data analysis, I have learned two ways of running regression ( that incorporate fixed effects)

Way 1: logit DepVar IndepVar1 IndepVar2 i.fixed_effect1 i.fixed_effect2, vce (robust)----this gives me Psudo R^2

Way 2: Two steps xtset fixed_effect1 and then xtlogit DepVar IndepVar1 IndepVar2 i.fixed_effect2, fe vce (bootstrap)--- it does not give Psudo R^2


Comparing Way 1 and Way 2, I can not use vce (robust) in Way 2. Moreover, Way 2 does not give Psudo R^2.

Is there any advantage in going for Way 2 instead of going for Way 1? Or is there any advantage in going for Way 1 instead of Way 2?