proc iml问题

1. proc iml;
   
2. use X1;
   
3. read all;
   
4. new_ac=J(500,1,0);
   
5. new_d_ar=J(500,1,0);
   
6. near=J(500,3,10000);
   
7. 
   
8. do i=1 to 500;
   
9. new_ac[i,1]=ac[i,1];       
   
10. new_d_ar[i,1]=d_ar[i,1];
    
11. diffni=abs(ni-ni[i,1]);       
    
12. diffni[i,1]=100;       
    
13. 
    
14. if revman[i,1]=1 | exman[i,1]=1 then do;
    
15. inds=indust[i,1];               
    
16. ye=year[i,1];       
    
17. do k=1 to 500;       
    
18. if indust[k,1]=inds & year[k,1]=ye then near[k,1]=diffni[k,1];         
    
19. if indust[k,1]=inds & year[k,1]=ye then near[k,2]=ac[k,1];        
    
20. if indust[k,1]=inds & year[k,1]=ye then near[k,3]=d_ar[k,1];       
    
21. 
    
22. end;
    
23. 
    
24. 
    
25. a=near[>:<,1];
    
26. new_ac[i,1]=near[a,2];
    
27. new_d_ar[i,1]=near[a,3];
    
28. end;
    
29. end;
    
30. 
    
31. result=new_ac || new_d_ar;
    
32. 
    
33. Create new  from result[colname={new_ac new_d_ar}];
    
34. Append from result;
    
35. Close new;       
    
36. 
    
37. quit;

复制代码

上面的程序中想要将“a=near[>:<,1]; ”相应的ID另存为数据集，请问要怎么修改？

不太懂iml语句，如果a在矩阵模块里是个矩阵的话，我用宏把它每个都导了出来，然后set到了一起



data a;
input a;
run;

%macro a;
proc iml;
use X1;
read all;
new_ac=J(500,1,0);
new_d_ar=J(500,1,0);
near=J(500,3,10000);

%do i=1 %to 500;
new_ac[i,1]=ac[i,1];        
new_d_ar[i,1]=d_ar[i,1];
diffni=abs(ni-ni[i,1]);        
diffni[i,1]=100;        

if revman[i,1]=1 | exman[i,1]=1 then do;
inds=indust[i,1];               
ye=year[i,1];        
do k=1 to 500;        
if indust[k,1]=inds & year[k,1]=ye then near[k,1]=diffni[k,1];         
if indust[k,1]=inds & year[k,1]=ye then near[k,2]=ac[k,1];         
if indust[k,1]=inds & year[k,1]=ye then near[k,3]=d_ar[k,1];        

end;


a=near[>:<,1];

create a_&i from a[colname=a];
append from a;
close a;

new_ac[i,1]=near[a,2];
new_d_ar[i,1]=near[a,3];
end;
%end;
result=new_ac || new_d_ar;

Create new  from result[colname={new_ac new_d_ar}];
Append from result;
Close new;        
quit;

%do i=1 %to 500;
data a;
set a a_&i;
run;
%end;

%mend;

%a;

intheangel 发表于 2014-3-12 13:23
不太懂iml语句，如果a在矩阵模块里是个矩阵的话，我用宏把它每个都导了出来，然后set到了一起

宏运行没有成功，不过感谢你的热心帮助