Re: 大家來討論統計作業吧...
%macro work1(work);
%do i=1 %to &work;
%let k=1500;
data all;
set hs01a hs02a hs03a hs04a hs05a hs06a hs07a hs08a hs09a
hs10a hs11a hs12a hs13a hs14a hs15a hs16a hs17a hs18a hs19a
hs20a hs21a hs63a hs64a hs71a hs72a
;
random=ranuni(0);
run;
PROC SORT DATA=all;
BY random;
run;
DATA work2;
SET all (DROP=RANDOM);
IF _N_ LE &K;
RUN;
上面這個好像是在做抽樣的動作吧?!
proc means noprint;
output out=mean&i mean=mean;
var age;
run;
%end;
%mend;
%work1(100000);
這個是在做十萬次的平均數,值得注意的是out=mean&i地方!
data qq;
%macro merge;
set
%do i = 1 %to 100000;
mean&i
%end;
%mend;
%merge;
這個地方是在將mean1到mean100000合併到qq裡面。
proc gchart data=qq;
var mean;
run;
proc means data=qq;
run;
最後是畫圖和求平均。
這應該是可以跑的,雖然我沒有實際run十萬次,
我跑一千次就花了五分鐘,也就是十萬次要花五百分鐘,
老師上次開應該沒有開九小時吧?!我覺得很多地方都可以在簡化的樣子!
但其實我還有個小問題想問大家,
就是畫圖的時候,如何調X軸的組數阿?
都跑了一千筆了,但它都只給顯十幾組,看起來都好鬆散,實在不夠大器!!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.166.25.185
推
10/23 14:48, , 1F
10/23 14:48, 1F
推
10/23 18:51, , 2F
10/23 18:51, 2F
→
10/23 21:23, , 3F
10/23 21:23, 3F
推
10/23 23:01, , 4F
10/23 23:01, 4F
推
10/23 23:12, , 5F
10/23 23:12, 5F
→
10/23 23:12, , 6F
10/23 23:12, 6F
→
10/23 23:13, , 7F
10/23 23:13, 7F
→
10/23 23:28, , 8F
10/23 23:28, 8F
推
10/23 23:32, , 9F
10/23 23:32, 9F
推
10/23 23:33, , 10F
10/23 23:33, 10F
推
10/23 23:33, , 11F
10/23 23:33, 11F
→
10/23 23:40, , 12F
10/23 23:40, 12F
推
10/23 23:42, , 13F
10/23 23:42, 13F
討論串 (同標題文章)
以下文章回應了本文:
完整討論串 (本文為第 4 之 6 篇):