如果没有缺失值，也可以简单地根据by组变量，分组统计一下观测数；

或者用ods：
data ex;
input
y x;
cards;
33.35 1854
34.78 1857
36.22 1948
37.66 2017
39.09 2092
40.53 2197
41.76 2608
42.81 2937
43.9 3119
44.94 3654
;
run;
ods output  NObs=obs;
proc reg data=ex;
model y=x;
run;
quit;

proc sort data=sashelp.class out=class;
   by sex;
run;


proc reg data=class;
   by sex;
   model weight=height;
   ods output nobs=n;
run;
quit;

模型输出的第一部分就会给出结果。读了多少观测。用了多少观测。

哦，可能我有点没表达清楚，我不是需要把观测输出到output中，我需要通过gvkey yrmth做几千个回归，再计算因变量预测值的平均值，并且把观测数小于15的删去。我是需要生成一张新表，这个表保存了以gvkey yrmth分组的每个回归的样本个数以及预测值，这样可以实现吗？

data ex;
input
y x category;
cards;
33.35 1854 1
34.78 1857 1
36.22 1948 1
37.66 2017 1
39.09 2092 2
40.53 2197 2
41.76 2608 2
42.81 2937 2
43.9 3119  2
44.94 3654 2
;
run;
proc sql;
create table ex1 as
select *,count(1) as count from ex
group by category
order by category;
quit;


proc reg data=ex1;
where count>=4;
by category;
model y=x;
output out=res(keep=category pred count) predicted=pred;
run;
quit;