你们好，我在做一个回归的时候，probit模型中有i.area，加了这个地区控制变量以后，stata的mfx估计边际效应不能用了，请问下应该怎么做。