data test;
   set quextion(keep=site);
   set quextion(keep=site rename=(site=nextsite) firstobs=2);
run;

proc sql;
create table test1 as
  select distinct site,nextsite,count(*) as num
    from test
        group by site,nextsite;
quit;

proc transpose data=test1 out=test2(drop=_name_) prefix=_;
   by site;
   var num;
   id nextsite;
run;

Dwyane, I don't know much about Macro.

没太看懂下面的那个矩阵的感觉

邓大，跑了一遍，貌似不太对，把log内容给贴上来了，你在帮忙看看。。。谢谢哈
NOTE: PROCEDURE PRINTTO used (Total process time):
      real time           0.01 seconds
      cpu time            0.00 seconds


208  data quextion;
209  infile datalines end=eof;
WARNING: The value of the INFILE END= option cannot be set for CARDS or DATALINES input.
210  input year site;
211  array x[5,5] _temporary_;
212  array a[5];
213  site_1=lag(site);
214  if _n_=1 then call pokelong(repeat(put(0,rb8.), 25), addrlong(x[1,1]), 25*8);
215  else x[site_1, site] = x[site_1, site]+1;
216  if eof then do i=1 to 5;
217  do j=1 to 5;
218   a[j] = x[i,j];
219  end;
220  output;
221  end;
222  keep a1-a5;
223  cards;

NOTE: The data set WORK.QUEXTION has 0 observations and 5 variables.
NOTE: DATA statement used (Total process time):
      real time           0.01 seconds
      cpu time            0.01 seconds

应该是我木有说清楚。。。在解释一下哈
其实，就是一个转移次数的频率计算，相邻两年间的site从一个值到另一个值记一次，譬如：
数据集a，如下：
year site
2000 1
2001 3
2002 2
2003 2
2004 3
那么该数据集的理论转移类型有9种，期频数分别记为：
1-1
1-2 2-1
1-3 3-1
2-2
2-3 3-2
3-3
然后，一次统计上述六种情况的频数，在列成表格，如下：
      1       2       3  
1   1-1   1-2    1-3
2   2-1   2-2    2-3
3   3-1   3-2    3-3目的就是想通过宏实现最后一张表格的统计结果。