**** random distribution;
data random;
retain x 0;
do i = 1 to 100 ;
y= ranuni(0);
z= ranuni(225);
x= x + 1;
output;
end;
run;
proc plot data= random;
plot x*y = 'a'
x*z = 'b'/overlay;
quit;
***normal disribution;
data normal;
retain x 0;
do i = 1 to 100 ;
y= rannor(0);
z= rannor(225);
x= x + 1;
output;
end;
run;
proc plot data= normal;
plot x*y = 'a'
x*z = 'b'/overlay;
quit;
这里有两段code,第一个是用ranuni()函数生成了两个函数值。再显示其分布,第二个是用rannor函数生成另外两个函数值。
下面是输出
第一个是随机函数产生的随机值分布
The SAS System 08:22 Sunday, November 25, 2006 4
Plot of x*y. Symbol used is 'a'.
Plot of x*z. Symbol used is 'b'.
x ?
?
100 ? b a a a b b
? a b b a b a b a
? b b a a b b a
? b b a b b a a a
? a a a b b b a b
80 ? a b a a b b b a
? a b b b a a b a
? a b a a b a b
? a b b b a a a b
? b a aa b b b a
60 ? b b b a a aa
? b a a b a a
? a b a b a b a b
? b b a b a a b a
? bb b a aa a
40 ? a b a b a b b a
? b a b a b baa
? a b a b a b a b
? ba a a b b ba
? b a b aa b b a
20 ? a bba ab a b
? a a b b a b a
? b b b a a a a
? a a b b b aa b
? b a b a bb a a
0 ? b a
?
妰儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍儍儍儍儍儓儍
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
y
这一个是正态分布的分布图
The SAS System 08:22 Sunday, November 25, 2006 5
Plot of x*y. Symbol used is 'a'.
Plot of x*z. Symbol used is 'b'.
x ?
?
100 ? b a a a b b
? b b b ab a a a
? a b a b b a a
? b b a b b a a
? a a a b a b bb
80 ? a a b b a b ba
? b a b a a b a
? a b a b b a b
? b a a a a bb
? b b a a a b a
60 ? aa b b b a a
? ba a b b b a
? a a b b b a b a
? a a b a ab b b
? b a b a b b a a
40 ? a b b ba b a a
? b a a a a b b b
? a a b a b b a
? a b b a bb a a
? ab b a b a b a
20 ? aa a b b b a
? a ba b a b ab
? b a bb a b a a
? a a a a bb b b
? a b b b b a a a
0 ? a b
?
妰儓儍儍儍儍垉儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍垉?
-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
y
从这里得出的结论:
1。随机函数中的seed 参数不影响抽样的随机性。换句话说,用不同seed 抽出来的两组样本集同方差的概率是>95%的,不会影响统计结果的,大家可以放心下seed
2。sas 提供了好几种随机函数,这些函数有各自的用途,大家要进行抽样的时候要仔细考虑一下用什么方式来进行随机抽样。