请问在stata中如何编写这样一个函数,如附件所示。非常感激，自己研究了很久着实写不出来。