求助：选分组配对用循环语句OR PSMATCH2？

数据为附件match2

现对数据稍作修改，把roa变成绝对值，其余不变(数据为附件的match2)，要求为“针对其中sp=1的公司，寻找sp=0，且同ind，同daq，roa最接近的样本（配对样本不能重复）”。

*原先的程序基本不变，只去掉变量year即可。另外，若配对样本观测值有多个，只取stkcd最小的那个。
use match2,clear
bys year ind daq sp (stkcd): g exp=_n if sp
bys year ind daq: egen s=max(exp)
drop if s==.
expandcl s if !sp,cl(year stkcd) gen(g)
bys year stkcd: replace exp=_n if !sp
bys year ind daq exp: egen r=sum(roa*sp)
g q=abs((roa-r)/r) if !sp
bys year ind daq exp: egen u=min(q) if !sp
drop if u!=q&!sp
egen mat=group(year ind daq exp)

bys mat sp (stkcd): replace exp=_n
drop if exp>1
drop exp-u

配对样本（即sp=0)不能重复的话，还要"expand cl s,cl(stkcd) gen(g)"干吗呢？这个命令不就是针对sp=1的样本比sp=0的样本还多的情形吗？

还要"expand cl s,cl(stkcd) gen(g)"干吗呢？

注意这里变量s的意义。

谢谢sungmoo。是这样的，输入命令“ duplicates report stkcd”，会发现配对样本有很多重复的：

Duplicates in terms of stkcd

--------------------------------------
   copies | observations       surplus
----------+---------------------------
        1 |          151             0
        2 |           38            19
        3 |            3             2
        4 |           12             9
        5 |            5             4
--------------------------------------
我前面指出条件有“配对样本（即sp=0）的不能重复”，如果重复，则删除sp=1的样本。

无尘第二 发表于 2012-3-15 22:36
谢谢sungmoo。是这样的，输入命令“ duplicates report stkcd”，会发现配对样本有很多重复的：

Duplica ...

“duplicates report stkcd”
这步是在什么时候执行？是在打开数据库之后立即执行吗？

不是打开数据库match1后立即执行。是跑完您给出的程序，完成配对之后。查看有无sp=0的配对样本重复，以及重复了多少个时用的这个命令。

不是打开数据库match1后立即执行。是跑完您给出的程序，完成配对之后。查看有无sp=0的配对样本重复，以及重复了多少个时用的这个命令。

*个人以为，看是否重复。应用(看mat的各值的频数是否超过2)：
ta mat

非常感谢sungmoo的耐心、细心解答！

sungmoo好！当year不是唯一时，上述程序运行遇到下列问题：
“expandcl s,cl(stkcd) gen(g)
expression is not constant within clusters
r(198);”
数据见附件match3.dta。
捣腾了好久，解决不了。。。。。。。