这个操作还真不清楚如何直接实现。

但是，我想可以采用两步法来实现：

假设模型为 y = a + b1*x1 + b2* x2 + c1*z1 + u

其中，b1，b2 为常数系数，而 c1 为可变系数。

第一步，采用 xtreg 命令估计那些你认为是常数的变量：
  xtreg y x1 x2, fe
    predict e_y, e      // 这一步相当于从 y 里剔除了 x1 和 x2 的影响，仅留下 z1 的影响
  xtreg z1 x1 x2
    predict e_z, e      //  这一步相当于从 z1 里剔除 x1 和 x2 的影响   
第二步，估计变系数模型
  xtsc e_y e_z        // 由此可以得到系数 c1
第三步，从新估计 b1 和 b2
    gen y2 = y - z1*_b[e_z]
    xtreg y2 x1 x2      // 由此得到 b1 和 b2 的无偏估计量

十分感谢连老师的指导