The multiple estimation syntax fails if combinations of variables are provided (f1+f2 in the example). This also fails in the formula parser refactor #1118, so will provide a fix there.
Example
import pyfixest as pf
data = pf.get_data(N=1_000, seed=0, model="Feols")
fml = "Y ~ X1 | sw0(f1, f2, f1 + f2)"
pf.feols(fml=fml, data=data)
My use case is running a specifications with all possible combinations of fixed effects (e.g., no fixed effects, unit fixed effects, time fixed effects, and unit and time fixed effects). fixest implements a special "multiverse stepwise" syntax for this, see #1136.