Proc model data=original;
  Param b0 b1;
  y=b0+b1*x1;
  fit y/gmm kernel=(bart,laglength+1, lnM/lnT);
  instruments x1;
by co;
run;
lnM/lnT (M=bandwidth,T=no.of the observations)
本文来自: 人大经济论坛 SAS专版 版，详细出处参考： https://bbs.pinggu.org/forum.php?mod=viewthread&tid=1357777&page=1&from^^uid=859611

以上是我找的用于计算Newey-West修正方差的程序，现在由于参数不一样，所以不能用by co直接计算，需要单独指定参数。
能不能用类似于 when co=1这样的语句（直接这样写不行），或者是谁能给个自动将数据集拆分的程序。

再请问怎么才能让SAS自动根据BIC准则计算laglength？
bandwidth有Newey-West automatic选项吗？