我之前是用R寫過大數法則以及***極限定理的小程式,但,在SAS上遇到困難,
我目前寫的程式如下:
DM LOG 'CLEAR';
DM OUTPUT 'CLEAR';

DATA one;
ARRAY k {100} k1-k100;
   DO n=1 TO 20;
   i=1;
       DO i=1 to 100;
       k(i)=RANNOR(i);
       END;
   OUTPUT;
   END;
RUN;
PROC PRINT NOOBS DATA=one;
RUN;

上半部似乎是沒有問題,因為可以print出來;但,下半部分,我想要把之前建立的one 中的k1 到 k100 combine 成一個k 矩陣,row=20,column=100 可以利用迴圈來處理嗎?

PROC IML;
/*read data file into iml*/
use one;
read all;
k=k1|| k2 || k3 || k4;  (希望可以一直combine到k100,但我只會手動,不會利用迴圈combine)
xmean=k[:,];  /*此部分是想要計算平均值*/
之後我還想要畫一個直條圖,但,在PROC IML下要如何畫呢?

RUN;


希望高手可以幫幫忙教我.


以上的部分為law of large  number....至於central limit theorem我還沒開始寫程式碼,
我想先把這部分弄懂後,也許就會寫出來了...
麻煩大家了...
p.s.
附檔類似的R程式碼,希望能幫助你們了解我要寫的SAS程式碼是什麼

T<-50
# T is sample size
N<-5000
# N is Replication number
Xmean<-matrix(nr=N,nc=1)
for (i in 1:N)
{

X<-rnorm(T,mean=0,sd=1)
Xmean[i,1]<-mean(X)
}
hist(Xmean)