proc iml;                                                                                                                              
m=                                                                                                                                      
{1  1  2   1,                                                                                                                           
2  4  2   1,                                                                                                                           
3  2  1   1,                                                                                                                           
3  2  3   0,                                                                                                                           
2  1  3   0,                                                                                                                           
1  1  2   2,                                                                                                                           
2  4  2   2,                                                                                                                           
3  2  1   2};                                                                                                                          
flag=unique(m[,4]);                                                                                                                     
sum=j(ncol(flag),4,0);                                                                                                                  
do i=1 to ncol(flag);                                                                                                                  
   x=m#repeat(m[,4]=flag[i],1,4);                                                                                                      
   sum[i,1:3]=x[+,1:3];                                                                                                                 
   sum[i,4]=flag[i];                                                                                                                    
end;                                                                                                                                    
prod=sum[#,1:3];                                                                                                                        
print sum;                                                                                                                              
print prod;                                                                                                                             
quit;

结果：
      SUM

5         3         6         0
6         7         5         1
6         7         5         2

             PROD

180       147       150

你是哪个 泉州-彬仔 吧？跑这来了 昨天那么问你也不说