Re: [程式] SAS選擇最近的值後轉置

[west1996 ()]

注意：沒有data可以測，憑想像空打，程式碼可能有錯，請自行debug一下


這支macro可以一次轉一個變數一個時間點，餵給他三個東西
1. 資料集名稱
2. 變數名稱
3. 月數(要給一個數字)
macro會產生該名稱該變數的結果在work裡，程式碼最後一行的範例是把yourdata這個
資料集的GOT變數的3個月資訊轉出來存成work.out_GOT_3M


P.S.要一次轉多個變數多個時間點需要另外加工包裝macro，比較麻煩，既然你的欄位
    不多，時間點也不多，data也不大，所以乾脆多跑幾次macro再自行merge
    資料就好，速度上不會差很多

程式碼如下：

%macro aggregatedata(dsn=, var=, month=);
data work.out_&var._&month.M;
        set &dsn.;
        by id;
        retain temp_value temp_month over_value over_month over_flag;
        length &var._&month.M $ 12;
        if first.id=1 then do; temp_value=.; temp_month=.; over_value=.; over_month=.; over_flag=0; end;
        if &month.-1<= month <= &month.+1 and &var. ^= . and temp_month ^= &month. and over_flag ^= 1 then do;
                if month <= &month. then do;
                        temp_value=&var.;
                        temp_month=month;
                end;
                else if month> &month. then do;
                        over_value=&var.;
                        over_month=month;
                        over_flag=1;
                end;
        end;
        if last.id=1 then do;
                if temp_value = . and over_value= . then &var._&month.M='.';
                else if temp_value=. and over_value ^=. then &var._&month.M=cats(over_value,'(',over_month,')');
                else if temp_month = &month. then &var._&month.M=cats(temp_value,'(',temp_month,')');
                else if &month.-temp_month <= over_month-&month. then &var._&month.M=cats(temp_value,'(',temp_month,')');
                else if &month.-temp_month > over_month-&month. then &var._&month.M=cats(over_value,'(',over_month,')');
                end;
                output;
        end;
        keep id &var._&month.M;
run;
%mend;

%aggregatedata(dsn=yourdata, var=GOT, month=3)


※ 引述《yonny (悠逆)》之銘言：
: [軟體程式類別]: SAS (R也可以 但是我跟R比較不熟)
: [程式問題]: 資料處理
: [軟體熟悉度]:熟悉
: [問題敘述]:
: 我的資料大概長像這樣
: No是每次檢驗的流水號 ID相同是同一個人
: GOT和GPT是在不同時間點做的檢驗數值
: 但不是每個時間點 都有兩個檢查都做
: 而且不同人做檢查的時間點也都不相同
: No ID Month GOT GPT
: 1   1  0.2  41  34
: 2   1  2.8  42  .
: 3   1  3.5  43  36
: 4   1  3.7  44  37
: 5   1  4.9  45  38
: 6   1  5.5  51  39
: 7   1  5.7  52  40
: 8   1  6.0  53  .
: 9   1  6.2  54  .
: 10  2  0.9  20  21
: 11  2  4.0  22  .
: 12  2  4.1  24  25
: 13  2  5.7  .   36
: 14  2  6.0  30  .
: 以下是我期望的資料型態
: ID GOT_0M   GOT_3M  GOT_6M   GPT_0M  GPT_3M  GPT_6M
: 1  41(0.2)  42(2.8) 53(6.0)  34(0.2) 36(3.5) 40(5.7)
: 2  20(0.9)  22(4.0) 30(6.0)  21(0.9) .       36(5.7)
: 我想要把資料轉置成橫的,以GOT為例,
: 我要創造三個變項 GOT_0M GOT_3M GOT_6M (分別為0,3,6個月的GOT)
: 這三個時間點若無完全符合的值, 則抓正負一個月內的值來代替
: 若正負一個月內都無值, 則missing
: 值後面括號(是用哪個時間點的值)
: ---
: 以上大概是我的資料型態
: 實際的資料 有好幾百人的資料
: 每個人大概都有30~100個檢驗的時間點
: 檢驗數值(如GOT,GPT之類的) 大概有10幾個
: 而我們要取的時間點實際上也不只有三個(0,3,6,12,18,24,36,48M)
: 感覺有點複雜(囧)
: 我目前大概查一下 可能用PROC SQL先寫出來
: 再看改成MACRO
: (我今天掙扎了一下午 只寫了一點點簡單的 完全離目標很遠XD)
: 想請問看板上的高手有無建議要用什麼statement寫比較適合?
: 我再研究一下
: 或是若有高手可以幫忙寫出來
: 小小3000P表達感謝!
: 非常謝謝!!

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 220.136.120.240
※ 文章網址: https://www.ptt.cc/bbs/Statistics/M.1458928188.A.1CC.html
推 yonny: 我試過大概沒問題了(倒數第二個end是多的) 非常謝謝!!
→ yonny: 其實我實際資料欄位*時間點 可能會有100多個 哈哈
→ yonny: 不過或許可以用MACRO+do loop合併看看
→ yonny: 我自己是用PROC SQL寫 可以同時轉多個時間點和變項
→ yonny: 只是寫法比較耗資源 或許資料大一點 時間要比較多
→ yonny: 差最後一些小地方還沒完成 之後寫完再跟您分享
→ yonny: 非常感謝你!!!!!!
→ west1996: 100多個對我來說算是小拉，因為ctrl+C, ctrl+V一下就結
→ west1996: 束了，要包成可以處理多個時間點多個變數難其實不是難
→ west1996: 在寫程式，而是需要規範input的變數清單與時間點清單的
→ west1996: 撰寫格式，這需要溝通，所以我當初才沒包這一段，否則
→ west1996: 包成那個版本其實比較有利，因為無論多大的資料量都只需
→ west1996: 要掃過一遍資料就處理完了~
推 yonny: 謝謝你啦! 下面我也有回我的作法 給您參考~