模拟数据输入：
set obs 10
gen year = 1990 +_n
expand 100
gen size = 100+ 10*invnormal(uniform())


根据题意我最先想到的方法：
sort year
gen group = 1
forval i = 10(10)90 {
 by year: egen p`i' = pctile(size), p(`i')
 replace group = (`i'/10 + 1) if size > p`i'
}
drop p*

可能比较简单的方法：
sort year size
by year: gen group2 = ceil(10*_n / _N)

验证通过：
assert group == group2

不知是否满足你的要求？也请论坛达人支招！

[此贴子已经被作者于2008-11-17 11:13:16编辑过]

不好意思，第二个方法确实有错——我原来的思路是先定义并计算rank和total两个新变量，然后再计算group2，但最终给出的是简化程序，无需定义rank和total，一忙乱却忘记对程序作相应调整了。现已重新编辑原帖子，第二种方法可以用了。

[此贴子已经被作者于2008-11-17 11:22:53编辑过]

可以直接用Pctile siz1=size, nq(10)

之后根据得到的10个分位点gen数据