求问大家，我用stata导出 probit的结果时，要导出边际系数和pr2等结果，就用margins, dydx(_all) post 语句将原本的回归系数替换为边际效应，但是这样导出结里的pseudo R2 等就是空的，应该是被margins全部替代了。请问大家代码该如何修改才能既使系数为边际系数，又有pseudo R2、adj. R2等结果呢？
代码如下

probit y1 x1 x2 x3
margins, dydx(_all) post
est store m1

probit y2 x1 x2 x3
margins, dydx(_all) post
est store m2


probit y2 x1 x2 x3
margins, dydx(_all) post
est store m3


esttab m1 m2 m3 using myfile.rtf,r2 ar2 pr2 se(5) b(5) star(* 0.1 ** 0.05 *** 0.01) replace