使用SAS CALL RANUNI生成指定范围随机数的技术求助

核心问题修正点 #

1. 固定种子导致重复值：原代码每次调用CALL RANUNI都用固定种子2020，会生成完全相同的随机数，所有变量值一致。需改用变量存储种子，让系统自动更新种子值。
2. 范围转换逻辑缺失：RANUNI默认生成0-1的均匀分布随机数，需通过公式映射到10-20区间。
3. 循环匹配需求：调整循环结构，生成10组样本，每组包含15个随机变量。

修改后的代码（简化版） #

DATA seven;
  /* 初始化种子变量，调用RANUNI后会自动更新 */
  seed = 2020;
  /* 循环生成10组样本 */
  DO group = 1 TO 10;
    /* 数组定义15个变量，简化重复代码 */
    array x[15];
    /* 每组生成15个10-20的随机变量 */
    DO j = 1 TO 15;
      CALL RANUNI(seed, rand_val);
      /* 范围转换：下限 + (上限-下限)*随机数 */
      x[j] = 10 + 10 * rand_val;
    END;
    OUTPUT; /* 每组完成后输出一条观测 */
  END;
  DROP seed j;
RUN;

PROC PRINT DATA=SEVEN NOOBS;
RUN;

代码说明 #

• 种子处理：用变量seed存储初始种子，每次调用RANUNI后种子自动更新，确保生成不同的随机数。
• 范围转换：通过10 + 10 * rand_val将0-1的随机数映射到10-20区间，其中10是区间下限，10是上下限差值（20-10）。
• 数组简化：使用array x[15]一次性定义15个变量，避免重复写15次CALL RANUNI，代码更简洁。
• 循环逻辑：外层循环控制10组样本，内层循环生成每组的15个变量，每组完成后输出一条观测，完美匹配需求。

保留原变量名的版本 #

如果需要沿用类似原代码的变量命名方式（如x1-x15），可以使用以下代码：

DATA seven;
  seed = 2020;
  DO group = 1 TO 10;
    CALL RANUNI(seed, x1); x1 = 10 + 10 * x1;
    CALL RANUNI(seed, x2); x2 = 10 + 10 * x2;
    CALL RANUNI(seed, x3); x3 = 10 + 10 * x3;
    CALL RANUNI(seed, x4); x4 = 10 + 10 * x4;
    CALL RANUNI(seed, x5); x5 = 10 + 10 * x5;
    CALL RANUNI(seed, x6); x6 = 10 + 10 * x6;
    CALL RANUNI(seed, x7); x7 = 10 + 10 * x7;
    CALL RANUNI(seed, x8); x8 = 10 + 10 * x8;
    CALL RANUNI(seed, x9); x9 = 10 + 10 * x9;
    CALL RANUNI(seed, x10); x10 = 10 + 10 * x10;
    CALL RANUNI(seed, x11); x11 = 10 + 10 * x11;
    CALL RANUNI(seed, x12); x12 = 10 + 10 * x12;
    CALL RANUNI(seed, x13); x13 = 10 + 10 * x13;
    CALL RANUNI(seed, x14); x14 = 10 + 10 * x14;
    CALL RANUNI(seed, x15); x15 = 10 + 10 * x15;
    OUTPUT;
  END;
  DROP seed;
RUN;

PROC PRINT DATA=SEVEN NOOBS;
RUN;

内容的提问来源于stack exchange，提问作者zannatus saba