SAS中文论坛

标题: 新手请教sas程序~ [打印本页]

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

作者: shiyiming 时间: 2008-10-30 17:09
标题: Re: 新手请教sas程序~
[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 时间: 2008-10-30 18:08
标题: Re: 新手请教sas程序~
十分感谢~~~

作者: shiyiming 时间: 2008-10-30 19:37
标题: Re: 新手请教sas程序~
再问一下~
在上面程序产生的数据集中要1、分别求取前5、25、100个变量的平均值，
2、给出 x1、m5、m25、m100的频数分布图；3、求出x1、m5、m25、m100的均值和方差
程序应该怎么写？先谢谢啦~~

作者: shiyiming 时间: 2008-10-31 11:40
标题: Re: 新手请教sas程序~
哪个好人愿意再帮帮我啊？？作业快交了。。。。。 <img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" />

作者: shiyiming 时间: 2008-10-31 16:31
标题: Re: 新手请教sas程序~
不明白你所说的频数表是对哪些值的频数
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 时间: 2008-10-31 18:54
标题: Re: 新手请教sas程序~
谢谢horse1~~~
我没有写清楚，那个m5,m25,m100就是前5，前25，前100个变量的平均值。。。。。

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)