Dear Stata enthusiasts,

Is there any command that can be used to estimate two or more logit models with seemingly unrelated regression while controlling for a long list of indicators (e.g., firm fixed effects)? I came across a paper that adopts this method but have no idea how to implement it in Stata. Thank you!