先bysort indid，然后计算，因为是分行业计算的

谢谢，但是我试了还是不行哦，因为每个行业里还有3个地区值。。。:
. bysort  indid: gen DIFF=  sumind -  sumind[_n-1]
结果
year        indid        prvcid        prod        sumind        DIFF
1999        11        120000        56423        68779       
1999        11        110000        12356        68779        0
2000        11        120000        66812        83709        14930
2000        11        110000        16897        83709        0
1999        12        120000        7842        13540       
1999        12        110000        5698        13540        0
2000        12        120000        5996        16231        2691
2000        12        110000        10235        16231        0
1999        13        110000        99632        109886       
2000        13        110000        78985        111985        2099
1999        13        120000        10254        109886        -2099
2000        13        120000        33000        111985        2099

bysort indid: gen re=(sumind-sumind[_n-1])/sumind[_n-1]
sort year

不是这样吗？

ps：你的数据也不是真正的面板数据

把数据collapse成每年每行业数据 然后计算
collapse (sum) prod, by(year indid)

所有程序为
collapse (sum) prod, by(year indid)
xtset indid year
gen base=L.prod
gen did=(prod-base)/base