求助SAS高手~~

sissi_baiye

1150

收藏 2011-03-16

万分感谢啦先~~

简介，用下面的公式设计一个SAS程序，生成n个随机数列

use something called a mixed congruential generator to manufacture our random numbers. It works

on the basis of recursion.

Xi+1 = (aXi + c)mod m, i=1,... n (1)

where X1 is the seed, a is the multiplier, c is the increment, and m =2^32-1.

Create sequence of uniform pseudo-random variates of user-speci ed length (likely 106) and parameters

of the congruential generator

扫码加我拉你入群

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

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

全部回复

sissi_baiye

2011-3-17 00:32:45

proc iml;
start ran_numbers(n,seed,multiplier,increment);
  modulars=(2^=32-1);
  thin=1;
  Store;
  do i = 1 to n by thin;
x1 = seed;
value(i)= x(i)*multiplier+increment
x(i+1)= mod((value(i),modulars);
Return ran_numbers;
output;
end;
finish ran_numbers;
  reset print;
  run;
quit;

这是我自己写的程序，但是do loop
                  value(i)= x(i)*multiplier+increment
x(i+1)= mod((value(i),modulars);
有问题，麻烦高手帮我改改

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群