如何用SAS宏计算这个公式？

数据集是一个Panel，结构如下：
country   year    vd     qw   qwe    fc     ......      sxb (50个变量)
china     2004     5     322   33      86    ......      322
china     2005     4     22     83      16    ......      75
china     2006     6     352   39      76    ......      32
.
.
.
china     2011     5     652   53      66    ......      82
us
us
us
us.
.
.
.
.
new_zealand
.
.
.
（100个国家）（2004-2011，共8年）     总计100*8=800个观测值

我要做的是：每年，对所有国家的每一个变量按公式换算成0-10的指数（连续的scale）。请大侠帮忙在SAS中实现。。。
比如说，对于2004年，v1变量有100个国家，即100个不同的观测值。将这些值通过公式计算成0-10的指数；同样，2004年v2也有指数。同理，2005-2011年。
这样变换之后仍然是一个panel，公式为：10*(max-obs)/(max-min)

请大侠将SAS代码贴在下面。。。

1. **create a test dataset;
   
2. data test;
   
3.   do country = 'c1','c2','c3';
   
4.    do i = 1 to 100;
   
5.      x = 2 * rannor(123) + 1;
   
6.          y = 3 * rannor(123) + 2;
   
7.      output;
   
8.    end;
   
9.   end;
   
10.   drop i;
    
11. run;
    
12. 
    
13. **create a macro variable for each variable;
    
14. proc sql;
    
15.   create table col as
    
16.   select name from dictionary.columns
    
17.   where libname='WORK' and memname='TEST' and name not in ('country');
    
18. quit;
    
19. 
    
20. data _null_;
    
21.   set col end=final;
    
22.   count+1;
    
23.   var = cats('v',count);
    
24.   call symputx(var,name);
    
25.   if final then call symputx('var_num',_n_);
    
26. run;
    
27. 
    
28. %macro rescale;
    
29. proc means data=test noprint nway;
    
30. class country;
    
31. var %do i = 1 %to &var_num;
    
32.       &&v&i
    
33.     %end;;
    
34. output out=out max= min=/autoname;
    
35. run;
    
36. 
    
37. proc sort data=test;
    
38. by country;
    
39. proc sort data=out;
    
40. by country;
    
41. run;
    
42. 
    
43. data new;
    
44.   merge test out;
    
45.   by country;
    
46.   %do i = 1 %to &var_num;
    
47.     n&&v&i = 10 * (&&v&i.._max - &&v&i)/(&&v&i.._max - &&v&i.._min);
    
48.   %end;
    
49.   keep country %do i = 1 %to &var_num;
    
50.                  n&&v&i
    
51.                %end;;
    
52. run;
    
53. %mend;
    
54. 
    
55. %rescale
    

复制代码

诸位大侠，揭榜啊~~~~~

一定要用宏？我觉得转个置就行了啊
下面是程序
data test;
infile 'E:\abc.txt';
input country $ year vd qw qwe fc ...... sxb
run;
proc sort data=test out=test1;
by year;
run;
proc transpose data=test1 out=test2;
by year;
id country;
run;
data test3;
set test2;
array raw(100) china us ...... new_zealand......;
maxz=max(china,us......new_zealand.......);
minz=min(china,us......new_zealand......);
do i=1 to 100
  raw(i)=10*(maxz-raw(i))/(maxz-minz);
end;
run;
proc Transpose data=test3 out=test4;
by         year;
id         _name_;
run;

啊，上面的程序少个冒号。另外国家里面应该有字符长于8个的的吧，所以还得加个length声明，我不知道你的国家最长是多少，暂时定个50吧，不够你自己改
data test;
infile 'E:\abc.txt';
length country $ 50;
input country $ year vd qw qwe fc ...... sxb;
run;
proc sort data=test out=test1;
by year;
run;
proc transpose data=test1 out=test2;
by year;
id country;
run;
data test3;
set test2;
array raw(100) china us ...... new_zealand......;
maxz=max(china,us......new_zealand.......);
minz=min(china,us......new_zealand......);
do i=1 to 100
  raw(i)=10*(maxz-raw(i))/(maxz-minz);
end;
run;
proc Transpose data=test3 out=test4;
by         year;
id         _name_;
run;

非常感谢二位！数据集太大的时候转置可能会比较麻烦，因此还是用宏比较稳妥。

3# guoluo

万分感谢！！！！

不过从我运行的结果来看，这个宏中计算0-10指数的时候，似乎是以每个国家来做的，即，v1变量的观测值在同一国家（c1国家）中转为0-10；我需要的是能将v1变量在同一年（如2004年）将所有国家的观测值变为0-10的指数。麻烦大侠改一下吧。另外，请设置一个txt文件，标价20论坛币，我来购买。。

改了一下

1. **create a test dataset;
   
2. 
   
3. data test;
   
4.   do country = 'c1','c2','c3','c4','c5';
   
5.    do year = 2004 to 2011;
   
6.      x = 2 * rannor(123) + 1;
   
7.      y = 3 * rannor(123) + 2;
   
8.      output;
   
9.    end;
   
10.   end;
    
11. run;
    
12. 
    
13. **create a macro variable for each variable;
    
14. 
    
15. proc sql;
    
16.   create table col as
    
17.   select name from dictionary.columns
    
18.   where libname='WORK' and memname='TEST' and name not in ('country','year');
    
19. quit;
    
20. 
    
21. data _null_;
    
22.   set col end=final;
    
23.   count+1;
    
24.   var = cats('v',count);
    
25.   call symputx(var,name);
    
26.   if final then call symputx('var_num',_n_);
    
27. run;
    
28. 
    
29. %macro rescale;
    
30. proc means data=test noprint nway;
    
31. class year;
    
32. var %do i = 1 %to &var_num;
    
33.       &&v&i
    
34.     %end;;
    
35. output out=out max= min=/autoname;
    
36. run;
    
37. 
    
38. proc sort data=test;
    
39. by year country;
    
40. proc sort data=out;
    
41. by year;
    
42. run;
    
43. 
    
44. data new;
    
45.   merge test out;
    
46.   by year;
    
47.   %do i = 1 %to &var_num;
    
48.     n&&v&i = 10 * (&&v&i.._max - &&v&i)/(&&v&i.._max - &&v&i.._min);
    
49.   %end;
    
50.   keep country year %do i = 1 %to &var_num;
    
51.                       n&&v&i
    
52.                     %end;;
    
53. run;
    
54. %mend;
    
55. 
    
56. %rescale

复制代码

也可以用stdize过程

1. %macro rescale2;
   
2. 
   
3. proc stdize data=test out=new2 method=range mult=10;
   
4. by year;
   
5. var %do i = 1 %to &var_num;
   
6.       &&v&i
   
7.     %end;;
   
8. run;
   
9. 
   
10. data new2;
    
11.   set new2;
    
12.   array newvar(*) %do i = 1 %to &var_num;
    
13.                  &&v&i
    
14.                %end;;
    
15.   do i = 1 to dim(newvar);
    
16.     newvar(i) = 10 - newvar(i);
    
17.   end;
    
18. run;
    
19. %mend;
    
20. 
    
21. %rescale2

复制代码

论坛币就不用了

8# guoluo
太感谢了！！！！！！！！！！！