问问题还卖论坛币，我也是醉了

问了别人，说是在logit后者是probit回归的时候控制了年份跟行业，psmatch2就不用再控制了

还好只是醉了，乱点的哈哈

最后还是不知道如何用pamatch2来解决问题，helo psmatch2里有如下一段Matching within strata

    The following code illustrates how to match within exact cells and then calculate the average effect for the whole population.

        g att = .
        egen g = group(groupvars)
        levels g, local(gr)
        qui foreach j of local gr {
                psmatch2 treatvar varlist if g==`j', out(outvar)
                replace att = r(att) if  g==`j'
        }
        sum att
但我还是没改好，应该是能实现分组匹配的。由于我想实现的是1:1匹配，在论坛上看到不错的帖子，改了变量就能用了。
bys year ind rank (stkcd):g exp=_n if rank
bys year ind :egen s=max(exp)
drop if s==.
expandcl s if !rank,cl(year stkcd) gen(g)
bys year stkcd:replace exp=_n if !rank
bys year ind exp:egen r=sum(pscore*rank)
g q=abs((pscore-r)/r) if !rank

su exp
loc m=r(max)
forv i=1/`m'{
bys year ind rank (q exp stkcd): g exps=exp[`i'] if !rank
bys year ind rank (q exp stkcd): g stkcds=stkcd[`i'] if !rank
bys year ind rank (q exp stkcd): replace q=. if exp==exps&_n>`i'
bys year ind rank (q exp stkcd): replace q=. if stkcd==stkcds&_n>`i'
drop exps stkcds
}
drop if q==.&!rank

bys year ind exp:egen u=min(q) if !rank
drop if u!=q&!rank
egen mat=group(year ind exp)
drop exp-u
sort mat rank
详情见https://bbs.pinggu.org/thread-498689-2-1.html。希望stata大神能把分组匹配的语句给举例说明下。