现有一个模拟回归分析，数据为全国29个省份的600多个县的个体数据，每县约几百个个体，其中有一个变量是表明个体代码的，十位数，前三位是省份代码，中间四位是市县代码，后三位是个体代码，例如“1432456001”，前三位“143”代表某省，中间四位“2456”代表某市某县，最后三位“001”代表某个人。总体样本量约为十万个，变量有四十多个，有些个体的变量里面有缺失值。我想做如下处理，如果一个县的某个人的某个变量有缺失，就用该县的平均值予以替代，如果某县全部缺失，则用该省的平均值予以替代，以此类推，这个怎么用stata命令实现呢？谢谢大家