全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
1185 1
2012-03-30
一共两个问题
1.用SAS模拟正太分布 T分布等 过程 模拟正太分布的表(没听太清楚 就是一般课本后面都有的附表 让把正太的给模拟出来 完全瞎了 不知道怎么动手 求高手
2.在SAS中产生服从指数为3的100个随机数,用T检验有多少落在95%置信区间
跪求高手指导程序
过两天就要交了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2015-4-2 16:21:39
一、SAS随机数函数和CALL子程序
SAS系统产生随机数,两种方式,利用SAS函数(Functions)或者CALL子程序(CALL Routines),它们的语法格式是:
          方式           代码                         说明
函数 var=name(seed,) var为存储随机数列的变量,name为特定的分布函数形式,seed为随机数种子,为特定分布要求的参数(可选)
CALL子程序 call name(seed,,var) 同上,记得seed=0, ±1,±2, , ± (2**31-2)
SAS可用的随机数函数列表如下(可以参见SAS Help and Documentation-SAS Products-Base SAS-SAS Language Dictionary-Functions and CALL Routines-Functions and CALL Routines by Category):
        分布     SAS函数         说明
二项分布(Binomial) ranBin(seed,n,p) n为独立实验的次数,p为成功概率
指数分布(Exponential) ranExp(seed)   
正态分布(Normal) ranNor(seed),or normal(seed) ranNor和normal是同质的,但normal没有相对应的CALL子程序
泊松分布(Poisson) ranPoi(seed,m) m为均值
均匀分布(Uniform) ranUni(seed),or uniform(seed) ranUni和uniform是同质的,但uniform没有相对应的CALL子程序
柯西分布(Cauchy) ranCau(seed)   
伽玛分布(Gamma) ranGam(seed,a) a>0为形状参数
由分布律表格决定的离散分布(tabled probability distribution) ranTbl(seed,p1,p2,...pn)  ∑p=1
三角分布(Triangular) ranTri(seed,h)  h为斜边(最可能值)  
上面的随机函数,除了normal和uniform,都可以由直接相应的CALL子程序调用。

例子
创建两个随机数变量,各包含3个记录,其中x1的种子为123,x2的种子为456:
data ranuni(drop=i);
    retain seed1 123 seed2 456;
    do i=1 to 3;
        x1=ranuni(seed1);
        x2=ranuni(seed2);        
        output;
    end;
run;
proc print data=ranuni;run;
结果为:
Obs    seed1    seed2       x1         x2
1      123      456     0.75040    0.32091
2      123      456     0.17839    0.90603
3      123      456     0.35712    0.22111
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群