我的数据在同一年、同一城市都有一个密度数据，我想要每个城市的密度按年份滞后一期。注意同一年份同一城市的数据有很多，不是只有一个，且因为数据库的原因，只有每条数据在样本中的家庭编码和在家庭中的编码（就是在样本中先以家庭为单位给个id，再在家庭中区分1、2、3）。求stata命令，要批量处理的那种（主要是前面的 xtset/tsset 不会写）。
年份       城市        密度   密度（-1）
2021        a            30        33
...            ...           ...         ...
2021        a            30        33
2021        b            20        23
...            ...           ...         ...
2021        b            20        23
2020        a            33
...            ...           ...
2020        a            33
2020        b            23