求助SAS高手~~

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

万分感谢啦 先~~

简介，用下面的公式设计一个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

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：MANUFACTURE multiplier parameters Generator Parameter SAS

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);
有问题，麻烦高手帮我改改