[討論] 作業四的第三部份 (SAS版)

[FFWillie (致敬吧！)]

本來以為很簡單結果也是搞了很久...
因為這次若是使用SAS的同學應該會不得不使用到MACRO這個語法
我覺得算是還滿有代表性的練習,所以po上來大家如果有什麼建議的話請大力鞭!


=======================程式碼分隔線=====================
data pop;
infile 'c:\2000\**\';
input age 45-47;
/*以上是讀取檔案*/

/*開始MARCO*/
%macro datago;
%do i=101 %to 1100;

proc surveyselect data=pop method=SRS rep=1
                  sampsize=1500  seed=&i out=srspop;
RUN;

proc means data=srspop;
output out=result mean=ma std=s;
run;

data out&i;
set result(keep=ma s);

x1=ma+1.96*(s/sqrt(1500));
x2=ma-1.96*(s/sqrt(1500));
if x1>32.9745429>x2 then i=1;else i=0;

data outall;
set out101~out1100;
/*這裡我是用excel直接把out後面的數字產出然後貼到程式碼裡,
非常的土法煉鋼,如果有人有好方法麻煩提供一下吧*/

%end;
%mend datago;
/*結束MARCO*/

%datago;/*執行MARCO*/

proc print;
run;

proc freq;
run;

===================================END=======================

以上就是這次作業的第三部份
我大概跑了10分鐘左右吧(沒有去量)
用SPSS真的要跑一兩個小時嗎?
那這樣如果上課前寫完不就跑不出結果orz
也不能在記中寫....

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.119.202.133
→ FFWillie:補充一下...我跑了兩次次數只差一次,都在95%內
→ FFWillie:我開始相信信賴區間是真的了orz
推 happyyaya:哈   我SPSS跑了5分鐘而已  我還以為會很久^^"
推 FFWillie:稍微改進一下....現在竟然只要跑2分鐘orz
→ shenaringo:是我自己電腦的問題才跑那麼久  它老了 ......orz
→ farlight:現在作業那麼難喔~需要寫到巨集?? 嘖嘖 辛苦了~~