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

[consequence (consequence)]

我果真沒出息，凡事都求人，這是求人之後得到的結果.....

/*讀入檔案*/
DATA raw;
INFILE "G:/2000/data/**";
INPUT age 45-47;
run;

/*跑母體平均數*/
PROC MEANS DATA=raw;
VAR age;
run;

%MACRO main;

/*讓以下指令重複執行*/
%DO i=1 %TO 1000;
        /*隨機抽樣1500個樣本*/
PROC SURVEYSELECT DATA=raw METHOD=SRS REP=1 SAMPSIZE=1500 OUT=sample&i;

        RUN;

        /*計算信賴區間*/
PROC MEANS DATA= sample&i CLM ALPHA=0.05 MAXDEC=3;
        VAR age;
OUTPUT OUT=sampleCI&i UCLM=x1 LCLM=x2;
        RUN;
%END;



DATA result;
%DO i=1 %TO 1000;
        /*判斷信賴區間是否包含母體平均數*/
        SET sampleCI&i;
        IF x1>32.9745429>x2 THEN inCI=1;
        ELSE inCI=0;
        OUTPUT;
%END;
        RUN;
%MEND main;
/*結束MARCO*/

%main;/*執行MARCO*/



/*計算抽樣1000次，有多少次的信賴區間包含母體平均數*/
PROC FREQ DATA=result;
TABLE inCI;
RUN;

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.119.27.51
推 yiyi311:ㄜ...沒出息之人＋1