範 例 七： 電泳膠片內的環帶頻率排序
 
 
程式檔名稱 : linEAR7.htm

data a;
input  id $ 1-7 primer $ 10-13 @ 19 bands $ 19-58;    讀入 primer 不同樣品電泳膠片內的條帶
l=length(bands);                                      判斷電泳膠片內的條帶數
marker_j=0;
   do until (marker_j = l);
      marker_j=marker_j+1;
      bandi=substr(bands,marker_j,1);
      if bandi='+' then bandiv=100;                   將字串 "+" "-" 值轉為 100, 0  的數值(bandiv) 
      if bandi='-' then bandiv=0;
      drop gel;
   output;
  end;

cards;
7971156  ak02  1  ----+-+-+----------+--------------------
7971157  ak02  2  ----+---+-------------------------------
7971158  ak02  3  ----+-+-+-------------+-----------------
7971162  ak02  5  ------+-+-------------------------------
7971167  ak02  6  ------+-+----------+-----------------+--
7971168  ak02  7  ------+-+-------------------------------
7971169  ak02  8  ----+-+-+-------------------------------
7971171  ak02  9  ----+-+-+----------------------------+--
7971174  ak02 10  ----+-+-+-------------------------------
7971178  ak02 11  ----+-+-+-------------------------------
7971179  ak02 12  ----+---+-------------------------------
7971184  ak02 13  ----+-+-+-------------------------------
7971196  ak02  1  ----+-+-+----------------------------+--
7971198  ak02  2  ----+-+-+----------------------------+--
7971296  ak02  1  --------+-------------+--+-+------------
7971315  ak02  2  --------+------+------+--+--------------
7971317  ak02  3  --------+------+------+--+--------------
7971319  ak02  4  --------+------+------+--+----+---------
7971322  ak02  5  --------+----------------+--------------
7971324  ak02  6  --------+------+---------+----+---------
7971329  ak02  8  --------+------+----+-+--+----+---------
7971332  ak02  9  --------+-------------+--+--------------
7971334  ak02 10  --------+------+---------+----+---------
7971336  ak02 11  --------+-----+----------+----+---------
7971338  ak02 12  ------+++-------------+--+--------------
7971340  ak02 13  ----+---+-------------+--+-+------------
;
run;

proc summary NWAY;                                電泳膠片內的環帶頻率計算
class primer marker_j;
var bandiv;
output out=bands mean=FQ n=n;
id primer marker_j;
run;

proc sort data=bands out=bands;                   電泳膠片內的環帶頻率排序
by descending FQ;
run;

data bands;
set bands;
run;

proc print noobs;
run;



輸出 7-1  電泳膠片內的條帶頻率排序(降序排列)


 PRIMER    MARKER_J    _TYPE_    _FREQ_         FQ     N

  ak02         9          3        26      100.000    26
  ak02         7          3        26       50.000    26
  ak02         5          3        26       46.154    26
  ak02        26          3        26       46.154    26
  ak02        23          3        26       34.615    26
  ak02        16          3        26       23.077    26
  ak02        31          3        26       19.231    26
  ak02        38          3        26       15.385    26
  ak02        20          3        26        7.692    26
  ak02        28          3        26        7.692    26
  ak02         8          3        26        3.846    26
  ak02        15          3        26        3.846    26
  ak02        21          3        26        3.846    26
  ak02         1          3        26        0.000    26
  ak02         2          3        26        0.000    26
  ak02         3          3        26        0.000    26
  ak02         4          3        26        0.000    26
  ak02         6          3        26        0.000    26
  ak02        10          3        26        0.000    26
  ak02        11          3        26        0.000    26
  ak02        12          3        26        0.000    26
  ak02        13          3        26        0.000    26
  ak02        14          3        26        0.000    26
  ak02        17          3        26        0.000    26
  ak02        18          3        26        0.000    26
  ak02        19          3        26        0.000    26
  ak02        22          3        26        0.000    26
  ak02        24          3        26        0.000    26
  ak02        25          3        26        0.000    26
  ak02        27          3        26        0.000    26
  ak02        29          3        26        0.000    26
  ak02        30          3        26        0.000    26
  ak02        32          3        26        0.000    26
  ak02        33          3        26        0.000    26
  ak02        34          3        26        0.000    26
  ak02        35          3        26        0.000    26
  ak02        36          3        26        0.000    26
  ak02        37          3        26        0.000    26
  ak02        39          3        26        0.000    26
  ak02        40          3        26        0.000    26