新手请教sas程序~

shiyiming

刚学sas不久，碰到一个题目向大家请教：编写SAS程序，采用随机函数，产生1000行100列（表示100个变量，表示为x1,…,x100，各观察了1000次）符合正态分布N(0,1)的随机数，存入永久数据集。
看到有个帖子和这个内容差不多，但是这个题目要求产生100个变量，这是怎么做的？还有直接用rannor(seed)函数可以么？
先谢谢大家的帮助了~~

shiyiming

[code:4iunuwfp]data aa;
     array nor{100} x1-x100;
         do i=1 to 1000;
            do j=1 to 100;
                   nor(j)=normal(1);
                end;
                output;
         end;
run;[/code:4iunuwfp]

shiyiming

十分感谢~~~

shiyiming

再问一下~
在上面程序产生的数据集中要1、分别求取前5、25、100个变量的平均值，
2、给出 x1、m5、m25、m100的频数分布图；3、求出x1、m5、m25、m100的均值和方差
程序应该怎么写？先谢谢啦~~

shiyiming

哪个好人愿意再帮帮我啊？？作业快交了。。。。。 <!-- s:( --><img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" /><!-- s:( -->

shiyiming

不明白你所说的频数表是对哪些值的频数
data aa;
     array nor{100} x1-x100;
    do i=1 to 1000;
       do j=1 to 100;
         nor(j)=normal(1);
      end;
      output;
    end;
run;

data bb;               /*球m5 m25 m100*/
     set aa;
         m5=mean(of x1-x5);
         m25=mean(of x1-x25);
         m100=mean(of x1-x100);
         keep x1 m5 m25 m100;
run;

/*以下是求x1的均值和方差， 对m5 m25 m100 可同样做*/
proc means data=bb noprint;
     var x1;
         output out=mx1 mean=mn1 std=std1;
run;

data mx1;
     set mx1;
         variance=std1**2;
         drop _freq_ _type_;
run;

shiyiming

谢谢horse1~~~
我没有写清楚，那个m5,m25,m100就是前5，前25，前100个变量的平均值。。。。。